[RFEDA原创]matlab改写的Electromagnetic Simulation Using  the FDTD  Method中的fdtd1_1.c

因为不想用C就改了个，随便看看，大家如果也有改的发上给偶学习学习，偶也刚开始学。 i47LX;}  
clc; ^sf,mM~D  
clear all; F@?-^ E@  
close all; %6 *c40  
<L]Gk]k_R  
cc = 3.0e8; N6._Jb  
muz = 4.0*pi*1e-7; 8%K{lg"  
epsz = 1/(cc*cc*muz); w1tM !4r  
_Ay^v#a  
ie = 200;                   5mZ2CDV  
ic = ie/2; MKfK9>a  
ib = ie-1; q}t]lD %C  
_^& q
,S  
spread = 12.0;          %refernce with width of Gauss pulse }&vD(hX  
t0 = 40; T_D3WHp  
t = 0; %>&ex
0j]  
n = 0; V[(fE=cIN~  
pulse = 0.0; p-iFe\+  
GURiW42  
ex(1:ie) = 0.0; {IgH0+z  
hy(1:ie) = 0.0; U 3< 3T  
?>B?*IK!  
x = linspace(1,1*ie,ie); Dz0D ^(;V  
subplot(2,1,1),plot(x,ex(1:ie),'r'),axis([0 ie -1 1]); <YW)8J  
ylabel('EX'); I(*3n"  
%axis([xmin xmax ymin ymax]) ':(AiD-}  
5f7id7SI  
subplot(2,1,2),plot(x,hy(1:ie),'b'),axis([0 ie -1 1]); zpwoK&T+  
ylabel('HY'); qKD  

m#UQ,
EM  
rect=get(gcf,'Position'); Q` &#u#  
rect(1:2)=[0 0]; 6=/sEzS'  
f0H.$UAL  
M=moviein(n/2,gcf,rect); )c4tGT<  
7Vy_Cec1  
for n =1:200 <W#G)c0  
    t = t+1; ;y,g%uqE  
    pulse = exp(-0.5*(((t0-t)/spread)^2)); <J<"`xKL  
    ex(ic) = pulse; 9I/l+IS"X  
    ex(2:ie) = ex(2:ie)+0.5*(hy(1:ib)-hy(2:ie)); i9L]h69r  
    hy(1:ib) = hy(1:ib)+0.5*(ex(1:ib)-ex(2:ie)); Nud =K'P=  
Ss%Cf6qdWL  
    if mod(n,10)==0; vcFR Td  
       subplot(2,1,1),plot(x,ex(1:ie),'r'),axis([0 ie -1 1]); c K\   
       ylabel('EX'); Wo "s;Z  
^UKAD'_#%O  
       subplot(2,1,2),plot(x,hy(1:ie),'b'),axis([0 ie -1 1]); [P+kQBLpL  
   & .. tM]Gu?6  
F- -g?Q^  

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

终于找到了，呵呵，请问楼主有没有写这本书的matlab版FDTD_2.2.C啊，我写的时候，用一个频率来做，结果就正确，换成三个频率，Fourier变换后Figure里就啥显示都没有了，我也不知道哪里出问题了，请教一下

学习一下~ 我也在学习用matlab                                        

haohaoaho!

谢谢,学习以下!!

交流学习下

好东西啊~！~~

好   啊   好   啊

看看～～