1543阅读
2回复

[求助]2D_Y型光波导形成图像问题（MATLAB程序）

UID ：29143

注册：2009-04-04
登录：2009-05-02
发帖：2
等级：仿真新人

0楼发表于: 2009-05-02 09:07:41

大家帮我看下这段程序  为什么会形成不规则的Y型  有什么地方需要修改吗
%tic
%Initial(原始的) parameters and other things.
W=0.36; %Normalized（规范的） frequency
%The following parameters are control parameters.
WaveGuide=1; %If this program is for Wave Guide? If so, please specify 1.
IsMovie=0; %If you want to play movie, please use 1.
IsFigure=1; %If it will plot the figures? If so, please specify 1.
WantToSeeEp=1; %Do you want to see the distrubution of Ep? If so, please specify 1.
%End of defining control parameters

MLatx=11; %11How many Lattice（格子） cell in x direction.
MLaty=11; %11How many Lattice cell in y direction.

NMlat=21; %21The grid number in each Lattice Cell.
          %SHOULD BE ODD INTEGER（奇整数）!!!!
if mod(NMlat,2)==0
   NMlat=NMlat+1;
end       %Force it to be a odd(奇数) integer!

NTx=MLatx*NMlat+1; %It is the number of the Grid along x axis.
NTy=MLaty*NMlat+1; %It is the number of the Grid along y axis.

if WaveGuide==1
   Nrow=5; %4位置问题The row number of columns between the PML boundary and the waveguide.
end

NPML=12; %12大小就是时间长短问题How many PML layers will be used in our computation.

NTimeSteps=2500; %2500Total number of Time Steps

Meach=10; %20Define the interval(间隔) for plot figures if IsFigure==1.
           %This also works for saving intervals.

R=0.2; %刚好全反射The radius of dielectric columns
ea=11; %11.4The dielectric constant of these columns.

Zmax=0.6;    %0.6The maximum value for z axis when plotting figures.
Colormax=0.6; %0.6越大越暗The maximum value for colormap when plotting figures.

%Some constants
mu0=4*pi*1.0e-7; %Epsilon Zero, if using Gauss Unit, it equals to 1.
e0=8.85*1e-12;   %Mu Zero, if using Gauss Unit, it equals to 1.
c=1/sqrt(mu0*e0); %The light speed.
factor=mu0/e0; %The factor between conductivity(导电率) and permeability（浸透性）.
               %Permeability=Conductivity*factor, in PML.

a=1;%e-6;   % 1 The lattice constant.
W=W*(2*pi*c/a); %frequency

Dx=a/NMlat; %Delta x.
Dy=Dx; %Delta y.
Dt=1/sqrt(1/(Dx*Dx)+1/(Dy*Dy))/c; %Time interval

%tic
%In the following partm we define the dielectric constants:
Ep=ones(NTx-1,NTy-1)*e0;
Ep_cell=ones(NMlat,NMlat)*e0;
x=-(NMlat-1)/2*Dx:Dx:(NMlat-1)/2*Dx;
[X,Y]=meshgrid(x);
X=X';
Y=Y';
flag=find(sqrt(X.^2+Y.^2)<R);
Ep_cell(flag)=e0*ea;
Ep=repmat(Ep_cell,MLatx,MLaty);

   if WaveGuide==1
   eb=1; %The dielectric constant in the waveguide
   Ep(  1:(MLatx-Nrow)*NMlat,                        Nrow*NMlat+1:(Nrow+1)*NMlat)=e0*ea;
   Ep( (MLatx-Nrow-1)*NMlat+1:(MLatx-Nrow+1)*NMlat,(Nrow+1)*NMlat+1:(Nrow+2)*NMlat)=e0*ea;%e0*eb
   Ep( (MLatx-Nrow)*NMlat+1:(MLatx-Nrow+2)*NMlat,(Nrow+2)*NMlat+1:(Nrow+3)*NMlat)=e0*ea;%e0*eb
   Ep( (MLatx-Nrow+1)*NMlat+1:(MLatx-Nrow+3)*NMlat,(Nrow+3)*NMlat+1:(Nrow+4)*NMlat)=e0*ea;%e0*eb
   Ep( (MLatx-Nrow+2)*NMlat+1:(MLatx-Nrow+4)*NMlat,(Nrow+4)*NMlat+1:(Nrow+5)*NMlat)=e0*ea;%e0*eb
   Ep( (MLatx-Nrow+3)*NMlat+1:(MLatx-Nrow+5)*NMlat,(Nrow+5)*NMlat+1:(Nrow+6)*NMlat)=e0*ea;%e0*eb
   Ep( (MLatx-Nrow+3)*NMlat+1:(MLatx-Nrow+5)*NMlat,(Nrow-5)*NMlat+1:(Nrow-4)*NMlat)=e0*ea;%e0*eb
   Ep( (MLatx-Nrow+2)*NMlat+1:(MLatx-Nrow+4)*NMlat,(Nrow-4)*NMlat+1:(Nrow-3)*NMlat)=e0*ea;%e0*eb
   Ep( (MLatx-Nrow+1)*NMlat+1:(MLatx-Nrow+3)*NMlat,(Nrow-3)*NMlat+1:(Nrow-2)*NMlat)=e0*ea;%e0*eb
   Ep( (MLatx-Nrow)*NMlat+1:(MLatx-Nrow+2)*NMlat,(Nrow-2)*NMlat+1:(Nrow-1)*NMlat)=e0*ea;%e0*eb
   Ep( (MLatx-Nrow-1)*NMlat+1:(MLatx-Nrow+1)*NMlat,(Nrow-1)*NMlat+1:Nrow*NMlat)=e0*ea;

   end
%toc

if WantToSeeEp==1
   x=0:Dx:(NMlat*MLatx-1)*Dx;
   x=x-(NMlat*MLatx-1)* ..

未注册仅能浏览部分内容，查看全部内容及附件请先登录或注册

附件: