| 查看: 496 | 回复: 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人回复
-
散金币，求好运，祝面上顺利！
已经有80人回复
-
物理学I论文润色/翻译怎么收费?
已经有193人回复
-
华东师范大学芯片设计徐珑真实水平咋样
已经有1人回复
-
调剂
已经有0人回复
-
【新加坡】纳米电子器件项目组有“联合培养博士生”名额
已经有0人回复
-
中德博士后2026项目5月30日截止- 国内博士后入站、德国两年全职
已经有0人回复
-
面向全球申请人的中德博士后2026项目5月30日截止- 国内博后入站德国两年全职
已经有0人回复
-
【新加坡】纳米电子器件项目组有“联合培养博士生”名额
已经有2人回复
-
重庆交大材料学院26年硕士招生复试实施细则已出。欢迎加入光子学微结构与器件课题组。
已经有0人回复
-
北工大宋晓艳教授团队招收2026年9月入学博士【学博和专博名额均有】！
已经有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人回复
点击这里搜索更多相关资源科研从小木虫开始，人人为我，我为人人
