数字信号处理实验报告
波器输出波形 %
带
通
滤
波
器
设
计
与
实
现
==================================================== 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); %滤波器软件实现 % 带通滤波器设计与实现绘图部分（省略） %
高
通
滤
波
器
设
计
与
实
现
================================================ fp=890;fs=600;
wp=2*fp/Fs;ws=2*fs/Fs;rp=0.1;rs=60; ?指标（低通滤波器的通、阻带边界频）
[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); %滤波器软件实现
- 26 -
数字信号处理实验报告
% 高低通滤波器设计与实现绘图部分（省略）
七、 实验程序运行结果
实验4程序exp4.m运行结果如图104.2所示。由图可见，三个分离滤波器指标参数选取正确，算耗函数曲线达到所给指标。分离出的三路信号y1(n)，y2(n)和y3(n)的波形是抑制载波的单频调幅波。
(a) 低通滤波器损耗函数及其分离出的调幅信号y1(t)
(b) 带通滤波器损耗函数及其分离出的调幅信号y2(t)
- 27 -
数字信号处理实验报告
(c)高通滤波器损耗函数及其分离出的调幅信号y3(t)
图104. 实验4程序exp4.m运行结果
八、思考题
（1）请阅读信号产生函数mstg，确定三路调幅信号的载波频率和调制信号频率。
（2）信号产生函数mstg中采样点数N=800，对st进行N点FFT可以得到6根理想谱线。如果取N=1000，可否得到6根理想谱线？为什么？N=2000呢？请改变函数mstg中采样点数N的值，观察频谱图验证您的判断是否正确。
（3）修改信号产生函数mstg，给每路调幅信号加入载波成分，产生调幅（AM）信号，重复本实验，观察AM信号与抑制载波调幅信号的时域波形及其频谱的差别。 提示：AM信号表示式：s(t)?[1?cos(2?f0t)]cos(2?fct)。
答：（1）三路调幅信号的载波频率分别为250Hz、500Hz、1000Hz。调制信号分别为50Hz、100Hz、200Hz。
（2） 因为信号st是周期序列，谱分析时要求观察时间为整
- 28 -
数字信号处理实验报告
数倍周期。所以，本题的一般解答方法是，先确定信号st的周期，在判断所给采样点数N对应的观察时间Tp=NT是否为st的整数个周期。但信号产生函数mstg产生的信号st共有6个频率成分，求其周期比较麻烦，故采用下面的方法解答。
分析发现，st的每个频率成分都是25Hz的整数倍。采样频率Fs=10kHz=25×400Hz，即在25Hz的正弦波的1个周期中采样400点。所以，当N为400的整数倍时一定为st的整数个周期。因此，采样点数N=800和N=2000时，对st进行N点FFT可以得到6根理想谱线。如果取N=1000，不是400的整数倍，不能得到6根理想谱线。
九、实验心得
通过本次试验，我对数字滤波器有了更深入的了解，也对matlab的应用更加熟练，这次试验提高了我的工程素质。
- 29 -
数字信号处理实验报告
实验四、FIR数字滤波器设计与软件实现
一、实验目的
（1）掌握用窗函数法设计FIR数字滤波器的原理和方法。 （2）掌握用等波纹最佳逼近法设计FIR数字滤波器的原理和方法。
（3）掌握FIR滤波器的快速卷积实现原理。 （4）学会调用MATLAB函数设计与实现FIR滤波器。
二、实验内容及步骤
（1）、认真复习第七章中用窗函数法和等波纹最佳逼近法设计FIR数字滤波器的原理；
（2）、调用信号产生函数xtg产生具有加性噪声的信号xt，并自动显示xt及其频谱，如图10.5.1所示；
图10.5.1 具有加性噪声的信号x(t)及其频谱如图
（3）、请设计低通滤波器，从高频噪声中提取xt中的单频调幅信号，要求信号幅频失真小于0.1dB，将噪声频谱衰减60dB。先观察xt的频谱，确定滤波器指标参数。
- 30 -