声振论坛

 找回密码
 我要加入

QQ登录

只需一步,快速开始

查看: 1240|回复: 0

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

[复制链接]
发表于 2007-8-21 17:58 | 显示全部楼层 |阅读模式

马上注册,结交更多好友,享用更多功能,让你轻松玩转社区。

您需要 登录 才可以下载或查看,没有账号?我要加入

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

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)

想请问大侠,侠女我该如何用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 编辑 ]
回复
分享到:

使用道具 举报

您需要登录后才可以回帖 登录 | 我要加入

本版积分规则

QQ|小黑屋|Archiver|手机版|联系我们|声振论坛

GMT+8, 2024-12-12 05:32 , Processed in 0.063944 second(s), 22 queries , Gzip On.

Powered by Discuz! X3.4

Copyright © 2001-2021, Tencent Cloud.

快速回复 返回顶部 返回列表