Vt=12;
V=120;
K=4;
r0=3000;
(我随便改了下,不知道对了没,图形是条直线)
q0=7*pi/18;
eta0=-pi/90;
syms r q psaiV eta t;
d1=diff(r);
d2=diff(q);
d3=diff(psaiV);
e1=sym('d1=-Vt*cos(q)-V*cos(psaiV)');
e2=sym('r*d2=Vt*sin(q)-V*sin(psaiV)');
e3=sym('psaiV=q-eta');
e4=sym('d3=K*d2');
r=solve(e1);
q=solve(e2);
eta=solve(e3);
psaiV=solve(e4);
d1=-Vt*cos(q)-V*cos(psaiV);
figure;
fplot('r',[0,5000])
xlabel('t-时间');
ylabel('r-距离');