数字的信号处理实验的2离散时间傅立叶变换

《数字的信号处理实验的2离散时间傅立叶变换》由会员上传分享，免费在线阅读，更多相关内容在工程资料-天天文库。

1、实用标准文案信号处理实验实验二：离散时间傅立叶变换精彩文档实用标准文案一、实验题目：离散时间傅里叶变换二、实验原理经由正、逆离散时间傅里叶变换表达的傅里叶表示式是信号分析的一个关键部分，下面是分析方程与综合方程。由以上公式知，离散时间傅里叶变换是w的周期复值函数，周期是2π，并且周期常选为【-π,π】.对离散时间傅里叶变换有两个问题：（1）DTFT的定义对无限长信号是有效的。（2）DTFT是连续变量的w函数。第二个问题是频率抽样问题。Matlab擅长在有线网格点上计算DTFT。通常选择足够多的频率以使绘

2、出的图平滑，逼近真实的DTFT。对计算有利的最好选择是在（-π，π）区间上一组均匀的隔开的频率，或者共轭对称变换选择【0，π】，采用上述抽样方法，DTFT式变为X（ejw）=X(ej2πk/N)=n=0L-1e-j(2πk/N)n,k=0,1,……N-1在对DTFT进行抽样时，并不要求N=L，尽管通常由DFT进行计算时，如果N=L计算很方便。通常，不可能计算一个无限长想信号的DTFT。但有一个重要的类型，其计算式容易的。这一类型的信号就是指数信号，其DTFT是e-jw有理函数。H（ejw）=B(ejw)

3、A(ejw)=l=0Qble-jwlk=0Pake-jwk指数信号h[n]=anu[n]是这类信号的一员，但是对它不能使用前面的dtft函数来处理。另一方面，很容易推导出它的dtft的表达式：若

4、a

5、<1,有h[n]=anu[n]H(ejw)=n=0∞anu[n]e-jwn=11-ae-jw三、实验内容（1）脉冲信号的DTFT设矩形脉冲r[n]由下式定义r[n]=10≤n≤L0其他a．证明r[n]的dtft可有下面的数学表达式得出R(ejw)=sin⁡(12wL)sin⁡(12w)e-jw(L-1)/2

6、该变换的第一项具有与dtft相关的特殊形式，称为混叠sinc函数：asinc(w,L)=精彩文档实用标准文案b．使用dtft函数计算12点脉冲信号的dtft。绘出在区间-pi≤w≤pi上对w的dtft。把实部和虚部分开绘出，但要注意这些图不是很有用。另绘出dtft的幅度。选择频率样本的数量是脉冲长度的5~10倍，以使绘出的图看上去平滑。程序：首先定义一个dtft函数：function[H,W]=dtft(h,N)N=fix(N);L=length(N);h=h(:);if(N

7、:#datasamplescannotexceed#freqsamples')endW=(2*pi/N)*[0:(N-1)]';mid=ceil(N/2)+1;W(mid:N)=W(mid:N)-2*pi;W=fftshift(W);H=fftshift(fft(h,N));计算12点脉冲的dtft：formatcompact,subplot(111)xn=ones(12,1);[X,W]=dtft(xn,120);subplot(211),plot(W,real(X));grid,title('DTF

8、T的实部')精彩文档实用标准文案xlabel('w'),ylabel('Re')subplot(212),plot(W,imag(X));gridxlabel('w'),ylabel('Im')title('DTFT的虚部')绘出在区间-pi≤w≤pi上对w的dtft（实部和虚部分开绘出）dtft的幅度formatcompactxn=ones(12,1);[X,W]=dtft(xn,120);plot(W,abs(X));grid,title('DTFT的幅度')xlabel('w'),ylabel('

9、abs')精彩文档实用标准文案c.注意asinc函数零点的位置是规则分布的。对奇数长脉冲，比如L=15的脉冲重复进行dtft计算并绘出幅度。程序：formatcompactxn=ones(15,1);[X,W]=dtft(xn,150);plot(W,abs(X));grid,title('DTFT的幅度')xlabel('w'),ylabel('abs')精彩文档实用标准文案d.对asinc函数零点的间距与asin函数的直流值，确定出通用规则。由图可知，asinc函数零点的间距一定。(2)asinc的

10、m文件编写一个matlab文件如asic（w，L），之间从式中计算在频格上的asinc（w，L），该函数有两个输入：L和W，函数必须检查被0除的情形。直接计算混叠sinc函数得到脉冲信号DTFT绘出幅度，保存该图以便与dtft得到的结果比较。程序：首先定义asinc函数functionq=asinc(w,L)ifw==0q=L;%检查被零除的情形elseq=sin(w*L/2)./sin(w/2);end运用混叠函数asinc计算12点脉冲