求助

hGU  m7  
有耗无色散条件下 模拟TM波    正弦激励源 FDTD  ,V^2Oa  
最后画出的图本来应该是 同心圆 可我模拟出来的却总是 1.回字形的波形 2.而且波形在很小的范围内 无论介质参数和波源强度怎么改 波形还是很小 总传不出去 ZLDO&}  
QmgO00{  
检查了好长时间 都不知咋回事 E] 6]c!2:  
请帮下忙 指点一下  UuS6y9@v  
我刚入 .. MO^Q 8v  
Ztl?*zL  

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

% 成功 ?  !!!! HVdB*QEH  
% 2D TM Uw]o9 e0S  
%  有耗介质 时间步进FDTD，强迫源，无边界和角点处理，须保证传播不到 O2":)zU.  
r $[{sW  
clear; Dm%Q96*VAq  
!UR3`Xk  
%  物理常数 s",
G w]8  
ip0=10^(-9)/36/pi; j115:f  
u0=4*pi*10^(-7); qm<-(Qc(W  
ita0=120*pi; Pv=]7>e  
c=3*10^(8); XqU0AbQ  
m70AWG  
%  波源参数 Jz4;7/  
e0=1000.0; k0&FUO  
lamda=1.0; B}P!WRNmln  
f=c/lamda; pHowioFx  
omiga=2*pi*f; )gz]F_  
rmi&{o:  
%介质系数 P1z:L
 
ip=30*ip0; IAWs}xIly  
uu=8*u0; 37biRXqLH  
sita=4; uJ=d!Kn  
H2xDC_Fs  
% 步长 DxD\o+:r  
dels=lamda/20; 39m8iI%w[  
delt=dels/2/c; Lxn-M5RPQ  
;s#I b_  
% 区域参数，以（1，1）为坐标原点，计算时坐标不加 1，作图或计算标准偏差时平移一下 l*(L"]  
% i 表列，j 表行，数组引用 i 在后，j 在前  ！！！ :fDzMD  
imax=321; mN l[D  
jmax=321; tSY4'  
k{'<J(Hb  
%  源坐标 Q@|"xKa  
is=151; 7Le-f  
js=151; d04gmc&*  
{3SK|J`  
%  时间步数 m^
zD']  
nmax=150; qW*k|;S  
t<"`gM^|  
% 实际区域 MR: H3  
imax2=21; |z8_]o+|r1  
jmax2=21; 1;?w#/&t  
~.6% %1?  
mE=Tj%+x  
% 数组初始化 EuA352x  
ez=zeros(jmax,imax); 'K}2m  
hx=zeros(jmax,imax); `
gF]  
hy=zeros(jmax,imax); V6+:g=@U-l  
aK(e%Ed t"  
>l=jJTJ;q  
ez2=zeros(jmax2,imax2); 3g:P>(  
*N
M*
 
xia=zeros(1,imax2); zlB[Eg^X  
zuo=zeros(jmax2,1); 4uh~@Lv  
FjI1'Ah\  
Ted tmX$  
%  迭代系数 Pf;RJeD  
ca=(1-sita*delt/2/ip)/(1+sita*delt/2/ip); z =\ENG|x#  
cb=(delt/ip)/(1+sita*delt/2/ip); VD $PoP  
cp=1; Z%b1B<u$  
cq=delt/uu; |MN2v[y  
[S-#}C?~  
+ rM]RFi  
% 开始迭代 @$CPTv3e  
AFeFH.G6Jr  
for n=0:nmax (w+SmD  
     IGo5b-ds  
%  计算 ez     X{[$4\di{  
    for i=2:(imax-1) +;*4.}  
        for j=2:(jmax-1)             Z^
Wv(:Nr  
|Bv,*7i&  
           if ( abs(i-is)<0.0001 ) & ( abs(j-js)<0.0001 )  ;5  
              lin1=omiga*n*delt; +
+ 5!8Nv  
              ez(j,i)=e0*sin(lin1); EMxMJ=  
           else 5E 9
R+N  
              lin1=hy(j,i)-hy(j,i-1); 3[?;s}61  
              lin2=hx(j,i)-hx(j-1,i); DwK$c^2q{.  
             {9) 
HB:  
              lin3=(lin1-lin2); h`eHoKJ#w  
                           &os:h] C  
              ez(j,i)=ca*ez(j,i)+cb*lin3/dels; (=Oo=8\  
           end 2|(lKFkQ  
           0bD\`Jiv,  
       end bYX.4(R  
   end snNB;hkj  
A;6ew
4  
%  计算 hx     E70o nR!i  
    for i=1:imax _ /Eg_dQ~@  
        for j=1:(jmax-1) %sPq*w.  
             8A/rkoht*  
            lin1=( ez(j+1,i)-ez(j,i) ); .81 ~ K[  
            hx(j,i)=cp*hx(j,i)-cq*lin1/dels; hBifn\dFr  
           s$
lJJL  
       end |+JC'b?,  
   end JuOCOl\  
$" =3e]<  
%  计算 hy   /%F,
  
    for i=1:(imax-1) T.B7QAI. H  
        for j=1:jmax <xrya_R?  
             5
{IbKj|  
            lin1=( ez(j,i+1)-ez(j,i) ); *+00  
            hy(j,i)=cp*hy(j,i)+cq*lin1/dels; W59xe&l  
           l<(jm{q?u  
       end 7j(gW  
   end fP8bWZ{  
(\wV)c9  
end   B5R7geC  
^&c &5S}  
%  绘制曲线 ttt4h  
P?jI:'u!R.  
    for i=141:161 F}@]Lq+  
        for j=141:161 W@%g_V}C*  
             G,1g~h%I$  
            ez2(j-140,i-140)=ez(j,i); t1$pl6&,  
           9[ o$/x}  
       end ?iamo.0zN  
   end   +CXq41g"c  
nQg_1+  
    for i=1:imax2         &IsPqO  
        xia(1,i)=(i-1)*dels;         o,?h}@  
    end     }D3hP|.X  
5Q8s{WQ  
    for j=1:jmax2         RX^Xtc"  
        zuo(j,1)=(j-1)*dels;         =+[`9  
    end ~at:\h4:  
     0bSnD|#I  
    contour( xia,zuo,ez2,40);   q]s_hWWv  
    xlabel('m'); I= cayR  
    ylabel('m'); t8.
3  
   ez(31,31) jz>b>;  
%  程序结束！！！

我觉得可能是你的递推公式里面的系数错了吧

我检查了好几遍   没错啊  

按照你的描述，无论源的幅度多大，都会很快衰减。这个很容易调试程序看出来的啊。