matlab实验程序

liweiyehot

运行结果.pdf (35.68 KB, 下载次数: 7933)

2011-12-26 10:56 上传

点击文件名下载附件
如下三个m文件用来求函数的DFT（傅里叶变换）和IDFT（傅里叶反变换）

***********主程序为求函数xn=3*exp(n)的傅里叶变换，及移位后的序列fn=xn(mod((n+m),N)+1)的傅里叶反变换。**************

testdecDfIDf.m（主程序）
clear all;
N=15;               %有限长序列的长度
n=0:N-1;            
xn=3*exp(n);        %移位前序列
m=5;                %移动长度，左移为正数，右移为负数
fn=xn(mod((n+m),N)+1);%移位后序列
XK=DFT(xn);         %将前面的DFT程序变换成子程序，这里直接调用
FK=IDFT(fn);        %移位后序列f(n)的IDFT
subplot(2,2,1);stem(1:N,xn);
xlabel('(a) 时域离散序列x(n)');
subplot(2,2,2);stem(1:N,XK);
xlabel('(b) x(n)的DFT变换结果');
subplot(2,2,3);stem(1:N,fn);
xlabel('(a) f(n)');
subplot(2,2,4);stem(real(FK));
xlabel('(b) f(n)的IDFT变换结果');


**********将傅里叶变换编写成m文件，方便以后直接调用*************

DFT.m（傅里叶变换）
function Xk=DFT(xn)
N=length(xn);
WNnK=dftmtx(N);                             %求旋转因子矩阵
Xk=xn*WNnK;                                 %计算x(n)的DFT

**********将傅里叶反变换编写成m文件，方便以后直接调用*************

IDFT.m（傅里叶反变换）
function FK=IDFT(fn)
N=length(fn);
k=0:N-1;
n=0:N-1;
WN=exp(-1j*2*pi/N);
WNnK=WN.^(-n'*k);           %计算旋转因子矩阵
FK=fn*WNnK/N;               %计算F(K)的IDFT

liweiyehot

请各位高手多多指教哦！

EnSafe

最好是有简单的说明文字，说明这个例子的目的、结果的含义。图应有坐标。

liweiyehot

EnSafe 发表于 2011-12-27 09:32
最好是有简单的说明文字，说明这个例子的目的、结果的含义。图应有坐标。

谢谢楼上的提醒，！