原方程为m*(d^2x/dt^2)=-c(Vw+dx/dt)^2,dx/dt(t=0)=v0*cosA
x(0)=0
令x1=x,x2=dx1/dt,则systematic form为：
dx2/dt=-(c/m)*(Vw+x2)^2，x1（0）=0， x2（0）=vocosA
2>,y（t）是位移函数，g是重力加速度，沿用上述常量，二阶方程为：d^y/dt^2=-g, dy/dt(t=0)=v0*sinA ,y(0)=0 令y=y1，y2=dy1/dt，则systematic form为：
dy2/dt=-g，y1（0）=0，y2（0）=vo*sinA
------------------------------------------------------
options=optimset('display','off');
[t,x]=ode45('f1',[t0,t1],[0,v0*cosA ],options,c,Vw,v0,A,m) %t0是初始时间，t1是终止时间。
--------------------------------------------
function doty=f1(t,x,flag,c,Vw,v0,A,m)
doty=zeros(2,1);
doty=[x(2);-c.*(Vw+x(2)).^2./m];
第二个方程解法类似，你有什么疑问可以看看help
|