你的系统是非线性的。毫无疑问,“解决”做不到这一点,因为它高度非线性。如果你能计算公式解决方案为什么使用解决?解决的另一个limitaion甚至是它不能解决高阶方程与象征性的线性系统,它停止。如果你还是想解决数值,使用“vpasolve”。
Matlab解决找不到明确solutionsfor象征性的方程。万博 尤文图斯
2视图(30天)
显示旧的评论
我有一个12方程组,我想找到定点当所有方程= 0。这是我的代码:
清晰的(x1, x2) x3 x4 x5 x6 r1, r2 r3 r4 r5 r6 f1 f2 f3 f4 f5 f6 c pαβb
信谊(x1, x2) x3 x4 x5 x6 r1, r2 r3 r4 r5 r6 f1 f2 f3 f4 f5 f6 c pαβb%所有符号
%定义常微分方程
x1dot x1 - p = f1 * * x1 * (r1 + r2β* +β* r4);
x2dot x2 = f2 * x2 - p * * (r2);
x3dot = f3 * x3 - p * x3 * (r3 +β* r2);
x4dot = f4 * x4 - p * x4 *(β* r1 + r4 +β* r5);
x5dot = f5 * x5 - p * * * (r5);
x6dot = f6 * x6 - p * x6 *(β* r5 + r6);
r1dot = c * ((x1 * r1 /αα(r1 + r2 * + * r4)) +(α* x4 * r1 / (r4 +α* r1 +α* r5))) - b * r1;
r2dot = c *((αr2 * x1 * /αα(r1 + r2 * + * r4)) + x2 + r2(α* x3 * / (r3 +αr2 *))) - b * r2;
r3 r3dot = c * ((x3 * /(α* r2 + r3))) - b * r3;
r4dot = c * (x1(α* * r4 /αα(r1 + r2 * + * r4)) + (x4 * r4 / (r4 +α* r1 +α* r5))) - b * r4;
r5dot = c *((α* x4 * r5 / (r4 +α* r1 +α* r5)) + x5 +(α* x6 * r5 /(α* r5 + r6))) - b * r5;
r6dot = c * (x6 * r6 /(α* r5 + r6)) - b * r6;
%定义系统常微分方程的函数
有趣= [x1dot x2dot、x3dot x4dot, x5dot, x6dot, r1dot, r2dot, r3dot, r4dot, r5dot, r6dot);
%计算固定的点,解决函数= = 0
S =解决(x1dot = = 0, x2dot = = 0, x3dot = = 0, x4dot = = 0, x5dot = = 0, x6dot = = 0, r1dot = = 0, r2dot = = 0, r3dot = = 0, r4dot = = 0, r5dot = = 0, r6dot = = 0, (x1, x2) x3 x4 x5 x6 r1, r2 r3 r4 r5 r6));
我能找到的解决方案。万博 尤文图斯我也使用了相同的代码找到固定的类似的方程组,所以我不确定什么是错的。谁能帮我算出什么问题呢?
有趣的是,这段代码适用于下列方程组:
x1dot x1 - p = f1 * * x1 * (r1 + r2β* +β* r5);
x2dot x2 = f2 * x2 - p * * (r2);
x3dot = f3 * x3 - p * x3 * (r3 +β* r2);
x4dot = f4 * x4 - p * x4 *(β* r2 + r4 +β* r5);
x5dot = f5 * x5 - p * * * (r5);
x6dot = f6 * x6 - p * x6 *(β* r5 + r6);
r1dot = c * (x1 * r1 /αα(r1 + r2 * + * r5)) - b * r1;
r2dot = c *((αr2 * x1 * /αα(r1 + r2 * + * r5)) + x2 + r2(α* x3 * / (r3 +αr2 *)) +(αr2 * x4 * / * r2(α+ r4 +α* r5))) - b * r2;
r3 r3dot = c * ((x3 * /(α* r2 + r3))) - b * r3;
r4dot = c * (x4 * r4 / (r4 +α* r2 +α* r5)) - b * r4;
r5dot = c * (x1(α* * r5 /αα(r1 + r2 * + * r5)) +(α* x4 * r5 /(αr2 * + r4 +α* r5)) + x5 +(α* x6 * r5 /(α* r5 + r6))) - b * r5;
r6dot = c * (x6 * r6 /(α* r5 + r6)) - b * r6;
这是非常相似,有相同数量的变量。事实上这段代码给我100多个解决方案和conidtions与他们有关。万博 尤文图斯
谢谢你！
