Fast Fourier transform _FFT
2016-08-23
0 0 0
no vote
Other
Earn points
int main(void)
{
FILE *fp_in,*fp_out;
int i,bitn,SLIP;
char asc[2222];
double xxx;
Fp_in=fopen ("sample FFT.txt", "RT");
fp_out=fopen("power.txt","wt"); bitn=10;
num=(1<<bitn);
if (num>LENGTH_MAX) return 2; fgets(asc,1111,fp_in);
sscanf(asc,"%lf ",&DT);
SLIP=0;
for (i=0;i<SLIP;i++) fgets(asc,1111,fp_in);
for (i=0;i<num;i++)
{
fgets(asc,1111,fp_in);
sscanf(asc,"%lf %lf %lf ",&xxx,&My[i],&Mz[i]);
}
for (i=0;i<num;i++)
{
Mz[i]=0.0;
x_r[i]=My[i];
x_i[i]=Mz[i];
}
for (i=0;i<num;i++) {x_r[i]=x_i[i];x_i[i]=0.0;};/ /for (i=0;i<num;i++) printf("%10.2lf %10.2lfn",x_r[i],x_i[i]);/* test
num=8;
for (i=0;i<num;i++)
{
x_r[i]=i;
x_i[i]=0;//i+10;
}
test */ FFT_Base2(x_r,x_i,num);
FFT_Magnitude_Phase(x_r,x_i,x_M,x_P,num);
DF=1.0/
{
FILE *fp_in,*fp_out;
int i,bitn,SLIP;
char asc[2222];
double xxx;
Fp_in=fopen ("sample FFT.txt", "RT");
fp_out=fopen("power.txt","wt"); bitn=10;
num=(1<<bitn);
if (num>LENGTH_MAX) return 2; fgets(asc,1111,fp_in);
sscanf(asc,"%lf ",&DT);
SLIP=0;
for (i=0;i<SLIP;i++) fgets(asc,1111,fp_in);
for (i=0;i<num;i++)
{
fgets(asc,1111,fp_in);
sscanf(asc,"%lf %lf %lf ",&xxx,&My[i],&Mz[i]);
}
for (i=0;i<num;i++)
{
Mz[i]=0.0;
x_r[i]=My[i];
x_i[i]=Mz[i];
}
for (i=0;i<num;i++) {x_r[i]=x_i[i];x_i[i]=0.0;};/ /for (i=0;i<num;i++) printf("%10.2lf %10.2lfn",x_r[i],x_i[i]);/* test
num=8;
for (i=0;i<num;i++)
{
x_r[i]=i;
x_i[i]=0;//i+10;
}
test */ FFT_Base2(x_r,x_i,num);
FFT_Magnitude_Phase(x_r,x_i,x_M,x_P,num);
DF=1.0/
c
fft
变换
快速
傅里叶
Related Source Codes
PClite
0
0
no vote
generic embeded FTP Client
0
0
no vote
UDP Test Demo
0
0
no vote
TMS320F28335 DSP transmits data to EEPROM through
0
0
no vote
TMS320F28335 CAN communication source code
0
0
no vote
No comment