gekko 错误：矩阵在结构上是奇异的怎么解决呢？

如何解决gekko 错误：矩阵在结构上是奇异的怎么解决呢？

我正在尝试遵循以下建议： How to apply m.connection from Gekko using arrays?

现在我收到这些错误：从 MC23A 返回错误，因为矩阵在结构上是奇异的，秩 = 9453 从 MA28A/AD 返回错误，因为从 mc23a/ad 返回错误错误：异常：访问冲突回溯：不可用，使用 -ftrace=frame 或 -ftrace=full 编译错误：找不到“results.json”

你知道我做错了什么吗？提前致谢。

更新后的代码：

<div id="container">
        <canvas id="board" width="1000px" height="1000px"></canvas>
    </div>

解决方法

为了运行脚本，我取消了这部分代码的注释：

hocl=m.Array(m.Var,2)
nh3=m.Array(m.Var,2)
nh2cl=m.Array(m.Var,2,lb=1e-10)
nh2cl_meas=m.Array(m.Param,2)
nhcl2=m.Array(m.Var,lb=1e-10)
h=m.Array(m.Var,2)
oh=m.Array(m.Var,2)
I=m.Array(m.Var,lb=1e-10)
ocl=m.Array(m.Var,2)
nh4=m.Array(m.Var,2)
h2co3=m.Array(m.Var,2)
hco3=m.Array(m.Var,2)
co32=m.Array(m.Var,2)
alk=m.Array(m.Var,2)
DOC1=m.Array(m.Var,2)
DOC2=m.Array(m.Var,2)
cnh3=m.Array(m.Var,2)
cnh2cl=m.Array(m.Var,2)

代码运行但给出了可疑的结果，浓度或比率为负。

错误Structurally Singular表示在IPOPT中使用线性求解器求解Ax=b方程时在寻找搜索方向时存在问题。奇异矩阵是不可逆的。如果矩阵的行列式为零或者对于数值解，如果 A 矩阵具有阻止 Ax=b 解的特定属性，则矩阵是奇异矩阵。有更多关于如何在 Interior Point algorithm 中使用 Ax=b 线性求解的信息。对于您的问题，很可能是因为问题在数值上难以解决，或者自由度的变量交换没有关联的方程或目标来确定该变量。