清华大学电路原理19非线性电路简介.ppt

19.6 非线性电阻电路方程的数值求解方法—— 牛顿-拉夫逊法（Newton-Raphson Algorithm） 一、具有一个未知量的非线性代数方程求解 0 x f(x) 设方程 f(x) = 0 解为x*， 则f(x *) = 0 。x*为 f(x) 与 x 轴交点。 基本思路：线性化。 具体做法：有指导的猜试、迭代。 设 x0 f(x0) = 0 x0 就是解 设 x1 f(x1) = 0 x1 就是解 ? ? ? f(xk) ? 停止 实际处理: ? xk － xk-1 ? ? xk 就是解。 利用牛顿-拉夫逊法求x* 步骤如下。 （1） 选取一个合理值x0，称为迭代初值。此时x0 一 般与 x* 不等。 （2） 迭代 。 取x1 =x0+?x0 作为第一次修正值，?x0 充分小。 将 f ( x0+?x0 ) 在? x0 附近展开成台劳级数： （3）若 ? xk+1 － xk ? ? 则 xk+1 就是方程的解 x* 。 ? ? ? 迭代公式 几何解释： 将 f(x) 在 x0 处线性化。 取线性部分，并令 ? xk+1 - xk ? ? k=k+1 k=0 x0 No Yes x* =xk+1 程序流程 0 x f(x) 叠代过程的几何解释： 收敛性：与函数本身有关，与初值有关。 切线 解：列节点方程 例 + ? IS1 Un I3 U3 R2 取 ，迭代结果如下表。 k 0 1 2 3 4 0 -2 0.85714 0.73469 0.67563 0.03295 0.66669 0.00009 0.66667 0.00001 经四次迭代后，得 注意：初估值选择不好会产生振荡 （迭代不收敛）。 二、具有多个未知量的非线性方程组的求解 设 n 个未知量 一般表示为 设第 k 次迭代时 若 ，则 即为所求的一组解答。 下面分析每次修正值?xj ( j=1，2，?，n)的计算 对 先选一组初估值 进 行第一次计算，然后不断修正，进行迭代运算。 若 ，则进行修正，寻找 将 在 xik 附近展成台劳级数，取线性部分，并令其 等于零， 得 简记为： J 称为雅可比矩阵 得方程组的解 X k +1 写成矩阵形式为： 返回目录 19.7 用友网络模型求解非线性电阻电路 基本思想： 找出由k次迭代值求第k+1次迭代值的线性化模型。 i = f ( u ) 令 uk， uk+1 分别为第 k 次和第 k+1 次的电压初值，其 对应的电流分别为 + ? u i 与牛顿法的区别： 牛顿法：对非线性方程不断线性化。 友网络：对非线性电阻伏安特性不断线性化。 把 在 处展开成台劳级数， 取线性部分，即将非线性电阻在 处线性化， 式中 为非线性电阻在 点处的动态电导。 得 得 将电路中所有非线性电阻分别用各自的线性化模型 代替，就可得到和原电路对应的“友网络模型”。逐次 迭代计算，即可得到所要求的解。 非线性电阻在第 k+1 次迭代时的“线性 化”模型。 + - 在进行第 k+1 次迭代时， 是已知的。 上述关系可用如下等效电路来描述： 解 画出友网络模型。 + - unk +1 例 + - IS1 Un I3 U3 R2 列节点方程（ k+1次迭代时）： 将上述关系代入 前 式，得 已知 u2 i2 2 几何解释 u i 0 u0 i0 0 u1 i1 1 P 切线 小结