在MATLAB中求解约束线性和非线性方程组
我正在使用直接刚度法在MATLAB中解决FEM问题。现在将问题表述为非线性方程组：
KU=F or CF=U with C = K^{-1}.
问题是在3-Dd中制定的，并且有200多个节点。该模型用作变形估计器，力总是输入和变形输出。所以有很多线性方程。现在需要扩展这个问题，需要实现2种类型的约束。一些节点必须受到约束，以便它们只能围绕给定点以圆周运动移动。并且一些节点必须以最大变形量绑定。
总共有数百个形式方程式的方程式：
a_1 X(1) + a_2 X(2) + ... - F(1) = 0
非线性形式的约10个等式约束（也可能是正弦曲线）：
X(1)^2 + X(2)^2 = L
和形式的10个不等式约束：
X(1) < 30\pi/180
同时解决具有此类约束的非线性方程组的最佳方法是什么？是否可以使用fsolve？还是fmincon？我一直在研究求解器，但是我无法找到一个可靠的结论，即解算器能够同时处理所有这些方程类型，以及在计算时间方面最有效的方法。