数字信号处理(IIR数字滤波器设计及软件实现)实验报告

计算机科学与工程学院

《数字信号处理》实验报告[4]

《数字信号处理》实验报告

2

计算机科学与工程学院 含载频成分。 (3)编程序调用 MATLAB 滤波器设计函数 ellipord 和 ellip 分别设计这三个椭圆滤波 器，并绘图显示其幅频响应特性曲线。 (4)调用滤波器实现函数 filter,用三个滤波器分别对信号产生函数 mstg 产生的信 号 st 进行滤波，分离出 st 中的三路不同载波频率的调幅信号 y1(n)、y2(n)和 y3(n), 并绘 图显示 y1(n)、y2(n)和 y3(n)的时域波形，观察分离效果。

三、信号产生函数 mstg 清单function st=mstg %产生信号序列向量 st,并显示 st 的时域波形和频谱 %st=mstg 返回三路调幅信号相加形成的混合信号，长度 N=1600 N=1600 %N 为信号 st 的长度。 Fs=10000;T=1/Fs;Tp=N*T; %采样频率 Fs=10kHz,Tp 为采样时间 t=0:T:(N-1)*T;k=0:N-1;f=k/Tp; fc1=Fs/10; %第 1 路调幅信号的载波频率 fc1=1000Hz, fm1=fc1/10; %第 1 路调幅信号的调制信号频率 fm1=100Hz fc2=Fs/20; %第 2 路调幅信号的载波频率 fc2=500Hz fm2=fc2/10; %第 2 路调幅信号的调制信号频率 fm2=50Hz fc3=Fs/40; %第 3 路调幅信号的载波频率 fc3=250Hz, fm3=fc3/10; %第 3 路调幅信号的调制信号频率 fm3=25Hz xt1=cos(2*pi*fm1*t).*cos(2*pi*fc1*t); %产生第 1 路调幅信号 xt2=cos(2*pi*fm2*t).*cos(2*pi*fc2*t); %产生第 2 路调幅信号 xt3=cos(2*pi*fm3*t).*cos(2*pi*fc3*t); %产生第 3 路调幅信号 st=xt1+xt2+xt3; %三路调幅信号相加 fxt=fft(st,N); %计算信号 st 的频谱 %====以下为绘图部分，绘制 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('幅度')

四、实验程序框图如图

调用函数 mstg 产生 st,自动绘图 显示 st 的时域波形和幅频特性曲 线

调用 ellipord 和 ellip 分别设计三个椭圆滤 波器，并绘图显示其幅频响应特性曲线。 《数字信号处理》实验报告

《数字信号处理》实验报告 4

《数字信号处理》实验报告 5

《数字信号处理》实验报告

6

《数字信号处理》实验报告

7

《数字信号处理》实验报告 8

《数字信号处理》实验报告

9

《数字信号处理》实验报告 10