维波导传播常数计算.docVIP

　　　　　　　　二维波导　传播常数计算(MatLab语言) clear c=2.9979E8; u0=1.256637E-6;e0=8-12; a=7.112e-3;b=3.556e-3; %a*b nt=7500;nf=300; %时间步数与频率点数 nx1=21;ny1=41; %fields array for Transverse section nt1=500;nt2=1000;nt3=5500; %取样时刻 pi2=pi*2;pia=pi/a; dpaj=-j*pi2; nx=nx1-1;ny=ny1-1; js=ny/2+7;jc=5;ic=nx1-4; %? dx=b/nx;dy=a/ny; %步长 sbeta=0 dt=0.99/c/sqrt(1/dx^2+1/dy.^2);