实验8IIR数字滤波器设计综合实验(预习报告).doc

实验8 IIR数字滤波器设计综合实验 1．实验程序及运行结果 信号产生函数：mstg.m 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(2,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(2,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）该函数生成的三路时域信号分别为 （2）合成信号为 （3）采样频率为，采样点数为点，采样时间。 （4）抑制载波单频调幅信号的数学表示式为 源程序：shzxhchlshiyan4_1 %shzxhchlshiyan4_1 % IIR数字滤波器设计及软件实现 clear all,close all,clc,clf; Fs=10000;T=1/Fs; %采样频率 %调用信号产生函数mstg产生由三路抑制载波调幅信号相加构成的复合信号st st=mstg; %低通滤波器设计与实现=========================================== fp=280;fs=450; wp=2*fp/Fs;ws=2*fs/Fs;rp=0.1;rs=60; %椭圆低通滤波器技术指标（低通滤波器的通、阻带边界频） [N,wp]=ellipord(wp,ws,rp,rs); %调用ellipord计算椭圆低通滤波器阶数N和通带截止频率wp [B,A]=ellip(N,rp,rs,wp); %调用ellip计算椭圆低通滤波器系统函数系数向量B和A y1t=filter(B,A,st); %求低通滤波器的零状态响应 [H,w]=freqz(B,A); %求低通滤波器的频响特性 % 绘制椭圆低通滤波器衰减特性曲线 figure(2); mag=abs(H); db=20*log10((mag+eps)/max(mag)); %eps浮点数的相对误差 subplot(2,1,1);plot(w/pi,db);%绘制损耗函数曲线 axis([0,1,-100,5]);grid on; title(低通滤波损耗函数曲线);xlabel(\omega/\pi );ylabel(幅度(dB)); set(gca,Xtick,[0,0.1,0.2,0.3,0.4,0.5,0.6,0.7,0.8,0.9,1]); set(gca,Ytick,[-100,-80,-60,-40,-20,0]); % 绘制225-275Hz范围内时域信号波形 N=1600;t=0:T:(N-1)*T; subplot(2,1,2);plot(t,y1t); title(低通滤波分离225-275Hz范围内时域信号波形); xlabel(t/s);ylabel(y_1(t));axis([0,0.16,-1.1,1.1]); %带通滤波