#include "stdio.h"
#include "math.h"
#define dt 0.001   //时间采样间隔
#define NP 2001    //函数长度（样点数）
#define df 0.5    //频率采样间隔
#define fc 20    // 方波的高频
#define Hf 40  //最大频率
#define PI 3.1415926
main()
{ FILE *fp;
  float ft[NP],t,Xr,Xi,Xamp;
  int i,k,Nf;
  Nf=Hf/df+1;
  printf("%d\n",Nf);
  fp=fopen("Sinc_AMP2.xls","w");
  for(i=-NP/2;i<=NP/2;i++)
 { t=i*dt;
  if(t==0)ft[i+NP/2]=1.0;
  else ft[i+NP/2]=sin(2*PI*fc*t)/(2*PI*fc*t);
 //       fprintf(fp,".4f  \t  .4f\n",i*dt,ft[i+NP/2]);
 }
  for(k=-Nf;k<=Nf;k++)
 { Xr=0.0;
  Xi=0.0;
  for(i=-NP/2;i<=NP/2;i++)
   { t=i*dt;
    Xr=Xr+ft[i+NP/2]*cos(2*PI*k*df*t);
    Xi=Xi+ft[i+NP/2]*sin(2*PI*k*df*t);
   }
  Xamp=sqrt(Xr*Xr+Xi*Xi)*dt;
  fprintf(fp,".4f  \t  .4f\n",k*df,Xamp);
    }
  fclose(fp);
}