MATLAB实现快速卷积算法函数

  • p5_190635
    了解作者
  • 711B
    文件大小
  • rar
    文件格式
  • 0
    收藏次数
  • VIP专享
    资源类型
  • 0
    下载次数
  • 2022-06-08 05:10
    上传日期
编写一个MATLAB函数,实现快速卷积算法 用一个N点离散傅立叶变换同时计算两个N点实序列的离散傅立叶变换
ksjj.rar
  • 2.txt
    1.3KB
内容介绍
1.编写一个MATLAB函数,实现快速卷积算法。 函数实现: function y=clconvt(x1,x2,N) if length(x1)>N errordlg('N必须>=x1的长度!’,’错误!’); elseif length(x2)>N errordlg('N必须>=x2的长度!’,’错误!’); else x1=[x1 zeros(1,N-length(x1))]; x2=[x2 zeros(1,N-length(x2))]; y1=fft(x1,N); y2=fft(x2,N); Y=y1.*y2; y=ifft(Y); end 例 x=[1,2,2,1]; >> h=[1,2,3,4,5]; >> N=9; >> y=clconvt(x,h,9) y = Columns 1 through 8 1.0000 4.0000 9.0000 15.0000 21.0000 21.0000 14.0000 5.0000 Column 9 0.0000 >> stem(y) 2.编写一个MATLAB函数,用一个N点离散傅立叶变换同时计算两个N点实序列的离散傅立叶变换,并将结果与直接使用两个N点离散傅立叶变换计算出来的结果进行比较。 函数实现如下: function [X1,X2]=qwdft(x1,x2,N) x=x1+j*x2; X=fft(x,N); Y=conj(X); K=0:N-1; X1=0.5*(X+Y(mod(-K,N)+1)); X2=-j*0.5*(X-Y(mod(-K,N)+1)); 源代码: x1=[1,3,2,4]; >> x2=[3,2,6,1]; >> [X1,X2]=qwdft(x1,x2,4) X1 = 10.0000 -1.0000 + 1.0000i -4.0000 -1.0000 - 1.0000i X2 = 12.0000 -3.0000 - 1.0000i 6.0000 -3.0000 + 1.0000i >> x1=[1,3,2,4]; >> X1=fft(x1) X1 = 10.0000 -1.0000 + 1.0000i -4.0000 -1.0000 - 1.0000i >> x2=[3,2,6,1]; >> X2=fft(x2) X2 = 12.0000 -3.0000 - 1.0000i 6.0000 -3.0000 + 1.0000i
评论
    相关推荐
    • matlab算法经典程序
      一些常用的算法matlab实现,包括各种智能算法,图论算法,以及规划算法
    • 电力系统状态估计MATLAB算法
      状态估计算法 MATLAB 内附readme 详细说明了使用方法和步骤 有专门的txt文件 可以输入自己的bus阵 line阵等 即可进行状态估计
    • matlab算法大全
      matlab 中常用的程序 函数使用示例
    • 粒子滤波matlab算法
      粒子滤波正在得到重视,越来越的人开始注意,在故障诊断和预测领域里,用来估计状态参数。附件中matlab程序介绍了基本粒子滤波算法
    • matlab算法大全
      本文档涵盖了很多数学算法,利用matlab实现,在工程技术上有很高的应用价值
    • Matlab算法大全
      Matlab算法 PDF 分章阅读 高清PDF 很好很全面的阅读材料
    • 各种MATLAB算法
      各种MATLAB算法,代码,可供初学者参考应用,
    • Matlab算法大全
      Matlab算法大全 第01章线性规划 第02章整数规划 第03章非线性规划 第04章动态规划 第05章图与网络 第06章排队论 第07章对策论 第08章层次分析法 第09章插值与拟合 第10章数据的统计描述和分析 第11章方差分析 第12章...
    • matlab算法大全
      matlab算法大全,很全,提供插值,数值微分、积分等功能代码的实现
    • matlab 算法程序
      matlab 算法程序,包括了插值、函数逼近、数值积分、非线性方程求解、统计分析、偏微分方程数值解法等17个部分,每个部分针对各种函数有m文件代码和相关解释说明。