| 查看: 535 | 回复: 0 | ||
[求助]
求助，微分代数方程组
|
|
function dx = yuanpan4(t,x) m=10;r=0.5;k=100;g=9.8;u=0.4; dx=zeros(5,1); if x(4)>0 dx(1)=x(2); dx(2)=(-k*x(1)+x(5))/m; dx(3)= (-2*x(5))/(m*r); dx(4)=x(2)-r*x(3)-x(4); dx(5)= x(5)+m*g*u; elseif x(4)==0 dx(1)=x(2); dx(2)=(-k*x(1)+x(5))/m; dx(3)=(-2*x(5))/(m*r); dx(4)=x(2)-r*x(3)-x(4); dx(5)=0; else dx(1)=x(2); dx(2)=(-k*x(1)+x(5))/m; dx(3)=(-2*x(5))/(m*r); dx(4)=x(2)-r*x(3)-x(4); dx(5)=x(5)-m*g*u; end end 这是函数。 tspan=[0,16]; x0=[3;3;5;0.5;-39.2]; M=diag([1,1,1,0,0]); options=odeset('Mass',M,'RelTol',1e-6,'AbsTol',[1e-10,1e-10,1e-10,1e-10,1e-10]); [t,x]=ode15s(@yuanpan4,tspan,x0,options);plot(t,x(:,1),'-b',t,x(:,4),'-r',t,x(:,5),'-g'); grid on; 这是主程序。 Warning: Failure at t=1.192445e-02. Unable to meet integration tolerances without reducing the step size below the smallest value allowed (2.775558e-17) at time t. 这是错误提示 |
回复此楼» 猜你喜欢
如何从铁电相到顺电相。
已经有1人回复
-
各位大佬，求一份最新Wien2k版本作为学习用途
已经有0人回复
-
物理学I论文润色/翻译怎么收费?
已经有146人回复
-
南科大活性流体和软物质课题组诚招2027级博士、硕士研究生和博后
已经有42人回复
-
瑞典林雪平大学博士后招聘｜PEC Water Splitting 方向
已经有24人回复
-
PRB投稿，Acknowledgment sent to author状态十天了
已经有3人回复
-
PVA溶解
已经有2人回复
-
基金申请
已经有44人回复
-
CSC与新西兰维多利亚大学PhD奖学金项目
已经有0人回复
-
新西兰Robinson研究所 招聘CSC公派访问人员
已经有0人回复
-
帮我的英语口语老师找学生
已经有0人回复
找到一些相关的精华帖子，希望有用哦~
帮忙：求解一组微分代数方程，奖励100金币
已经有7人回复
- matlab 解非线性偏微分方程组
已经有6人回复
- 请问这个偏微分方程组有一般的解法么？
已经有12人回复
- 求解线性代数的问题
已经有5人回复
- matlab求三元4次符号代数方程组，出现warning且答案不正确，该怎么办？
已经有12人回复
- matlab 非线性微分方程求解
已经有3人回复
- 求推荐求解偏微分方程组的书
已经有3人回复
- 搞传热的虫子求助：遇到有变量相乘的偏微分方程组如何离散？
已经有16人回复
- 求助-帮我用mathematics计算下下面的方程组
已经有13人回复
- matlab解偏微分方程求助
已经有12人回复
- 【求助】matlab求解偏微分方程
已经有4人回复
- 求助matlab拟合方程组
已经有12人回复
- 偏微分方程数值计算书籍推荐
已经有10人回复
- 第5章_解线性方程组的直接方法.ppt
已经有13人回复
- matlab求解非线性方程组
已经有16人回复
- matlab数值求解边界条件微分方程组
已经有7人回复
- 【求助】有限元法求解二阶偏微分方程组
已经有6人回复
点击这里搜索更多相关资源科研从小木虫开始，人人为我，我为人人
20