实验三离散傅立叶变换.ppt

由题目可得 x=0.5*sin(2*pi*15*t)+2*sin(2*pi*40*t) fs=100 N=128/1024 例：已知信号由15Hz幅值0.5的正弦信号和40Hz幅值2的正弦信号组成，数据采样频率为100Hz,试绘制N=128点DFT的幅频图。 第三十页 fs=100; N=128; n=0:N-1; t=n/fs; x=0.5*sin(2*pi*15*t)+2*sin(2*pi*40*t); y=fft(x,N); f=(0:length(y)-1)*fs/length(y); mag=abs(y); stem(f,mag); title(‘N=128点’) 第三十一页 第三十二页 实验三离散傅立叶变换 第一页 一、实验目的 加深对离散傅立叶变换(DFT)的理解。 掌握利用MATLAB语言进行离散傅立叶变换和逆变换的方法。 加深对离散傅立叶变换基本性质的理解。 掌握离散傅立叶变换快速算法的应用。 第二页 二、实验原理及方法 建立以时间t为自变量的“信号”与以频率f为 自变量的“频率函数”(频谱)之间的某种变换关系。 所以“时间”或“频率”取连续还是离散值, 就形成 各种不同形式的傅里叶变换对。 傅里叶变换 第三页 四种不同傅里叶变换对 傅里叶级数(FS)：连续时间, 离散频率的傅里叶变换。周期连续时间信号傅里叶级数(FS)得到非周期离散频谱密度函数。 傅里叶变换(FT)：连续时间, 连续频率的傅里叶变换。非周期连续时间信号通过连续付里叶变换(FT)得到非周期连续频谱密度函数。 离散时间的傅里叶变换(DTFT):离散时间,连续频率的傅里叶变换。非周期离散的时间信号(单位园上的Z变换(DTFT))得到周期性连续的频率函数。 离散傅里叶变换(DFT):离散时间, 离散频率的傅里叶变换。 第四页 上面讨论的前三种傅里叶变换对,都不适用在计算机上运算, 因为至少在一个域( 时域或频域)中, 函数是连续的。因为从数字计算角度我们感兴趣的是时域及频域都是离散的情况, 这就是第四种离散傅里叶变换。 第五页 离散傅里叶级数(DFS) 离散时间序列x(n)满足x(n)=x(n+rN)，称为离散周期序列，其中N为周期，x(n)为主值序列。 由傅立叶分析知道周期函数可由复指数的线性组合叠加得到。其频率为基本频率的倍数。从离散时间傅立叶变换的频率周期性，我们知道谐波次数是有限的，其频率为 周期序列可表示成： 第六页 其中 叫做离散傅立叶级数系数，也称为周期序列的频谱，可由下式表示 注意 也是一个基本周期为N的周期序列。 上面两式称为周期序列的傅立叶级数变换对。 令 表示复指数，可以得到以下： 第七页 例：求出下面周期序列的DFS x(n)={……,0,1,2,3,0,1,2,3,0,1,2,3,……} 基本周期为N=4，WN=W4=-j， 因而 第八页 MATLAB实现 矩阵-向量相乘运算来实现。 由于 和 均为周期函数，周期为N，可设 和 代表序列 和 的主值区间序列，则前面的两个表达式可写成： 式中，矩阵WN为方阵——DFS矩阵。 第九页 利用MATLAB实现傅立叶级数计算 编写函数实现DFS计算 function xk=dfs(xn,N) n=[0:1:N-1]; %n的行向量 k=n; %k的行向量 WN=exp(-j*2*pi/N); %WN因子 nk=n’*k; %产生一个含nk值的N乘N维矩阵 WNnk=WN.^nk; %DFS矩阵 xk=xn* WNnk; %DFS系数行向量 第十页 例：xn=[0,1,2,3]，N=4 xn=[0,1,2,3]; N=4; xk=dfs(xn,N)’ 第十一页 逆运算IDFS function xn=idfs(xk,N) n=[0:1:N-1]; k=n; WN=exp(-j*2*pi/N); nk=n’*k; WNnk=WN.^(-nk); xn=（xk*WNnk）/N; 第十二页 xn=idfs(xk,4) x=xn 第十三页 周期重复次数对序列频谱的影响 理论上讲 ，周期序列不满足绝对可积条件，要对周期序列进行分析，可以先取K个周期进行处理，然后让K无限增大，