[求助]自己写的求细直导线内感应电流的程序及结果，请高手帮忙分析一下 。

[attachment=9002] 参考论坛中的一些程序，分别用hallen方程和harrington矩量法书上第79页方法编写了处于外电场中的细直导线的感应电流。（程序及结果请见附件）。现遇到的问题是： x2 s%qZ#  
（1）两种方法的大小是不一样的，而且差很多； _;8aiZt|u  
（2）在程序中的波长是105000米。而导线长60米。根本不满足半波的条件，但是得到的电路的形状确是半波的？对半波全波的问题一直不是很理解，希望有高手能指点一二。 urJ>dw?FI  
非常感谢！ +Xb )bfN  
FW?zJ  
现把程序贴出来，请高手指教： sn|q EH  
1、hallen程序： 41pk )8~pt  
u0=4*pi*10^(-7);epth0=8.85*10^(-12); %真空中的介电常数 &) '5_#S  
tlamda=350*10^(-6);   %波长时间 H|E{n/g  
c0=3*10^8; |2!!>1k  
dtl=H/nn; jafIKSD]%  
w0=2*pi/tlamda; P>*g'OK^!G  
lamda=tlamda*c0;  %波长 lkj^<%N"r  
kdl=w0/c0;      %常数k 
qbrf;`  
eta=120*pi;%波阻抗 r6B\yH2  
%************************************************************************* F4!,8)}  
%                       hallen方程 AnB]f~Yjl  
%************************************************************************* Qv3g 4iJ  
for mm=1:(nn+1) eSV_.uvsb  
    Zb(mm)=H-(mm-1)*dtl;%竖直方向上将天线离散化 0SoU\/kUi  
end c5>'1L  
z=zeros(nn+1,nn+1);%生成阻抗矩阵 -EwtO4vLJ  
for mm=1:(nn+1) cfb8kNn~+  
    for kk=1:(nn+1) @;T>*_Yhn  
        if mm==kk         %对应格林函数奇点 <tT*.n
M\  
           zz(mm,kk)=log(dtl/adl)/(2*pi)-(j*kdl*dtl)/(4*pi); -3YsrcJi  
        else Z*/*P4\  
            R1=sqrt((Zb(mm)-Zb(kk))^2+adl^2); f87>ul!*  
            zz(mm,kk)=exp(-j*kdl*R1)*dtl/(4*pi*R1);%基函数和权函数均取狄拉克函数 ?DAW~+,!7o  
        end; P'4oI0Bw  
    end; o6@Hj+,,  
end; b!>w4MPe  
for mm=1:(nn-1) Ihe/P {t]J  
    Z(:,mm)=zz(:,mm+1); /0w?"2-  
end;
Yl65|=ne  
for mm=1:(nn+1) ?*I _'2  
    Z(mm,nn)=cos(kdl*Zb(mm));%把Z矩阵最后两列补足 8J&K_JC^  
    Z(mm,nn+1)=sin(kdl*Zb(mm)); U}c[oA  
end a59l"b  
for i=1:(nn+1) ;%Z%]n
IS  
     temp=zeros(1,n+1); EYwDv4H,g  
     for mm=1:nn {gSR49!Q  
         temp=temp+dEkk(2*mm-1,:)*sin(kdl*abs(Zb(i)-Zb(mm)))+dEkk(2*mm+1,:)*sin(kdl*abs(Zb(i)-Zb(mm+1))); owYSR?aG  
     end #MZ0Sd8]&  
     V(i,:)=temp*dtl/(j*2*eta); N9W\>hKaeh  
end T>"GH M  
I=inv(Z)*V;%最后两个元素不是电流项，而是系数c1和c2! 9!XW):  
current=zeros(nn+1,n+1); ^':Az6Z  
for i=2:nn 1' #%UA  
    current(i,:)=I(i-1,:); ELF,T(  
end n&?]GyQ  
z=H:-dtl:0; Z19d Ted33  
plot(z,abs(current(:,2)),'b') UOWOOdWSB  
hold on .x$!
Rc}  
GN36:>VWb  
2、harrington程序  （《计算电磁场的矩量法》p79页的方法） /o<tmK_m  
%********************************************************************* vt;<+"eps  
u0=4*pi*10^(-7);epth0=8.85*10^(-12); %真空中的介电常数 )~hsd+ 0t
 
tlamda=350*10^(-6);   %波长时间 uUAib<wdPL  
c0=3*10^8; 'o2V}L'nG  
dtl=H/nn; G$ l>By  
w0=2*pi/tlamda; mp_
(ke  
lamda=tlamda*c0;  %波长 X@,xwsM%tb  
kdl=w0/c0;      %常数k ]rEFWA  
eta=120*pi;%波阻抗 R/H?/  
fzz=zeros(nn-1,nn-1);fzf=zeros(nn-1,nn-1);fff=zeros(nn-1,nn-1);ffz=zeros(nn-1,nn-1); m"B)%?C#  
for m=1:(nn-1)      %左正右正组合  b;!oPT  
    for k=1:(nn-1)
&Wj %`T{  
        rzz(m,k)=sqrt((abs(m-k)*dtl)^2+adl^2); %4Nq T  
        rff(m,k)=rzz(m,k); >$gWeFu  
        if m<k %ZV a{Nc  
            rzf(m,k)=sqrt((abs(k-m+1)*dtl)^2+adl^2); G0sg\]  
            rfz(m,k)=sqrt((abs(k-m-1)*dtl)^2+adl^2); @@U'I^iG  
        elseif m==k E'p5  
            rzf(m,k)=sqrt(dtl^2+adl^2); >
+v)^7c  
            rfz(m,k)=rzf(m,k); &hmyfH&S  
        else 9v(&3,)a  
            rzf(m,k)=sqrt((abs(m-k-1)*dtl)^2+adl^2); %eoO3"//  
            rfz(m,k)=sqrt((abs(m-k+1)*dtl)^2+adl^2); opz.kP[e,  
        end urT!?*g,  
        if m-k==1 _7SOl.5ZE  
            fzf(m,k)=log10(dtl/adl)/(2*pi*dtl)-j*kdl/(4*pi); "FC;k >m  
            ffz(m,k)=exp(-j*kdl*rfz(m,k))/rfz(m,k); UQl3Tq4QM  
            fzz(m,k)=exp(-j*kdl*rzz(m,k))/rzz(m,k); 2?LZW1
4$d  
            fff(m,k)=exp(-j*kdl*rff(m,k))/rff(m,k); j)Z3m @Ii5  
        elseif k-m==1 .tB[8Y=J  
            ffz(m,k)=log10(dtl/adl)/(2*pi*dtl)-j*kdl/(4*pi); U8?%Dq%i  
            fzz(m,k)=exp(-j*kdl*rzz(m,k))/rzz(m,k); l?_Iu_Qp  
            fff(m,k)=exp(-j*kdl*rff(m,k))/rff(m,k); B\2<r5|QG  
            fzf(m,k)=exp(-j*kdl*rzf(m,k))/rzf(m,k); mw<LNnT{8  
        elseif m==k GgH=w`;_  
            fzz(m,m)==log10(dtl/adl)/(2*pi*dtl)-j*kdl/(4*pi); Sqw.p#  
            fff(m,m)==log10(dtl/adl)/(2*pi*dtl)-j*kdl/(4*pi);
UzVnC:  
            fzf(m,k)=exp(-j*kdl*rzf(m,k))/rzf(m,k); zosJ=$L  
            ffz(m,k)=exp(-j*kdl*rfz(m,k))/rfz(m,k); =K#D^c~  
        else CT6Ca,  
            fzz(m,k)=exp(-j*kdl*rzz(m,k))/rzz(m,k); +{ QyB  
            ffz(m,k)=exp(-j*kdl*rfz(m,k))/rfz(m,k); rj"oz"  
            fzf(m,k)=exp(-j*kdl*rzf(m,k))/rzf(m,k); Y_&D W4  
            fff(m,k)=exp(-j*kdl*rff(m,k))/rff(m,k); #uNQ+US0  
        end Mn(:qQo^&`  
        if m==k a4FvQH#j  
           fdl(m,m)=log10(dtl/adl)/(2*pi*dtl)-j*kdl/(4*pi); .63:
G<  
           Zdl(m,m)=j*w0*u0*dtl^2*fdl(m,m)+(fzz(m,m)-ffz(m,m)-fzf(m,m)+fff(m,m))/(j*w0*epth0); #_{3W-35*  
        else ]Y;EIn  
&n .. h^ ex?  
^-T!(P:  

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

导线长 / 波长  << 1，在做计算仿真前，就可以知道电流大概分布是什么样的了，电流分布肯定都是差不多的，都是底部最大，然后减少，顶端为零的曲线。 但这不是半波了，半波的话，你要确信在底部的电流梯度为零才行。

还是不太明白。我以为半个正弦波就是半波了呢。看来我还需要再好好学习一下。谢谢。

半个正弦波是半波，但你看到的不是半个了。你看一下幅度分布是不是和我画的差不多。 *fs'%"w-  
比如简单情况下，电流分布是 I(z) = I0 * (exp(i*kz) - exp(i*k(2*l - z))) ，幅度等于I0 * sin(k*(l - z))，因为l << 波长，你可以写个程序，把你的实际参数代入到这个表达式中，画一下结果看看，你看样子是不是差不多。

引用第3楼gwzhao于2008-09-02 16:40发表的  : [
K(|V  
半个正弦波是半波，但你看到的不是半个了。你看一下幅度分布是不是和我画的差不多。 p i;,?p-  
比如简单情况下，电流分布是 I(z) = I0 * (exp(i*kz) - exp(i*k(2*l - z))) ，幅度等于I0 * sin(k*(l - z))，因为l << 波长，你可以写个程序，把你的实际参数代入到这个表达式中，画一下结果看看，你看样子是不是差不多。

tbx* }uy2  
非常感谢您！ %/b3G*$W  
我照您写的式子把参数输进去的确像您画的那样。 g;N)K3\2
 
但是我写程序得到的图却是正弦波的一半，不知是哪错了。我再好好看看，也希望您多指教。谢谢！

好帖子，不过显然是程序有问题。

引用第5楼plantice于2008-09-03 00:06发表的  : :krdG%r  
好帖子，不过显然是程序有问题。

ej^pFo  
谢谢您！我因为是高电压专业的，对天线这些问题也是因为要用才临时补充了一点，所以理解肯定只是表面的。还请您多指教。 *D_pFS^l  
程序中的电场dEdl是在我的程序中先算出来的，然后用来计算处于其中的直导线上感应的电流。我抓紧时间再好好学习一下，希望各位高手不辞吝教，谢谢！

啊，真的是正弦波的一半啊，波长差那么多，长度这么小，想不通啊。要不，你把代码贴上来看看？ {0t-Q k  

引用第4楼xiaoyudida于2008-09-02 20:24发表的  :
非常感谢您！ v.~uJ.T  
我照您写的式子把参数输进去的确像您画的那样。 Z<2j#rd  
但是我写程序得到的图却是正弦波的一半，不知是哪错了。我再好好看看，也希望您多指教。谢谢！

引用第7楼gwzhao于2008-09-03 09:19发表的  : 
G/w&yd4  
啊，真的是正弦波的一半啊，波长差那么多，长度这么小，想不通啊。要不，你把代码贴上来看看？ *A>I)a<:  

rTgCmr'&  
谢谢您的帮助。代码我也贴上去了。

我觉得出现半个波是因为在求解时设定的边界条件，即导线两端的电流为0。这样中间的电流大于0所以就出现这样的波形了。但是我看得文献上求参数的时候不都是令两端的电流为0吗？