Traducción del programa Matlab

k = 1; asignar un valor a k

b =[0 25 40]* k; multiplicar cada elemento del array por k

a=conv([1 9+15j ] , [1 9-15j]); Encuentra la convolución de dos matrices

sys=tf(b, a); Encuentra la función de transferencia del dominio de transformación z. El coeficiente del numerador es B y. el coeficiente del denominador es a.

t = 0:0.1:10; genera una serie de tiempo, de 0 a 10, con un intervalo de 0,1.

h = impulso (sys, t); genera una respuesta de impulso en el dominio del tiempo con sys como función de transferencia.

Subtrama (221); dividida en cuatro subtramas, dibuja la primera subtrama.

Plot (h); dibujar h

Grid

xlabel(' t '); nombre de la coordenada del eje x

ylabel( ' h(t)'); nombre de la coordenada del eje y

Título (' h(t)'); , 100); Encuentre el espectro

h 1 = ABS(h); Obtenga la amplitud del espectro, es decir, las características de amplitud-frecuencia.

h2 =Ángulo (h); toma el ángulo de fase del espectro, es decir, las características fase-frecuencia.

Gráfico de rama (222); dibuja la segunda imagen

Gráfico (w, h 1 gráfico de características de velocidad de fotogramas

Cuadrícula

< p); >Xlabel('Frecuencia angular (w)'); nombre de las coordenadas del eje x

Ylabel('amplitud'); nombre de las coordenadas del eje y

Título ('h(jw) ' características de amplitud-frecuencia); nombre del gráfico

Gráfico de rama (223) tercera imagen

gráfico (w, H2 * 180/pi);

Cuadrícula

Xlabel('Frecuencia angular (w)');

Ylabel('Fase (grados)');

Título (Características de frecuencia de fase de 'h(jw)');