Gráficas y simualción
Enviado por Hopepr • 29 de Octubre de 2024 • Práctica o problema • 668 Palabras (3 Páginas) • 17 Visitas
% Definir la función por partes
t1 = 0:0.01:3;
y1 = 5 * t1 / 3;
t2 = 3:0.01:6;
y2 = 2 * t2 - 6;
t3 = 6:0.01:9;
y3 = 7 * t3 / 3 - 14;
% Combinar las partes de la función
t = [t1, t2, t3];
y = [y1, y2, y3];
% Graficar la función
plot(t, y, 'LineWidth', 2);
% Agregar etiquetas y título
xlabel('t');
ylabel('f(t)');
title('Gráfico de la función por partes');
% Mostrar la cuadrícula
grid on;
% Agregar leyenda para cada parte
legend('5t/3, 0<=t<3', '2t-6, 3<=t<6', '7t/3-14, 6<=t<9');
[pic 1]
% Parámetros dados
T = 9;
w0 = 2 * pi / T;
a0 = 6;
% Definición de la función an y bn
n_values = 1:50;
an = -1.5 * cos(2 * pi * n_values / 3) ./ (pi^2 * n_values.^2) + 5 * sin(2 * pi * n_values / 3) ./ (pi * n_values) - 7.5 ./ (pi^2 * n_values.^2) - 1.5 * cos(4 * pi * n_values / 3) ./ (pi^2 * n_values.^2) + 6 * sin(4 * pi * n_values / 3) ./ (pi * n_values) + 10.5 ./ (pi^2 * n_values.^2);
bn = -1.5 * sin(2 * pi * n_values / 3) ./ (pi^2 * n_values.^2) - 5 * cos(2 * pi * n_values / 3) ./ (pi * n_values) - 1.5 * sin(4 * pi * n_values / 3) ./ (pi^2 * n_values.^2) - 6 * cos(4 * pi * n_values / 3) ./ (pi * n_values) - 7 ./ (pi * n_values);
% Definición de la función f(t)
f_t = @(t) a0 / 2 + sum(an .* cos(n_values * w0 * t) + bn .* sin(n_values * w0 * t));
c_m =
% Generar puntos de tiempo para graficar (por ejemplo, 1000 puntos en un periodo)
t = linspace(0, T, 1000);
% Calcular los valores de la función f(t) para cada punto de tiempo
y = arrayfun(f_t, t);
% Graficar la función
figure;
plot(t, y);
title('Gráfica de f(t)');
xlabel('Tiempo (s)');
ylabel('f(t)');
grid on;
[pic 2]
% Parámetros dados
T = 9;
w0 = 2 * pi / T;
a0 = 6;
% Definición de la función an y bn
...