数字低通滤波器的指标.doc

《数字低通滤波器的指标.doc》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。

1、数字低通滤波器的指标如下:信号的取样频率是1000Hz，截止频率是1000Hz，截止频率是200Hz，滤波器的阶数是81，用矩形窗、汗宁窗、凯赛窗、三角窗分别设计低通滤波器。程序如下：wpass=0.4*pi;%通带截止频率为0.4πw1=boxcar(81);%矩形窗时域函数w2=hamming(81);%汉宁窗时域函数w3=kaiser(81,7.865);%凯赛窗时域函数w4=triang(81);%三角窗时域函数n=1:1:81;hd=sin(wpass*(n-41))./(pi*(n-41));hd(41)=wpass/pi;h

2、1=hd.*rot90(w1);h2=hd.*rot90(w2);h3=hd.*rot90(w3);h4=hd.*rot90(w4);[a1,rad]=freqz(h1);%频率响应函数求其幅度和频率[a2,rad]=freqz(h2);[a3,rad]=freqz(h3);[a4,rad]=freqz(h4);subplot(2,2,1);plot(rad,20*log10(abs(a1)));xlabel('归一化频率');ylabel('幅度(dB)');title('利用矩形窗设计的滤波器');gridon;subplot(2,2

3、,2);plot(rad,20*log10(abs(a2)));xlabel('归一化频率');ylabel('幅度(dB)');title('利用汗宁窗设计的滤波器');gridon;subplot(2,2,3);plot(rad,20*log10(abs(a3)));xlabel('归一化频率');ylabel('幅度(dB)');title('利用凯赛窗设计的滤波器');gridon;subplot(2,2,4);plot(rad,20*log10(abs(a4)));xlabel('归一化频率');ylabel('幅度(dB)')

4、;title('利用三角窗设计的滤波器');gridon;FFT(三个)书上P140第一个FFT已知信号由频率15Hz、幅值为0.5的正xuan信号和频率为40Hz、幅值为2的正XUAN信号组成，数据取样频率为FS=100Hz，试分别绘制N=128点DFT的幅频图和N=1024点DFT幅频图由题可知信号可写为x(t)=0.5sin(2πf1t)+2sin(2πf2t)%MatLABExampleclffs=100;N=128;n=0:N;t=n/fs;x=0.5*sin(2*pi*15*t)+2*sin(2*pi*40*t);y=fft(

5、x,N);mag=abs(y);f=(0:length(y)-1)'*fs/length(y);subplot(221)plot(f,mag);xlabel(Frequency(Hz)');ylabel(Magnitude');title('N=128')reidsubplot(222)plot(f(1:N/2),mag(1:N/2));xlabel('Frequency(Hz)');ylabel('Magnitude');title('N=128')gridfs=100;N=1024n=0:N-1;t=n/fs;x=0.5*sin(2*

6、pi*15*t)+2*sin(2*pi*40*t);y=fft(x,N);mag=abs(y);f=(0:length(y)-1)'*fs/length(y);subplot(223)plot(f,mag);xlabel('Frequency(Hz)');ylabel('Magnitude');title('N=1024')grodsubplot(224)plt(f(1:N/2),mag(1:N/2));xlabel('Frequency(Hz)');ylabel('Magnitude');title('N=1024')grid要改