数字信号处理实验10.4重点.pptx

IIR数字滤波器设计及软件实现 实验四 实验分为： 设计IIR数字滤波器一般采用间接法（脉冲响应不变法和双线性变换法），应用最广泛的是双线性变换法。基本设计过程是： ①先将给定的数字滤波器的指标转换成过渡模拟滤波器的指标； ②设计过渡模拟滤波器； ③将过渡模拟滤波器系统函数转换成数字滤波器的系统函数。MATLAB信号处理工具箱中的各种IIR数字滤波器设计函数都是采用双线性变换法。 ④在编译实验程序之前，建立本实验所需的三个库函数： mstg.m; myplot.m; tplot.m 。 第六章介绍的滤波器设计函数butter、cheby1 、cheby2 和ellip可以分别被调用来直接设计巴特沃斯、切比雪夫1、切比雪夫2和椭圆模拟和数字滤波器。 实验理论 st=mstg; %====以下为绘图部分，绘制st的时域波形和幅频特性曲线subplot(3,1,1) plot(t,st);grid;xlabel(t/s);ylabel(s(t)); axis([0,Tp/8,min(st),max(st)]);title((a) s(t)的波形) subplot(3,1,2) stem(f,abs(fxt)/max(abs(fxt)),.);grid;title((b) s(t)的频谱) axis([0,Fs/5,0,1.2]); xlabel(f/Hz);ylabel(幅度) 实验程序 实验结果图 实验图 1 fp=280;fs=450; wp=2*fp/Fs;ws=2*fs/Fs;rp=0.1;rs=60; %DF指标（低通滤波器的通、阻带边界频） [N,wp]=ellipord(wp,ws,rp,rs); %调用ellipord计算椭圆DF阶数N和通带截止频率wp [B,A]=ellip(N,rp,rs,wp); %调用ellip计算椭圆带通DF系统函数系数向量B和A 实验程序 y1t=filter(B,A,st);%滤波器软件实现 % 低通滤波器设计与实现绘图部分 figure(2);subplot(3,1,1); myplot(B,A); %调用绘图函数myplot绘制损耗函数曲线 yt=y_1(t); subplot(3,1,2);tplot(y1t,T,yt); %调用绘图函数tplot绘制滤波器输出波形 实验结果图 实验图 2 fpl=440;fpu=560;fsl=275;fsu=900; wp=[2*fpl/Fs,2*fpu/Fs];ws=[2*fsl/Fs,2*fsu/Fs];rp=0.1;rs=60; [N,wp]=ellipord(wp,ws,rp,rs); %调用ellipord计算椭圆DF阶数N和通带截止频率wp [B,A]=ellip(N,rp,rs,wp); %调用ellip计算椭圆带通DF系统函数系数向量B和A 实验程序 y2t=filter(B,A,st); %滤波器软件实现 figure(3);subplot(3,1,1); myplot(B,A); %调用绘图函数myplot绘制损耗函数曲线 yt=y_2(t); subplot(3,1,2);tplot(y2t,T,yt); %调用绘图函数tplot绘制滤波器输出波形 实验结果图 实验图 3 fp=890;fs=600; wp=2*fp/Fs;ws=2*fs/Fs;rp=0.1;rs=60; %DF指标（低通滤波器的通、阻带边界频） [N,wp]=ellipord(wp,ws,rp,rs); %调用ellipord计算椭圆DF阶数N和通带截止频率wp [B,A]=ellip(N,rp,rs,wp,high); %调用ellip计算椭圆带通DF系统函数系数向量B和A 实验程序 y3t=filter(B,A,st); %滤波器软件实现 figure(4);subplot(3,1,1); myplot(B,A); %调用绘图函数myplot绘制损耗函数曲线 yt=y_3(t); subplot(3,1,2);tplot(y3t,T,yt); %调用绘图函数tplot绘制滤波器输出波形 实验结果图 实验图 4 （1）请阅读信号产生函数mstg，确定三路调幅信号的载波频率和调制信号频率。 答：三路调幅信号的载波频率分别为250Hz、500Hz、1000Hz. 调制信号频率分别为100Hz、50Hz、25Hz （2）信号产生函数mstg中采样点数N=800，对st进行N点FFT可以得到6根理想谱线。如果取N=1000，可否得到6根理想谱线？为什么？N=2000呢？请改变函数mstg中采样点数N的值，观察频谱图验证您的判断是否正确。 答：分析发现，st的每个频率成分都是25Hz的整数倍。采样频率Fs=10kH