c语言与matlab生成FFT结果对比

马剑楠

于 2019-03-25 08:54:52 发布

阅读量1.8k

点赞数 1

CC 4.0 BY-SA版权

本文链接：https://blog.youkuaiyun.com/Jas817ma/article/details/88649836

本文探讨了C语言和MATLAB实现快速傅里叶变换(FFT)的过程。C语言实现包括了复数类型定义、FFT算法、初始化变换核、变址计算等步骤。MATLAB中，FFT函数直接调用，适用于长度为2的整数次幂的序列。通过对比，发现当输入序列长度与MATLAB指定的N不匹配时，MATLAB会通过补零或截断处理，导致结果差异。

摘要生成于 C知道，由 DeepSeek-R1 满血版支持，前往体验 >

首先用C语言实现FFT的输出
一、程序内容
#include <stdio.h>
#include <math.h>
#include <stdlib.h>
#define N 1000
//定义复数类型
typedef struct
{
double real;
double img;
}complex;

void fft(); // 快速傅里叶变换
void initW(); // 初始化变换核
void change(); // 变址

void add(complex ,complex ,complex *); // 复数加法
void mul(complex ,complex ,complex *); // 复数乘法
void sub(complex ,complex ,complex *); // 复数减法
void divi(complex ,complex ,complex *); // 复数除法

void output(); // 输出快速傅里叶变换的结果
complex x[N], *W; // 输入序列,变换核
int size_x=0; // 输入序列的大小，在本程序中仅限2的次幂
double PI; // 圆周率

int main()
{
int i; // 输出结果
system(“cls”);
PI=atan(1)*4;
printf(“Please input the size of x:\n”); // 输入序列的大小
scanf("%d",&size_x);
printf(“Please input the data in x[N]:\n”); // 输入序列的实部和虚部

printf("Enter the sequence\n",i);
for (i = 0; i < size_x; i++)
{
	scanf("%lf%lf",&am

最低0.47元/天解锁文章

200万优质内容无限畅学