| 查看: 514 | 回复: 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. 这是错误提示 |
回复此楼» 猜你喜欢
大豆异黄酮分离
已经有0人回复
-
湖南大学材料学院急招2026年博士生，临时增加一名博士联培指标
已经有10人回复
-
物理学I论文润色/翻译怎么收费?
已经有69人回复
-
天津理工大学晶体材料全国重点实验室刘红军教授课题组招收博士生1-2名
已经有1人回复
-
中国科学院物理研究所谌志国研究员团队招收2027年博士研究生
已经有5人回复
-
2026年中德博士后交流项目 - 新型量子和磁性材料：材料制备表征和中子散射研究
已经有12人回复
-
26申博推荐：南京航空航天大学国际前沿院光学方向招收博士生！
已经有1人回复
-
如何从铁电相到顺电相。
已经有1人回复
找到一些相关的精华帖子，希望有用哦~
帮忙：求解一组微分代数方程，奖励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人回复
点击这里搜索更多相关资源科研从小木虫开始，人人为我，我为人人
