7、实验程序清单 8、%实验1程序exp1.m 9、% IIR数字滤波器设计及软件实现 10、clear all;close all 11、Fs=10000;T=1/Fs; %采样频率 12、%调用信号产生函数mstg产生由三路抑制载波调幅信号相加构成的复合信号st 13、st=mstg; 14、%低通滤波器设计与实现========================================= 15、fp=280;fs=450; 16、wp=2*fp/Fs;ws=2*fs/Fs;rp=0.1;rs=60; ?指标（低通滤波器的通、阻带边界频） 17、[N,wp]=ellipord(wp,ws,rp,rs); %调用ellipord计算椭圆DF阶数N和通带截止频率wp 18、[B,A]=ellip(N,rp,rs,wp); %调用ellip计算椭圆带通DF系统函数系数向量B和A 19、y1t=filter(B,A,st); %滤波器软件实现 20、% 低通滤波器设计与实现绘图部分 21、figure(5); 22、subplot(2,1,1); 23、myplot(B,A); %调用绘图函数myplot绘制损耗函数曲线 24、yt='y_1(t)'; 25、subplot(2,1,2); 26、tplot(y1t,T,yt); %调用绘图函数tplot绘制滤波器输出波形 27、%带通滤波器设计与实现==================================================== 28、fpl=440;fpu=560;fsl=275;fsu=900; 29、wp=[2*fpl/Fs,2*fpu/Fs];ws=[2*fsl/Fs,2*fsu/Fs];rp=0.1;rs=60; 30、[N,wp]=ellipord(wp,ws,rp,rs); %调用ellipord计算椭圆DF阶数N和通带截止频率wp 31、[B,A]=ellip(N,rp,rs,wp); %调用ellip计算椭圆带通DF系统函数系数向量B和A 32、y2t=filter(B,A,st); %滤波器软件实现 33、% 带通滤波器设计与实现绘图部分 34、figure(3); 35、subplot(2,1,1); 36、myplot(B,A); %调用绘图函数myplot绘制损耗函数曲线 37、yt='y_2(t)'; 38、subplot(2,1,2); 39、tplot(y2t,T,yt); %调用绘图函数tplot绘制滤波器输出波形 40、%高通滤波器设计与实现================================================ 41、fp=890;fs=600; 42、wp=2*fp/Fs;ws=2*fs/Fs;rp=0.1;rs=60; ?指标（低通滤波器的通、阻带边界频） 43、[N,wp]=ellipord(wp,ws,rp,rs); %调用ellipord计算椭圆DF阶数N和通带截止频率wp 44、[B,A]=ellip(N,rp,rs,wp,'high'); %调用ellip计算椭圆带通DF系统函数系数向量B和A 45、y3t=filter(B,A,st); %滤波器软件实现 46、% 高低通滤波器设计与实现绘图部分 47、figure(4); 48、subplot(2,1,1); 49、myplot(B,A); %调用绘图函数myplot绘制损耗函数曲线 50、yt='y_3(t)'; 51、subplot(2,1,2); 52、tplot(y3t,T,yt); %调用绘图函数tplot绘制滤波器输出波形实验程序运行结果 实验4程序exp4.m运行结果如图104.2所示。由图可见，三个分离滤波器指标参数选取正确，算耗函数曲线达到所给指标。分离出的三路信号y1(n)，y2(n)和y3(n)的波形是抑制载波的单频调幅波。 (a) 低通滤波器损耗函数及其分离出的调幅信号y1(t) (b) 带通滤波器损耗函数及其分离出的调幅信号y2(t) (c)高通滤波器损耗函数及其分离出的调幅信号y3(t) 图104. 实验4程序exp4.m运行结果 53、思考题 （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)。