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

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

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

)7<JGzBZ1  
Y值为负的那几个点是异常值 nmn$$=~)  
Q1mz~r  
clc tQ< ou,  
clear K 4j'e6  
close all kG[u$[B  
nx = 0:0.6:150; 9w[7X"#n  
n = length(nx)-1; AFGWlC#`  
figure t/yGMR=  
xz = []; A-aukJg9  
yz = []; ;hA>?o_i(  
for ii = 1:n H2 5Mx>|d  
    f=300; / *Z(;-  
    h=150-nx(ii); K%P$#a  
    b=atan(h/f); 1"RO)&  
    d=(b/pi)*180 9n]zh-  
    syms c dms R>Q  
   Y=solve(sin(c-b)-1.494*sin(3*c-pi)==0,c); C|5eV=f)P  
   C=double(Y); ` :eXXE  
   e=(C/pi)*180; 1Y$ gt  
    x=0:0.1:0.6; 6AKH0t|4  
    y=(0.6-x)*tan(C); *F1!=:&s  
    xz = [xz x+nx(ii)] (s4w0z  
    yz = [yz y]; zGHP{a1O7  
    plot(xz,yz,'k-') c
x_FtD  
    hold on oJK1~;:  
    pause(0.01) {1li3K&0s  
end 8G; t[9  
L(XGD  
'e_^s+l)a  
biKom|<nm  
lZ.x@hDS