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

数字低通滤波器的指标如下:信号的取样频率是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; h1=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( 利用矩形窗设计的滤波器); grid on; subplot(2,2,2); plot(rad,20*log10(abs(a2))); xlabel(归一化频率); ylabel(幅度(dB)); title(利用汗宁窗设计的滤波器); grid on; subplot(2,2,3); plot(rad,20*log10(abs(a3))); xlabel(归一化频率); ylabel(幅度(dB)); title(利用凯赛窗设计的滤波器); grid on; subplot(2,2,4); plot(rad,20*log10(abs(a4))); xlabel(归一化频率); ylabel(幅度(dB)); title(利用三角窗设计的滤波器); grid on; 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) %MatLAB Example clf fs=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(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) reid subplot(222) plot(f(1:N/2),mag(1:N/2)); xlabel(Frequency(Hz)); ylabel(Magnitude); title(N=128) grid fs=100; N=1024 n=0:N-1; t=n/fs; x=0.5*sin(2*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) grod subplot(224) plt(f(1:N/2),mag(1:N/2)); xlabel(Frequency(Hz)); ylabel(Magnitude); title(N=1024) grid 要改