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

"=9)|{=m  
y Xi$w.gr  
%%  电场平均 97VS xhr  
for i=1:length g0RfvR  
    ez_outave(:,i)=(ez_out(:,i)+ez_out(:,i+1))/2; :D=y<n;S+  
end R6G%_,p$7  
[b@9V
_  
%初始化 h/#s\>)T  
current=zeros(4,length);    %  电流z方向 P Z+Rz1x  
magnetic_x=zeros(4,length); %  磁流x方向 0 oHnam  
magnetic_y=zeros(4,length); %  磁流y方向 a]R1Fi0n  
$
C6O<A  
%%  电流z 6XGqZ!2  
current(1,:)=hx_out(1,:);   %  下 `~ R%}ID  
current(2,:)=-hx_out(2,:);  %  上 ldWr-  
current(3,:)=-hy_out(3,:);  %  左 " G0HsXi  
current(4,:)=hy_out(4,:);   %  右 p87VJ}  
ro@`S:  
%%  磁流x
I~7eu&QZ  
magnetic_x(1,:)=ez_outave(1,:);   %  下 Jgf=yri  
magnetic_x(2,:)=-ez_outave(2,:);  %  上 +?bOGU
ik  
magnetic_x(3,:)=0;                %  左 Dad*6;+N  
magnetic_x(4,:)=0;                %  右 v iM6q<Ht  
q,_ 1?A)
 
%%  磁流y /;rk-I  
magnetic_y(1,:)=0;                  %  下 y Ide]  
magnetic_y(2,:)=0;                  %  上 FtxmCIVIV~  
magnetic_y(3,:)=-ez_outave(3,:);    %  左 e/
V8lo  
magnetic_y(4,:)=ez_outave(4,:);     %  右 ~ny4Ay$#  
%!Ak]|[7  
%%  外推边界坐标X,Y  描述一个外推边界 2sXX0kq~V  
posx=zeros(4,length); %初始化 /'>#1J|TlK  
posy=zeros(4,length); cvQMZ,p  
4w0Y(y  
posx(1,:)=linspace(pa*ddx,pb*ddx,length)+ddx/2;   %  下 P/hIJV[  
posx(2,:)=linspace(pa*ddx,pb*ddx,length)+ddx/2;   %  上 KdkZ-.  
posx(3,:)=pa*ddx;  %  左 VtC1TZ3-7  
posx(4,:)=pb*ddx;  %  右 x-ShY&k  
,;-55|o\V  
posy(1,:)=qa*ddx;  %  下 C<\O;-nHH  
posy(2,:)=qb*ddx;  %  上 F /% 5 r{  
posy(3,:)=linspace(qa*ddx,qb*ddx,length)+ddx/2;  %  左 POc-`]6<F  
posy(4,:)=linspace(qa*ddx,qb*ddx,length)+ddx/2;  %  右 EX^j^#N  
/YH5s=  
%%  归一化fz,fmx,fmy求解 % qUF}rlS=r  
fz=ddx*current/abs(inh);   iMr/i?`i  
fmx=ddx*magnetic_x/abs(inh); lzStJ,NPqn  
fmy=ddx*magnetic_y/abs(inh); bcZ s+FOPd  
L;k9}HWpP  
%%  RCS求解; [#6Eax,j  
co=0; 3~^}R  
for phi=pi/2:pi/180:pi+pi/2; [1gWc`#  
    co=co+1; A!j6JY.w  
    ss=0;
<eG8xC  
    for t=1:4; R} aHo0r  
        for s=1:length uWKc .  
            kr=posx(t,s)*cos(phi)+posy(t,s)*sin(phi); kkE1
CHY  
            ss=ss+(-fz(t,s)-fmx(t,s)*sin(phi)+fmy(t,s)*cos(phi))*exp(sqrt(-1)*k*kr);  %%  外推公式 (E(k
w="  
        end; a).bk!G  
    end; <xQHb^:  
    rcs(co)=k/4*(abs(ss))^2; !<xeAo%8  
end; y%x2  
rcs=10*log10(rcs/(lamda)); hM @F|t3  
f_GqJ7Gk]  
figure; `?R~iLIAq  
plot(rcs,'g'); ?Wm.'S'to  
xlabel('Angel'),ylabel('RCS'); 5//.q;z  
1]vrpJw  
toc; wWR9dsB.;  
Gdf*x<T1  
N!&$fhY)  
VXtW{*{"  
这是我在论坛上看到的有外推的程序。。。想问下在RCS计算那里，其中的角度的问题。。。 6?v)Hb}J%d  
（1） F2$Z4%x#  
kr=posx(t,s)*cos(phi)+posy(t,s)*sin(phi); _e3kO6X  
            ss=ss+(-fz(t,s)-fmx(t,s)*sin(phi)+fmy(t,s)*cos(phi))* .. U[UjL)U  
L0O},O  

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

没人回答么？  帮帮忙啊。。。。 U!L<v!$  
# 4AyA$t  
for phi=pi/2:pi/180:pi+pi/2; 那么这里的phi应该是葛书中的远区观察点的方位角。 ugtb`d{ Sl  
可这角度变化，不是应该是从0至pi么？这里却是90度到270度？

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