fdtd_1d.m

fdtd 例程

clc;clear;
low1=0;
low2=0;
high1=0;
high2=0;

figure(1);

line('erasemode','xor');

for i=1:201;  %  电场初始化
    e(i)=0;
end;

for i=1:200;  %  磁场初始化
    h(i)=0;
end;

for t=0:1000;
    
    for k=2:200;
        e(k)=e(k)+0.5*(h(k-1)-h(k));  %  电场迭代 
    end;
    
    e(101)=e(101)+exp((-0.5)*((40-t)/12)^2); %  软源
   
    
%              e(101)=exp((-0.5)*((40-t)/12)^2);  %  硬源
    %         e(101)=e(101)+cos(0.2*t); %  正弦源
    
            e(1)=low2;  %  吸收边界
            low2=low1;
            low1=e(2);
            
            e(201)=high2;
            high2=high1;
            high1=e(200);
%     
%         e(1)=0;  金属边界
%         e(201)=0;

    
    for k=1:200;
        h(k)=h(k)+0.5*(e(k)-e(k+1)); %  磁场迭代
    end;
    
   plot(e);
   axis([1 201 -1.5 1.5]); %  吸收边界显示需要  
   drawnow;
   pause(0.01);
    
end;