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

因为不想用C就改了个，随便看看，大家如果也有改的发上给偶学习学习，偶也刚开始学。 tSe[*V4{'  
clc; )2Bb,p<Wr  
clear all; C-\3,  
close all; BS ]:w(}[  
B
<ue}t  
cc = 3.0e8; ZQ>Q=eCs 1  
muz = 4.0*pi*1e-7; C80< L5\  
epsz = 1/(cc*cc*muz);  /PTq.  
Km,%p@`m  
ie = 200;                   ZL:SJ,C  
ic = ie/2; ra2{8 x  
ib = ie-1; I E{:{b\  
X`J86G)  
spread = 12.0;          %refernce with width of Gauss pulse |*DkriYY  
t0 = 40; ym+Ezb#o  
t = 0; P&Uj?et"  
n = 0; K>C@oE[W  
pulse = 0.0; FmhAUe  
10)RLh|+  
ex(1:ie) = 0.0; T0~~0G)k  
hy(1:ie) = 0.0; 'rTJ*1i  

`F]  
x = linspace(1,1*ie,ie); 0 vYG#S  
subplot(2,1,1),plot(x,ex(1:ie),'r'),axis([0 ie -1 1]); &;~?\>?I  
ylabel('EX'); [jksOC)@4  
%axis([xmin xmax ymin ymax]) r-*j"1 e  
b 0qA  
subplot(2,1,2),plot(x,hy(1:ie),'b'),axis([0 ie -1 1]); fzA Fn$[  
ylabel('HY'); c3 )jsf  
_zJ /z  
rect=get(gcf,'Position'); CI-1>= "OE  
rect(1:2)=[0 0]; OL%}C*Zq  
)%PMDG|  
M=moviein(n/2,gcf,rect); rNurzag  
ioEjbqD<  
for n =1:200 n/x((d%"E  
    t = t+1; &er,Wyc(  
    pulse = exp(-0.5*(((t0-t)/spread)^2)); v*&jA8D  
    ex(ic) = pulse; 3y,2RernK  
    ex(2:ie) = ex(2:ie)+0.5*(hy(1:ib)-hy(2:ie)); CO9PQ`9+  
    hy(1:ib) = hy(1:ib)+0.5*(ex(1:ib)-ex(2:ie)); {3.n!7+  
wa~zb!y<  
    if mod(n,10)==0; =XK}eQ_d  
       subplot(2,1,1),plot(x,ex(1:ie),'r'),axis([0 ie -1 1]); f'w`<  
       ylabel('EX'); ,FXc_BCx4  
d3Y;BxEz  
       subplot(2,1,2),plot(x,hy(1:ie),'b'),axis([0 ie -1 1]); KMK8jJ  
   & .. D/JSIDd  
E-($Xc  

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

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

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

haohaoaho!

谢谢,学习以下!!

交流学习下

好东西啊~！~~

好   啊   好   啊

看看～～