数字信号处理上机报告
实验三、信号的频域与Z域分析
班级： 学号： 姓名： 成绩：
1实验目的
（1）理解序列离散傅里叶变换（DTFT）的定义，熟悉序列DTFT的计算及其主要性质； （2）掌握Z变换的计算和主要性质，熟悉Z变换的收敛域及其与序列特性的关系，以及Z变换与DTFT的关系；
（3）掌握时域离散线性时不变系统的频域分析方法，深刻理解系统的频率响应。了解系统的稳态响应和暂态响应、相位延迟和群延迟等概念；
（4）掌握时域离散线性时不变系统的z域分析方法，深刻理解离散系统的系统函数及其零极点分布，熟悉零极点分布与系统的因果性和稳定性关系、零极点分布对系统频率特性的影响、差分方程的Z变换解法等；
2 实验内容
（1） 设计计算机程序，产生序列并计算序列的DTFT，绘制其幅频特性和相频特性曲
线；
（2） 根据系统的单位脉冲响应和差分方程，计算系统的频率响应，绘制系统频率响应
的幅频特性和相频特性曲线；
（3） 根据系统的单位脉冲响应和差分方程，计算系统的系统函数、零极点分布；改变
系统的零极点分布，观察系统频率响应的变化。
3实验步骤
（1）设计有限长序列Rn；计算序列的DTFT，绘制幅频特性和相频特性曲线
（2）改变系统的系统函数的零点分布，绘制系统改变前和改变后的频率响应的幅频特性和相频特性曲线
4 程序设计
x=[1,1,1,1];nx=[0:3];%x(n)=R(n)
w=linspace(-2.8*pi,2.8*pi,100000);%取100000个点
数字信号处理上机报告
X=x*exp(-j*nx'*w);%DTFT figure(1);
subplot(3,2,1),plot(w/pi,abs(X));xlabel('\\omega/\\pi');ylabel('|X(e^j^\\omega)|') subplot(3,2,2),plot(w/pi,angle(X));xlabel('\\omega/\\pi');ylabel('\\phi(\\omega)/\\pi') %差分方程求解 a=[1,-0.4];b=[1]; [H,w]=freqz(b,a,'whole');
subplot(3,2,3),plot(w/pi,abs(H));xlabel('\\omega/\\pi');ylabel('|X(e^j^\\omega)|') subplot(3,2,4),plot(w/pi,angle(H));xlabel('\\omega/\\pi');ylabel('\\phi(\\omega)/\\pi') %零极点分布
a=[1,-1.6,0.9425];%分母 b1=[1,-0.3];b2=[1,-0.8];%分子 [F,w]=freqz(b1,a,'whole'); figure(2);
subplot(2,2,1),zplane(b1,a);%零极点分布图
subplot(2,2,3),plot(w/pi,abs(F));xlabel('\\omega/\\pi');ylabel('|X(e^j^\\omega)|') subplot(2,2,4),plot(w/pi,angle(F));xlabel('\\omega/\\pi');ylabel('\\phi(\\omega)/\\pi') figure(3);%改变零极点分布，观察频率响应变化 [F,w]=freqz(b2,a,'whole'); subplot(2,2,1),zplane(b2,a);
subplot(2,2,3),plot(w/pi,abs(F));xlabel('\\omega/\\pi');ylabel('|X(e^j^\\omega)|') subplot(2,2,4),plot(w/pi,angle(F));xlabel('\\omega/\\pi');ylabel('\\phi(\\omega)/\\pi')
数字信号处理上机报告
5实验结果及分析
零点改变前系统频率响应：