程序有问题,在140.4、139.8、132等点处有异常值,请教一下大家怎么改,应该是Y=…那一步的式子有问题,谢谢
BvsSrse
e6>G8d Y值为负的那几个点是异常值
OxYAM,F DnFzCJ clc
/IG{j} clear
,f2oO?L} close all
Q"ZpT nx = 0:0.6:150;
4~&3.1 n = length(nx)-1;
a_V\[V{R= figure
.wD
$Bsm`t xz = [];
>x JzV yz = [];
":8\2Qp for ii = 1:n
;aV3j/ f=300;
9xO@_pkX h=150-nx(ii);
@<{%r b=atan(h/f);
kqm(D# d=(b/pi)*180
DH
yv^ syms c
Q(]m1\a Y=solve(sin(c-b)-1.494*sin(3*c-pi)==0,c);
k9f|R*LM C=double(Y);
h@Ea5x e=(C/pi)*180;
jwT` Z x=0:0.1:0.6;
zQ|2D*W y=(0.6-x)*tan(C);
`VKFA<T xz = [xz x+nx(ii)]
R\5,H!V9n yz = [yz y];
?a3wBy plot(xz,yz,'k-')
\^*:1=|7u] hold on
&J&'J~N pause(0.01)
*:@KpYWx" end
o 2Nu@^+ :31_WJ^ e{@RBYX@+c eO <N/?t m2\\!C]f