国内做GPT做的最好的是科大讯飞，之前的版本跟傻子一样，最近又使用，发现提升了不少，基本上接近了chatGPT,大家可以注册试用。网页版如下，只需要通过手机验证码即可注册登录使用。
https://xinghuo.xfyun.cn/
下面跟大家看看效果
#include <complex.h>
#include <math.h>
#include <stdio.h>
#ifndef M_PI
#define M_PI 3.14159265358979323846
#endif
void fft(complex double *a, int n) {
if (n <= 1) return;
complex double a0[n / 2], a1[n / 2];
for (int i = 0; i < n / 2; i++) {
a0[i] = a[i * 2];
a1[i] = a[i * 2 + 1];
}
fft(a0, n / 2);
fft(a1, n / 2);
for (int i = 0; i < n / 2; i++) {
complex double t = cexp(-2.0 * I * M_PI * i / n) * a1[i];
a[i] = a0[i] + t;
a[i + n / 2] = a0[i] - t;
}
}
int main() {
int n = 8;
complex double a[] = {1, 1, 1, 1, 0, 0, 0, 0};
fft(a, n);
for (int i = 0; i < n; i++) {
printf("%f%+fi
", creal(a[i]), cimag(a[i]));
}
return 0;
}
（内容由讯飞星火AI生成）
继续提问
继续提问：如果使用Arduino来实现FFT，该如何编写代码？
讯飞星火：
使用Arduino实现FFT需要对C语言和数字信号处理有一定的了解。以下是一个简化的Arduino FFT实现示例：
#include <complex.h>
#include <math.h>
#ifndef M_PI
#define M_PI 3.14159265358979323846
#endif
void fft(complex double *a, int n) {
if (n <= 1) return;
complex double a0[n / 2], a1[n / 2];
for (int i = 0; i < n / 2; i++) {
a0[i] = a[i * 2];
a1[i] = a[i * 2 + 1];
}
fft(a0, n / 2);
fft(a1, n / 2);
for (int i = 0; i < n / 2; i++) {
complex double t = cexp(-2.0 * I * M_PI * i / n) * a1[i];
a[i] = a0[i] + t;
a[i + n / 2] = a0[i] - t;
}
}
void setup() {
complex double a[] = {1, 1, 1, 1, 0, 0, 0, 0};
int n = sizeof(a) / sizeof(a[0]);
fft(a, n);
for (int i = 0; i < n; i++) {
Serial.print(creal(a[i]));
Serial.print(" + ");
Serial.print(cimag(a[i]));
Serial.println("i");
}
}
void loop() {
}
这个实现中，我们使用了C语言中的复数类型complex double来表示复数，其中实部和虚部都是浮点数。我们还使用了cexp函数来计算复数的指数，以及creal和cimag函数来提取复数的实部和虚部。
在setup函数中，我们初始化了一个长度为8的复数数组a，并调用fft函数计算其DFT。然后，我们将结果输出到串口。
在loop函数中，我们没有执行任何操作，因为FFT计算只需要在程序开始时执行一次。如果你需要在循环中不断更新FFT结果，可以将FFT计算部分的代码移动到loop函数中，并在每次循环时重新计算FFT。
（内容由讯飞星火AI生成）
大家觉得这个GPT工具如何呢？想不想赶紧试用一下。
| 
楼主
标题置顶
标题高亮
