1. 课程导入与问题界定

在工程计算与科学研究的绝大多数场景中，我们最终要回答的问题往往可以归结为一句话：当某个量取什么值时，系统恰好满足给定的平衡条件或设计指标？把这句话写成数学形式，通常就得到一个方程 f(x) = 0。若 f 关于 x 是线性的，问题可以用消元法或矩阵分解一次求解；但真实世界中的 f 几乎总是非线性的——它可能包含三角函数、指数、对数、幂次、分段判断，甚至本身就是一个只能通过数值仿真才能给出取值的"黑箱"。这时，解析求根公式往往不存在，我们只能借助迭代方法逐步逼近真解。本课程讨论的正是这套方法体系：如何构造迭代、如何判断收敛、如何控制误差、如何在多解中锁定工程上真正需要的那个根。

为什么非线性求根是"绕不开"的基本功

先看几个不同领域的典型问题，它们表面上差异极大，数学结构却高度一致。

• 机构设计：四杆机构中，给定曲柄转角 θ，求连杆与摇杆的夹角 φ，满足闭环矢量方程 acosθ + bcosφ − c − dcosψ = 0 一类的关系，含三角函数，必须迭代求解。
• 化工与热力学：用范德华方程求给定温度压力下的摩尔体积 v，方程是 v 的三次多项式，且存在气液两相多个实根，需要按物理意义挑选。
• 电路分析：含二极管的非线性电路，支路电流满足 I = I_s(e^V/nV_T − 1)，与线性电阻网络联立后形成超越方程组。
• 控制与优化：求闭环系统的临界增益，等价于求特征方程在虚轴上的根，或求某个增益函数等于 1 且相角为 −180° 的交点。
• 反问题与标定：已知实测响应，反推材料参数或模型系数，本质是求残差函数零点或极小值点。

这些问题的共同点是：方程本身不难写，难在"解不出来"。工程上真正消耗时间的，往往不是建模，而是让模型"算出数来"并确认这个数可信。因此，非线性方程求根不是一门孤立的数学技巧，而是连接模型与结论的关键一环。

问题的数学界定：我们到底在求什么

设 f: ℝ → ℝ 连续。求根即寻找 x* 使 f(x*) = 0。看似简单，实际包含四个必须分别回答的子问题，它们决定了方法选择与结果可信度：

1. 存在性：区间 [a, b] 上是否有根？连续函数若 f(a)·f(b) < 0，由介值定理至少存在一个根，这是二分法的理论基石。
2. 唯一性与多解：单调区间内根唯一；非单调时可能出现多个根，甚至重根（f 与 f′ 同时为零），后者会使收敛速度显著退化。
3. 收敛性：迭代序列 x_k 是否收敛？收敛速度是线性、超线性还是二次？这直接决定迭代次数。
4. 误差与停止准则：残差 |f(x_k)| 小，是否等于 x_k 接近真解？在 f 的斜率很小时，残差小但解可能偏差很大，因此需要同时监控步长 |x_k+1 − x_k| 与残差。

本节核心关系的可视化

下图把"实际问题 → 数学方程 → 方法选择 → 收敛判据 → 工程可用解"这条主线，与四类典型非线性形态并列呈现，帮助建立全局认识：不同函数形态（单根、多根、重根、黑箱）对应不同的方法策略与风险点。

一个可动手的示例：从问题到第一次迭代

为了让"问题界定"落到可操作层面，考虑一个教学示例（参数为合理假设，非实测）：某散热器设计要求在给定工况下满足 f(x) = x³ − 2x − 5 = 0，其中 x 为归一化设计变量。先做界定工作：

1. 存在性：f(2) = 8 − 4 − 5 = −1 < 0，f(3) = 27 − 6 − 5 = 16 > 0，故 [2, 3] 内至少一根。
2. 唯一性：f′(x) = 3x² − 2，在 [2, 3] 上恒大于 0，函数严格单调，根唯一。
3. 方法选择：f 光滑且导数易得，取 x₀ = 2，用牛顿法 x_k+1 = x_k − f(x_k)/f′(x_k)。
4. 首步计算：f(2) = −1，f′(2) = 10，得 x₁ = 2 − (−1)/10 = 2.1；继续迭代会快速逼近约 2.0946。
5. 停止准则：同时要求 |x_k+1 − x_k| < 10⁻⁸ 且 |f(x_k)| < 10⁻¹⁰，避免"残差小但解偏"的假收敛。

这个例子虽小，却完整体现了本课程的工作方式：先界定（存在、唯一、形态），再选法，再迭代，最后用双判据验收。后续章节将逐一展开每一步的理论依据、收敛证明、失效情形与工程实现细节，并扩展到非线性方程组、多解筛选与鲁棒算法设计。

学习本课程后，你应当能够：面对一个工程问题，独立完成方程建立与根的定性分析；根据函数可导性、初值信息与精度要求选择合适方法；写出可运行的迭代程序并设置合理的停止准则；对不收敛或多解情形做出诊断与修正；最终给出带有误差说明、可供工程决策使用的数值解。