MATLAB desenho filtro digital (IIR e FIR) (contendo o código Matlab)

MATLAB desenho filtro digital (IIR e FIR) (contendo o código Matlab)

As vantagens e desvantagens dos filtros de FIR
vantagens: boa linearidade fase, rápida velocidade de processamento, a estabilidade do circuito fechado de realimentação não é mais forte do que o IIR.
Desvantagens: FIR amplitude de frequência IIR característica baixo de precisão em comparação com
os filtros IIR
vantagens: a filtragem de ordem melhor IIR ao mesmo
inconveniente: um alto nível de fase n linear, necessidades de filtragem bidireccionais para ser corrigida quando a correcção não é fácil de controlar.
A chave é que a função de janela de filtro FIR: Aqui estão
janela Hamming:
FIR - a janela de Hamming

1. A amplitude dos lóbulos laterais inferiores, em particular, o primeiro lóbulo lateral;
2. Sidelobe amplitude para diminuir mais rapidamente, a fim de facilitar o aumento de atenuação de banda atenuada;
3. A largura do lóbulo principal é mais estreita, de modo que a banda de transição de filtro é mais estreito.
   N = 40
   factor de vazamento de 0,04%, de lóbulo lateral relativa supressão -42.1dB, 0,0625 largura do lóbulo principal (a frequência normalizada)
   
   de N, maior o tempo de cálculo é, o atraso de grupo mais. Tomar N = 40 (de acordo com a diferentes situações, a ordem do filtro podem ter valores diferentes, aqui filtrada do sinal de áudio de um período específico, tendo, assim, 40)
   N maior, quanto maior o tempo de cálculo é, quanto maior for o atraso de grupo. Tomar N = 40

resultados filtrados mostram:

sinal frente azul filtrada; o vermelho é FIR sinal filtrado
com o filtro IIR elíptica filtragem do sinal de áudio:

filtrar resultado mostra:
canções directos IIR_ Eliptic (oval) filtro de
filtragem FIR (amarelo), IIR filtrada (roxo)
para diferentes sinais de áudio, a análise original espectro, só muda para os parâmetros de projeto do filtro, para implementar um filtro digital.
Veja o código abaixo:

[x,fs]=audioread('C:\Users\a\Desktop\大三上\精密仪器\讨论课题目\新闻(加噪声).wav');
n=length(x);
  y=fft(x,n);       %做FFT变换
  f=fs*(0:n/2-1)/n;
  subplot(2,1,1);
  plot(x);
 axis([0 1400000 -1.5 1.5]);
  title('原始新闻信号采样后的时域波形'); 
xlabel('时间轴')
ylabel('幅值A')
subplot(2,1,2); 
plot(f,abs(y(1:n/2))); 
axis([0 1000 0 5000]);
title('原始新闻信号采样后的频谱图'); 
xlabel('频率Hz'); 
ylabel('频率幅值');
% n=0:N-1;
% w=2*n*pi/N;
FIRsignal=IIR_BANDPASSCheby;
outsignal=filter(FIRsignal,x);
%  IIRsignal=IIR_BANDPASSEliptic;
%  outsignal2=filter(IIRsignal,outsignal);
% set(gca,'fontsize',15,'tickdir','out')
% axis([0 5.0E-4 -2 2])
% xlabel('Time');ylabel('Amplitude')
% legend('FIR低通数字滤波器后的信号')  
 m=fft(outsignal,n);
 f=48000*(0:n/2-1)/n;
 figure
 subplot(211);
  plot(outsignal);
   axis([0 1400000 -1.5 1.5]);
  title('新闻滤波信号时域波形');
xlabel('时间轴')
ylabel('幅值A')
 subplot(212);
 plot(f,abs(m(1:n/2)))
 axis([200 2700 0 5000]);
  title('新闻滤波信号频谱图'); 
xlabel('频率Hz'); 
ylabel('频率幅值');
    sound(1.5*outsignal,44100);
%    audiowrite('C:\Users\a\Desktop\精密仪器\讨论课题目\新闻(IIR滤波后).wav',outsignal2,44100);
%    sound(x,44100);

A primeira questão, congratulamo-nos com o intercâmbio de correção, não gosto não spray! ! ! !