DFT.rar_DFT_Minimum phase DFT_dft c++_离散傅里叶变换

#include <stdio.h> #include <math.h> #include "dft.h" #define N 64 #define L 2*N-1 int main() { int i,j; double h[N],x[N],y[N],yr[N],yi[N],temp[N],ay[N],fai[N]; //x[]---抽样数据cydata.txt h[]---自定义数据脉冲hdata.txt for(i=0;i<N;i++) {yr[i]=0.0; yi[i]=0.0; temp[i]=0.0; } //init(xr,xi,temp); readfile(h); j=cy(30,x); dispx(j,x); //输出文件在抽样后数据cydata.txt disph(j,h); //输出文件在自定义数据脉冲hdata.txt juanji(x,h,y,N,N,L); dispy(j,y); //卷积结果存在ydata.txt dft(y,temp,yr,yi,N,1); dispyr(j,yr); //频谱实部在yrdata.txt dispyi(j,yi); //频谱虚部在yidata.txt zf(ay,yr,yi); dispzf(j,ay); //输出文件在振幅谱aydata.txt xw(fai,yr,yi); dispxw(j,fai); return 0; }