Otros Obtener la señal resultante de pasar una señal periódica a través de un sistema lineal con SCILAB

[Katerina]

Tomar la señaltemporal correspondiente al desarrollo en serie de Fourier de un trende pulsos rectangulares definido en unintervalo de tiempo de 100 s centrado en t = 0. Utilizar 1000 armónicos para componer la señal:

Utilizando el teorema derespuesta lineal y el principio de superposición, obtener la señalresultante de pasar esta señal periódica a través del sistemalineal e invariante con el tiempo representado por la siguientefunción de transferencia, donde ω₀ es un parámetro quedetermina el ancho de banda del mismo:

Obtener resultados paralos tres anchos de banda del sistema ω₀ = 2π/5, 2π y10π rad/s, presentando superpuestas la señal temporal antes ydespués de pasar por el sistema para cada ancho de banda, con el finde observar y comentar mejor las distorsiones que produce el sistema.Explicar detalladamente cómo se ha implementado en el código elteorema de respuesta lineal y el principio de superposición.
NOTA: Se recuerda que elteorema de respuesta lineal expone que una señal con un únicoarmónico que pasa a través de un sistema lineal estacionario vemodificada su amplitud y la fase en base al valor de la función detransferencia del sistema en la frecuencia del armónico. La amplitudse modifica por un factor que se corresponde con el valor del módulode la función de transferencia en esa frecuencia, y la fase semodifica una cantidad correspondiente al valor de la fase de lafunción de transferencia también en esa frecuencia. La frecuenciadel armónico no se modifica por ser lineal. El principio desuperposición dice que el resultado de pasar una señal por unsistema lineal estacionario equivale a la suma de las modificacionesque se produzcan en cada uno de sus armónicos individualmente (esdecir, equivale a pasar cada uno de los armónicos por el sistema yluego sumarlos todos).

Este es mi código:


function y = transferencia(w, w0)
y = 1 ./ (1 + (%i * (w / w0)))
endfunction


t = -50:0.01:50
T = 2
sumatorio = 0
_{[FONT=Times New Roman, serif]w01[/FONT]}[FONT=Times New Roman, serif]= 2 * %pi / 5[/FONT]
[FONT=Times New Roman, serif]w[/FONT]_{[FONT=Times New Roman, serif]02[/FONT]}[FONT=Times New Roman, serif]= 2 * %pi[/FONT]
[FONT=Times New Roman, serif]w[/FONT]_{[FONT=Times New Roman, serif]03[/FONT]}[FONT=Times New Roman, serif]= 10 * %pi[/FONT]
for n = 1:1000
_wn = 2 * %pi* n / T
sumatorio = sumatorio +(10 * 1 * sin(w_n * 1 / 2) * cos(w_n * t) / (w_n* 1 / 2))
[FONT=Times New Roman, serif]y[/FONT]_{[FONT=Times New Roman, serif]1[/FONT]}[FONT=Times New Roman, serif]= transferencia(w[/FONT]_{[FONT=Times New Roman, serif]n[/FONT]}[FONT=Times New Roman, serif],w[/FONT]_{[FONT=Times New Roman, serif]01[/FONT]}[FONT=Times New Roman, serif])[/FONT]
[FONT=Times New Roman, serif]y[/FONT]_{[FONT=Times New Roman, serif]2[/FONT]}[FONT=Times New Roman, serif]= transferencia(w[/FONT]_{[FONT=Times New Roman, serif]n[/FONT]}[FONT=Times New Roman, serif],w[/FONT]_{[FONT=Times New Roman, serif]02[/FONT]}[FONT=Times New Roman, serif])[/FONT]
[FONT=Times New Roman, serif]y[/FONT]_{[FONT=Times New Roman, serif]3[/FONT]}[FONT=Times New Roman, serif]= transferencia(w[/FONT]_{[FONT=Times New Roman, serif]n[/FONT]}[FONT=Times New Roman, serif],w[/FONT]_{[FONT=Times New Roman, serif]03[/FONT]}[FONT=Times New Roman, serif])[/FONT]
; end
f = (10 * 1 / T) + (2 *sumatorio / T)

No sé está bien así. ¿Cómo podría representar los resultados con plot?