|
马上注册,结交更多好友,享用更多功能,让你轻松玩转社区。
您需要 登录 才可以下载或查看,没有账号?我要加入
x
clear
global m b g ;
m=10;b=10;g=9.8;
vx=100*cos(pi*35/180);
vy=100*sin(pi*35/180);
tspan=[0 7];
y0=[0 vx 0 vy]';
[t,y]=ode23('ptyd',tspan,y0);
plot(y(:,1),y(:,3),'b');
xlabel('距离/s');ylabel('高度/m');
grid
h2=figure
plot(t,y(:,1),'b');
xlabel('时间/s');ylabel('距离/m');
grid
h3=figure
plot(t,y(:,3),'b');
xlabel('时间/s');ylabel('高度/m');
grid
h4=figure
plot(t,y(:,2),'b');
xlabel('时间/s');ylabel('水平速度/m');
grid
function yp=ptyd(t,y)
global m b g;
yp=[y(2)-b*y(2)/my(4)-g-b*y(4)/m]';
Line: 24 Column: 1
Function definitions are not permitted at the prompt or in scripts. |
|