非线性方程式有自变量的分段函数该如何用ode函数求解

eatnche594

这是我已经降阶好的运动方程式

function  yp=ivp13(t,y)

yp=zeros(1,16);  % note: column vector
yp(1)=y(2);
yp(2)=((6.15*10.^6*(y(3)-y(1)))-(6.15*10.^6*y(1)))/(13.2)-y(16);
yp(3)=y(4);
yp(4)=((8.1*10.^7*(y(5)-y(3)))-(6.15*10.^6*(y(3)-y(1))))/(1.3)-y(16);
yp(5)=y(6);
yp(6)=(6.03*10.^6*(y(7)-y(5))-8.1*10.^7*(y(5)-y(3)))/(2.28)-y(16);
yp(7)=y(8);
yp(8)=(8.1*10.^7*(y(9)-y(7))-6.03*10.^6*(y(7)-y(5)))/(1.12)-y(16);
yp(9)=y(10);
yp(10)=((V+0.5*A*y(9)+0.5*J*y(9).^2)*8.1*10.^7.*(D+y(11)-(D+V*y(9)+0.5*A*y(9).^2)))-2.01*10.^6*(y(9)-y(7))/(3.66)-y(16);
yp(11)=y(12);
yp(12)=(2.01*10.^6*(y(13)-y(11))-2.83*10.^8.*(D+y(11)-(D+V*y(9)+0.5*A*y(9).^2)))/(0.35)-y(16)*V-A*y(15).^2;
yp(13)=y(14);
yp(14)=(1000-2.01*10.^6*(y(13)-y(11)))/(5)-y(16)*V-A*y(15).^2;
yp(15)=y(16);
y(16)=(107.257-D*(0.25+0.104+5)*(V*y(16).^2+1000))/(0.096+13.2+1.3+2.28+1.12+3.66+((0.25+0.104+5)*D.^2));
yp=yp';

其中D,V,A,J是有自变量的分段函数,如附件 Appendix_A.doc (56.5 KB, 下载次数: 0)

2015-9-20 03:26 上传

点击文件名下载附件

想请问大侠，侠女我该如何用ode函数求解

給定条件为
tspan=[0 0.4];
x0=[0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0];
求解yp(16),yp(14)

麻煩了 感激不尽
下载次数: 1
2007-8-21 16:16

[ 本帖最后由 花如月 于 2007-8-22 10:27 编辑 ]