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

因为不想用C就改了个，随便看看，大家如果也有改的发上给偶学习学习，偶也刚开始学。 :TRhk.  
clc; |y DaFv  
clear all; Wu@v%!0  
close all; (X/dP ~  
t2hI^J0y  
cc = 3.0e8; A|<i7QVY  
muz = 4.0*pi*1e-7; r}M2t$nv  
epsz = 1/(cc*cc*muz); (_nU}<y_i  
b~Un=-@5a  
ie = 200;                   /;;$9O9  
ic = ie/2; &'"dYZj{  
ib = ie-1; *T-v^ndJh  
Z@gEJ^"yA"  
spread = 12.0;          %refernce with width of Gauss pulse M-zqD8D  
t0 = 40; ?F^O7\rw  
t = 0; $0,lE+7*  
n = 0; u!TMt8+c  
pulse = 0.0; tk66Ggi[K  
xz*MFoE  
ex(1:ie) = 0.0; "VgPaz
#  
hy(1:ie) = 0.0; q vGP$g  
+L-(Lz[p  
x = linspace(1,1*ie,ie); s0/y> ok  
subplot(2,1,1),plot(x,ex(1:ie),'r'),axis([0 ie -1 1]); kn= fW1  
ylabel('EX');
Ocf:73t  
%axis([xmin xmax ymin ymax]) _r?.%]\.  
WKB K)=  
subplot(2,1,2),plot(x,hy(1:ie),'b'),axis([0 ie -1 1]); W0\ n?$ZC~  
ylabel('HY'); b@ OF  
4O'%$6KR(  
rect=get(gcf,'Position'); RCY}JH>}  
rect(1:2)=[0 0]; NxrfRhaU3  
.mvpFdn  
M=moviein(n/2,gcf,rect); LNOz.2fr>  
;t\C!A6  
for n =1:200 4Ixu%  
    t = t+1; .q5WK#^  
    pulse = exp(-0.5*(((t0-t)/spread)^2)); UE 1tm  
    ex(ic) = pulse; &H:2TL!  
    ex(2:ie) = ex(2:ie)+0.5*(hy(1:ib)-hy(2:ie)); +P|Z1a -jB  
    hy(1:ib) = hy(1:ib)+0.5*(ex(1:ib)-ex(2:ie)); {O5(O oDa  
Rd,5&X$  
    if mod(n,10)==0; CH] +S>$  
       subplot(2,1,1),plot(x,ex(1:ie),'r'),axis([0 ie -1 1]); r#\Lq;+-B  
       ylabel('EX'); a]t| /Mq  
@2/xu  
       subplot(2,1,2),plot(x,hy(1:ie),'b'),axis([0 ie -1 1]); UUR` m  
   & .. OY,iz  
6I-Qq?L[H  

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

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

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

haohaoaho!

谢谢,学习以下!!

交流学习下

好东西啊~！~~

好   啊   好   啊

看看～～