2.m

来自「滤波器的设计代码」· M 代码 · 共 39 行

fs=2000;   %采样频率
N=200;    % 采样点
F=fs/N;    %频率分辨率
n=0:199;   %采样点序列
Tp=1/F;    %采样时间
 t=n/200*Tp;


w=linspace(0,2*pi,250);

X=sin(20*2*pi*t)+sin(220*2*pi*t)+sin(750*2*pi*t) ;   %对原时域的信号进行离散化，采样200点
plot(w,X);
XX=fft(X,250);
plot(w,XX);

wp=[60 350];
ws=[40 400];
rp=3;
rs=30;
wp1=wp/fs*2;
ws1=ws/fs*2;
[N1,Wn]=buttord(wp1,ws1,rp,rs); 
[b,a]=butter(N1,Wn);
% freqz(b,a,512,fs)
% Hs=20*log10(freqz(b,a,512,fs));
% plot(abs(Hs))

HH=freqz(b,a,250,fs/2)
YY=XX.*HH';
% plot(abs(YY))
y=ifft(YY,250)


plot(w,y)
shg