它可能会不利地改变条件数，是的，这将使其更难解决。你如何强迫这个矩阵是对称的？应该说，如果你解决 10 阶线性问题，你可能做得很好。

正如 whpowell 在他们的评论中指出的那样，matlab 将为对称与非对称使用不同的求解程序。它将在对称矩阵的两个求解器之一和非对称矩阵的不同两个求解器之间进行选择。因此，您不仅要更改矩阵条件数，还要更改求解器。

不是真正的答案，但评论太长了。

根据我的经验，保持对称性和处理边界条件的最佳方法是仅使用内部点。“惩罚方法”（或大数字方法，我知道它有几个名称）是另一种可能性。

如果您的枚举是$x_1,\ldots,x_{n+1}$，然后继续工作$x_2, \ldots, x_{n}$，作为的值$u$在$x_1,x_{n+1}$给出（在你的情况下，只是$0$）。ig他们不是$0$，它们将在 rhs-vector 上移动

您可以通过以下方式减少刚度矩阵的自由度$1$，并保持对称性。在 lin 之后，您将获得一个解向量。系统求解过程，仅由$n-1$内部价值观。然后，您可以在 sol.txt 的末尾添加已知值。如果你愿意，矢量。