警告:在t = 9.057068 e-01失败。无法满足集成公差没有减少步长低于最小值允许在时间t (1.776357 e15汽油)。

2视图(30天)

显示旧的评论

贾丝廷娜Ranceva 2021年10月3日

0
链接

这个问题直接联系

//www.tatmou.com/matlabcentral/answers/1465784 -警告-失败- - t - 9 - 057068 e - 01 -没- -满足集成公差-不-减少步骤

评论道: 沃尔特·罗伯森 2021年10月3日

                          tspan = (10 0);
                         
                          y0 = [0;0;0;0;0;0;0;0;0;0);
                         
                          选择= odeset (“质量”@ (t, q)质量(t, q, P));
                         
                          (t, q) =数值(@ (t, q) f (t, q, P), tspan, y0,选择);

我得到了这个消息:

警告:在t = 9.057068 e-01失败。无法满足集成公差没有减少步长低于允许的最小值

在时间t (1.776357 e15汽油)。

我必须做什么?

0评论
显示隐藏1以上的评论

登录置评。

在回答这个问题。

答案(1)

沃尔特·罗伯森 2021年10月3日

0
链接

直接链接到这个答案

//www.tatmou.com/matlabcentral/answers/1465784 -警告-失败- - t - 9 - 057068 e - 01 -没- -满足集成公差-不-减少answer_800629步

你的f函数在时间t = 9.057068 e-01奇点

你没有把你的代码质量或者你的代码 f 你没有把你的 P 变量,所以我们不能测试或检验方程更具体。

4评论
显示隐藏 3年长的评论

贾丝廷娜Ranceva 2021年10月3日

f函数:

                                   函数dydt = f (t, q, P)
                                  
                                   m1 = P.m1;
                                  
                                   m2 = P.m2;
                                  
                                   m3 = P.m3;
                                  
                                   m4 = P.m4;
                                  
                                   I1 = P.I1;
                                  
                                   k1 = P.k1;
                                  
                                   k2 = P.k2;
                                  
                                   k3 = P.k3;
                                  
                                   k4 = P.k4;
                                  
                                   c1 = P.c1;
                                  
                                   c2 = P.c2;
                                  
                                   c3 = P.c3;
                                  
                                   c4 = P.c4;
                                  
                                   b1 = P.b1;
                                  
                                   a1 = P.a1;
                                  
                                   a2 = P.a2;
                                  
                                   L = P.L;
                                  
                                   fi0 = P.fi0;
                                  
                                   g =供货商;
                                  
                                   L1 = P.L1;
                                  
                                   z1 = P.z1;
                                  
                                   z2 = P.z2;
                                  
                                   dydt = [q (2)
                                  
                                   m1 * g-k1 * (q (1) z1) + k2 * (q (5) - q (1) a1 * (q(7)))往上平移* q (2) + c2 * (q (6) - q (2) a1 * (q (8)))
                                  
                                   问(4)
                                  
                                   m2 * g-k3 * (q (3) z2) + k4 * (q (5) - q (3) a2 * (q (7))) c3 * q (4) + c4 * (q (6) - q (4) + a2 * (q (8)))
                                  
                                   问(6)
                                  
                                   m3 * g-k2 * (q (5) - q (1) a1 * (q (7))) k4 * (q (5) - q (3) a2 * (q (7))) c2 * (q (6) - q (2) a1 * (q (8))) c4 * (q (6) - q (4) + a2 * (q (8)))
                                  
                                   问(8)
                                  
                                   a1 * k2 * (q (5) - q (1) a1 * (q (7))) a2 * k4 * (q (5) - q (3) a2 * (q (7))) + a1 * c2 * (q (6) - q (2) a1 * (q (8))) a2 * c4 * (q (6) - q (4) + a2 * (q (8)))
                                  
                                   问(10)
                                  
                                   m4 * g * L * cos (q (9)));
                                  
                                   结束

贾丝廷娜Ranceva 2021年10月3日

质量功能:

                                   函数M =质量(t, q, P)%质量矩阵函数
                                  
                                   %提取参数
                                  
                                   m1 = P.m1;
                                  
                                   m2 = P.m2;
                                  
                                   m3 = P.m3;
                                  
                                   m4 = P.m4;
                                  
                                   I1 = P.I1;
                                  
                                   k1 = P.k1;
                                  
                                   k2 = P.k2;
                                  
                                   k3 = P.k3;
                                  
                                   k4 = P.k4;
                                  
                                   c1 = P.c1;
                                  
                                   c2 = P.c2;
                                  
                                   c3 = P.c3;
                                  
                                   c4 = P.c4;
                                  
                                   b1 = P.b1;
                                  
                                   a1 = P.a1;
                                  
                                   a2 = P.a2;
                                  
                                   L = P.L;
                                  
                                   fi0 = P.fi0;
                                  
                                   g =供货商;
                                  
                                   L1 = P.L1;
                                  
                                   z1 = P.z1;
                                  
                                   z2 = P.z2;
                                  
                                   %设置质量矩阵的非零元素
                                  
                                   M = 0 (10,10);
                                  
                                   M (1,1) = 1;
                                  
                                   米(2,2)= 1;
                                  
                                   米(3、3)= 1;
                                  
                                   米(4,4)=平方米;
                                  
                                   米(5,5)= 1;
                                  
                                   米(6,6)= m3 + m4;
                                  
                                   米(6、8)= m4 * L * cos (fi0 + q (7));
                                  
                                   10米(6)= m4 * L * sin (q (9));
                                  
                                   (7)= 1;
                                  
                                   6米(8日)= m4 * L * cos (fi0 + q (7));
                                  
                                   8米(8日)= I1 + L1 * (cos (fi0 + q (7))) ^ 2 + L1 * (sin (fi0 + q (7))) ^ 2;
                                  
                                   10米(8日)= m4 * (L1 * cos (fi0 + q (7)) * L *罪(q (9)) + L1 * sin (fi0 + q (7)) * L * cos (q (9)));
                                  
                                   9米(9日)= 1;
                                  
                                   6米(10日)= m4 * L * sin (q (9));
                                  
                                   8米(10日)= m4 * (L * sin (q (9)) * L1 * cos (fi0 + q (7)) + L * cos (q (9)) * L1 * sin (fi0 + q (7)));
                                  
                                   米(10,10)= m4 * ((L * sin (q (9))) ^ 2 + (L * cos (q (9))) ^ 2);
                                  
                                   结束

沃尔特·罗伯森 2021年10月3日

目前,我没有看到任何明显的问题,只是阅读代码。我以后再试试当我回到我的桌子上。

登录置评。

在回答这个问题。

类别

工程机械工程计算流体动力学(CFD)

找到更多的在计算流体动力学(CFD)在帮助中心和文件交换

社区寻宝

找到宝藏在MATLAB中央,发现社区如何帮助你!

开始狩猎!

警告:在t = 9.057068 e-01失败。无法满足集成公差没有减少步长低于最小值允许在时间t (1.776357 e15汽油)。

0评论
显示隐藏1以上的评论

答案(1)

4评论
显示隐藏 3年长的评论

另请参阅

类别

标签

社区寻宝

警告:在t = 9.057068 e-01失败。无法满足集成公差没有减少步长低于最小值允许在时间t (1.776357 e15汽油)。

0评论 显示隐藏1以上的评论

答案(1)

4评论 显示隐藏 3年长的评论

另请参阅

类别

标签

社区寻宝

0评论
显示隐藏1以上的评论

4评论
显示隐藏 3年长的评论