| 查看: 489 | 回复: 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. 这是错误提示 |
回复此楼» 猜你喜欢
求标准粉末衍射卡号 ICDD 01-076-1802
已经有0人回复
-
新西兰Robinson研究所招收全奖PhD
已经有0人回复
-
物理学I论文润色/翻译怎么收费?
已经有226人回复
-
石墨烯转移--二氧化硅衬底石墨烯
已经有0人回复
-
笼目材料中量子自旋液体基态的证据
已经有0人回复
-
数学教学论硕士可以读数学物理博士吗？
已经有0人回复
-
德国亥姆霍兹Hereon中心汉堡分部招镁合金腐蚀裂变SCC课题方向2026公派博士生
已经有4人回复
-
澳门大学 应用物理及材料工程研究院 潘晖教授课题组诚招博士后
已经有11人回复
-
求助NH4V4O10晶体的CIF文件
已经有0人回复
-
英国全奖博士招聘-深度学习与量子物理
已经有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人回复
点击这里搜索更多相关资源科研从小木虫开始，人人为我，我为人人