[求助]咨询一下关于外推的程序问题

F"_SCA?9?  
"&An9H'  
%%  电场平均 U_+>4zdm  
for i=1:length *
5 5yF`  
    ez_outave(:,i)=(ez_out(:,i)+ez_out(:,i+1))/2; Gg_i:4F  
end 3a =KgOvp  
Gk<h_1WW
K  
%初始化 z7}@8F  
current=zeros(4,length);    %  电流z方向 nlK
WZYv  
magnetic_x=zeros(4,length); %  磁流x方向 MD&Ebq5V  
magnetic_y=zeros(4,length); %  磁流y方向 ,+f'%)s_x  
3]T2Zp&;  
%%  电流z 3*
ZE``  
current(1,:)=hx_out(1,:);   %  下 mwBOhEefNJ  
current(2,:)=-hx_out(2,:);  %  上 yB&s2J  
current(3,:)=-hy_out(3,:);  %  左 C{4[7  
current(4,:)=hy_out(4,:);   %  右 8[eH8m#~$  
`6'fX[j5  
%%  磁流x }I!D65-#'  
magnetic_x(1,:)=ez_outave(1,:);   %  下 J?V8u
Ely  
magnetic_x(2,:)=-ez_outave(2,:);  %  上 /Ue
~W,|  
magnetic_x(3,:)=0;                %  左 ~G"5!,J  
magnetic_x(4,:)=0;                %  右 R{/nlS5  
r'"H8>UZ%  
%%  磁流y 4Xlq Ym  
magnetic_y(1,:)=0;                  %  下 `Q,03W#GJ%  
magnetic_y(2,:)=0;                  %  上 4]B(2FR[8  
magnetic_y(3,:)=-ez_outave(3,:);    %  左 ,?728pfw  
magnetic_y(4,:)=ez_outave(4,:);     %  右 -(dc1?COi  
zUWu5JI  
%%  外推边界坐标X,Y  描述一个外推边界 wTG6>l]H  
posx=zeros(4,length); %初始化 w18kTa!4@  
posy=zeros(4,length); kd2+k4@#  
C"K(-/  
posx(1,:)=linspace(pa*ddx,pb*ddx,length)+ddx/2;   %  下 EP*"=_  
posx(2,:)=linspace(pa*ddx,pb*ddx,length)+ddx/2;   %  上 R'He(x  
posx(3,:)=pa*ddx;  %  左 PKk_9Xd  
posx(4,:)=pb*ddx;  %  右 -B'<*Y  
Que-  
posy(1,:)=qa*ddx;  %  下 y<MXd,eE  
posy(2,:)=qb*ddx;  %  上 zl$'W=[rFs  
posy(3,:)=linspace(qa*ddx,qb*ddx,length)+ddx/2;  %  左 nF]lSg&]X  
posy(4,:)=linspace(qa*ddx,qb*ddx,length)+ddx/2;  %  右 kLK}N>
v}X  
=98@MX%P  
%%  归一化fz,fmx,fmy求解 % n"*A.  
fz=ddx*current/abs(inh);   +eQg+@u  
fmx=ddx*magnetic_x/abs(inh); "4Cb dD//  
fmy=ddx*magnetic_y/abs(inh); "??$yMW  
Ahm*_E2E  
%%  RCS求解; YjAwt;%-D  
co=0; f/aSqhAW  
for phi=pi/2:pi/180:pi+pi/2; ;BsyN[bF  
    co=co+1; n1X7T0'  
    ss=0; EHmw(%a|+  
    for t=1:4; h?jKq2`  
        for s=1:length z^b
v)u  
            kr=posx(t,s)*cos(phi)+posy(t,s)*sin(phi); Zgd| J T7  
            ss=ss+(-fz(t,s)-fmx(t,s)*sin(phi)+fmy(t,s)*cos(phi))*exp(sqrt(-1)*k*kr);  %%  外推公式 (X;D.
s  
        end; >yiK&LW^?  
    end; qq G24**9v  
    rcs(co)=k/4*(abs(ss))^2; um~U_&>  
end; @uApm~}  
rcs=10*log10(rcs/(lamda)); U*b1yxt  
)cxLpTr  
figure; qXcHf6  
plot(rcs,'g'); f5'Cq)Vw_  
xlabel('Angel'),ylabel('RCS'); B*W)e$  
_XvSe]`f`  
toc; <'f+nC=2  
[f 4Nq \i  
+G)L8{FY(  
E
VX*YGxx6  
这是我在论坛上看到的有外推的程序。。。想问下在RCS计算那里，其中的角度的问题。。。 yL{X}:;}  
（1） as6a)t.^  
kr=posx(t,s)*cos(phi)+posy(t,s)*sin(phi);  [,JUC<  
            ss=ss+(-fz(t,s)-fmx(t,s)*sin(phi)+fmy(t,s)*cos(phi))* .. 7,X5]U&A<x  
ZjLuqo  

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

没人回答么？  帮帮忙啊。。。。 09_5niaz[  
rb_FBa%  
for phi=pi/2:pi/180:pi+pi/2; 那么这里的phi应该是葛书中的远区观察点的方位角。 Q |1-j  
可这角度变化，不是应该是从0至pi么？这里却是90度到270度？

观察方位角是你设定的，也可以是0～360，这由你的需要决定。另外，还要考虑入射方向。