数字信号处理实验标准报告.docVIP

南昌航空大学数字信号实验 年 月 日 课程名称： 数字信号处理 实验名称： 系统及系统响应 班级： 姓名： 同组人： 指导老师评定： 签名： 一、实验目的 (1) 观察离散系统的频率响应。 (2) 熟悉连续信号经理想采样前后的频谱变化关系，加深对时域采样定理的理解。 (3) 利用序列的FT对连续信号、离散信号及系统响应进行频域分析。 (4) 利用卷积方法观察分析系统的时域特性。 二、实验内容 (1) 给定一因果系统，求出并绘制H(z)的幅频响应与相频响应。 (2) 对信号进行采样可得到采样序列 其中A=444.128，，， 实现下列实验内容： a.取采样频率＝1kHz，即T＝1ms观察所得采样xa(n)的幅频特性和图中的||在折叠频率附近有无明显差别。 b.改变采样频率＝300Hz，观察的变化，并做记录；进一步降低采样频率，＝200Hz，观察频谱混叠是否明显存在，说明原因，并记录这时的曲线。 （4）求通过系统的响应y(n)。 三、程序清单 1、 b=[1,sqrt(2),1]; %系统函数的分母系数 a=[1,-0.67,0.9]; %系统函数的分子系数 [h,w]=freqz(b,a); %函数freqz基于FFT算法计算数字滤波器的Z变换频率响应 am=abs(h); %得到幅频响应 subplot(2,1,1);plot(w,am);xlabel(w); ph=angle(h); %得到相频响应 subplot(2,1,2);plot(w,ph) 图2-1 系统函数的幅频响应和相频响应 2、 fs=10000; t=0:1/fs:0.1; %采集信号长度为0.1秒,时间间隔为1/fs A=444.128;a=50*sqrt(2)*pi;b=a; xa=A*exp(-a*t).*sin(b*t); %近似模拟信号 k=0:511; w=2*pi*k/512; Xa=dtft(xa,w); %近似模拟信号频谱 figure(2);subplot(3,2,1); plot(t,xa); %画出原始波形 axis([0,max(t),min(xa),max(xa)]);title(模拟信号); xlabel(t(s));ylabel(xa(t)); line([0,max(t)],[0,0]); f=fs*k/512; %模拟频率f subplot(3,2,2); plot(f,abs(Xa)); title(模拟信号的幅度频谱); axis([0,500,0,max(abs(Xa))]); xlabel(f(Hz));ylabel(|Xa(jf)|); fs1=1000;t1=0:1/fs1:0.1; %采集信号长度为0.1秒 x1=A*exp(-a.*t1).*sin(b*t1); %1kHZ采样序列x1(n) X1=dtft(x1,w); %x1(n)的512点dtft fs2=300;t2=0:1/fs2:0.1; %采集信号长度为0.1秒 x2=A*exp(-a.*t2).*sin(b*t2);%1kHZ采样序列x2(n) X2=dtft(x2,w); %x1(n)的512点dtft subplot(3,2,3);stem(t1,x1,.); axis([0,max(t1),min(x1),max(x1)]); title(采样序列x1(n)); xlabel(n);ylabel(x1(n));line([0,max(t1)],[0,0]); subplot(3,2,4);plot(w/pi,abs(X1)); title(x1(n)的幅度频谱); axis([0,2,0,max(abs(X1))]); xlabel(w/pi);ylabel(|X1(jw)|); subplot(3,2,5);stem(t2,x2,.); axis([0,max(t2),min(x2),max(x2)]); title(采样序列x2(n)); xlabel(n);ylabel(x2(n)); line([0,max(t2)],[0,0]); subplot(3,2,6);plot(w/pi,abs(X2)); title(x2(n)的幅度频谱); axis([0,2,0,max(abs(X2))]); xlabel(w/pi);ylabel(|X2(jw)|); function X=dtft(x,w) X=x*exp(-j*[