cbb.h

采用相互垂直的两个振子构成圆极化天线的FDTD程序．可以计算其圆极化波瓣宽度等天线指标

		 }
		  else
		  {sgmz=0;sgmmz=0;}


		 }
     if(loop1==Nx/2&&loop2==Ny/2&&((loop3>=Nz/2-DD1/2-1&&loop3<Nz/2-1)||(loop3>Nz/2&&loop3<=Nz/2+DD1/2)))
	 {       Hyz[loop1][loop2][loop3]=((2*u-sgmmz*dt)/(2*u+sgmmz*dt))*Hyz[loop1][loop2][loop3]-(2*dt)/(2*u+sgmmz*dt)*1/((dx-r0)*dz)*((Exy[loop1][loop2][loop3+1]*(dx-r0)+Exz[loop1][loop2][loop3+1]*(dx-r0)-Exy[loop1][loop2][loop3]*(dx-r0)-Exz[loop1][loop2][loop3]*(dx-r0)));
             Hyx[loop1][loop2][loop3]=((2*u-sgmmx*dt)/(2*u+sgmmx*dt))*Hyx[loop1][loop2][loop3]+(2*dt)/(2*u+sgmmx*dt)*1/((dx-r0)*dz)*((Ezx[loop1+1][loop2][loop3]*dz+Ezy[loop1+1][loop2][loop3]*dz));
	 }
	 else if(loop1==Nx/2-1&&loop2==Ny/2&&((loop3>=Nz/2-DD1/2-1&&loop3<Nz/2-1)||(loop3>Nz/2&&loop3<=Nz/2+DD1/2)))
	 {       Hyz[loop1][loop2][loop3]=((2*u-sgmmz*dt)/(2*u+sgmmz*dt))*Hyz[loop1][loop2][loop3]-(2*dt)/(2*u+sgmmz*dt)*1/((dx-r0)*dz)*((Exy[loop1][loop2][loop3+1]*(dx-r0)+Exz[loop1][loop2][loop3+1]*(dx-r0)-Exy[loop1][loop2][loop3]*(dx-r0)-Exz[loop1][loop2][loop3]*(dx-r0)));
             Hyx[loop1][loop2][loop3]=((2*u-sgmmx*dt)/(2*u+sgmmx*dt))*Hyx[loop1][loop2][loop3]-(2*dt)/(2*u+sgmmx*dt)*1/((dx-r0)*dz)*((Ezx[loop1][loop2][loop3]*dz+Ezy[loop1][loop2][loop3]*dz));
	 }
	 else if(((loop1>=Nx/2-DD2/2-1&&loop1<Nx/2-1)||(loop1>Nx/2&&loop1<=Nx/2+DD2/2))&&loop2==Ny/2+movedis&&loop3==Nz/2)
	 {
	         Hyz[loop1][loop2][loop3]=((2*u-sgmmz*dt)/(2*u+sgmmz*dt))*Hyz[loop1][loop2][loop3]+(2*dt)/(2*u+sgmmz*dt)*1/(dx*(dz-r1))*((Ezx[loop1+1][loop2][loop3]*(dz-r1)+Ezy[loop1+1][loop2][loop3]*(dz-r1)-Ezx[loop1][loop2][loop3]*(dz-r1)-Ezy[loop1][loop2][loop3]*(dz-r1)));
             Hyx[loop1][loop2][loop3]=((2*u-sgmmx*dt)/(2*u+sgmmx*dt))*Hyx[loop1][loop2][loop3]-(2*dt)/(2*u+sgmmx*dt)*1/(dx*(dz-r1))*((Exy[loop1][loop2][loop3+1]*dx+Exz[loop1][loop2][loop3+1]*dx));
	 }
	 else if(((loop1>=Nx/2-DD2/2-1&&loop1<Nx/2-1)||(loop1>Nx/2&&loop1<=Nx/2+DD2/2))&&loop2==Ny/2+movedis&&loop3==Nz/2-1)
	 {
	         Hyz[loop1][loop2][loop3]=((2*u-sgmmz*dt)/(2*u+sgmmz*dt))*Hyz[loop1][loop2][loop3]+(2*dt)/(2*u+sgmmz*dt)*1/(dx*(dz-r1))*((Ezx[loop1+1][loop2][loop3]*(dz-r1)+Ezy[loop1+1][loop2][loop3]*(dz-r1)-Ezx[loop1][loop2][loop3]*(dz-r1)-Ezy[loop1][loop2][loop3]*(dz-r1)));
             Hyx[loop1][loop2][loop3]=((2*u-sgmmx*dt)/(2*u+sgmmx*dt))*Hyx[loop1][loop2][loop3]+(2*dt)/(2*u+sgmmx*dt)*1/(dx*(dz-r1))*((Exy[loop1][loop2][loop3]*dx+Exz[loop1][loop2][loop3]*dx));
	 }
     else
	 {   
          Hyz[loop1][loop2][loop3]=((2*u-sgmmz*dt)/(2*u+sgmmz*dt))*Hyz[loop1][loop2][loop3]-(2*dt)/(2*u+sgmmz*dt)*((Exy[loop1][loop2][loop3+1]+Exz[loop1][loop2][loop3+1]-Exy[loop1][loop2][loop3]-Exz[loop1][loop2][loop3])/dz);
          Hyx[loop1][loop2][loop3]=((2*u-sgmmx*dt)/(2*u+sgmmx*dt))*Hyx[loop1][loop2][loop3]+(2*dt)/(2*u+sgmmx*dt)*((Ezx[loop1+1][loop2][loop3]+Ezy[loop1+1][loop2][loop3]-Ezx[loop1][loop2][loop3]-Ezy[loop1][loop2][loop3])/dx);
	 }
		  if(loop1>NNN-1&&loop1<Nx-NNN&&loop2>NNN-1&&loop2<Ny-NNN&&loop3>NNN-1&&loop3<Nz-NNN)
		  {Hyz[loop1][loop2][loop3]+=Hyx[loop1][loop2][loop3];
		   Hyx[loop1][loop2][loop3]=0;}

	   } 
	  }
	  }
/****************************************************************************************************************************************************************************************************************************************/       
   /*z向磁场*/
	  for(loop1=0;loop1<Nx;loop1++)
	  {for(loop2=0;loop2<Ny;loop2++)
	  {for(loop3=0;loop3<Nz+1;loop3++)
	 {   if(loop1>NNN-1&&loop1<Nx-NNN&&loop2>NNN-1&&loop2<Ny-NNN&&loop3>NNN-1&&loop3<Nz-NNN)
		{sgmx=0;sgmy=0;sgmz=0;
		 sgmmx=0;sgmmy=0;sgmmz=0;
		}
	     else
		 { 
		  if(loop1>=0&&loop1<=NNN-1)
		 {  sgmx=sgmmax*(NNN-loop1-1)/NNN*(NNN-loop1-1)/NNN;
		    sgmmx=sgmx*142122.3034;
		 }
		  else if(loop1>=Nx-NNN&&loop1<Nx)
		 {  sgmx=sgmmax*(NNN+loop1-Nx)/NNN*(NNN+loop1-Nx)/NNN;
		    sgmmx=sgmx*142122.3034;
		 }
		  else
		  {sgmx=0;sgmmx=0;}


           if(loop2>=0&&loop2<=NNN-1)
		 {  sgmy=sgmmax*(NNN-loop2-1)/NNN*(NNN-loop2-1)/NNN;
		    sgmmy=sgmy*142122.3034;
		 }
		  else if(loop2>=Ny-NNN&&loop2<Ny)
		 {  sgmy=sgmmax*(NNN+loop2-Ny)/NNN*(NNN+loop2-Ny)/NNN;
		    sgmmy=sgmy*142122.3034;
		 }
		  else
		  {sgmy=0;sgmmy=0;}


		   if(loop3>=0&&loop3<=NNN-1)
		 {  sgmz=sgmmax*(NNN-loop3-1)/NNN*(NNN-loop3-1)/NNN;
		    sgmmz=sgmz*142122.3034;
		 }
		  else if(loop3>=Nz-NNN&&loop3<Nz)
		 {  sgmz=sgmmax*(NNN+loop3-Nz)/NNN*(NNN+loop3-Nz)/NNN;
		    sgmmz=sgmz*142122.3034;
		 }
		  else
		  {sgmz=0;sgmmz=0;}


		 }
           if(((loop1>=Nx/2-DD2/2-1&&loop1<Nx/2-1)||(loop1>Nx/2&&loop1<=Nx/2+DD2/2))&&loop2==Ny/2+movedis&&loop3==Nz/2)
			{
	               Hzx[loop1][loop2][loop3]=((2*u-sgmmx*dt)/(2*u+sgmmx*dt))*Hzx[loop1][loop2][loop3]-(2*dt)/(2*u+sgmmx*dt)*1/(dx*(dy-r1))*((Eyz[loop1+1][loop2][loop3]*(dy-r1)+Eyx[loop1+1][loop2][loop3]*(dy-r1)-Eyz[loop1][loop2][loop3]*(dy-r1)-Eyx[loop1][loop2][loop3]*(dy-r1)));
                   Hzy[loop1][loop2][loop3]=((2*u-sgmmy*dt)/(2*u+sgmmy*dt))*Hzy[loop1][loop2][loop3]+(2*dt)/(2*u+sgmmy*dt)*1/(dx*(dy-r1))*((Exy[loop1][loop2+1][loop3]*dx+Exz[loop1][loop2+1][loop3]*dx));
			}
	        else if(((loop1>=Nx/2-DD2/2-1&&loop1<Nx/2-1)||(loop1>Nx/2&&loop1<=Nx/2+DD2/2))&&loop2==Ny/2-1+movedis&&loop3==Nz/2)
			{
	               Hzx[loop1][loop2][loop3]=((2*u-sgmmx*dt)/(2*u+sgmmx*dt))*Hzx[loop1][loop2][loop3]-(2*dt)/(2*u+sgmmx*dt)*1/(dx*(dy-r1))*((Eyz[loop1+1][loop2][loop3]*(dy-r1)+Eyx[loop1+1][loop2][loop3]*(dy-r1)-Eyz[loop1][loop2][loop3]*(dy-r1)-Eyx[loop1][loop2][loop3]*(dy-r1)));
                   Hzy[loop1][loop2][loop3]=((2*u-sgmmy*dt)/(2*u+sgmmy*dt))*Hzy[loop1][loop2][loop3]-(2*dt)/(2*u+sgmmy*dt)*1/(dx*(dy-r1))*((Exy[loop1][loop2][loop3]*dx+Exz[loop1][loop2][loop3]*dx));
			}
            else  if(loop1==Nx/2&&loop2==Ny/2&&((loop3>=Nz/2-DD1/2-1&&loop3<Nz/2)||(loop3>Nz/2&&loop3<=Nz/2+DD1/2+1)))
			{
	               Hzx[loop1][loop2][loop3]=((2*u-sgmmx*dt)/(2*u+sgmmx*dt))*Hzx[loop1][loop2][loop3]-(2*dt)/(2*u+sgmmx*dt)*1.0/(dx*dy-Pi*r0*r0/4)*(Eyz[loop1+1][loop2][loop3]*dy+Eyx[loop1+1][loop2][loop3]*dy-Eyz[loop1][loop2][loop3]*(dy-r0)-Eyx[loop1][loop2][loop3]*(dy-r0));
                   Hzy[loop1][loop2][loop3]=((2*u-sgmmy*dt)/(2*u+sgmmy*dt))*Hzy[loop1][loop2][loop3]+(2*dt)/(2*u+sgmmy*dt)*1.0/(dx*dy-Pi*r0*r0/4)*(Exy[loop1][loop2+1][loop3]*dx+Exz[loop1][loop2+1][loop3]*dx-Exy[loop1][loop2][loop3]*(dx-r0)-Exz[loop1][loop2][loop3]*(dx-r0));	       
			}
		    else  if(loop1==Nx/2&&loop2==Ny/2-1&&((loop3>=Nz/2-DD1/2-1&&loop3<Nz/2)||(loop3>Nz/2&&loop3<=Nz/2+DD1/2+1)))
			{
	                Hzx[loop1][loop2][loop3]=((2*u-sgmmx*dt)/(2*u+sgmmx*dt))*Hzx[loop1][loop2][loop3]-(2*dt)/(2*u+sgmmx*dt)*1.0/(dx*dy-Pi*r0*r0/4)*(Eyz[loop1+1][loop2][loop3]*dy+Eyx[loop1+1][loop2][loop3]*dy-Eyz[loop1][loop2][loop3]*(dy-r0)-Eyx[loop1][loop2][loop3]*(dy-r0));
                    Hzy[loop1][loop2][loop3]=((2*u-sgmmy*dt)/(2*u+sgmmy*dt))*Hzy[loop1][loop2][loop3]+(2*dt)/(2*u+sgmmy*dt)*1.0/(dx*dy-Pi*r0*r0/4)*(Exy[loop1][loop2+1][loop3]*(dx-r0)+Exz[loop1][loop2+1][loop3]*(dx-r0)-Exy[loop1][loop2][loop3]*dx-Exz[loop1][loop2][loop3]*dx);
			}
            else  if(loop1==Nx/2-1&&loop2==Ny/2&&((loop3>=Nz/2-DD1/2-1&&loop3<Nz/2)||(loop3>Nz/2&&loop3<=Nz/2+DD1/2+1)))
		   {
	                Hzx[loop1][loop2][loop3]=((2*u-sgmmx*dt)/(2*u+sgmmx*dt))*Hzx[loop1][loop2][loop3]-(2*dt)/(2*u+sgmmx*dt)*1.0/(dx*dy-Pi*r0*r0/4)*(Eyz[loop1+1][loop2][loop3]*(dy-r0)+Eyx[loop1+1][loop2][loop3]*(dy-r0)-Eyz[loop1][loop2][loop3]*dy-Eyx[loop1][loop2][loop3]*dy);
                    Hzy[loop1][loop2][loop3]=((2*u-sgmmy*dt)/(2*u+sgmmy*dt))*Hzy[loop1][loop2][loop3]+(2*dt)/(2*u+sgmmy*dt)*1.0/(dx*dy-Pi*r0*r0/4)*(Exy[loop1][loop2+1][loop3]*dx+Exz[loop1][loop2+1][loop3]*dx-Exy[loop1][loop2][loop3]*(dx-r0)-Exz[loop1][loop2][loop3]*(dx-r0));
    	   }
            else  if(loop1==Nx/2-1&&loop2==Ny/2-1&&((loop3>=Nz/2-DD1/2-1&&loop3<Nz/2)||(loop3>Nz/2&&loop3<=Nz/2+DD1/2+1)))
		   {
	                Hzx[loop1][loop2][loop3]=((2*u-sgmmx*dt)/(2*u+sgmmx*dt))*Hzx[loop1][loop2][loop3]-(2*dt)/(2*u+sgmmx*dt)*1.0/(dx*dy-Pi*r0*r0/4)*(Eyz[loop1+1][loop2][loop3]*(dy-r0)+Eyx[loop1+1][loop2][loop3]*(dy-r0)-Eyz[loop1][loop2][loop3]*dy-Eyx[loop1][loop2][loop3]*dy);
                    Hzy[loop1][loop2][loop3]=((2*u-sgmmy*dt)/(2*u+sgmmy*dt))*Hzy[loop1][loop2][loop3]+(2*dt)/(2*u+sgmmy*dt)*1.0/(dx*dy-Pi*r0*r0/4)*(Exy[loop1][loop2+1][loop3]*(dx-r0)+Exz[loop1][loop2+1][loop3]*(dx-r0)-Exy[loop1][loop2][loop3]*dx-Exz[loop1][loop2][loop3]*dx);
    	   }

 
           else
		   { 
          Hzx[loop1][loop2][loop3]=((2*u-sgmmx*dt)/(2*u+sgmmx*dt))*Hzx[loop1][loop2][loop3]-(2*dt)/(2*u+sgmmx*dt)*((Eyz[loop1+1][loop2][loop3]+Eyx[loop1+1][loop2][loop3]-Eyz[loop1][loop2][loop3]-Eyx[loop1][loop2][loop3])/dx);
          Hzy[loop1][loop2][loop3]=((2*u-sgmmy*dt)/(2*u+sgmmy*dt))*Hzy[loop1][loop2][loop3]+(2*dt)/(2*u+sgmmy*dt)*((Exy[loop1][loop2+1][loop3]+Exz[loop1][loop2+1][loop3]-Exy[loop1][loop2][loop3]-Exz[loop1][loop2][loop3])/dy);
		   }
		  if(loop1>NNN-1&&loop1<Nx-NNN&&loop2>NNN-1&&loop2<Ny-NNN&&loop3>NNN-1&&loop3<Nz-NNN)
		  {Hzx[loop1][loop2][loop3]+=Hzy[loop1][loop2][loop3];
		   Hzy[loop1][loop2][loop3]=0;}

	   }  
	  }
	  }          /*设置波动源，向前波动1个时间周期*/


       return(1);    /*此步仅为便于调试*/
       
}