[求助]matlab编程绘图，三角函数，菲涅尔 [复制链接]

程序有问题，在140.4、139.8、132等点处有异常值，请教一下大家怎么改，应该是Y=…那一步的式子有问题，谢谢 tUW^dGo.  

图片:})WX7H~ND~3OE[$Y(~WEX_R.png

}sJ}c}b  
Y值为负的那几个点是异常值 @MoCEtt  
Oiqc]4TL  
clc *b!.9pK  
clear PR AP~P&^  
close all 7q 5 \]J[  
nx = 0:0.6:150; uZ@qlq8  
n = length(nx)-1; [} d39  
figure lPC{R k.\C  
xz = []; ^^24a_+2  
yz = []; -UAMHd}4  
for ii = 1:n DHyQ:0q  
    f=300; \d:Uq5d)0  
    h=150-nx(ii); wlh%{l  
    b=atan(h/f); jq_4x[  
    d=(b/pi)*180 *lu*h&Y  
    syms c UF\k0oLz  
   Y=solve(sin(c-b)-1.494*sin(3*c-pi)==0,c); ,J<+Wxz  
   C=double(Y); 0B2f[A  
   e=(C/pi)*180; tu@-+<*  
    x=0:0.1:0.6; aI}htb{m`  
    y=(0.6-x)*tan(C); M1k_ldP  
    xz = [xz x+nx(ii)]  BPKrRex  
    yz = [yz y]; OwgPgrV  
    plot(xz,yz,'k-') /qq*"R  
    hold on - K"L6m|  
    pause(0.01) _5p]Arg?}&  
end V`i(vC(  
"jq6FT)O  
X4Uy3TV>  
Se\iMs  
q[Ey!h)xq