1.m

来自「由于待提取的分量为直流分量」· M 代码 · 共 55 行

close all;
clear all;
clc
while(1)
    clc
N=input('请输入采样点数');
n=1:N;
r=input('请输入躁声均方差');
noise=r*randn(1,N);
signal=5;
temp=input('   1——IIR 滤波     2——FIR 滤波   其他—退出')
if(temp==1) 
    clc
    disp('IIR 滤波器');
a=input('请输入a(0--1)的值');
b=1-a;
bb=[b];
aa=[1 -a];
elseif(temp==2)
  clc
    disp('IIR 滤波器');
    h0=input('请输入 h0=');
    h1=input('请输入 h1=');
    h2=input('请输入 h2=');
    h3=input('请输入 h3=');
    
    bb=[ h0 h1 h2  h3];
    aa=[1];
else
    break;
end
xn=noise+signal;
yn=filter(bb,aa,xn);
subplot(311);
plot(n,xn);
title('   x(n)');
xlabel('n');
ylabel('幅度');
subplot(312);
plot(n,yn);
title('   y(n)');
xlabel('n');
ylabel('幅度');
subplot(313);
[h,w]=freqz(bb,aa);
plot(w,abs(h));
title('滤波器幅频响应');
xlabel('W/pi')
ylabel('幅度');
temp=input('是否要继续:  y|Y----是    n|N  ------否 :','s');
if(temp==['n']|temp==['N'])  
    break;
end
end