Gráficas y simualción

% 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

n_values = -15:15; % Cambiando el rango a -15 a 15

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);

% Calcular el espectro de amplitud c(n)

c = sqrt(an.^2 + bn.^2) / 2;

% Graficar el espectro de amplitud

figure;

stem(n_values, c, 'o');

title('Espectro de Amplitud (Simétrico)');

xlabel('n');

ylabel('c(n)');

grid on;

[pic 3]

% Par% Parámetros dados

T = 9;

w0 = 2 * pi / T;

a0 = 6;

% Definición de la función an y bn con valor absoluto de n

n_values = -15:15; % Cambiando el rango a -15 a 15

abs_n_values = abs(n_values);

an = -1.5 * cos(2 * pi * abs_n_values / 3) ./ (pi^2 * abs_n_values.^2) + 5 * sin(2 * pi * abs_n_values / 3) ./ (pi * abs_n_values) - 7.5 ./ (pi^2 * abs_n_values.^2) - 1.5 * cos(4 * pi * abs_n_values / 3) ./ (pi^2 * abs_n_values.^2) + 6 * sin(4 * pi * abs_n_values / 3) ./ (pi * abs_n_values) + 10.5 ./ (pi^2 * abs_n_values.^2);

...