实验二离散时间傅里叶变换.
实验二 离散时间傅里叶变换
实验原理
1、经由正、逆离散时间傅里叶变换表达的信号傅里叶表示式是信号分析的一个关键部分。
)=(3.9)
(3.10)
类似地,当LTI系统用于滤波时,作为冲击响应离散时间傅里叶变换的频率响应,提供了LTI系统简介的描述。离散时间傅里叶变换X()是的周期复值函数,周期总是2π,并且基周期通常选在区间[-π,π)上。对离散时间傅里叶变换DTFT来说有两个问题:
(1)DTFT的定义对无限长信号是有效的。
(2)DTFT是连续变量的函数。
在MATLAB中,任何信号(向量)必须是有限长度的,仅此就是第一点成为问题。因此,不可能使用MATLAB计算无限长信号的DTFT。有一个值得注意的例外情形,当能从变换定义式推导出解析式并只是计算它时,可以使用MATLAB计算无限长信号的DTFT。
2、对于频率抽样问题。MATLAB擅长在有限网格点上计算DTFT。通常选择足够多的频率以使绘出的图平滑,逼近真实的DTFT。对计算有利的最好选择是在(-π,π)区间上一组均匀地隔开的频率,或者对共轭对称变换选择[0,π]区间。采用上述抽样办法,DTFT式变成X()=
DTFT的周期性意味着在-π≤<0区间上的数值是那些对k>N/2的数值。因为上市是在有限数量的频率点=2πk/N处计算,并在有限范围内求和,因此它是可计算的。由于信号长度必须是有限的(0≤n
在对DTFT进行抽样时,并不要求N=L,尽管通常经由DFT进行计算。在正确应用FFT计算N点DFT前,需要对x[n]进行时间混叠。
计算DTFT需要两个函数,MATLAB的freqz函数计算无限长信号,dtft(h,H)函数计算有限长信号的DTFT。
实验要求
理解数值计算在离散时间傅里叶变换中的作用。
实验内容
脉冲信号的DTFT
(1)要求:设矩形脉冲r[n]= 1 0≤n
0 其他
a.证明r[n]的DTFT可由(3.13)得出,记asinc(,L)(3.14)
b.使用dtft函数计算12点脉冲信号的DTFT。绘出在区间-π≤
c.注意asinc函数零点的位置是规律分布的。对奇数长脉冲,比如L=15的脉冲重复进行DTFT计算并绘出幅度,同样再次检验零点位置,注意峰值高度。
d.对于asinc函数零点的间距与asinc函数的直流值,确定出通用规则。
(2)程序
M文件
function [H,W] = dtft(h,N)
N=fix(N);
L=length(h);
h=h(:);
if(N
error
end
W=(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));
%b
nn=0:11;
u=ones(1,12);
[X,W]=dtft(u,72);
subplot(221),plot(W,real(X));grid,title('REAL RESPONSE')
xlabel('FREQUENCY W'),ylabel('REAL A')
subplot(222),plot(W,imag(X));
grid,title('IMAGE RESPONSE')
xlabel('FREQUENCY W'),ylabel('IMAGE A')
subplot(223),plot(W,abs(X));
grid,title('MAGNITUDE RESPONSE')
xlabel('FREQUENCY W'),ylabel('|H(w)|')
subplot(224),plot(W,angle(X));
grid,title('PHASE RESPONSE')
xlabel('FREQUENCY W'),ylabel('DEGREES')
运行结果
%c
nn=0:14;
u=ones(1,15);
[X,W]=dtft(u,90);
[Y,W]=dtft(X,90);
subplot(111),plot(W,abs(Y));
grid,title('MAGNITUDE RESPONSE')
xlabel('FREQUENCY W'),ylabel('|H(w)|')
运行结果
%d
如图L=12时由R(e^jw)=0得sin(wL/2)=0
即wL/2=k*pi 则w=k*pi/36
所以零点间距为pi/6
直流值:1