Traducción del programa Matlab
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 xYlabel('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)');