实验二离散时间傅里叶变换.docx

实验二 离散时间傅里叶变换一．实验原理1、经由正、逆离散时间傅里叶变换表达的信号傅里叶表示式是信号分析的一个关键部分。X()=(3.9)（3.10）类似地，当LTI系统用于滤波时，作为冲击响应离散时间傅里叶变换的频率响应，提供了LTI系统简介的描述。离散时间傅里叶变换X()是的周期复值函数，周期总是2π，并且基周期通常选在区间[-π，π）上。对离散时间傅里叶变换DTFT来说有两个问题：DTFT的定义对无限长信号是有效的。DTFT是连续变量的函数。在MATLAB中，任何信号（向量）必须是有限长度的，仅此就是第一点成为问题。因此，不可能使用MATLAB计算无限长信号的DTFT。有一个值得注意的例外情形，当能从变换定义式推导出解析式并只是计算它时，可以使用MATLAB计算无限长信号的DTFT。2、对于频率抽样问题。MATLAB擅长在有限网格点上计算DTFT。通常选择足够多的频率以使绘出的图平滑，逼近真实的DTFT。对计算有利的最好选择是在（-π，π）区间上一组均匀地隔开的频率，或者对共轭对称变换选择[0,π]区间。采用上述抽样办法，DTFT式变成X()= DTFT的周期性意味着在-π≤0区间上的数值是那些对kN/2的数值。因为上市是在有限数量的频率点=2πk/N处计算，并在有限范围内求和，因此它是可计算的。由于信号长度必须是有限的（0≤nL），这个求和式不适用于x[n]=u[n]的情形。在对DTFT进行抽样时，并不要求N=L，尽管通常经由DFT进行计算。在正确应用FFT计算N点DFT前，需要对x[n]进行时间混叠。计算DTFT需要两个函数，MATLAB的freqz函数计算无限长信号，dtft（h，H）函数计算有限长信号的DTFT。二．实验内容：1.脉冲信号的DTFT设矩形脉冲r[n]由下式定义（1）要求：a.证明r[n]的DTFT可由下面的数学表达式得出该变换的第一项时常具有与DTFT相关的特殊形式，称为混叠sinc函数：b.使用dtft函数计算12点脉冲信号的DTFT。绘出在区间上对的DTFT。把实部和虚部分开绘出。另绘出DTFT的幅度。选择频率样本的数量是脉冲长度的5到10倍，以使绘出的图看上去平滑，用不同数量的频率样本做实验。绘图时注意要正确标出频率坐标的变量。c.注意asinc函数零点的位置是规则分布的。对奇数脉冲，比如L=15的脉冲重复进行DTFT计算并绘出幅度；同样再次检验零点位置，注意幅值高度。d.对asinc函数零点的间距与asinc函数的直流值，确定出通用规则。（2）程序：B-2N=96;L=12;h=ones(1,L);n=0:(L-1);[H,W]=dtft(h,N);subplot(312)plot(W/2/pi,real(H));grid;title(DTFT实部)ylabel(real(H1))xlabel(frequency)subplot(3,1,1)plot(W/2/pi,abs(H));grid;title(DTFT幅度)ylabel(abs(H1))xlabel(frequency)subplot(313)plot(W/2/pi,imag(H));grid;title(DTFT虚部)ylabel(real(H1))xlabel(frequency)function [H,W]=dtft(h,N)N=fix(N);L=length(h);h=h(:);if(NL) error(DTFT:#data samples cannot exceed# freq samples)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));endB-1N=60;L=12;h=ones(1,L);n=0:(L-1);[H,W]=dtft(h,N);subplot(312)plot(W/2/pi,real(H));grid;title(DTFT实部)ylabel(real(H1))xlabel(frequency)subplot(3,1,1)plot(W/2/pi,abs(H));grid;title(DTFT幅度)ylabel(abs(H1))xlabel(frequency)subplot(313)plot(W/2/pi,imag(H));grid;title(DTFT虚部)ylabel(real(H1))xlabel(frequency)B-3N=120;L=12;h=ones(1,L);n=0:(L-1);[H,W]=dtft(h,N);subplot(312)plot(W/2/pi,real(H));grid;title(DTFT实部)ylabel(real(H1))