信号与系统实验六..docxVIP

离散时间信号与系统频域分析实验报告信号与系统实验报告六 学院：计算机与信息工程学院专业：通信工程班级：2012级计算机与信息技术学院设计性实验报告 专业：通信工程 年级/班级：2012级 2013—2014学年第2学期课程名称信号与系统指导教师学号姓名实验地点计科楼216实验时间2014-5-21项目名称离散时间信号与系统的频域分析实验类型设计性一、实验目的1、掌握离散时间信号与系统的频域分析方法，从频域的角度对信号与系统的特性进行分析。2、掌握离散时间信号傅里叶变换与傅里叶逆变换的实现方法。3、掌握离散时间傅里叶变换的特点及应用4、掌握离散时间傅里叶变换的数值计算方法及绘制信号频谱的方法二.实验仪器或设备装有MATLAB软件设备的计算机一台三 总体设计（设计原理、设计方案及流程等）1. 离散时间系统的频率特性系统单位冲激响应h[n]的傅里叶变换H () 成为LTI 系统的频率响应。与连续时间LTI 系统类似，通过系统频率响应可以分析出系统频率特性。与系统单位冲激响应h[n]一样，系统的频率响应H () 反映了系统内在的固有特性，它取决于系统自身的结构及组成系统元件的参数，与外部激励无关，是描述系统特性的一个重要参数， H () 是频率的复函数可以表示为：H () =|H ()| 其中，|H ()| 随频率变化的规律称为幅频特性；?() 随频率变化的规律称为相频特性。2. 离散时间信号傅里叶变换的数值计算方法算法原理，由傅里叶变换原理可知：F()=序列f [n]的离散时间傅里叶变换F()是的连续函数。由于数据在matlab中以向量的形式存在,F()只能在一个给定的离散频率的集合中计算。然而，只有类似F()=形式的的有理函数，才能计算其离散时间傅里叶变换。3.涉及到的Matlab函数3.1 freqz函数：实现离散时间系统频率响应特性的求解调用格式:[H,w]=freqz(B,A,N)B 和A 分别为离散系统的系统函数分子、分母多项式的系数向量，返回量H 则包含了离散系统频响在0~??范围内N 个频率等分点的值（其中N 为正整数），w 则包含了范围内N 个频率等分点。调用默认的N 时，其值是512。由于F () 是的连续函数，需要尽可能大地选取N 的值，以使得产生的图形和真实离散傅里叶变换的图形尽可能一致。为更加方便快速地运算，应将N 的值选为2 的幂，如256 或512.3.2 real 函数：求复数的实部调用格式:real_f=real (f)；3.3 imag函数：求复数的虚部调用格式:imag_f=imag (f)；3.4 abs 函数：求复数的模调用格式:abs_f=abs (f)；3.5 angle 函数：求复数的相位调用格式:angle_f=angle (f)；3.6 fft函数：实现离散信号f[n]的傅里叶变换值F[k]调用格式：F = fft(f)，计算序列f[n]的离散傅里叶变换值F[k]，其中F[k]长度与f[n]长度相同F = fft(f,L)，计算序列f[n]的L点离散傅里叶变换值F[k],其中L不小于N。若L大于N，则需要在计算离散傅里叶变换之前，对f[n]尾部补足L-N 个零。3.7 ifft函数：实现离散信号的傅里叶逆变换调用格式：f = ifft(F)四 实验步骤（包括主要步骤、代码分析等）1 离散时间傅里叶变换(1)下面参考程序是如下序列在范围?4??≤??≤ 4??的离散时间傅里叶变换F()=%计算离散时间傅里叶变换的频率样本clear all;w=-4*pi:8*pi/511:4*pi;num=[2 1]; den=[1 -0.6];h=freqz(num,den,w);subplot(2,1,1)plot(w/pi,real(h)); grid;title(‘实部’)xlabel(‘omega/\pi’);ylabel(‘振幅’);subplot(2,1,2)plot(w/pi, imag(h)); grid;title(‘虚部’)xlabel(‘omega/\pi’);ylabel(‘振幅’);figure;subplot(2,1,1)plot(w/pi, abs(h)); grid;title(‘幅度谱’)xlabel(‘omega/\pi’);ylabel(‘振幅’);subplot(2,1,2)plot(w/pi, angle (h)); grid;title(‘相位谱’)xlabel(‘omega/\pi’);ylabel(‘以弧度为单位的相位’);修改程序，在范围0≤≤??内计算如下有限长序列的离散时间傅里叶变换h[n]=[1 2 3 4 5 6 7 8 9](2)利用(1)的程序，通过比较结果的幅度谱和相位谱，验证离散