修改代码,通过arduino中的fft库实现幅频转换

时间: 2024-04-16 17:28:24 浏览: 17
好的,我已经修改代码,使用 Arduino 的 FFT 库来实现幅频转换。请注意,这里假设你已经安装了 Arduino FFT 库,并将其包含在代码中。 ```c++ #include <Wire.h> #include <Adafruit_GFX.h> #include <Adafruit_SSD1306.h> #include <arduinoFFT.h> #define OLED_RESET 4 Adafruit_SSD1306 display(OLED_RESET); const int wavePin = A0; const int numSamples = 128; const int screenWidth = 128; const int screenHeight = 64; const int fftHeight = 32; const int waveOffsetY = 16; float waveData[numSamples]; float fftData[numSamples / 2]; arduinoFFT fft; void setup() { Serial.begin(9600); display.begin(SSD1306_SWITCHCAPVCC, 0x3C); display.clearDisplay(); display.setTextColor(WHITE); display.setTextSize(1); pinMode(wavePin, INPUT); for (int i = 0; i < numSamples; i++) { waveData[i] = 0; fftData[i] = 0; } fft.windowFunction(FFT_BLACKMAN_HARRIS); fft.reorder(1); } void loop() { getWaveData(); calculateFFT(); display.clearDisplay(); // 右侧波形图 drawWaveform(); // 右侧 FFT 转换图 drawFFT(); // 左侧信息 drawInfo(); display.display(); } void getWaveData() { for (int i = numSamples - 1; i > 0; i--) { waveData[i] = waveData[i - 1]; } int sensorValue = analogRead(wavePin); float voltage = map(sensorValue, 0, 1023, 0, 5) - 2.5; waveData[0] = voltage; } void calculateFFT() { fft.run(waveData, fftData, FFT_FORWARD); } void drawWaveform() { int waveOffsetX = screenWidth / numSamples; for (int i = 0; i < numSamples; i++) { int x = i * waveOffsetX; int y = map(waveData[i], -2.5, 2.5, screenHeight - waveOffsetY, waveOffsetY); display.drawPixel(x, y, WHITE); } } void drawFFT() { int fftOffsetX = screenWidth / (numSamples / 2); for (int i = 0; i < numSamples / 2; i++) { int x = i * fftOffsetX; int h = map(fftData[i], -100, 100, fftHeight, 0); display.drawLine(x, screenHeight - h, x, screenHeight, WHITE); } } void drawInfo() { display.setCursor(0, 0); // 左侧幅值 display.print("Amplitude: "); display.println(getAmplitude()); // 左侧中值 display.print("Median: "); display.println(getMedian()); // 左侧频率 display.print("Frequency: "); display.println(getFrequency()); // 左侧采样率 display.print("Sampling Rate: "); display.println(getSamplingRate()); } float getAmplitude() { float amplitude = 0; for (int i = 0; i < numSamples; i++) { amplitude += abs(waveData[i]); } return amplitude / numSamples; } float getMedian() { float sortedData[numSamples]; for (int i = 0; i < numSamples; i++) { sortedData[i] = waveData[i]; } // 使用冒泡排序对数据进行排序 for (int i = 0; i < numSamples - 1; i++) { for (int j = 0; j < numSamples - i - 1; j++) { if (sortedData[j] > sortedData[j + 1]) { float temp = sortedData[j]; sortedData[j] = sortedData[j + 1]; sortedData[j + 1] = temp; } } } return sortedData[numSamples / 2]; } float getFrequency() { float maxVal = -1; int maxIndex = -1; for (int i = 0; i < numSamples / 2; i++) { if (fftData[i] > maxVal) { maxVal = fftData[i]; maxIndex = i; } } // 计算频率的方法根据你的具体应用而定 // 这里只是一个示例,可能需要根据你的传感器和信号特性进行调整 float frequency = maxIndex * getSamplingRate() / (numSamples / 2); return frequency; } float getSamplingRate() { float samplingRate = 0; // 计算采样率的方法根据你的具体应用而定 // 这里只是一个示例,可能需要根据你的传感器和信号特性进行调整 return samplingRate; } ``` 在这个修改后的代码中,我们使用了 Arduino FFT 库的 `arduinoFFT` 类来计算 FFT 转换。在 `setup` 函数中,我们设置了窗函数为 Blackman-Harris,并对 FFT 数据进行重新排序。在 `calculateFFT` 函数中,我们调用 `fft.run` 来计算 FFT 转换数据。 另外,由于 Arduino FFT 库没有提供排序函数,我们使用了一个简单的冒泡排序算法 `getMedian` 函数来计算中位数。 请注意,在 `getFrequency` 和 `getSamplingRate` 函数中,计算频率和采样率的方法是示例代码,你可能需要根据你的传感器和信号特性进行调整。 希望这个修改后的代码能够满足你的需求!如果你还有其他问题,请随时提问。

相关推荐

最新推荐

recommend-type

PID库(Arduino PID Library)使用说明

Arduino PID Library - Version 1.2.1 by Brett Beauregard &lt;br3ttb@gmail.com&gt; brettbeauregard.com This Library is licensed under the MIT License
recommend-type

####这是一篇对python的详细解析

python
recommend-type

zigbee-cluster-library-specification

最新的zigbee-cluster-library-specification说明文档。
recommend-type

管理建模和仿真的文件

管理Boualem Benatallah引用此版本:布阿利姆·贝纳塔拉。管理建模和仿真。约瑟夫-傅立叶大学-格勒诺布尔第一大学,1996年。法语。NNT:电话:00345357HAL ID:电话:00345357https://theses.hal.science/tel-003453572008年12月9日提交HAL是一个多学科的开放存取档案馆,用于存放和传播科学研究论文,无论它们是否被公开。论文可以来自法国或国外的教学和研究机构,也可以来自公共或私人研究中心。L’archive ouverte pluridisciplinaire
recommend-type

MATLAB柱状图在信号处理中的应用:可视化信号特征和频谱分析

![matlab画柱状图](https://img-blog.csdnimg.cn/3f32348f1c9c4481a6f5931993732f97.png) # 1. MATLAB柱状图概述** MATLAB柱状图是一种图形化工具,用于可视化数据中不同类别或组的分布情况。它通过绘制垂直条形来表示每个类别或组中的数据值。柱状图在信号处理中广泛用于可视化信号特征和进行频谱分析。 柱状图的优点在于其简单易懂,能够直观地展示数据分布。在信号处理中,柱状图可以帮助工程师识别信号中的模式、趋势和异常情况,从而为信号分析和处理提供有价值的见解。 # 2. 柱状图在信号处理中的应用 柱状图在信号处理
recommend-type

用Spring boot和vue写一个登录注册界面

好的,这是一个比较常见的技术问题,我可以为您提供一些参考资料和思路。首先,您可以使用Spring Boot作为后端框架来处理用户认证和注册的逻辑,比如使用Spring Security实现用户登录认证。同时,您还需要设计相应的数据模型和数据库表结构来存储用户信息。在前端方面,您可以使用Vue.js作为框架来构建登录注册页面,使用Axios来发起API请求并和后端进行交互。当然,在实现过程中,还需要考虑一些具体细节,比如数据校验、安全性和用户体验等方面。希望这些信息能够帮助到您。
recommend-type

JSBSim Reference Manual

JSBSim参考手册,其中包含JSBSim简介,JSBSim配置文件xml的编写语法,编程手册以及一些应用实例等。其中有部分内容还没有写完,估计有生之年很难看到完整版了,但是内容还是很有参考价值的。
recommend-type

"互动学习:行动中的多样性与论文攻读经历"

多样性她- 事实上SCI NCES你的时间表ECOLEDO C Tora SC和NCESPOUR l’Ingén学习互动,互动学习以行动为中心的强化学习学会互动,互动学习,以行动为中心的强化学习计算机科学博士论文于2021年9月28日在Villeneuve d'Asq公开支持马修·瑟林评审团主席法布里斯·勒菲弗尔阿维尼翁大学教授论文指导奥利维尔·皮耶昆谷歌研究教授:智囊团论文联合主任菲利普·普雷教授,大学。里尔/CRISTAL/因里亚报告员奥利维耶·西格德索邦大学报告员卢多维奇·德诺耶教授,Facebook /索邦大学审查员越南圣迈IMT Atlantic高级讲师邀请弗洛里安·斯特鲁布博士,Deepmind对于那些及时看到自己错误的人...3谢谢你首先,我要感谢我的两位博士生导师Olivier和Philippe。奥利维尔,"站在巨人的肩膀上"这句话对你来说完全有意义了。从科学上讲,你知道在这篇论文的(许多)错误中,你是我可以依
recommend-type

MATLAB柱状图在数据分析中的作用:从可视化到洞察

![MATLAB柱状图在数据分析中的作用:从可视化到洞察](https://img-blog.csdnimg.cn/img_convert/1a36558cefc0339f7836cca7680c0aef.png) # 1. MATLAB柱状图概述** 柱状图是一种广泛用于数据可视化的图表类型,它使用垂直条形来表示数据中不同类别或组别的值。在MATLAB中,柱状图通过`bar`函数创建,该函数接受数据向量或矩阵作为输入,并生成相应的高度条形。 柱状图的优点在于其简单性和易于理解性。它们可以快速有效地传达数据分布和组别之间的比较。此外,MATLAB提供了广泛的定制选项,允许用户调整条形颜色、
recommend-type

命名ACL和拓展ACL标准ACL的具体区别

命名ACL和标准ACL的主要区别在于匹配条件和作用范围。命名ACL可以基于协议、端口和其他条件进行匹配,并可以应用到接口、VLAN和其他范围。而标准ACL只能基于源地址进行匹配,并只能应用到接口。拓展ACL则可以基于源地址、目的地址、协议、端口和其他条件进行匹配,并可以应用到接口、VLAN和其他范围。