问题描述
我正在尝试遵循以下建议: 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
线性求解的信息。对于您的问题,很可能是因为问题在数值上难以解决,或者自由度的变量交换没有关联的方程或目标来确定该变量。