求助

w0tlF:Eg  
有耗无色散条件下 模拟TM波    正弦激励源 FDTD  M_$pqVm  
最后画出的图本来应该是 同心圆 可我模拟出来的却总是 1.回字形的波形 2.而且波形在很小的范围内 无论介质参数和波源强度怎么改 波形还是很小 总传不出去 tCtR(mG=A  
Zdj~B1  
检查了好长时间 都不知咋回事 ?i)-K?4Sb  
请帮下忙 指点一下  NZ9`8&93  
我刚入 .. >Mrz$ z{x  
XC0G5rtB  

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

% 成功 ?  !!!! YPf&y"E&H  
% 2D TM o=w&&B  
%  有耗介质 时间步进FDTD，强迫源，无边界和角点处理，须保证传播不到 Gv?3}8Wp  
qNC.|R  
clear; |~8iNcIS  
L a0H  
%  物理常数 wgkh}b   
ip0=10^(-9)/36/pi; n KDX=73  
u0=4*pi*10^(-7); +3]@0VM26;  
ita0=120*pi; p
jIXZ=  
c=3*10^(8); -t:yy:4  
_BHb0zeot  
%  波源参数 X<,sc;"b`k  
e0=1000.0; ToDNBt.u{+  
lamda=1.0; `"A\8)6-  
f=c/lamda; @6h=O`X>  
omiga=2*pi*f; 7!<cU  
 UZmzk  
%介质系数 z/6kxV89  
ip=30*ip0; )-2OraUm<  
uu=8*u0; c6E@+xU  
sita=4; q,0o:
nI  
Fg5>CppH  
% 步长 ]*[S#Jk  
dels=lamda/20; G?'L1g[lc  
delt=dels/2/c; ,Z&"@g  
QOy+T6en  
% 区域参数，以（1，1）为坐标原点，计算时坐标不加 1，作图或计算标准偏差时平移一下 t+q:8HNh  
% i 表列，j 表行，数组引用 i 在后，j 在前  ！！！ -.ha\t0J  
imax=321; $LFzpg  
jmax=321; %$!}MxUM  
jP@H$$-=wH  
%  源坐标 k2k/v[60  
is=151; 'JKFEUzM  
js=151; VAsaJ`vcb  
5a5I+* c  
%  时间步数 %CD}A%~  
nmax=150; ;u!>( QQ  
ZfL\3Mn  
% 实际区域 J3S@1"   
imax2=21; K}&|lCsb  
jmax2=21; ?J%$;"q  
9_5Fl,u z  
+\\,FO_  
% 数组初始化 `<b 3e(A  
ez=zeros(jmax,imax); cDm_QYQ  
hx=zeros(jmax,imax); x_2 [+Ol  
hy=zeros(jmax,imax); )z2Tm4>iql  
)i.p[  
E
G@*J*|S  
ez2=zeros(jmax2,imax2); :&&Ps4\Sq  
wrac\.  
xia=zeros(1,imax2); ?`,Xb.NA$K  
zuo=zeros(jmax2,1); 0Evmq3,9  
T AwA)Zg  
3oE *86  
%  迭代系数 6_^u}me  
ca=(1-sita*delt/2/ip)/(1+sita*delt/2/ip); a}hpcr({?  
cb=(delt/ip)/(1+sita*delt/2/ip); R`!x<J  
cp=1; #wk'&XsC#z  
cq=delt/uu; "'^#I_*Mf  
[@uL)*o_#  
q 1Rk'
k4+  
% 开始迭代 `GqS.O}C  
1EyM,$On  
for n=0:nmax u"?cmg<.1  
     mI@E>VCV[  
%  计算 ez     ]y2(ZTNTs  
    for i=2:(imax-1) ,*@m<{DX)  
        for j=2:(jmax-1)             bF|j%If%  
" F~uTo  
           if ( abs(i-is)<0.0001 ) & ( abs(j-js)<0.0001 ) vd9l1"S  
              lin1=omiga*n*delt; N&-d8[~  
              ez(j,i)=e0*sin(lin1); w2@ `0  
           else LXJ"ct  
              lin1=hy(j,i)-hy(j,i-1); 6^if%62l&  
              lin2=hx(j,i)-hx(j-1,i); UyEyk$6SU  
             " <m)Fh;  
              lin3=(lin1-lin2); CLdLOu"  
                           ]uWx<aDB  
              ez(j,i)=ca*ez(j,i)+cb*lin3/dels; #J\
2/~  
           end O)Nj'Hcu  
           [NQ`S ~_:  
       end %h4pIA  
   end g(;t,Vy,I  
)DI/y1  
%  计算 hx     ])7t!<  
    for i=1:imax OzA"i
y  
        for j=1:(jmax-1) ,fT5I6l  
             H9xxId?3u  
            lin1=( ez(j+1,i)-ez(j,i) ); f:KKOLm  
            hx(j,i)=cp*hx(j,i)-cq*lin1/dels; BYqDC<Fq  
           1Ty{k^%  
       end kK16+`\+  
   end 0>)('Kv  
)67Kd]  
%  计算 hy   1YTnOiYS1  
    for i=1:(imax-1) (9*=d_=  
        for j=1:jmax A%"mySW  
             : N>5{  
            lin1=( ez(j,i+1)-ez(j,i) ); +s V$s]U  
            hy(j,i)=cp*hy(j,i)+cq*lin1/dels; V2^(qpM!  
           d-#MRl$rtK  
       end lw\OsB$  
   end H
d U1gV>  
ZcYxH|Gn  
end   "c\ZUx_i6  
#]_S{sO  
%  绘制曲线 w5bD  
Aq^1(-g  
    for i=141:161 E@}t1!E<  
        for j=141:161 94 H\,}i8  
             WID4{>G2  
            ez2(j-140,i-140)=ez(j,i);  G;Q)A$-  
           LrX7WI  
       end A%Pjg1(uX  
   end   _Su? VxU  
G)?VC^Q  
    for i=1:imax2         C
KNC"Y*X  
        xia(1,i)=(i-1)*dels;        
LR@rn2Z  
    end     $*Njvr7  
yxf|Njo0  
    for j=1:jmax2         J-:\^uP  
        zuo(j,1)=(j-1)*dels;         y"<nx3  
    end w-Y-;*S  
     5 UpN/\He  
    contour( xia,zuo,ez2,40);   a-MDZT<xA+  
    xlabel('m'); ]uI#4t~  
    ylabel('m'); W~$YKBW  
   ez(31,31) j="{^b  
%  程序结束！！！

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

我检查了好几遍   没错啊  

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