my.h

一个用于天线系统中的近远场转换VC++程序

#include <stdio.h>
#include <math.h>



#define  NNN	6      		
#define  Nx		60      
#define  Ny		60    
#define  Nz		60     
#define  NN     1000              	
#define  Np     4000     

#define  mx     Nx/2
#define  my     Ny/2
#define  mz     Nz/2
#define  tmax   1000          
#define  Gap    0.02

const double    f=1.5e8;
const double	pi=3.1415926;
const double    miu=4*pi*0.0000001;
const double    ee=8.851*1e-12;
const double    c=3.0*1e8;
const double	dx=0.10;
const double    dy=0.10;
const double    dz=0.10;
const double    dt=dx/(2*c);
const double    r=dx/5;   
const double    R0=0.00001;      
const double    sgmmax=-log(R0)/(2*NNN*dx)*3*ee*c;  
const double    sgmmay=-log(R0)/(2*NNN*dy)*3*ee*c;  
const double    sgmmaz=-log(R0)/(2*NNN*dz)*3*ee*c;  
const double    sgmmmax=sgmmax*miu/ee; 
const double    sgmmmay=sgmmay*miu/ee;  
const double    sgmmmaz=sgmmaz*miu/ee; 

const double    l=1;
const double    mydl=Nx/2-l/dx/2;
const double    myur=Nx/2+l/dx/2;


float Ezx[Nx+1][Ny+1][Nz];
float Exy[Nx][Ny+1][Nz+1];
float Eyx[Nx+1][Ny][Nz+1];


float Hzx[Nx][Ny][Nz+1];
float Hxy[Nx+1][Ny][Nz];
float Hyx[Nx][Ny+1][Nz];

float Ezy[Nx+1][Ny+1][Nz];
float Exz[Nx][Ny+1][Nz+1];
float Eyz[Nx+1][Ny][Nz+1];

float Hzy[Nx][Ny][Nz+1];
float Hxz[Nx+1][Ny][Nz];
float Hyz[Nx][Ny+1][Nz];

double Ex[Nx][Ny+1][Nz+1]={0};
double Ey[Nx+1][Ny][Nz+1]={0};
double Ez[Nx+1][Ny+1][Nz]={0};

double Hx[Nx+1][Ny][Nz]={0};
double Hy[Nx][Ny+1][Nz]={0};
double Hz[Nx][Ny][Nz+1]={0};

const int DDD=3;  
void cal_J(int M);
void waitui(int M);

double sgmx[Nx],sgmy[Ny],sgmz[Nz];   
double sgmmx[Nx],sgmmy[Ny],sgmmz[Nz];   

float E1[NN+1],E2[NN+1];

float si[180],co[180];
void Initial();      
void fangxiangtu();       

 struct bb
 {
   float shi;
   float xu;
 }pp1[Np],pp2[Np],zk[Np/2],fxhs[180];

struct aa
 {
   struct bb sita;
   struct bb fai;
 }fxhs1;
struct cc
 {
     float zhengfu;
     float xiangwei;
 };

struct cc Jxu[Nx-2*NNN-2*DDD][Ny-2*NNN-2*DDD+1],Jxd[Nx-2*NNN-2*DDD][Ny-2*NNN-2*DDD+1],Jmxu[Nx-2*NNN-2*DDD][Ny-2*NNN-2*DDD+1],Jmxd[Nx-2*NNN-2*DDD][Ny-2*NNN-2*DDD+1];
struct cc Jxr[Nx-2*NNN-2*DDD][Nz-2*NNN-2*DDD+1],Jxl[Nx-2*NNN-2*DDD][Nz-2*NNN-2*DDD+1],Jmxr[Nx-2*NNN-2*DDD][Nz-2*NNN-2*DDD+1],Jmxl[Nx-2*NNN-2*DDD][Nz-2*NNN-2*DDD+1];
struct cc Jyu[Nx-2*NNN-2*DDD+1][Ny-2*NNN-2*DDD],Jyd[Nx-2*NNN-2*DDD+1][Ny-2*NNN-2*DDD],Jmyu[Nx-2*NNN-2*DDD+1][Ny-2*NNN-2*DDD],Jmyd[Nx-2*NNN-2*DDD+1][Ny-2*NNN-2*DDD];
struct cc Jyf[Ny-2*NNN-2*DDD+1][Nz-2*NNN-2*DDD],Jyb[Ny-2*NNN-2*DDD+1][Nz-2*NNN-2*DDD],Jmyf[Ny-2*NNN-2*DDD+1][Nz-2*NNN-2*DDD],Jmyb[Ny-2*NNN-2*DDD+1][Nz-2*NNN-2*DDD];
struct cc Jzr[Nx-2*NNN-2*DDD][Nz-2*NNN-2*DDD+1],Jzl[Nx-2*NNN-2*DDD][Nz-2*NNN-2*DDD+1],Jmzr[Nx-2*NNN-2*DDD][Nz-2*NNN-2*DDD+1],Jmzl[Nx-2*NNN-2*DDD][Nz-2*NNN-2*DDD+1];
struct cc Jzf[Ny-2*NNN-2*DDD+1][Nz-2*NNN-2*DDD],Jzb[Ny-2*NNN-2*DDD+1][Nz-2*NNN-2*DDD],Jmzf[Ny-2*NNN-2*DDD+1][Nz-2*NNN-2*DDD],Jmzb[Ny-2*NNN-2*DDD+1][Nz-2*NNN-2*DDD];