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

因为不想用C就改了个，随便看看，大家如果也有改的发上给偶学习学习，偶也刚开始学。 bY~@}gC**@  
clc; 5HS~op2n/  
clear all; 1t9.fEmT  
close all; tD$lNh^  
,LOQDIyn  
cc = 3.0e8; iK_c.b  
muz = 4.0*pi*1e-7; V,ZY
*f0  
epsz = 1/(cc*cc*muz); >\A8#@1  
JmpsQ,,  
ie = 200;                   {I{:GcS  
ic = ie/2; "gW7<ilw   
ib = ie-1; rp7W }P+uU  

V^il$'  
spread = 12.0;          %refernce with width of Gauss pulse K^t?gt@k}  
t0 = 40; j*;N\;iL!*  
t = 0; Vga-@  
n = 0; M*pRv  
pulse = 0.0; %}>dqUyQ  
a1N!mQ^  
ex(1:ie) = 0.0; u2(eaP8d  
hy(1:ie) = 0.0; o0z67(N&g  
SNT5Amz!  
x = linspace(1,1*ie,ie); DW(~Qdk  
subplot(2,1,1),plot(x,ex(1:ie),'r'),axis([0 ie -1 1]); _w ]4~V9  
ylabel('EX'); =wq;@'U  
%axis([xmin xmax ymin ymax]) B ?%L  
Z%gx%$  
subplot(2,1,2),plot(x,hy(1:ie),'b'),axis([0 ie -1 1]); m|svQ-/j  
ylabel('HY'); IJU0[EA]F  
AN9[G  
rect=get(gcf,'Position'); 1QJBb \  
rect(1:2)=[0 0]; pz doqAVI  

1q.(69M  
M=moviein(n/2,gcf,rect); Cef:tdk7  
F:37MUQi  
for n =1:200 T,JA#Rk|1N  
    t = t+1; ;Pb8YvG1$  
    pulse = exp(-0.5*(((t0-t)/spread)^2)); bZipm(e  
    ex(ic) = pulse; \.s`n2.w  
    ex(2:ie) = ex(2:ie)+0.5*(hy(1:ib)-hy(2:ie)); Ey&aBYR  
    hy(1:ib) = hy(1:ib)+0.5*(ex(1:ib)-ex(2:ie)); b*,
3< 9  
w>~M}Ahj  
    if mod(n,10)==0; 0NO1M)HQv  
       subplot(2,1,1),plot(x,ex(1:ie),'r'),axis([0 ie -1 1]); G\iyJSj[P  
       ylabel('EX'); t0m*PJcF  
~0-g%C?R  
       subplot(2,1,2),plot(x,hy(1:ie),'b'),axis([0 ie -1 1]); J7i+c];!<  
   & .. `3iQZui  
?n'OFpd  

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

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

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

haohaoaho!

谢谢,学习以下!!

交流学习下

好东西啊~！~~

好   啊   好   啊

看看～～