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

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

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

B"zB=Aw  
Y值为负的那几个点是异常值 ;xiwyfqgE  
!v(^wqna\  
clc dwQ1~  
clear <*WGvCh%w  
close all KVh#"]<WV  
nx = 0:0.6:150; 9X,iQ  
n = length(nx)-1; 1uppE|  
figure `6lOqH  
xz = []; i7Cuc+j8  
yz = []; /.=aA~|  
for ii = 1:n ![nL/  
    f=300; ];"40/X  
    h=150-nx(ii); .6LRg  
    b=atan(h/f); cTS.yN({G  
    d=(b/pi)*180 5IOGH*'U8  
    syms c 9D Nd} rXO  
   Y=solve(sin(c-b)-1.494*sin(3*c-pi)==0,c); D4C:%D  
   C=double(Y); Z4c'1-lh  
   e=(C/pi)*180; LokH4A17U  
    x=0:0.1:0.6; 9_nbMs  
    y=(0.6-x)*tan(C); K~2sX>l  
    xz = [xz x+nx(ii)] &3;"$P  
    yz = [yz y]; 1KbZ6Msy  
    plot(xz,yz,'k-') 93:oXyFjD  
    hold on \uZ|2WG`  
    pause(0.01) !icI Rqcf=  
end 2K{'F1"RM  
_ E-\aS{  
TnKv)%VF  
LMx/0  
Yh!=mW!OY