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

9@ 4]t6h[  
;GE6S{~-  
%%  电场平均 dU*$V7  
for i=1:length "n{';Q)  
    ez_outave(:,i)=(ez_out(:,i)+ez_out(:,i+1))/2; -Bq]E,Xf)  
end qzNXz_#+u  
!rs }83w!  
%初始化 n=)LB& m  
current=zeros(4,length);    %  电流z方向 <)pPq+  
magnetic_x=zeros(4,length); %  磁流x方向 pTaC$Ne  
magnetic_y=zeros(4,length); %  磁流y方向 15VvZ![$V  
lI%RdA[  
%%  电流z }.(DQwC}1k  
current(1,:)=hx_out(1,:);   %  下 dWi<U4  
current(2,:)=-hx_out(2,:);  %  上 ]#[4eaCg  
current(3,:)=-hy_out(3,:);  %  左 "@5qjLz]  
current(4,:)=hy_out(4,:);   %  右 eIy:5/s  
w:9`R<L  
%%  磁流x ju5o).!bg  
magnetic_x(1,:)=ez_outave(1,:);   %  下  _X  
magnetic_x(2,:)=-ez_outave(2,:);  %  上 f[vm]1#  
magnetic_x(3,:)=0;                %  左 Y4w
]jIv  
magnetic_x(4,:)=0;                %  右 KwV!smi2  
r@zs4N0WP  
%%  磁流y E=8GSl/Jx  
magnetic_y(1,:)=0;                  %  下 {\`y)k 7  
magnetic_y(2,:)=0;                  %  上 & kVa*O  
magnetic_y(3,:)=-ez_outave(3,:);    %  左 A9g/At_
  
magnetic_y(4,:)=ez_outave(4,:);     %  右 e|W;(@$<  
$tF\7.e@  
%%  外推边界坐标X,Y  描述一个外推边界 ]8R@2L3s  
posx=zeros(4,length); %初始化 T{2)d]Y  
posy=zeros(4,length); Q>$l
f.)  
C:t?HLY)fG  
posx(1,:)=linspace(pa*ddx,pb*ddx,length)+ddx/2;   %  下 U.V/JbX
X  
posx(2,:)=linspace(pa*ddx,pb*ddx,length)+ddx/2;   %  上 ps?su`  
posx(3,:)=pa*ddx;  %  左 g_n_Qlo  
posx(4,:)=pb*ddx;  %  右 *3P+K:2lNG  
Niu |M@  
posy(1,:)=qa*ddx;  %  下 V;g) P  
posy(2,:)=qb*ddx;  %  上 daP_Kz/2K  
posy(3,:)=linspace(qa*ddx,qb*ddx,length)+ddx/2;  %  左 .N>*+U>>P  
posy(4,:)=linspace(qa*ddx,qb*ddx,length)+ddx/2;  %  右 `\|@w@f|;  
oOc-1C y  
%%  归一化fz,fmx,fmy求解 % dl3;A_ 2  
fz=ddx*current/abs(inh);   q6*i/"mN*  
fmx=ddx*magnetic_x/abs(inh); $a8,C\me?  
fmy=ddx*magnetic_y/abs(inh); FZ}^)u}o  
K2e68GU  
%%  RCS求解; FZ RnIg  
co=0; N@Oe[X8  
for phi=pi/2:pi/180:pi+pi/2; 5
~Ek_B  
    co=co+1; AbhR*  
    ss=0; WO{7/h</  
    for t=1:4; ~(^pGL3<  
        for s=1:length U}4I29M  
            kr=posx(t,s)*cos(phi)+posy(t,s)*sin(phi); <KK.f9^o(  
            ss=ss+(-fz(t,s)-fmx(t,s)*sin(phi)+fmy(t,s)*cos(phi))*exp(sqrt(-1)*k*kr);  %%  外推公式 w
x`.  
        end; pEz^z9  
    end; qou\4YZ  
    rcs(co)=k/4*(abs(ss))^2; tBfmjxv  
end; `1n^~  
rcs=10*log10(rcs/(lamda)); ji>LBbnHdE  
_;0RW  
figure; pH.&C 5kA  
plot(rcs,'g');
vE8BB$D  
xlabel('Angel'),ylabel('RCS'); mn=b&{')e  
/F9Dg<#a  
toc; {}8C/4iP  
'^C *%"I]  
;~L,Aqn7  
si)>:e  
这是我在论坛上看到的有外推的程序。。。想问下在RCS计算那里，其中的角度的问题。。。 Nd"IW${Kg  
（1） T)%6"rPL3!  
kr=posx(t,s)*cos(phi)+posy(t,s)*sin(phi); <,0/BMz  
            ss=ss+(-fz(t,s)-fmx(t,s)*sin(phi)+fmy(t,s)*cos(phi))* .. A`JE(cIz3  
?%8%1d  

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

没人回答么？  帮帮忙啊。。。。 gl`J(  
*A8*FX>\F  
for phi=pi/2:pi/180:pi+pi/2; 那么这里的phi应该是葛书中的远区观察点的方位角。 &}Wi@;G]2
 
可这角度变化，不是应该是从0至pi么？这里却是90度到270度？

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