信号与系统实验六.docxVIP

信号鸟索傥实验 实验六离散时间信号与系统的频域分析 小组成员： 黄涛胡焰焰』084219— 洪燕一、 实验目的 1、 学握离散时间信号与系统的频域分析方法，从频域的角度对信号与系统的特性进行 分析。 2、 掌握离散时间信号傅里叶变换与傅里叶逆变换的实现方法。 3、 掌握离散吋间傅里叶变换的特点及应用 4、 掌握离散时间傅里叶变换的数值计算方法及绘制信号频谱的方法 二、 预习内容 1、 离散时间信号的傅里叶变换与逆变换。 2、 离散时间信号频谱的物理含义。 3、 离散吋间系统的频率特性。 4、 离散时间系统的频域分析方法。 三、 实验原理 离散时间系统的频率特性 在离散LTI系统时域分析中得到系统的单位冲激响应可以完全表征系统，进阳通过灿川 特性来分析系统的特性。系统单位冲激响应/卄］的傅里叶变换H（严）成为LTI系统的频率 响应。与连续时间LTI系统类似，通过系统频率响应可以分析出系统频率特性。 与系统单位冲激响应灿加一样，系统的频率响应H（严）氏映了系统内在的固有特性， 它取决于系统自身的结构及组成系统元件的参数，与外部激励无关，是描述系统特性的一个 重要参数，H（严是频率的复函数口J以表示为： 其屮，IH（RQ）I随频率变化的规律称为幅频特性：0（初随频率变化的规律称为相频特 性。 离散时间信号傅里叶变换的数值计算方法 算法原理，由傅里叶变换原理可知： F （严）=￡伽才创 序列/［刃的离散时间傅里叶变换F （严）是?的连续函数。由于数据在matlab中以向 量的行驶存在，F（e~只能在一个给定的离散频率的集合中计算。然而，只有类似 FD= “0 + P 严+ #2幺一八+ ...+ /切旷 I 丿—+ d{e^ + d2ej2a, +... + clNe~j0,N 形式的旷用的有理函数，才能计算其离散时间傅里叶变换。 涉及到的Matlab函数 3」freqz函数：实现离散时间系统频率响应特性的求解 调用格式:[H,w]=freqz(B, A,N) B和A分别为离散系统的系统函数分子、分母多项式的系数向量，返回量H则包含了 离散系统频响在()?兀范围内N个频率等分点的值(其屮N为正整数)，w则包含了范囤内 N个频率等分点。调用默认的N时，其值是512。由于FM)是的连续函数，需要尽可 能大地选取N的值，以使得产生的图形和真实离散傅里叶变换的图形尽可能一致。为更加 方便快速地运算，应将N的值选为2的幕，如256或512. 3.2 real函数：求复数的实部 调用格式:rcal_f=rcal (f); 3.3imag函数：求复数的虚部 调用格式:imag_f=imag (f)； abs函数：求复数的模 调用格式:abs_f=abs (f)； angle函数：求复数的相位 调用格式:angle_f=angle (f)； fft函数：实现离散信号f[n]的傅里叶变换值F[k] 调用格式： F = fh(f),计算序列/1加的离散傅里叶变换值F[k]f其中F[幻长度与伽]长度相同 F二fft(f,L),计算序列/[]的L点离散傅里叶变换值F[幻，其中L不小于N。若L大 于N,则需要在计算离散傅里叶变换Z前，对/[川尾部补足L?N个零。 ifft函数：实现离散信号的傅里叶逆变换 调用格式：f=ifft(F) 四、实验内容 1离散时间傅里叶变换 (1)下面参考程序是如下序列在范围4兀的离散时间傅里叶变换 F(e)= 2 + e I 7 1-0.6 严 %计算离散吋间傅里叶变换的频率样本 clear all; w=-4*pi;8*pi/511;4*pi; num=[2 1]; den=[l -0.6]; h=freqz(num,den,w); subplot(2J J) plot(w/pi,real(h)); grid; title。实部 J xlabelC omegaApi J; ylabclC 振幅)； subplot(2,l,2) plot(w/pi, imag(h)); grid; titled 虚部 J xlabelComegaApi); ylabel(振幅)； figure; subplot(2,1,l) plot(w/pi, abs(h)); grid; titled幅度谱J xlabel(*omegaApi,); ylabel(4 振幅 J; subplot(2,1,2) plot(w/pi, angle (h)); grid; title(湘位谱J xlabel(*omegaApi,); ylabcR4以弧度为单位的相位)； 修改程序，在范m0CO7T内计算如卜?有限长序列的离散时间傅里叶变换 hlnJ=Ll 2345678 9] clear all; w=():pi/51 l:p