做完了傅里叶变换及其谐波分量的提取,并进行两组数据的比较。

clear all

clc

clear

原始信号

I1=[55.3 56.6 57.4 59.6 57.4 56.5 54.3 52.6 51.6

51.0 49.3 47.8 46.1 44.2 42.7 43.1 46.0 48.4 52.7 55.2 57.6 55.5

54.4 53.0 52.9 52.2 51.2 49.4 47.1 45.0 43.6 42.1 39.4 37.6 36.0

38.8 41.7 43.3 45.0 46.6 47.1 51.7 53.7 56.7 60.6 57.4 54.3 51.0

49.3 47.8 45.5 43.3 39.7 37.7 34.7 35.9 37.7 39.3 41.8 43.6 46.2

50.6 54.5 58.2 55.3 54.4 53.0 51.2 49.4 45.5 43.1 39.5 36.2 35.5

32.5];

I2=[31.8 32.6 34.6 34.8 30.2 32.5 32.8 31.9 30.9

32.3 34.6 32.6 33.8 32.7 30.5 31.1 33.5 32.7 33.4 34.5 32.3 30.8

31.6 33.5 35.6 32.2 34.5 33.8 31.9 32.5 32.4 31.6 32.6 33.1 34.2

35.2 34.4 32.6 32.8 31.6 33.2 33.8 31.6 35.4 33.6 32.0 34.5 32.1

33.6 34.5 33.4 35.8 34.2 32.9 32.1 30.7 31.6 33.4 35.2 34.8 35.4

36.1 32.8 31.5 33.6 31.2 30.4 32.4 33.5 31.7 32.0 30.4 30.7 34.1

33.5];

y=fft(I1);

figure,plot(fftshift(abs(y))),title('fft');

xlabel('Frequency (Hz)');

ylabel('amplitude');

for k=1:1:75

yt=abs(2/75*y(k));%谐波分量提取

end

两组数据傅里叶变换

第一组:

直流:95.8427

谐波:

2.3215 1.3957 3.1908 7.5419 3.1123 1.4075 0.8825 1.1974 1.7597 0.4917 0.9282 1.2333 0.2993 0.2221

0.4623 0.4592 0.6114 0.3298 0.5716 0.3230 0.2118 0.3583 0.3568 0.4724 0.3335 0.4274 0.4121 0.3366 0.4144 0.4678 0.5667 0.3161 0.3299 0.4093 0.2250 0.2029 0.3736 0.3736 0.2029 0.2250

0.4093 0.3299 0.3161 0.5667 0.4678 0.4144 0.3361 0.4121 0.4274 0.3335 0.4724 0.3568 0.3583 0.2158 0.3230 0.5716 0.3298 0.6114 0.4592 0.4623 0.2221 0.2993 1.2333 0.9282 0.4917 1.7597 1.1974 0.8825 1.4075 3.1123 7.5419 3.1908 1.3957 2.3215

n =

0 8 14 14 10 6 0 0 4 18

第二组;

直流:65.8853

谐波:

0.5206 0.2146 0.2236 0.3109 0.3324 0.6096 0.1841 0.2449 1.0108 0.1162 0.3072 0.2688 0.3097 0.1874 0.1829 0.1890 0.2626 0.5247 0.1296 0.2039 0.2067 0.6307 0.2577 0.1903 0.2582 0.2841 0.3508 0.4312 0.1697 0.2627 0.3272 0.1979 0.1341 0.0277 0.0353 0.1793 0.1014 0.1014 0.1793 0.0353 0.0277 0.1341 0.1979 0.3272 0.2627 0.1697 0.4312 0.3508 0.2841 0.2582 0.1903 0.2577 0.6307 0.2067 0.2039 0.1296 0.5247 0.2626 0.1890 0.1829 0.1874 0.3097 0.2688 0.3072 0.1162 1.0108 0.2449 0.1841 0.6096 0.3324 0.3109 0.2236 0.2146 0.5206

n =

12 26 22 4 4 4 0 0 0 2

还要继续往下做,提取出两组数据明显的区别。

python傅里叶谐波分析_傅里叶变换及其谐波提取相关推荐

  1. python音频 降噪_从视频中提取音频数据,然后应用傅里叶对音频降噪(python)...

    视频准备 QQ有热键 然后随便打开一个视频网站进行录屏 我选择B站 从视频中提取音频 需要安装包moviepy pip install moviepy 提取代码 from moviepy.editor ...

  2. python傅里叶逆变换_C# 傅里叶变换 逆变换 调用MathNet包|简明python教程|python入门|python教程...

    https://www.xin3721.com/eschool/pythonxin3721/ 上司给我讲采样率,我觉得自己有点钻牛角尖了,想来想去还是不明白. 因为说的是日语,サンプリング周波数,我就 ...

  3. python手写汉字识别_中文简历表格提取,手写汉字识别(Python+OpenCV)

    原标题:中文简历表格提取,手写汉字识别(Python+OpenCV) 所有代码获取: 简历 网上对表格框的提取的相关资料较少,尤其是Python+OpenCV的实现方面. 整体流程 如今OpenCV官 ...

  4. python timedelta函数_通过日期字段提取年月日、timedelta提取时分秒并进行小时汇总...

    原博地址 知识梳理不易,请尊重劳动成果,文章仅发布在CSDN网站上,在其他网站看到该博文均属于未经作者授权的恶意爬取信息 如若转载,请标明出处,谢谢! 1 业务需求 在进行数据分析,经常会遇到时间处理 ...

  5. python低通滤波器_[Python图像处理] 二十三.傅里叶变换之高通滤波和低通滤波

    [Python图像处理] 二十三.傅里叶变换之高通滤波和低通滤波 发布时间:2019-04-28 20:35, 浏览次数:2739 , 标签: Python 该系列文章是讲解Python OpenCV ...

  6. 《OpenCv视觉之眼》Python图像处理十二 :Opencv图像轮廓提取之基于一阶导数的Roberts算法、Prewitt算法及Sobel算法

    本专栏主要介绍如果通过OpenCv-Python进行图像处理,通过原理理解OpenCv-Python的函数处理原型,在具体情况中,针对不同的图像进行不同等级的.不同方法的处理,以达到对图像进行去噪.锐 ...

  7. 《OpenCv视觉之眼》Python图像处理十四 :Opencv图像轮廓提取之Scharr算法和Canny算法

    本专栏主要介绍如果通过OpenCv-Python进行图像处理,通过原理理解OpenCv-Python的函数处理原型,在具体情况中,针对不同的图像进行不同等级的.不同方法的处理,以达到对图像进行去噪.锐 ...

  8. python 时间序列预测_使用Python进行动手时间序列预测

    python 时间序列预测 Time series analysis is the endeavor of extracting meaningful summary and statistical ...

  9. python机器学习预测_使用Python和机器学习预测未来的股市趋势

    python机器学习预测 Note from Towards Data Science's editors: While we allow independent authors to publish ...

最新文章

  1. 表达不同与构建不同: 对计算机的唯一真正的要求 (TrustNo.1 ) -- 待修改!!!...
  2. Hive远程模式安装(1.00)
  3. 13.2.8 用户认证
  4. Leet Code OJ 27. Remove Element [Difficulty: Easy]
  5. 不要重新发明轮子_是否重新发明轮子
  6. 数字游戏(水博客ing / csp-J T1 / luogu 5660)
  7. v-region 地址组件无法回显问题(选择地址保存后,打开编辑,发现刚才选择的地址【省市县乡】信息消失,信息不回显)
  8. “你都硕士博士了,竟然还不如我!”
  9. Spring Boot 项目的这些文件都是干啥用的?
  10. 【Vegas原创】重建Exchange 2007 OWA的虚拟目录
  11. 纯java程序jdbcjar包放在哪_oracle 11g jdbc jar包在哪个文件目录
  12. 字符集不同导致的ORA-00972
  13. java String 常用方法集合
  14. python线性链表
  15. matlab rsenc函数,RS码的MATLAB实现(附MATLAB仿真程序).pdf
  16. Cisco 路由器密码重置步骤
  17. Android 开发 -- 开发第一个安卓程序、Android UI开发(布局的创建:相对布局和线性布局、控件单位:px pt dp sp、常用控件 、常见对话框、ListView)
  18. 对话《深入理解Java虚拟机》作者周志明:电竞选手成为Java大神之路
  19. c语言洗牌发牌结构体,C语言实现洗牌发牌小程序
  20. android热更新bugly

热门文章

  1. 电镀废水处理工艺 离子交换树脂回收镍 硫酸镍经济效益越来越高 电镀废水如何获利 含重金属废水处理工艺
  2. 魅族flashfire_魅族Pro6s更新安卓7.1后完美root及xp框架
  3. 在通信设备商工作那几年技术上的得与失
  4. android 快盘接口,Android云盘试用体验第二期:金山快盘
  5. Science Robotics 封面论文 | 绝影狗:陌生环境自适应“摸爬滚打”「AI核心算法」...
  6. 【zz】免费的午餐——编程利用Google日历API发短信、Email
  7. 适用于多类移动场景的融合认证机制设计
  8. 老男孩mysql 百度云_英语语录:除了你,没人能掌控你的幸福
  9. 给自己八年前的“网红面经”写个续篇
  10. 微信IPAD协议(最新) 支持小程序授权、公众号授权、APP授权、抢红包