实验七20110599.docVIP

数字信号处理实验 姓名：刘书梅 学号班级：11信工2班 实验名称:用窗函数法设计数字滤波器 实验目的 掌握用窗函数法设计FIR数字滤波器的原理及具体方法； 深入理解吉布斯现象，理解不同窗函数的特点。 二、 习题1、代码与图形如下 clc;clear;clf; clf;N=1000;Fs=1000; %数据总数和采样频率 fc=200; n=[0:N-1];t=n/Fs; %时间序列 f1=50;f2=250; x=sin(2*pi*f1*t)+sin(2*pi*f2*t); %输入信号 b=fir1(40,fc*2/Fs); %设计40阶的低通滤波器，归一化截止频率 yfft=fftfilt(b,x,256); %对数据进行滤波 n1=81:241;t1=t(n1); %选择采样点间隔 x1=x(n1); %与采样点对应的输入信号 subplot(2,1,1); plot(t1,x1); grid on; %绘制输入信号 title(输入信号); n2=n1-40/2;t2=t(n2); %输出信号，扣除了相位延迟N/2 y2=yfft(n2); subplot(2,1,2); plot(t2,y2); %绘制输出信号 title(输出信号); grid on; xlabel(时间/s) 说明：此题首先通过fir1(40,fc*2/Fs)命令设计40阶的低通滤波器，再用此滤波器对输入信号x=sin(2*pi*f1*t)+sin(2*pi*f2*t)进行滤波。用subplot命令把图形分割成两部分，并分别通过plot命令绘出原输入信号的波形和经滤波后的信号的波形。 三、 习题2、代码与图形如下 clc;clear;clf; clf;N=300;Fs=1000; %数据总数和采样频率 Order=100; %滤波器阶数 n=[0:N-1];t=n/Fs; %时间序列 wn=[100 200]/(Fs/2); %边界频率转换为归一化频率，据6-19式 b=fir1(Order,wn,stop); % 设计 100 阶的带阻滤波器 [H,f]=freqz(b,1,512,Fs); figure(1) subplot(2,1,1) plot(f,20*log10(abs(H))) xlabel(频率/Hz);ylabel(振幅/dB);grid on; subplot(2,1,2),plot(f,180/pi*unwrap(angle(H))) xlabel(频率/Hz);ylabel(相位/^o);grid on; f1=50;f2=150; %输入信号频率 x=sin(2*pi*f1*t)+sin(2*pi*f2*t); %输入信号 y=fftfilt(b,x,256);% 对数据进行滤波 figure(2) subplot(2,1,1) plot(t,x); grid on; %绘制输入信号 title(输入信号); subplot(2,1,2); plot(t,y); %绘制输出信号 hold on; plot(Order/2/Fs*ones(1,2),ylim, r) %绘制延迟到的时刻 title(输出信号); xlabel(时间/s) grid on; 说明:此程序是来设计阶数为100，采样频率为1000的数字滤波器。fir1(Order,wn,stop)命令是用来设计带阻滤波器的，并用[H,f]=freqz(b,1,512,Fs)求出其频率响应。再把第一个图形分割为两部分，分别在两部分绘出滤波器的幅频响应和相频响应的图形。在第二个图形分割出的两部分分别绘出输入信号x=sin(2*pi*f1*t)+sin(2*pi*f2*t)的波形及经滤波器后的的图形，由图可见滤波后的输出信号平滑了很多，开始的一段不规则图形是由于延时造成的。 四、 习题3 用fdatool命令在MATLAB中设计一个IIR数字滤波器。滤波器的图形及解说如下 图1 图2 图3 图4 说明：图形1为阶数是30，截止频率为0.3的无限长切比雪夫I型数字低通滤波器的幅频特性图。图2为设计的滤波器的相频特性图，