| 查看: 493 | 回复: 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. 这是错误提示 |
回复此楼» 猜你喜欢
基于基元I统一理论的量子力学本源推导
已经有1人回复
-
推荐一款可以AI辅助写作的Latex编辑器SmartLatexEditor，超级好用，AI润色，全免费
已经有20人回复
-
物理学I论文润色/翻译怎么收费?
已经有66人回复
-
【EI|Scopus 双检索】第六届智能机器人系统国际会议（ISoIRS 2026）
已经有1人回复
-
2026年第四届电动车与车辆工程国际会议（CEVVE 2026）
已经有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人回复
点击这里搜索更多相关资源科研从小木虫开始，人人为我，我为人人