El programa de procesamiento de filtrado y análisis del espectro de señales se diseñó utilizando MATLAB.
%Escribir título
%Diseñar filtro de paso bajo:
[N, Wc]= pero ord()
Se estima que el orden mínimo n y la frecuencia de corte de 3 dB Wc del filtro de paso bajo Butterworth son %.
[a, b] = butter(N, Wc); % Diseño del filtro de paso bajo Butterworth
[h, f] = freqz() % Encuentra el número Respuesta de frecuencia del filtro de paso bajo
Figura (2); %Abrir ventana 2
Gráfico de rama (221) %Ventana dividida de visualización gráfica
plot( f, ABS(h));% Dibuje el diagrama de respuesta amplitud-frecuencia del filtro de paso bajo de Butterworth.
Título(Filtro de paso bajo bárbaro ' ');
Cuadrícula; % dibujar imagen con cuadrícula
sf = filter(a, b , s); nueva función de la función de superposición s después del filtro de paso bajo
Gráfico de rama (222);
plot(t, SF % dibuja la función de superposición El diagrama de dominio de tiempo de S después); pasando el filtro de paso bajo.
Xlabel('segundo');
Ylabel("tiempo de amplitud");
SF=fft(sf, 256 % en la función de superposición S); Después de pasar por el filtro de paso bajo, se realiza una transformada rápida de Fourier de base 2 de 256 puntos en la nueva función.
W=% nueva frecuencia angular de señal
Gráfico de rama (223);
plot() % traza la función superpuesta S después de pasar por el bajo; -pasar espectro de filtro.
Title('Espectrograma filtrado de paso bajo');
%Diseño de filtro de paso alto
[N, Wc]= pero ord()
Se estima que el orden más bajo n y la frecuencia de corte de 3 dB Wc del filtro de paso alto de Butterworth son %.
[a, b]=butter(N, Wc, 'high'); %Diseño de filtro de paso alto Butterworth
[h, f]= freqz(); la respuesta de frecuencia del filtro de paso alto digital
Figura (3);
Gráfico de rama (221);
plot()); Gráfico de respuesta amplitud-frecuencia del filtro de paso alto de Watts.
Title('Filtro de paso alto bárbaro');
Grid; % Dibuja una imagen con una cuadrícula
SF = filter(); % Función de superposición nueva función de s después de pasar el filtro de paso alto
Gráfico de rama (222);
plot(t, SF); % Dibuja el tiempo de la función de superposición S después de pasar el diagrama de dominio del filtro de paso alto.
xlabel('Tiempo (segundos)');
ylabel("Forma de onda de tiempo");
w % frecuencia angular de la nueva señal
Gráfico de rama (223);
plot()); % traza el espectro de la función superpuesta S después de pasar por el filtro de paso alto.
Title('Espectrograma filtrado de paso alto');
% Diseño de filtro de paso de banda
[N, Wc]= botón([) p>
El orden mínimo n y la frecuencia de corte de 3 dB Wc del filtro de paso de banda Butterworth se estiman en %.
[a, b] = butter(N, Wc); % Diseño del filtro de paso de banda Butterworth
[h, f] = freqz() % Encuentra el número de respuesta de frecuencia del filtro de paso de banda
Figura (4);
Gráfico de rama (221);
plot(f, ABS(h));% Dibuja la amplitud-frecuencia Gráfico de respuesta del filtro de paso de banda Butterworth.
Title('Filtro de paso de banda de mantequilla');
Cuadrícula; % dibujar imagen con cuadrícula
sf = filter(a, b, s); función de la función de superposición s después de pasar el filtro de paso de banda
Gráfico de rama (222);
plot(t, SF % dibuja la función de superposición S Diagrama de dominio de tiempo después del filtro de paso de banda); .
xlabel('Tiempo (segundos)');
ylabel("Forma de onda de tiempo");
SF = FFT() % en la función de superposición S. Después de pasar por el filtro de paso de banda, la nueva función se somete a una transformada rápida de Fourier de base 2 de 256 puntos.
W=(%Nueva frecuencia angular de señal
Gráfico de rama (223);
plot(') % dibuja la función de superposición S a través de banda- pasar el espectro de filtrado después del filtro.
Título('Espectrograma filtrado de paso de banda');