movon.m

二维光子晶体弯曲波导的fdtd程序,matlab语言 ,TE波, PML边界条件,平面波光源

%--------------------------------------------------------------------------------------
Ex(:,1:J-1,3)=C(:,1:J-1).*(-D(:,1:J-1).*Ex(:,1:J-1,2)...
                  +(Hz(:,2:J,2)-Hz(:,1:J-1,2))/dy);

%--------------------------------------------------------------------------------------
Ey(1:I-1,:,3)=C(1:I-1,:).*(-D(1:I-1,:).*Ey(1:I-1,:,2)...
                   -(Hz(2:I,:,2)-Hz(1:I-1,:,2))/dx);
               
% Ey(Is,Js,3)=Ey(Is,Js,3)+S/2; %   设置波源
%--------------------------------------------------------------------------------------
Hz(2:I-1,2:J-1,3)=Hz(2:I-1,2:J-1,2)...
                 -(Ey(2:I-1,2:J-1,3)-Ey(1:I-2,2:J-1,3))/2+(Ex(2:I-1,2:J-1,3)-Ex(2:I-1,1:J-2,3))/2;
             


Hz(1,2:J-1,3)=xcdt(1,2,Hz);
Hz(I,2:J-1,3)=xcdt(I,I-1,Hz);
Hz(2:I-1,1,3)=ycdt(1,2,Hz);
Hz(2:I-1,J,3)=ycdt(J,J-1,Hz);%  计算边界


Hz(1,1,3)=coner(1,Hz);
Hz(I,1,3)=coner(2,Hz);
Hz(I,J,3)=coner(3,Hz);
Hz(1,J,3)=coner(4,Hz);%  计算角顶点

Hz(Is,Js,3)=Hz(Is,Js,3)+S; %  设置波源

%--------------------------------------------------------------------------------------



%   继续计算下一个时间步的场强值
Hz(:,:,1)=Hz(:,:,2);
Hz(:,:,2)=Hz(:,:,3);

Ex(:,:,1)=Ex(:,:,2);
Ex(:,:,2)=Ex(:,:,3);

Ey(:,:,1)=Ey(:,:,2);
Ey(:,:,2)=Ey(:,:,3);