小球空中运动方程,只要稍微有一点高中物理和高等数学基础的朋友应该都可以看懂这个方程吧:
1.gif (7.12 KB, 下载次数: 28)
2011-3-16 10:17 上传
上面是相当简单的一个常微分方程组,要求解这个方程组也相当容易,MATLAB提供的ode45函数足以胜任。
编写主函数main.m
function main()
global m mu g
m=1; % 小球质量
mu=0.3; % 摩擦系数
V0=100; % 初始速度
theta=pi/6; % 初始角度
g=9.8; % 重力加速
k=0.8; % 撞击系数
H=10; % 初始高度
tspan=[0 10]; % 求解时间范围
x0=[0,V0*cos(theta),H,V0*sin(theta)]; % 微分方程初值
[t,x]=ode45(@odefun,tspan,x0);
X=x(:,1);
Vx=x(:,2);
Y=x(:,3);
Vy=x(:,4);
figure
plot(t,x)
legend('X','Vx','Y','Vy')
xlabel('时间')
ylabel('变量')
grid on
figure
plot(X,Y)
xlabel('水平方向')
ylabel('竖直方向')
title('小球运动轨迹')
grid on
继续编写微分方程组描述函数odefun.m
function dx=odefun(t,x)
% 微分方程描述函数
global m mu g
% X=x(1);
Vx=x(2);
% Y=x(3);
Vy=x