文章标题

[code=c] double aa=Nf+1-N;
double bb=1-aa;
double an[3500],bn[3500];
double xn[3500],xa[3500],xb[3500],xx[3500];
double pi=3.14159265;
double sxa,sxb;
for(int i5=0;i5<3500;i5++)
{
int J=2*i5+1;
xn[i5]=cd1d[J];
}
for( int n=N;n<3500;n++)
{
xa[n]=xn[n-N]+bb*(aa-cos((N-1)*2*pi/Nf))*xn[n-1]+aa*aa*xn[n];
xb[n]=-bb*xn[n-1]*sin((N-1)*2*pi/Nf);
sxa=0,sxb=0;
for(int k=1;k