syms x positive A=0; B=-0.052; C=-0.999; Q=[A^2 A*B A*C;A*B B^2 B*C; A*C B*C C^2]; P=[0 -C B;C 0 -A;-B A 0]; D=3*pi/180; H=[1 0 0;0 cos(D) -sin(D);0 sin(D) cos(D)]; I=[cos(x) -sin(x) 0;sin(x) cos(x) 0;0 0 1]; J=[1 0 0;0 cos(D) sin(D);0 0 1]; R=H*I*J; R=vpa(R); K=[300;0;0]; L=R*K; L=vpa(L); M=[204.811;1311.088;11.395]; N=[0;1350;0]; O=[1 0 0;0 1 0;0 0 1]; E=(L-N)'*(O-Q)*(M-N); F=(L-N)'*(P)*(M-N); G=(L-N)'*(Q)*(M-N)+1/2*((K-M)'*(K-M)-(L-N)'*(L-M)-(M-N)'*(M-N)); S=E^2+F^2-G^2; T=S^0.5; U=-F-T; V=U/(G-E); y=atan(V); x=0:pi/180:10*pi/180; plot(x,y) |
最佳答案
2 条回复