实验二 离散时间信号的频域分析.docxVIP

实验二 离散时间信号的频域分析 实验目的 信号的频域分析是信号处理中的一种有效的工具。在离散信号的时域分析中，通常将信号表示成单位采样序列δ(n)的线性组合，而在频域中，将信号表示成复变量e-jwn或e-j2πNn 本实验中，将学习利用MATLAB计算离散时间信号的DTFT和DFT，并加深对其相互关系的理解。 实验原理 （1） DTFT和DFT的定义及其相互关系。序列x(n)DTFT定义，X(ejw)=n=-∞∞x(n)e-jnw,它是关于自变量w的复变函数，且是以2π为周期的连续函数。 X(ejw)=Xre(ejw)+jXim(ejw) 其中，Xre(ejw)和Xim(ejw)分别是X(ejw)的实部和虚部；还可以表示为 X(ejw)=X(ejw 它们都是w的实函数，也是2π为周期的周期函数。 （2） 使用到的MATLAB命令有基于DTFT离散时间信号分析函数以及求解序列的DFT函数。 1）基于DTFT离散时间信号函数有：freqz,real,imag,abs,angle。函数freqz可以用来计算一个以ejw的有理分式形式给出的序列的DTFT值。Freqz的形式多样，常见的有H=freqz(num,den,w)，其中num表示序列有理分式DTFT的分子多项式系数，den表示分母多项式系数（均按z的降幂次排列），矢量w表示在0～2π中给定的一系列频率点集合。在求出DTFT直后，可以使用函数real,imag,abs和angle分别求出并绘出其实部，虚部，幅度和相位谱。 2）求解序列DFT的函数有：fft,ifft。函数fft(x)可以计算R点序列的R点DFT值；而fft(x,N)则计算R点序列的N点DFT，若RN，则直接截取R点DFT的前N点，若RN，则x先进行补零扩展为N点序列再求N点DFT。函数ifft(X)可以计算R点的谱序列的R点IDFT值；而ifft(X,N)同fft(x,N)的情况。 实验内容 编程计算并画出下面DTFT的实部，虚部，幅度和相位谱。 X(ejw)=0.0518-0.1553e 编制该函数的程序命令如下，绘出图形如图1-1所示。 num=[0.0518 -0.1553 0.1553 0.0518]; den=[1 1.2828 1.0388 0.3418]; w=0:0.001:2*pi; H=freqz(num,den,w); figure subplot(221) plot(w,real(H)) title(实部) grid subplot(222) plot(w,imag(H)) title(虚部) grid subplot(223) plot(w,abs(H)) title(幅度) grid subplot(224) plot(w,angle(H)) title(相位谱) 图1-1 计算32点序列x(n)=cos5π16n,0≤n≤31的32点和64点DFT，分别绘出幅度谱图形，并绘出该序列的DTFT 该序列的32点和64点DFT命令如下，幅度谱图如图2-1所示。 n1=0:31; n2=0:63; x1=cos((5*pi*n1)/16); x2=cos((5*pi*n2)/16); X1=fft(x1); X2=fft(x2); subplot(211) plot(n1,abs(X1)) title(32点的DFT) subplot(212) plot(n2,abs(X2)) title(64点的DFT) 图2-1 该序列32点和64点DTFT命令如下，图形如图2-2所示。 n=[0:1:199]; x=cos(5*pi*n/32); n1=[0:1:14];y1=[x(1:1:32),zeros(1,168)]; Y=fft(y1);y2=abs(Y(1:1:101));k2=0:1:100;w2=2*pi/100*k2; subplot(211) plot(w2/pi,y2); title(32点的DTFT) n1=[0:1:14];y1=[x(1:1:64),zeros(1,168)]; Y=fft(y1);y2=abs(Y(1:1:101));k2=0:1:100;w2=2*pi/100*k2; subplot(212) plot(w2/pi,y2); title(64点的DTFT) 图2-2 实验设备 计算机，MATLAB软件。 实验步骤 （1）对形式为X(eiw)=P(ejw)D(e （2）计算有限长序列的DFT和IDFT，绘出其幅度谱图，并解释DFT和DTFT的关系。 实验总结 通过本次实验，复习运用MATLAB软件对信号的运算及绘图功能，学会用MATLAB计算离