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

程序有问题，在140.4、139.8、132等点处有异常值，请教一下大家怎么改，应该是Y=…那一步的式子有问题，谢谢 >0$5H]1u  

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

|a3b2x,  
Y值为负的那几个点是异常值 ?!vW&KJZx  
XRin~wz|S  
clc HX[#tT|m~  
clear ?RyvM_(N6  
close all ymqhI\>y#  
nx = 0:0.6:150; CjCnh7tm  
n = length(nx)-1; ]sE^=;Pv?  
figure '?4[w]0J<  
xz = []; .!! yj,bQz  
yz = []; v<**GW]neD  
for ii = 1:n {[dqXG$v `  
    f=300; yK;I<8+>_  
    h=150-nx(ii); KcE=m\h  
    b=atan(h/f); Nrh`DyF0D!  
    d=(b/pi)*180 _l<"Qqt  
    syms c ~a Rq\fx{  
   Y=solve(sin(c-b)-1.494*sin(3*c-pi)==0,c); {:cA'6f.b  
   C=double(Y); ?,[w6O*  
   e=(C/pi)*180; wVB8PO8  
    x=0:0.1:0.6; x;/3_"$9>\  
    y=(0.6-x)*tan(C); B7C6Mau  
    xz = [xz x+nx(ii)] XO>Y*7rO  
    yz = [yz y]; 7Q&P4{hi0  
    plot(xz,yz,'k-') oKiu6=  
    hold on zyE yZc?  
    pause(0.01) "tIx$?I  
end FeJ5^Gh.  
#B
[>\D"*  
MvA_tRO  
#t*c*o  
@(L|