[求助]谁能帮我看看这个程序？FDTD  matlab

哪位高手帮我看看这个程序有什么错误，为什么出不了图形？请各位高手帮忙啊！是关于同轴电缆的FDTD的~~~大家帮帮忙~~~~~感激不尽 ?}y{tav=  
8cqH0{  
Z^AOV:|m  
qJY'"_Q{  
*>x~`  
Nqih LUv  
clear; >j [> 0D  
clc; [yzDa:%  
close all; |P{K\;-  
%*********************************************************************** GfEg][f  
%    基本常数 H-&Z+4 +Xs  
%*********************************************************************** Rj&V~or  
pi = 4*atan(1);  Vo@[  
cc=2.99792458e8;            %真空中的波速 ,2H@xji [  
mu0=4.0*pi*1.0e-7;          %真空中的磁导系数 %1E:rw@
 
eps0=1.0/(cc*cc*mu0);       %真空中的介电系数 FyZp,uD  
Fr=500000000;               %光波频率500MHz (DzV3/+p^  
Lamda=cc/Fr;                %波长 %nU8 Ca  
k0=2*pi/Lamda;              %波数 !+qy~h  
bat0=k0*sqrt(eps0);         %传播常数 'Xl[ y  
v=2*pi*Fr/bat0;             %波速 }yde9b?F  
%*********************************************************************** 7_,)"J2^  
%     同轴电缆参数，参考王均宏论文上的模型 Ep/kb-~-  
%*********************************************************************** ?.&]4z([  
b=8;                        %内导体半径 cVrses^yE  
a=20;                       %外导体半径 L B:wo.X  
c=22.5;                     %护套(绝缘）厚度2.5mm ich\`j[i  
L=100;                      %电缆长度 t@KN+ C  
%*********************************************************************** h^{D "  
%     网格参数 =yn|.%b  
%*********************************************************************** <I}O_:%  
dr=0.5;                    %r方向上的栅距 %
i!
&Fr  
dfi=10;                    %φ方向上的栅距 Z:Hk'|q}I  
dz=0.5;                    %z方向上的栅距 Y9h~ hD  
dt=0.5/(2*v);              %时间步 #b[B$  
ds=v*dt;                   %空间步 yQkj4v{  
tmax=100; 5ZG-3qj  
t0=4; P'^#I[G'  
Nrr=floor((a-b)/dr);        %就近取整，r方向上内导体的外边界到外导体的内边界所占网格数 c< ke)@  
Nr=floor((c-b)/dr);         %r方向上内导体的外边界到护套外层所占网格数 qla$}dnvc  
Nfi=2*pi/dfi;                %φ方向上的网格数 cTy;?
(E  
Nz=floor(L/dz);             %z方向上的网格数 ^4Ta0kDn  
Nrb=Nr+1; 4~< :Pj  
Nfib=Nfi+1; o$k1&hyH  
Nzb=Nz+1; p=T,JAIt  
%*********************************************************************** o<\CA[   
%     求相对介质系数epsr和r(i) ?n$;l-m[  
%*********************************************************************** 577H{;pW  
for i=1:1:Nrr             %内导体和外导体之间的介电常数为1.26 P I gbeP  
    for j=1:1:Nfi 9UD~$_<\  
        for k=1:1:Nz X@$x(Zc  
            epsr=1.26*eps0; }"?v=9.G  
        end ~HbZRDcJc  
    end ~(2G7x)  
end fD8A+aA  
for i=Nrr+1:1:Nr "Dbjp5_  
    for j=1:1:Nfi          %绝缘层的介电常数为2.3 *KF-q?PBb  
        for k=1:1:Nz 'je=.{[lWt  
            epsr=2.3*eps0; wyqXD.of  
        end J%ym1A9  
    end joa|5v'  
for i=1:Nrb Zqa
Ce>  
    r(i,:,:)=dr*i+a; %%%%%%%?????????????????????????????????????????????????????????????????????????? zY@|KV"^r  
end t/?x#X  
%*********************************************************************** HtWuZq;w  
%     入射场初始化 Z-(Vfp
4  
%*********************************************************************** E|!rapa  
Er=zeros(Nr,Nfib,Nzb); vi+k#KE  
Ef=zeros(Nrb,Nfi,Nzb); >_ji`/d{  
Ez=zeros(Nrb,Nfib,Nz); vM5u]u!  
Hr=zeros(Nrb,Nfi,Nz); L_Q1:nL-0  
Hf=zeros(Nr,Nfib,Nz); %T3
L-{s5  
Hz=zeros(Nr,Nfi,Nzb); F<wwuCbF  
%*********************************************************************** =Z ^=
 
%     FDTD迭代 T\g%.  
%***********************************************************************      Eeemy*U  
for Nt=1:1:tmax qq0bIfF\4  
    KsZXdM/  
%%%%%%%%%%%%迭代系数%%%%%%%%%%%% cwK6$Ax  
a0=dt/mu0; ^MPl wx  
a1=a0/dr; h(aF>a\Z  
a2=a0/dr/dfi; (uBevU\  
a3=a0/dz; Q_<CG[,6D1  
b0=dt/eps; X(m&  
b1=b0/dr; 8+5-7)  
b2=b0/dr/dfi; vO#4$,  
b3=b0/dz;  r"s <;  
%%%%%%第一步 求H%%%%%%%%%%% D'y/pv}!  
Hr(2:Nr,1:Nfi,1:Nz)=Hr(2:Nr,1:Nfi,1:Nz)-... 8$ dJh]\Y  
                    a2/r(2:Nr,1:Nfi,1:Nz)*(Ez(2:Nr,2:Nfib,1:Nz)-Ez(2:Nr,1:Nfi,1:Nz))+... `>^2MHF3LT  
                    a3*(Ef(2:Nr,1:Nfi,2:Nzb)-Ef(2:Nr,1:Nfi,1:Nz)); H]a;<V9[  
Hf(1:Nr,2:Nfi,1:Nz)=Hf(1:Nr,2:Nfi,1:Nz)-... !"\UT&  
                    a3*(Er(1:Nr,2:Nfi,1:Nz)-Er(1:Nr,2:Nfi,2:Nzb))+... F?]nPb|  
                    a1*(Ez(2:Nrb,2:Nfi,1:Nz)-Ez(1:Nr,2:Nfi,1:Nz)); Q<$I,C]  
Hz(1:Nr,1:Nfi,2:Nz)=Hz(1:Nr,1:Nfi,2:Nz)-... ve.rpF\  
                    a1/(r(1:Nr,1:Nfi,2:Nz)+ds/2)*((r(1:Nr,1:Nfi,2:Nz)+ds)*Ef(1:Nr,1:Nfi,2:Nz)-r(1:Nr,1:Nfi,2:Nz)*Ef(2:Nbr,1:Nfi,2:Nz))+... $E;`Y|r%WK  
                    a2/(r(1:Nr,1:Nfi,2:Nz)+ds/2)*(Er(1:Nr,2:Nfib,2:Nz)-Er(1:Nr,1:Nfi,2:Nz)); oi4tj.!J  
              G
q4~9Tm)*  
%%%%%%第一步 求E%%%%%%%%%%% 9OYsI
 
Er(1:Nr,2:Nfb,2:Nz)=Er(1:Nr,2:Nfi,2:Nz)+... AhOvI{  
                      b2/(r(1:Nr,1:Nfi,2:Nz)+ds/2)*(Hz(1:Nr,2:Nfi,2:Nz)-Hz(1:Nr,1:Nfi-1,2:Nz))-... p"Ki$.Y  
                      b3*(Hf(1:Nr,2:Nfi,2:Nz)-Hf(1:Nr,2:Nfi,1:Nz-1)); Rn={:u4  
Ef(2:Nr,1:Nfi,2:Nz)=Ef(2:Nr,1:Nfi,2:Nz)+... 2J;h}/!H  
                     b3*(Hr(2:Nr,1:Nfi,2:Nz)-Hr(2:Nr,1:Nfi,1:Nz-1))-... KD]8n]c
 
                     b1*(Hz(2:Nr,1:Nfi,2:Nz)-Hz(1:Nr-1,1:Nfi,2:Nz)); " $5J7  
                 {9wBb`.n^  
Ez(2:Nr,2:Nfi,1:Nz)=Ez(2:Nr,2:Nfi,1:Nz)+... Z/=x(I0  
                     b1/r(2:Nr,2:Nfi,1:Nz)*((r(2:Nr,2:Nfi,1:Nz)+ds/2)*Hf(2:Nr,2:Nfi,1:Nz)-(r(2:Nr,2:Nfi,1:Nz)-ds/2)*Hf(1:Nr-1,2:Nfi,1:Nz))-... Snx_NH#tA  
                     b2/r(2:Nr,2:Nfi,1:Nz)*(Hr(2:Nr,2:Nfi,1:Nz)-Hr(2:Nr,1:Nfi-1,1:Nz)); .VF4?~+M-  
%*********************************************************************** {5}UP@h  
%     导体处电场切向为0 ]<D9Q>  
%***********************************************************************  0*?~I;.2m$  
Ef(Nr,1:Nfi,1:Nz)=0; E]PHO\f-m}  
Ez(Nr,1:Nfi,1:Nz)=0; 9N^&~O|1  
%*********************************************************************** j'lfH6_')e  
%     激 .. YKbR#DC\  
^Xjh?+WM  

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

为什么大家都不回复？

大家都过来看看行不行哦？

程序中是不是还有定义错误啊？

不要着急，我先看下你的程序。你说出不了图像是什么意思，详细描述下。

程序还有运行的错误啊。

错误太多了，改不过来啊。

Ef(2:Nbr,1:Nfi,2:Nz))里面的Nbr是什么东西？？ iM;Btv[|  
Nfb又是什么东西，都没有定义？

谢谢楼上的帮助，感激你~~