www.pudn.com > fdtd_twomedium.rar > fdtd_singlemedium.m, change:2011-05-16,size:1593b


clear all; 
 
sigma=0.05;      %介质1中的电导率 
eps=8.854185*10^(-12); 
eps1=5*eps;           %介质1中的介电常数 
eps=8.854185*10^-12;  %真空介电常数 
deltt=2.5*10^(-12);%时间步长 
line('erasemode','xor'); 
D=deltt*sigma/eps1; 
high1=0; 
high2=0; 
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)-h(k-1));  % 规约化的电场迭代  
    end; 
     
  e(1)=exp((-0.5)*((40-t)/12)^2);  %  电场源:高斯源 
 % e(1)=cos(0.2*t); %  电场源:正弦源 
 
for k=1:200; 
          h(k)=h(k)-0.5*(e(k+1)-e(k)); %  磁场迭代 
      end; 
     subplot(211)   
   plot(e); 
     subplot(212) 
       plot(h);  
   axis([1 201 -1.5 1.5]);  
   drawnow; 
   pause(0.01); 
     
end;