xiangyang26.cpp

3D-FDTD程序计算平行双阵子

		  fmy.shi+=Jmyf[j-NNN-DDD][k-NNN-DDD].zhengfu*cos(-Jmyb[j-NNN-DDD][k-NNN-DDD].xiangwei+cansh1*(Nx/2-NNN-DDD)*si[loop1]*co[loop2]+cansh2*(j-Ny/2.0)*si[loop1]*si[loop2]+cansh3*(k+0.5-Nz/2.0)*co[loop1]);
		  fmy.xu+=Jmyf[j-NNN-DDD][k-NNN-DDD].zhengfu*sin(-Jmyb[j-NNN-DDD][k-NNN-DDD].xiangwei+cansh1*(Nx/2-NNN-DDD)*si[loop1]*co[loop2]+cansh2*(j-Ny/2.0)*si[loop1]*si[loop2]+cansh3*(k+0.5-Nz/2.0)*co[loop1]);
		}

	   for(j=NNN+DDD;j<=Ny-NNN-DDD;j++)                    //前面z向电磁场
	     for(k=NNN+DDD;k<Nz-NNN-DDD;k++)
		{
		  fy.shi+=Jyb[j-NNN-DDD][k-NNN-DDD].zhengfu*cos(-Jyf[j-NNN-DDD][k-NNN-DDD].xiangwei-cansh1*(Nx/2-NNN-DDD)*si[loop1]*co[loop2]+cansh2*(j+0.5-Ny/2.0)*si[loop1]*si[loop2]+cansh3*(k-Nz/2.0)*co[loop1]);
		  fy.xu+=Jyb[j-NNN-DDD][k-NNN-DDD].zhengfu*sin(-Jyf[j-NNN-DDD][k-NNN-DDD].xiangwei-cansh1*(Nx/2-NNN-DDD)*si[loop1]*co[loop2]+cansh2*(j+0.5-Ny/2.0)*si[loop1]*si[loop2]+cansh3*(k-Nz/2.0)*co[loop1]);
		  fmy.shi-=Jmyb[j-NNN-DDD][k-NNN-DDD].zhengfu*cos(-Jmyf[j-NNN-DDD][k-NNN-DDD].xiangwei-cansh1*(Nx/2-NNN-DDD)*si[loop1]*co[loop2]+cansh2*(j-Ny/2.0)*si[loop1]*si[loop2]+cansh3*(k+0.5-Nz/2.0)*co[loop1]);
		  fmy.xu-=Jmyb[j-NNN-DDD][k-NNN-DDD].zhengfu*sin(-Jmyf[j-NNN-DDD][k-NNN-DDD].xiangwei-cansh1*(Nx/2-NNN-DDD)*si[loop1]*co[loop2]+cansh2*(j-Ny/2.0)*si[loop1]*si[loop2]+cansh3*(k+0.5-Nz/2.0)*co[loop1]);
		}



	   for(i=NNN+DDD;i<Nx-NNN-DDD;i++)                    //右面x向电磁场
	     for(k=NNN+DDD;k<=Nz-NNN-DDD;k++)
		{
		  fz.shi-=Jzr[i-NNN-DDD][k-NNN-DDD].zhengfu*cos(-Jzr[i-NNN-DDD][k-NNN-DDD].xiangwei+cansh1*(i-Nx/2.0)*si[loop1]*co[loop2]+cansh2*(Ny/2-NNN-DDD)*si[loop1]*si[loop2]+cansh3*(k+0.5-Nz/2.0)*co[loop1]);
		  fz.xu-=Jzr[i-NNN-DDD][k-NNN-DDD].zhengfu*sin(-Jzr[i-NNN-DDD][k-NNN-DDD].xiangwei+cansh1*(i-Nx/2.0)*si[loop1]*co[loop2]+cansh2*(Ny/2-NNN-DDD)*si[loop1]*si[loop2]+cansh3*(k+0.5-Nz/2.0)*co[loop1]);
		  fmz.shi+=Jmzr[i-NNN-DDD][k-NNN-DDD].zhengfu*cos(-Jmzr[i-NNN-DDD][k-NNN-DDD].xiangwei+cansh1*(i+0.5-Nx/2.0)*si[loop1]*co[loop2]+cansh2*(Ny/2-NNN-DDD)*si[loop1]*si[loop2]+cansh3*(k-Nz/2.0)*co[loop1]);
		  fmz.xu+=Jmzr[i-NNN-DDD][k-NNN-DDD].zhengfu*sin(-Jmzr[i-NNN-DDD][k-NNN-DDD].xiangwei+cansh1*(i+0.5-Nx/2.0)*si[loop1]*co[loop2]+cansh2*(Ny/2-NNN-DDD)*si[loop1]*si[loop2]+cansh3*(k-Nz/2.0)*co[loop1]);
		 }

	   for(i=NNN+DDD;i<Nx-NNN-DDD;i++)                    //左面x向电磁场
	     for(k=NNN+DDD;k<=Nz-NNN-DDD;k++)
		{
		  fz.shi+=Jzl[i-NNN-DDD][k-NNN-DDD].zhengfu*cos(-Jzl[i-NNN-DDD][k-NNN-DDD].xiangwei+cansh1*(i-Nx/2.0)*si[loop1]*co[loop2]-cansh2*(Ny/2-NNN-DDD)*si[loop1]*si[loop2]+cansh3*(k+0.5-Nz/2.0)*co[loop1]);
		  fz.xu+=Jzl[i-NNN-DDD][k-NNN-DDD].zhengfu*sin(-Jzl[i-NNN-DDD][k-NNN-DDD].xiangwei+cansh1*(i-Nx/2.0)*si[loop1]*co[loop2]-cansh2*(Ny/2-NNN-DDD)*si[loop1]*si[loop2]+cansh3*(k+0.5-Nz/2.0)*co[loop1]);
		  fmz.shi-=Jmzl[i-NNN-DDD][k-NNN-DDD].zhengfu*cos(-Jmzl[i-NNN-DDD][k-NNN-DDD].xiangwei+cansh1*(i+0.5-Nx/2.0)*si[loop1]*co[loop2]-cansh2*(Ny/2-NNN-DDD)*si[loop1]*si[loop2]+cansh3*(k-Nz/2.0)*co[loop1]);
		  fmz.xu-=Jmzl[i-NNN-DDD][k-NNN-DDD].zhengfu*sin(-Jmzl[i-NNN-DDD][k-NNN-DDD].xiangwei+cansh1*(i+0.5-Nx/2.0)*si[loop1]*co[loop2]-cansh2*(Ny/2-NNN-DDD)*si[loop1]*si[loop2]+cansh3*(k-Nz/2.0)*co[loop1]);
		}

	   for(j=NNN+DDD;j<=Ny-NNN-DDD;j++)                    //后面z向电磁场
	     for(k=NNN+DDD;k<Nz-NNN-DDD;k++)
		{
		  fz.shi+=Jzf[j-NNN-DDD][k-NNN-DDD].zhengfu*cos(-Jzb[j-NNN-DDD][k-NNN-DDD].xiangwei+cansh1*(Nx/2-NNN-DDD)*si[loop1]*co[loop2]+cansh2*(j-Ny/2.0)*si[loop1]*si[loop2]+cansh3*(k+0.5-Nz/2.0)*co[loop1]);
		  fz.xu+=Jzf[j-NNN-DDD][k-NNN-DDD].zhengfu*sin(-Jzb[j-NNN-DDD][k-NNN-DDD].xiangwei+cansh1*(Nx/2-NNN-DDD)*si[loop1]*co[loop2]+cansh2*(j-Ny/2.0)*si[loop1]*si[loop2]+cansh3*(k+0.5-Nz/2.0)*co[loop1]);
		  fmz.shi-=Jmzf[j-NNN-DDD][k-NNN-DDD].zhengfu*cos(-Jmzb[j-NNN-DDD][k-NNN-DDD].xiangwei+cansh1*(Nx/2-NNN-DDD)*si[loop1]*co[loop2]+cansh2*(j+0.5-Ny/2.0)*si[loop1]*si[loop2]+cansh3*(k-Nz/2.0)*co[loop1]);
		  fmz.xu-=Jmzf[j-NNN-DDD][k-NNN-DDD].zhengfu*sin(-Jmzb[j-NNN-DDD][k-NNN-DDD].xiangwei+cansh1*(Nx/2-NNN-DDD)*si[loop1]*co[loop2]+cansh2*(j+0.5-Ny/2.0)*si[loop1]*si[loop2]+cansh3*(k-Nz/2.0)*co[loop1]);
		}

	   for(j=NNN+DDD;j<=Ny-NNN-DDD;j++)                    //前面z向电磁场
	     for(k=NNN+DDD;k<Nz-NNN-DDD;k++)
		{
		  fz.shi-=Jzb[j-NNN-DDD][k-NNN-DDD].zhengfu*cos(-Jzf[j-NNN-DDD][k-NNN-DDD].xiangwei-cansh1*(Nx/2-NNN-DDD)*si[loop1]*co[loop2]+cansh2*(j-Ny/2.0)*si[loop1]*si[loop2]+cansh3*(k+0.5-Nz/2.0)*co[loop1]);
		  fz.xu-=Jzb[j-NNN-DDD][k-NNN-DDD].zhengfu*sin(-Jzf[j-NNN-DDD][k-NNN-DDD].xiangwei-cansh1*(Nx/2-NNN-DDD)*si[loop1]*co[loop2]+cansh2*(j-Ny/2.0)*si[loop1]*si[loop2]+cansh3*(k+0.5-Nz/2.0)*co[loop1]);
		  fmz.shi+=Jmzb[j-NNN-DDD][k-NNN-DDD].zhengfu*cos(-Jmzf[j-NNN-DDD][k-NNN-DDD].xiangwei-cansh1*(Nx/2-NNN-DDD)*si[loop1]*co[loop2]+cansh2*(j+0.5-Ny/2.0)*si[loop1]*si[loop2]+cansh3*(k-Nz/2.0)*co[loop1]);
		  fmz.xu+=Jmzb[j-NNN-DDD][k-NNN-DDD].zhengfu*sin(-Jmzf[j-NNN-DDD][k-NNN-DDD].xiangwei-cansh1*(Nx/2-NNN-DDD)*si[loop1]*co[loop2]+cansh2*(j+0.5-Ny/2.0)*si[loop1]*si[loop2]+cansh3*(k-Nz/2.0)*co[loop1]);
		}
		  tt1=120*Pi*(fx.shi*co[loop1]*co[loop2]+fy.shi*co[loop1]*si[loop2]-fz.shi*si[loop1])+(-fmx.shi*si[loop2]+fmy.shi*co[loop2]);
		  tt2=120*Pi*(fx.xu*co[loop1]*co[loop2]+fy.xu*co[loop1]*si[loop2]-fz.xu*si[loop1])+(-fmx.xu*si[loop2]+fmy.xu*co[loop2]);
          fxhs.sita.shi=tt1;
		  fxhs.sita.xu=tt2;
		  tt1=120*Pi*(fx.shi*si[loop2]-fy.shi*co[loop2])+(fmx.shi*co[loop1]*co[loop2]+fmy.shi*co[loop1]*si[loop2]-fmz.shi*si[loop1]);
		  tt2=120*Pi*(fx.xu*si[loop2]-fy.xu*co[loop2])+(fmx.xu*co[loop1]*co[loop2]+fmy.xu*co[loop1]*si[loop2]-fmz.xu*si[loop1]);
		  Efai=sqrt(tt1*tt1+tt2*tt2);
		  fxhs.fai.shi=tt1;
		  fxhs.fai.xu=tt2;
		    fx.shi=0;fx.xu=0;
	        fy.shi=0;fy.xu=0;
	        fmx.shi=0;fmx.xu=0;
	        fmy.shi=0;fmy.xu=0;
			fz.shi=0;fz.xu=0;
			fmz.shi=0;fmz.xu=0;
			printf("%d\n",loop1);
}
  void distill(int M)
  { FILE  *out; 
   int i,j,loop1=0,loop2=0;
   if ((out=fopen("1","a"))==NULL)
   {printf("cannot open outfile\n");
   }
   else  
	   {fprintf(out,"%f   \n",(Hyz[NNN+DDD][NNN+DDD][Nz-NNN-DDD-1]+Hyz[NNN+DDD][NNN+DDD][Nz-NNN-DDD])/2);
       }


     fclose(out);

       for(i=NNN+DDD;i<Nx-NNN-DDD;i++)                    //上面y向电磁场
	   {
	     for(j=NNN+DDD;j<=Ny-NNN-DDD;j++)
		{
		   if((Hyz[i][j][Nz-NNN-DDD-1]+Hyz[i][j][Nz-NNN-DDD])/2.0>Jxu[loop1][loop2].zhengfu)
		  {Jxu[loop1][loop2].zhengfu=(Hyz[i][j][Nz-NNN-DDD-1]+Hyz[i][j][Nz-NNN-DDD])/2.0;Jxu[loop1][loop2].xiangwei=(M-NN+100)/100.0*2.0*pi;}
		   if((Hyz[i][j][NNN+DDD-1]+Hyz[i][j][NNN+DDD])/2.0>Jxd[loop1][loop2].zhengfu)
		  {Jxd[loop1][loop2].zhengfu=(Hyz[i][j][NNN+DDD-1]+Hyz[i][j][NNN+DDD])/2.0;Jxd[loop1][loop2].xiangwei=(M-NN+100)/100.0*2.0*pi;}

           if(Eyz[i][j][Nz-NNN-DDD]>Jmxu[loop1][loop2].zhengfu)
		  {Jmxu[loop1][loop2].zhengfu=Eyz[i][j][Nz-NNN-DDD];Jmxu[loop1][loop2].xiangwei=(M-NN+100)/100.0*2.0*pi;}
           if(Eyz[i][j][NNN+DDD]>Jmxd[loop1][loop2].zhengfu)
		  {Jmxd[loop1][loop2].zhengfu=Eyz[i][j][NNN+DDD];Jmxd[loop1][loop2].xiangwei=(M-NN+100)/100.0*2.0*pi;}

           if((Hzx[i][Ny-NNN-DDD-1][j]+Hzx[i][Ny-NNN-DDD][j])/2.0>Jxr[loop1][loop2].zhengfu)
		  {Jxr[loop1][loop2].zhengfu=(Hzx[i][Ny-NNN-DDD-1][j]+Hzx[i][Ny-NNN-DDD][j])/2.0;Jxr[loop1][loop2].xiangwei=(M-NN+100)/100.0*2.0*pi;}
           if(Ezx[i][Ny-NNN-DDD][j]>Jmxr[loop1][loop2].zhengfu)
		  {Jmxr[loop1][loop2].zhengfu=Ezx[i][Ny-NNN-DDD][j];Jmxr[loop1][loop2].xiangwei=(M-NN+100)/100.0*2.0*pi;}

		   if((Hzx[i][NNN+DDD-1][j]+Hzx[i][NNN+DDD][j])/2.0>Jxl[loop1][loop2].zhengfu)
		  {Jxl[loop1][loop2].zhengfu=(Hzx[i][NNN+DDD-1][j]+Hzx[i][NNN+DDD][j])/2.0;Jxl[loop1][loop2].xiangwei=(M-NN+100)/100.0*2.0*pi;}
           if(Ezx[i][NNN+DDD][j]>Jmxl[loop1][loop2].zhengfu)
		  {Jmxl[loop1][loop2].zhengfu=Ezx[i][NNN+DDD][j];Jmxl[loop1][loop2].xiangwei=(M-NN+100)/100.0*2.0*pi;}
           
		   if((Hxy[i][Ny-NNN-DDD-1][j]+Hxy[i][Ny-NNN-DDD][j])/2.0>Jzr[loop1][loop2].zhengfu)
		  {Jzr[loop1][loop2].zhengfu=(Hxy[i][Ny-NNN-DDD-1][j]+Hxy[i][Ny-NNN-DDD][j])/2.0;Jzr[loop1][loop2].xiangwei=(M-NN+100)/100.0*2.0*pi;}
           if(Exy[i][Ny-NNN-DDD][j]>Jmzr[loop1][loop2].zhengfu)
		  {Jmzr[loop1][loop2].zhengfu=Exy[i][Ny-NNN-DDD][j];Jmzr[loop1][loop2].xiangwei=(M-NN+100)/100.0*2.0*pi;}

		   if((Hxy[i][NNN+DDD-1][j]+Hxy[i][NNN+DDD][j])/2.0>Jzl[loop1][loop2].zhengfu)
		  {Jzl[loop1][loop2].zhengfu=(Hxy[i][NNN+DDD-1][j]+Hxy[i][NNN+DDD][j])/2.0;Jzl[loop1][loop2].xiangwei=(M-NN+100)/100.0*2.0*pi;}
		   if(Exy[i][NNN+DDD][j]>Jmzl[loop1][loop2].zhengfu)
		  {Jmzl[loop1][loop2].zhengfu=Exy[i][NNN+DDD][j];Jmzl[loop1][loop2].xiangwei=(M-NN+100)/100.0*2.0*pi;}


		  loop2++;
		}
		 loop1++;loop2=0;
       }
        loop1=0;loop2=0;

	    for(i=NNN+DDD;i<=Nx-NNN-DDD;i++)                    //上面y向电磁场
	   {
	     for(j=NNN+DDD;j<Ny-NNN-DDD;j++)
		{
			 if((Hxy[i][j][Nz-NNN-DDD-1]+Hxy[i][j][Nz-NNN-DDD])/2.0>Jyu[loop1][loop2].zhengfu)
			 {Jyu[loop1][loop2].zhengfu=(Hxy[i][j][Nz-NNN-DDD-1]+Hxy[i][j][Nz-NNN-DDD])/2.0;Jyu[loop1][loop2].xiangwei=(M-NN+100)/100.0*2.0*pi;}
			 if(Exy[i][j][Nz-NNN-DDD]>Jmyu[loop1][loop2].zhengfu)
			 {Jmyu[loop1][loop2].zhengfu=Exy[i][j][Nz-NNN-DDD];Jmyu[loop1][loop2].xiangwei=(M-NN+100)/100.0*2.0*pi;}

			 if((Hxy[i][j][NNN+DDD-1]+Hxy[i][j][NNN+DDD])/2.0>Jyd[loop1][loop2].zhengfu)
			 {Jyd[loop1][loop2].zhengfu=(Hxy[i][j][NNN+DDD-1]+Hxy[i][j][NNN+DDD])/2.0;Jyd[loop1][loop2].xiangwei=(M-NN+100)/100.0*2.0*pi;}
			 if(Exy[i][j][NNN+DDD]>Jmyd[loop1][loop2].zhengfu)
			 {Jmyd[loop1][loop2].zhengfu=Exy[i][j][NNN+DDD];Jmyd[loop1][loop2].xiangwei=(M-NN+100)/100.0*2.0*pi;}

			 if((Hzx[Nx-NNN-DDD-1][i][j]+Hzx[Nx-NNN-DDD][i][j])/2.0>Jyf[loop1][loop2].zhengfu)
			 {Jyf[loop1][loop2].zhengfu=(Hzx[Nx-NNN-DDD-1][i][j]+Hzx[Nx-NNN-DDD][i][j])/2.0;Jyf[loop1][loop2].xiangwei=(M-NN+100)/100.0*2.0*pi;}
			 if(Ezx[Nx-NNN-DDD][i][j]>Jmyf[loop1][loop2].zhengfu)
			 {Jmyf[loop1][loop2].zhengfu=Ezx[Nx-NNN-DDD][i][j];Jmyf[loop1][loop2].xiangwei=(M-NN+100)/100.0*2.0*pi;}

			 if((Hzx[NNN+DDD-1][i][j]+Hzx[NNN+DDD][i][j])/2.0>Jyb[loop1][loop2].zhengfu)
			 {Jyb[loop1][loop2].zhengfu=(Hzx[NNN+DDD-1][i][j]+Hzx[NNN+DDD][i][j])/2.0;Jyb[loop1][loop2].xiangwei=(M-NN+100)/100.0*2.0*pi;}
			 if(Ezx[NNN+DDD][i][j]>Jmyb[loop1][loop2].zhengfu)
			 {Jmyb[loop1][loop2].zhengfu=Ezx[NNN+DDD][i][j];Jmyb[loop1][loop2].xiangwei=(M-NN+100)/100.0*2.0*pi;}

			 if((Hyz[Nx-NNN-DDD-1][i][j]+Hyz[Nx-NNN-DDD][i][j])/2.0>Jzf[loop1][loop2].zhengfu)
			 {Jzf[loop1][loop2].zhengfu=(Hyz[Nx-NNN-DDD-1][i][j]+Hyz[Nx-NNN-DDD][i][j])/2.0;Jzf[loop1][loop2].xiangwei=(M-NN+100)/100.0*2.0*pi;}
			 if(Eyz[Nx-NNN-DDD][i][j]>Jmzf[loop1][loop2].zhengfu)
			 {Jmzf[loop1][loop2].zhengfu=Eyz[Nx-NNN-DDD][i][j];Jmzf[loop1][loop2].xiangwei=(M-NN+100)/100.0*2.0*pi;}

			 if((Hyz[NNN+DDD-1][i][j]+Hyz[NNN+DDD][i][j])/2.0>Jzb[loop1][loop2].zhengfu)
			 {Jzb[loop1][loop2].zhengfu=(Hyz[NNN+DDD-1][i][j]+Hyz[NNN+DDD][i][j])/2.0;Jzb[loop1][loop2].xiangwei=(M-NN+100)/100.0*2.0*pi;}
			 if(Eyz[NNN+DDD][i][j]>Jmzb[loop1][loop2].zhengfu)
			 {Jmzb[loop1][loop2].zhengfu=Eyz[NNN+DDD][i][j];Jmzb[loop1][loop2].xiangwei=(M-NN+100)/100.0*2.0*pi;}
      
		  loop2++;
		}
		 loop1++;loop2=0;
       }
        loop1=0;loop2=0;
  
  }





void zukang()
{
  int loop1;
  float tt,df=0;

  for(loop1=1;loop1<300;loop1++)
  {  df+=(1/dt)/Np;
	 pr[loop1-1]=chu(pp2[loop1-1],pp1[loop1-1]);
     zk[loop1-1].shi=pr[loop1-1].xu/(2*Pi*loop1/(dt*NN)*dt);
	 zk[loop1-1].xu=1/(2*Pi*loop1/(dt*Np)*ee*er*dx*dy/dz)-pr[loop1-1].shi/(2*Pi*loop1/(dt*Np)*dt);
	
  }

}
void zhubobi()
{
    int loop1;float tt1,tt2;
   for(loop1=1;loop1<300;loop1++)
  {  
     tt1=sqrt((zk[loop1-1].shi+50)*(zk[loop1-1].shi+50)+zk[loop1-1].xu*zk[loop1-1].xu);

	 tt2=sqrt((zk[loop1-1].shi-50)*(zk[loop1-1].shi-50)+zk[loop1-1].xu*zk[loop1-1].xu);
	 rou[loop1-1]=(tt1+tt2)/(tt1-tt2);
  }
}

/*
void zukang(int M)
{
   int i;
   float can;
   struct bb pr;
   can=f*dt*NN-floor(f*dt*NN);
   printf("%f\n",can);
   if(can>=0.5)
	   i=floor(f*dt*NN)+1;
   else
	   i=floor(f*dt*NN);

    printf("\n%d\n",i);
     pr=chu(pp2[i],pp1[i]);
	 zk[M].shi=pr.xu/(2*Pi*f*dt);
	 zk[M].xu=1/(2*Pi*f*ee*er*dx*dy/dz)-pr.shi/(2*Pi*f*dt);
	 printf("输入阻抗:%f+(%f)*i",zk[M].shi,zk[M].xu);
     if(M==12)
	   printf("狂爽!!!");
}
*/
 struct bb chu(struct bb p,struct bb q)
 {
     struct bb t1;
     t1.shi=(p.shi*q.shi+p.xu*q.xu)/(q.shi*q.shi+q.xu*q.xu);
     t1.xu=(p.xu*q.shi-p.shi*q.xu)/(q.shi*q.shi+q.xu*q.xu);
	 return t1;
 }
/*
void zhubobi(int M)
{
	float tt1,tt2;
    tt1=sqrt((zk[M].shi+50)*(zk[M].shi+50)+zk[M].xu*zk[M].xu);

	 tt2=sqrt((zk[M].shi-50)*(zk[M].shi-50)+zk[M].xu*zk[M].xu);
	 rou[M]=(tt1+tt2)/(tt1-tt2);

}
*/
 void ffft(float *p1,float *p2)
 {  int loop,loop1;
     for(loop=0;loop<NN;loop++)
	 {   for(loop1=1;loop1<Np;loop1++)      //NN点fft变换
	   {
	     pp2[loop1-1].shi+=*(p1+loop)*cos(2*Pi/Np*(loop-1)*(loop1-1));
		 pp2[loop1-1].xu-=*(p1+loop)*sin(2*Pi/Np*(loop-1)*(loop1-1));
	   
	   }
	   printf("%d\n",loop);
	 }
	
       for(loop=0;loop<NN;loop++)
	   {	   for(loop1=1;loop1<Np;loop1++)      //NN点fft变换
	   {
	     pp1[loop1-1].shi+=*(p2+loop)*cos(2*Pi/Np*(loop-1)*(loop1-1));
		 pp1[loop1-1].xu-=*(p2+loop)*sin(2*Pi/Np*(loop-1)*(loop1-1));
	   
	   }
	   printf("%d\n",loop);
	   }
 }
/*
void waitui(int M)
{
   int loop1,loop2;
   for (loop1=0;loop1<=1000;loop1++)
	   for (loop2=0;loop2<=DD/2;loop2++)
	   {fxhs[loop1].shi-=current1[loop2].zhengfu*dz/(2*(c/f))*120*Pi*sin((85+loop1*0.01)/180.0*Pi)*sin(current1[loop2].xiangwei+2*Pi/(c/f)*(D/2+loop2+0.5)*dz*cos((85+loop1*0.01)/180.0*Pi)-2*Pi/(c/f)*100)*1.0/100;
	    fxhs[loop1].xu+=current1[loop2].zhengfu*dz/(2*(c/f))*120*Pi*sin((85+loop1*0.01)/180.0*Pi)*cos(current1[loop2].xiangwei+2*Pi/(c/f)*(D/2+loop2+0.5)*dz*cos((85+loop1*0.01)/180.0*Pi)-2*Pi/(c/f)*100)*1.0/100;
        fxhs[loop1].shi-=current2[loop2].zhengfu*dz/(2*(c/f))*120*Pi*sin((85+loop1*0.01)/180.0*Pi)*sin(current2[loop2].xiangwei-2*Pi/(c/f)*(D/2+loop2+0.5)*dz*cos((85+loop1*0.01)/180.0*Pi)-2*Pi/(c/f)*100)*1.0/100;
		fxhs[loop1].xu+=current2[loop2].zhengfu*dz/(2*(c/f))*120*Pi*sin((85+loop1*0.01)/180.0*Pi)*cos(current2[loop2].xiangwei-2*Pi/(c/f)*(D/2+loop2+0.5)*dz*cos((85+loop1*0.01)/180.0*Pi)-2*Pi/(c/f)*100)*1.0/100;
       }

}
*/