仅供参考
syms A1 E1 A2 E2 %测站1、2测得的A角和E角
x1=400;y1=0;z1=1000;
x2=600;y2=0;z2=-1000;
yt=600;zt=500;%当xt是向量时,
i=0;
for xt=0:3
i=i+1;
eq1=x1-xt+((x1-x2).*tan(A2)-(z1-z2))./...
(tan(A1)-tan(A2));
eq3=z1-zt+((x1-x2).*tan(A2)-(z1-z2)).*tan((A1))./...
(tan((A1))-tan((A2)));
[A1s,A2s]=solve([eq1,eq3],[A1,A2]);
A1ss(i)=double(A1s);
A2ss(i)=double(A2s);
end
A1ss