Acerca del uso de MATLAB para diseñar un programa para análisis de espectro y filtrado de señales

Programa completo

Escribe un título

Diseña un filtro de paso bajo:

[N, Wc]=buttord()

Estime el orden mínimo N y la frecuencia de corte de 3 dB Wc del filtro de paso bajo Butterworth

[a, b]=butter(N, Wc);

[h, f]=freqz(); Encuentra la respuesta de frecuencia del filtro de paso bajo digital

figura(2) Abrir ventana

subplot(221); ventana dividida de visualización gráfica

plot(f, abs(h)); dibuja el diagrama de respuesta de amplitud-frecuencia del filtro de paso bajo de Butterworth

title(Butterworth low) -pass filter' ');

grid; Dibuja una imagen con una cuadrícula

sf=filter(a, b, s); filtro de paso bajo

sf=filter(a, b, s); p>

subplot(222);

plot(t, sf); gráfica de la función de superposición S después de pasar el filtro de paso bajo

xlabel('time( segundos)');

ylabel('Time by amplitude');

SF=fft(sf, 256); La nueva función después de que la función de superposición S pasa a través de un filtro de paso bajo Realiza una transformada rápida de Fourier de base 2 de 256 puntos

w= nueva señal frecuencia angular

subplot(223);

plot()); dibuja la función de superposición S Espectrograma después del filtro de paso bajo

title('Espectro después del bajo -pass filtering');

Diseñar filtro de paso alto

[N , Wc]=buttord()

Estimar el orden mínimo N y la frecuencia de corte de 3dB Wc del filtro de paso alto Butterworth

[a, b]=butter(N, Wc,' high'); Diseño del filtro de paso alto Butterworth

[h, f] =freqz(); Encuentra la respuesta de frecuencia del filtro de paso alto digital

figure(3);

subplot(221);

plot()) Dibuje el diagrama de respuesta de amplitud-frecuencia del filtro de paso alto de Butterworth

title('filtro de paso alto de Butterworth');

dibuje una imagen con una cuadrícula;

sf=filter(); superpone la nueva función S después de pasar el filtro de paso alto

subplot(222);

plot(t, sf) ; Dibujar el gráfico en el dominio del tiempo de la función de superposición S después de pasar el filtro de paso alto

xlabel('Tiempo(segundos)');

ylabel('Forma de onda de tiempo') ;

w; Nueva frecuencia angular de señal

subplot(223);

plot()

); Dibujar el espectro de la función de superposición S después de pasar el filtro de paso alto

title('Espectro después del filtrado de paso alto');

Diseñar el filtro de paso de banda<); /p>

[N, Wc]=buttord([)

Estimar el orden mínimo N y la frecuencia de corte de 3dB Wc del filtro de paso de banda Butterworth

[a, b] =butter (N, Wc); Diseño del filtro de paso de banda Butterworth

[h, f]=freqz(); Encuentre la respuesta de frecuencia del filtro de paso de banda digital

figura( 4);

subplot(221);

plot(f, abs(h)); Dibuje el diagrama de respuesta amplitud-frecuencia del filtro de paso de banda de Butterworth.

title(' butter bandpass filter');

grid; Dibuja una imagen con una cuadrícula

sf=filter(a, b, s); S después de pasar por el filtro de paso de banda Function

subplot(222);

plot(t, sf); Dibuje el gráfico en el dominio del tiempo de la función superpuesta S después de pasar por el filtro de paso de banda.

xlabel( 'Tiempo(segundos)');

ylabel('Forma de onda de tiempo');

SF=fft(); Realizar la nueva función después la función de superposición S pasa la transformada rápida de Fourier de 256 puntos base-2 del filtro de paso de banda

w=(nueva frecuencia angular de señal

subplot(223);

plot(')); superposición del gráfico Espectrograma de la función S después de pasar por el filtro de paso de banda

title('Espectro después del filtrado de paso de banda');