空信谊:0-by-1

16个视图(30天)
nune pratyusha
nune pratyusha 2022年11月26日
编辑: 约翰D 'Errico 2022年11月26日
clc
清除所有
关闭所有
信谊(x1, x2) x3 x4 t
tspan = [0 100];
[solx1 solx2 solx3 solx4] =解决([x1 + x2 = = 0, x1 + x2。* x3 + (0.0001 + 0.0001 * x4。^ 4)。* x2 = = 0, 1 x2。^ 2 = = 0 x2 = = 0], [x1, x2, x3, x4))
x1val = vpa (solx1)
x2val = vpa (solx2)
x3val = vpa (solx3)
x4val = vpa (solx4)
E1 = [x1val (1:), x2val (1:), x3val (1:), x4val (: 1)]
:E2 = [x1val (2), x2val (2:), x3val (2:), x4val (2:)
E3 = [x1val (3) x2val (3:), x3val (3:), x4val (3:)
J =雅可比矩阵([x1 + x2, x1 + x2。* x3 + (0.0001 + 0.0001 * x4。^ 4)。* x2, 1 x2。^ 2 x2), (x1, x2, x3, x4))
J - 1 =潜艇(J, {x1, x2, x3, x4}, {E1 (1:)});
eigE1 =圆(vpa (eig (j - 1)), 4)
J2 =潜艇(J, {x1, x2, x3, x4}, {E2 (1:)});
eigE2 =圆(vpa (eig (J2)), 4)
J3 =潜艇(J, {x1, x2, x3, x4}, {E3 (1:)});
:eigE3 =圆(vpa (eig (J3)), 4)

答案(1)

约翰D 'Errico
约翰D 'Errico 2022年11月26日
编辑:约翰D 'Errico 2022年11月26日
不要写整个混乱的代码,实际上withput思考是否virtially第一(重要)你写产生一个可行的结果。
信谊(x1, x2) x3 x4 t
tspan = [0 100];
[solx1 solx2 solx3 solx4] =解决([x1 + x2 = = 0, x1 + x2。* x3 + (0.0001 + 0.0001 * x4。^ 4)。* x2 = = 0, 1 x2。^ 2 = = 0 x2 = = 0], [x1, x2, x3, x4))
solx1 =空信谊:0-by-1 solx2 =空信谊:0-by-1 solx3 =空信谊:0-by-1 solx4 =空信谊:0-by-1
现在仔细看看你写的方程。我将把它们摊开成4行:
x1 + x2 = = 0 (1)
x1 + x2。* x3 + (0.0001 + 0.0001 * x4。^ 4)。* x2 = = 0 (2)
1 x2。^ 2 = = 0 (3)
x2 = = 0 (4)
考虑最后一个方程。告诉你,x2 = = 0。
哦。看第三个方程。从,如果x2 = = 0,有解决方案吗?这只会发生在1 = = 0,我记得数学,是一种罕见的事件。
Seriusly,不!永远不会。从来没有。没有什么结果。Nil。不会发生。
这个系统不存在解决方案。现在,回顾解决告诉你。空信谊。我会更仔细地看看你想解决的问题。重新审视如何生成这些方程。

类别

找到更多的在数学帮助中心文件交换

标签

下载188bet金宝搏


释放

R2021b

社区寻宝

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

开始狩猎!