实验三离散时间信号的频域分析韩梦欣.docVIP

实验三 离散时间信号的频域分析 实验目的 熟悉并应用信号的频域分析方法对信号进行处理； 掌握信号的复变量或的线性组合表示方式，能够将时域的离散序列映射到频域处理； 利用MATLAB计算离散时间信号的DTFT和DFT，并理解其相互关系。 二、实验原理 DTFT和DFT的定义及其相互关系。序列DTFT定义为。它是关于自变量的复函数，且是以为周期的连续函数。可以表示为 其中，分别是的实部和虚部；还可以表示为 其中，的幅度函数和相位函数；它们都是w的实函数，也是以为周期的周期函数。 序列的N点DFT定义为 是周期为N的序列。 与的关系：是对在一个周期中的谱的等间隔N点采样，即 而可以通过对内插获得，即 使用到的MATLAB命令有基于DTFT离散时间信号分析函数以及求解序列的DFT函数。 1）基于DTFT离散时间信号分析函数有。函数可以用来计算一个以的有理分式形式给出的序列的DTFT值。的形式多样，常见的有，其中表示序列有理分式DTFT的分子多项式系数，表示分母多项式系数（均按z的降幂排列），矢量w表示在0~中给定的一系列频率点集合。函数的其他形式参见帮助文件。在求出DTFT值后，可以使用函数分别求出并绘出其实部、虚部、幅度和相位谱。例如利用函数，然后利用函数分别求出幅频特性和相频特性，最后利用命令绘出曲线。 2） 求解序列DFT的函数有：fft，ifft。函数可以计算R点序列的R点DFT值；而则计算R点序列的N点DFT，若RN，则直接截取R点DFT的前N点，若RN，则x先进行补零扩展为N点序列再求N点DFT。函数可以计算R点的谱序列的R点IDFT值；而同的情况。 实验设备 计算机、MATLAB软件 四、实验内容 1.编程计算并画出下面DTFT的实部、虚部、幅度和相位谱。 程序如下 a=[0.0518 -0.1553 0.1553 0.0518];%分子多项式 b=[1 1.2828 1.0388 0.3418]; %分母多项式 [H w]=freqz(a,b,512); %一般默认值为512 subplot(2,2,1),plot(w,real(H)),title(实部);%绘出实部图像 subplot(2,2,2),plot(w,imag(H)),title(虚部);%绘出虚部图像 subplot(2,2,3),plot(w,abs(H)),title(幅度); %绘出幅度图像 subplot(2,2,4),plot(w,angle(H)),title(相位);%绘出相位图像 图像如下 2.计算32点序列，0n31的32点和64点DFT，分别绘出幅度谱图形，并绘出该序列的DTFT图形。 （1）程序如下 n=0:31; %n的范围：0=n=31 x=cos(((5*pi)/16)*n); xd=fft(x,32); %调用fft函数，计算32点序列的32点DFT subplot(2,1,1),plot(abs(xd)),title(幅度谱(32点)); %绘32点幅度谱 xd=fft(x,64); %调用fft函数，计算32点序列的64点DFT。（RN)先对x进行补零扩展为64点，再计算 subplot(2,1,2),plot(abs(xd)),title(幅度谱(64点)); %绘64点幅度谱 图像如下 对形式为的序列DTFT编程，绘出一个周期中实部、虚部、幅度及相角的图形 （1）程序如下 b=[ 0.5 -0.25] a=[0.5 -0.55 1]; [H w]=freqz(b,a,513); subplot(2,2,1),plot(w,real(H)),title(实部); subplot(2,2,2),plot(w,imag(H)),title(虚部); subplot(2,2,3),plot(w,abs(H)),title(幅度); subplot(2,2,4),plot(w,angle(H)),title(相位); （2）图像如下 计算有限长序列的DFT和IDFT，绘出其幅度谱图形 程序如下 A=1; n=0:15; %n的范围：0=n=15 x=cos(5*pi/16*n); %有限长序列 B=x; w=0:0.01:2*pi*2; [H]=freqz(B,A,w); %计算频率响应 magH=abs(H); %计算幅度 phaH=angle(H); %计算相位 subplot(2,1,1); plot(w,magH); %画幅度响应曲线 grid; %画网格 ylabel(Magnitude