• 设为首页
  • 点击收藏
  • 手机版
    手机扫一扫访问
    迪恩网络手机版
  • 关注官方公众号
    微信扫一扫关注
    公众号

用matlab脚本语言写M文件函数时用三种方法简单实现实现DFT(离散傅里叶变换) ...

原作者: [db:作者] 来自: [db:来源] 收藏 邀请
%用二重循环实现DFT:
function xk=dt_0(xn);     %define a function
N=length(xn);           %caculate the length of the variable
WN=exp(-j.*2.*pi./N);    
xk=zeros(1,N);          %define a non-zero 一维矩阵 
sum=zeros(1,N);          %define a non-zero 一维矩阵 
for k=1:N            %二重循环实现离散傅里叶变换DFT
    
    for n=1:N
        
       sum(n)=xn(n).*WN.^(k.*n);
       xk(k)=xk(k)+sum(n);
       
 
    end         
    
    
    
end
end



%用一重循环和内积实现DFT:
function xk=dt_1(xn);

 N=length(xn);
 WN=exp(-j.*2.*pi./N);    
 xk=zeros(1,N);      
  n=[0:N-1];
 for k=0:1:N-1;
   xk(k+1)=xn*WN.^(k.*n\');  %此处下标一定得从1开始,因为matlab的下标是从1开始的

 end    
    
end



%不用循环,仅有内积相乘实现DFT:
function xk=dt_2(xn);
N=length(xn);

WN=exp(-j*2*pi/N);   
n=0:1:N-1;     %定义一个一维矩阵,即行向量,从0到N-1
k=0:1:N-1;
nk=k\'*n;      %行向量k变换为列向量 乘上 行向量n ,得到一个N x N的矩阵
WNnk=WN.^(nk);      %做幂运算后的参数仍为一个 N x N的系数矩阵
xk=xn*WNnk;       %行向量 乘以  N x N的系数矩阵 即为DFT变换后的矩阵



end

以下是输入一个行向量xn=[1,2,3,4],MATLAB中用以上三种方法进行DFT的结果如下图所示:



 


鲜花

握手

雷人

路过

鸡蛋
该文章已有0人参与评论

请发表评论

全部评论

专题导读
上一篇:
Delphi对Excel表格的各种操作发布时间:2022-07-18
下一篇:
Delphi中字符串补齐方法发布时间:2022-07-18
热门推荐
阅读排行榜

扫描微信二维码

查看手机版网站

随时了解更新最新资讯

139-2527-9053

在线客服(服务时间 9:00~18:00)

在线QQ客服
地址:深圳市南山区西丽大学城创智工业园
电邮:jeky_zhao#qq.com
移动电话:139-2527-9053

Powered by 互联科技 X3.4© 2001-2213 极客世界.|Sitemap