function [xk] = my_fft(x,n)
longx=length(x);
if n==0%如果输入序列只有一位,fft是他本身xk=x(1)*exp(-2*1j*pi*0*0);
elsex1=[];for i=1:longx/2x1=[x1,x(2*i-1)];end% 提取离散序列下标为 0 2 4 6...的值,构成偶序列x2=[];for i=1:longx/2x2=[x2,x(2*i)];end% 提取离散序列下标为 1 3 5 7...的值,构成偶序列x1k=my_fft(x1,n-1);%嵌套自己x2k=my_fft(x2,n-1);wnk=[];for i=1:longx/2wnk=[wnk,exp(-2*1j*pi/(longx)*(i-1))];endx1kk=x1k+wnk.*x2k;%蝶形运算x2kk=x1k-wnk.*x2k;xk=[x1kk x2kk];
end
end

适用于x为长2^n长度的序列

用matlab自己编写代码实现基2—FFT相关推荐

  1. matlab rbe神经网络代码,径向基神经网络知识介绍

    RBF网络原理 RBF网络,即径向基神经网络,也是前馈型网络的一种.它的设计思想和BP网络完全不一样. Cover定理:将复杂的模式分类问题非线性的投射到高维空间将比投射到低维空间更可能是线性可分的. ...

  2. 基2FFT算法matlab程序编写,频率抽取(DIF)基2FFT算法的MATLAB实现

    频率抽取(DIF)基2FFT算法和时间抽取(DIT)基2FFT算法是两种等价的FFT算法,其相同之处: (1)DIF与DIT两种算法均为原位运算. (2)DIF与DIT运算量相同. 不同之处: (1) ...

  3. MATLAB编写代码来模拟消色差超透镜的行为

    您好! 在 MATLAB 中编写消色差超透镜的代码可以通过利用光学理论来模拟它的行为. 首先,您需要定义一个函数来描述超透镜的光学特性,包括折射率和反射率的分布.然后,您可以使用 MATLAB 中的光 ...

  4. 长理2017 f0803 求和:s=1-(1/2)+(1/3)-(1/4).........+(1/n),编写函数fun(char*s),实现字符串的逆置编写程序从基类派生圆柱,设计(均为可运行代码)

    求和:s=1-(1/2)+(1/3)-(1/4).........+(1/n), #include <stdio.h>int main() {int n,i;double s=1.0;sc ...

  5. java编写代码用什么_如何学习用Java编写代码:为什么要学习以及从哪里开始

    java编写代码用什么 by John Selawsky 约翰·塞劳斯基(John Selawsky) 如何学习用Java编写代码:为什么要学习以及从哪里开始 (How to learn to cod ...

  6. 基于matlab的64QAM,通信调制体制设计之64QAM性能分析MATLAB仿真及代码

    通信调制体制设计之64QAM性能分析MATLAB仿真及代码 通信调制体制设计之64QAM性能分析MATLAB仿真及代码 任务背景 弗雷泽岛旅游经理在审查您之前建立无线链路任务的解决方案时,正在研究使用 ...

  7. matlab安装MinGW-64代码C/C++编译器的方法

    一.说明 matlab中编写的m文件转成C/C++代码,以及matlab中执行C/C++带啊,需要用到mex 查看matlab是否支持可以输入mex -setup: 二.安装 如果不支持,没找到,则需 ...

  8. matlab矩阵指定行最大值,求Matlab程序:在2行矩阵中,如何求第1行最大值和第2行相应的最大,请问,如何在MATLAB中编写程序实现求两矩阵A*B,A.*...

    导航:网站首页 > 求Matlab程序:在2行矩阵中,如何求第1行最大值和第2行相应的最大,请问,如何在MATLAB中编写程序实现求两矩阵A*B,A.* 求Matlab程序:在2行矩阵中,如何求 ...

  9. 基于蝙蝠算法优化BP神经网络的数据分类算法及其MATLAB实现-附代码

    基于蝙蝠算法优化BP神经网络的数据分类算法及其MATLAB实现-附代码 文章目录 基于蝙蝠算法优化BP神经网络的数据分类算法及其MATLAB实现-附代码 1 蝙蝠算法与BP神经网络分类模型 1.1 蝙 ...

最新文章

  1. 必须了解的MySQL三大日志:binlog、redo log和undo log
  2. 世界首富太空争霸:从地上斗到天上,马斯克VS贝索斯,谁能赢
  3. Nutch URL过滤配置规则
  4. 测试 / 单元测试的重要性
  5. python作业:高级FTP程序
  6. 《活出生命的意义》:人生有何意义?
  7. NOIP2011提高组day2
  8. 多进程服务器(python 版)
  9. fasterrcnn论文_【每周CV论文推荐】 深度学习人脸检测入门必读文章
  10. 单元测试java小白_一位Java应届毕业生工作7个月后的小结
  11. 【肥朝】从一次问题讨论聊聊我对阅读源码的思考
  12. Software--Architecture--SOA 面向服务体系结构
  13. 音视频播放器—快进快退及逐帧播放
  14. 布线问题—分支限界—java实现
  15. 普通硬盘和固态硬盘的区别?
  16. mysql建数据库实例_MySQL数据库基础(四)——MySQL数据库创建实例
  17. 一起聊聊什么是P问题、NP问题、NPC问题
  18. EasyNVR网页Chrome无插件播放摄像机视频功能二次开发之云台控制接口示例代码
  19. c语言筛法求一个范围内的素数,运用埃拉托色尼筛法求解一定范围内的素数.docx...
  20. spring.main.allow-bean-definition-overriding = true

热门文章

  1. java koala_Koala应用程序:用于现代Web开发的跨平台应用程序
  2. 超级网管员——网络应用
  3. 腾讯云人脸核身相关问题
  4. vim编辑文档提示readonly问题解决方案
  5. swiper 点击分页器之后不能继续轮播
  6. 关于写blog这件事
  7. matlab 图像转视频教程,Matlab制作视频并转换成gif动态图的两种方法
  8. 琐记之零散知识点(java相关的零散知识都丢这里)
  9. java 拦截器HandlerInterceptor 自定义返回结果
  10. 记录一下在ROS里调用ORK的苦逼历程