单击此处编辑母版文本样式,第二级,第三级,第四级,第五级,,,,点击此处结束放映,单击此处编辑母版标题样式,,人民邮电出版社,第十章 电路的计算机辅助分析,10.1 概 论,10.2 元器件及特性,10.3 电路的拓扑描述,10.4 线性电路的直流和正弦交流稳态分析,10.5 非线性电阻电路的直流分析,10.6 电路的动态分析,10.7 灵敏度与容差分析,10.8 电路的最优化设计,第十章 电路的计算机辅助分析10.1 概 论,10.1 概 论,1.,CAA,的一般步骤,用,CAA,对电路性能进行分析的一般步骤如下1) 建立电路元器件的数学模型,(2) 选择算法，编制程序,(3) 输入和输出结果,10.1 概 论1. CAA的一般步骤,2.程序编制的基本要求,(1) 易读性,(2) 效率高,(3) 通用性,(4) 使用方便，界面友好,2.程序编制的基本要求,10.2 元器件及特性,1.电阻(或电导),(1) 线性电阻,(2),流控型非线性电阻,(3),压控型非线性电阻,10.2 元器件及特性1.电阻(或电导),2.电容,(1) 线性时不变电容,(2),荷控型非线性电容,(3),压控型非线性电容,2.电容,3.电感,(1) 线性时不变电感,(2),磁控型非线性电感,(3),流控型非线性电感,3.电感,4.独立电压源和独立电流源,(1) 独立电压源,(2),独立电流源,4.独立电压源和独立电流源,5.,受控源,(1) 电压控电流源,(2),电压控电压源,(3),电流控电流源,(4),电流控电压源,5.受控源,6.互感线圈,7.,理想变压器,8.,零泛器,(1) 零器(,Nullator,),(2),泛器(,Norator,),6.互感线圈,10.3 电路的拓扑描述,1.网络图论的基本概念,(1) 有向线段,任意二端元件(或二端支路)在图论中均由一条带有方向的线段表示，其方向表示该元件(或支路)的电流或电压的参考方向，这样带方向的线段称为有向线段(可为直线或弧线，但不能为折线)。
10.3 电路的拓扑描述1.网络图论的基本概念,(2) 节点及节点集合,(3) 支路、支路集合,(4) 线图(也称拓扑图)、子图,(5) 有向图与无向图,(6) 连通图,(7) 树、树支,(8) 回路及基本回路,(9) 割集与基本割集,(2) 节点及节点集合,2.电路的拓扑矩阵,前已指出，一个网络可以抽象成由节点和支路表示成的拓扑结构图的关联性有节点与支路，支路与回路，支路与割集等它们均可用矩阵的形式来表示，这些矩阵被称为网络的拓扑矩阵下面着重介绍有向图的三个拓扑矩阵，即关联矩阵、基本回路矩阵及基本割集矩阵2.电路的拓扑矩阵,(1) 关联矩阵,关联矩阵,A,a,可以用来说明拓扑图中每个节点连接哪几条支路，以及所连接的支路相对于该节点的方向2) 基本回路矩阵,基本回路为单连支回路，因此规定连支的方向为基本回路的方向这样，基本回路矩阵,B,f,=(,b,ij,),(,b,-,n,)×,b,(1) 关联矩阵,(3) 基本割集矩阵,基本割集为单树支割集，因此规定树支的方向为方割集的方向，基本割集矩阵用,Q,f,表示，定义为,Q,f,=(,q,ij,),n,×,b,n,和,b,的含义同上3) 基本割集矩阵,3.拓扑约束关系的矩阵形式,KCL,和,KVL,为电路的拓扑约束，只与电路的有向拓扑结构有关。
在已知电路拓扑矩阵的情况下，,KCL,和,KVL,可用矩阵形式表示1),KCL,的矩阵形式,①,用关联矩阵,A,表示：,AI,b,=0,②,用基本割集矩阵,Q,f,表示：,Q,f,I,b,=0,③,用基本回路矩阵,B,f,表示：,I,b,=,B,f,T,I,l,3.拓扑约束关系的矩阵形式,(2),KVL,的矩阵形式,① 用基本回路矩阵,B,表示：,B,f,U,b,=0,②,用关联矩阵,A,表示：,U,b,=,A,T,U,n,③,用基本割集矩阵,Q,f,表示：,U,b,=,Q,f,T,U,t,(2) KVL的矩阵形式,10.4 线性电路的直流和正弦交流稳态分析,10.4.1 线性电路方程组的建立,电路方程组的建立方法一般有拓扑矩阵法和元件贡献直接添加法10.4 线性电路的直流和正弦交流稳态分析10.4.1 线,1.建立线性电路方程组的拓 扑矩阵法,拓扑矩阵法的基本思路：由电路的两类约束出发，对不同的分析法，消去相应的非自变量1) 支路伏安关系,在电路的计算机分析中，为了减少电路中支路和节点的数目，同时将有源器件和无源元件合在一个支路中处理，引入了“组合支路”的概念，如图10-8所示1.建立线性电路方程组的拓 扑矩阵法,(2) 电路的拓扑方程,① 支路电流法,②,支路电压法,③ 基本回路方程,④,基本割集方程,⑤ 节点电位方程,(2) 电路的拓扑方程,⑥ 改进节点法,在改进节点法中，将网络元件分为可用导纳描述和不可用导纳描述两类，前者以节点电压作自变量，后者以支路电流为自变量。
⑥ 改进节点法,2.用元件贡献直接添加法建 立电路方程,元件贡献直接添加法的关键在于找出各种元件对节点方程的贡献为便于理解，先分别讨论电路中出现某一类元件的情况1) 导纳型支路,(2) 理想电压源,2.用元件贡献直接添加法建 立电路方程,(3) 理想变压器,(4) 互感线圈,(5),VCCS,(6),VCVS,(7),CCCS,(8),CCVS,(3) 理想变压器,10.4.2 线性方程组的求解,1.高斯(,Gauss),消元法,高斯消元法是一种古老的方法我们在中学学过消元法，高斯消元法就是它的标准化的适合在计算机上自动计算的一种方法1) 高斯消元法的基本思想,将方程,AX,＝,b,通过消元化为等价的三角方程组，然后回代解之10.4.2 线性方程组的求解,(2) 高斯消元法公式,记,AX,＝,b,为,A,,(1),X,＝,b,,(1),，,A,(1),,和,b,(1),,的元素记为,a,(1),,ij,和,b,,(1),,i,，,i,，,j,＝1，2…,n,2) 高斯消元法公式,2.列主元高斯消元法,(1) 基本思想,对第,k,次消元，从第,k,行到第,n,行选出第,k,列中绝对值为最大的元素,a,k,q,,k,，,作为第,k,次消元主元，交换第,k,行和第,q,行，然后按高斯消元法进行消元。
回代与高斯消元法一样2) 数学公式,2.列主元高斯消元法,3.,LU,分解法,对高斯消元法，当,AX,＝,b,中,A,不变而,b,向量不断改变时，则对应每一组向量,b,就要作一次高斯消元运算，这显然很不经济1) 基本思路,对矩阵方程,AX,＝,b,，,将,A,化为两个三角矩阵的乘积其中一个为单位下三角阵，另一个为上三角阵，即,A,＝,LU,3.LU分解法,10.5 非线性电阻电路的直流分析,10.5.1非线性电路,非线性电路分析在电路分析中是一个很重要的组成部分，因为实际电路中广泛地应用着各种非线性器件10.5 非线性电阻电路的直流分析10.5.1非线性电路,10.5.2非线性电阻电路方程 组的建立,1.节点电位方程组,若电路中仅含有线性电阻、独立电源和电压控制型非线性电阻，受线性或非线性电阻元件上电压控制的电流源(,VCCS,)，,受线性或非线性元件中电流控制的电流源(,CCCS,)，,且不含有零电阻支路，则这类电路可用节点法来分析10.5.2非线性电阻电路方程 组的建立,2.混合型方程组(改进节点法),方程组变量由节点电位向量、独立电压源支路电流向量及非线性电阻上的电流向量组成。
2.混合型方程组(改进节点法),(1) 先将电路中支路进行分类,仅讨论含以下4种元件的情况：,① 线性支路(复合支路)；,② 纯独立电压源支路；,③ 电流控制型非线性电阻(复合支路)；,④ 受上述支路电流控制的线性,CCCS,支路1) 先将电路中支路进行分类,10.5.3 非线性代数方程组的求解,1.单变量非线性方程根的数值解,(1) 一般迭代,(2),Newton,迭代,① 一般迭代法的加速,②,Newton,迭代法,,10.5.3 非线性代数方程组的求解,2.多变量非线性方程组的,Newton,迭代,当含,n,个变量时，非线性方程组的一般形式为,F,i,=(,x,1,，,x,2,，…，,x,n,)=0,i,=1，2，…，,n,2.多变量非线性方程组的 Newton迭代,10.6 电路的动态分析,10.6.1 概述,含电感、电容的电路称为动态电路由于电感、电容的,VAR,出现微分，电路方程不再像电阻电路那样是代数方程，而是微分方程一方面，如果微分方程的阶次很高，激励又是一般的激励的情况下，解微分方程是相当因难的，有的甚至没有解析解10.6 电路的动态分析10.6.1 概述,10.6.2 状态方程的建立,状态方程的建立方法可分为两大类型：直接法与间接法。
下面我们重点介绍拓扑矩阵法和常规的直观偏写法1.拓扑矩阵法,拓扑矩阵法通过拓扑矩阵利用计算机自动编写电路方程10.6.2 状态方程的建立,(1) 选规范树，支路分类，编号,(2),KCL，KVL,约束,(3),VCR,约束,(4) 状态方程,,(1) 选规范树，支路分类，编号,2.直观编写法,对不太复杂的电路，可用直观编写法建立电路的状态方程一般步骤：,① 选取独立电感电流和独立电容电压为状态变量；,② 对含有独立电感的独立回路列写,KVL,方程；,2.直观编写法,③ 对接有独立电容的独立节点列写KCL方程；,④ 将KVL和KCL方程中不是状态变量的量列写出由状态变量表示的约束关系式；,⑤ 消去非状态变量，把方程写为状态方程矩阵形式的标准表示式③ 对接有独立电容的独立节点列写KCL方程；,10.6.3 常微分方程初值 问题的数值解,常微分方程数值解法的基本原理：对求解区间进行剖分，然后把微分方程离散成在节点上的近似公式或近似方程，最后结合定解条件求出近似解下面介绍几种常用的数值解法10.6.3 常微分方程初值 问题的数值解,1.欧拉方法,(1) 前向欧拉法,(2) 欧拉隐式公式和欧拉中点公式,(3) 梯形公式和预估校正法,1.欧拉方法,2.龙格—库塔法(,R,-,K,法),欧拉法使用方便，但精度较低。
下面介绍一种比较适用的算法——龙格——库塔法由于算法导出较为复杂，只给出算法公式1) 二阶,R,-,K,法,(2) 四阶,R,—,K,法,2.龙格—库塔法(R-K法),10.7 灵敏度与容差分析,10.7.1 灵敏度分析,1.灵敏度的定义,10.7 灵敏度与容差分析10.7.1 灵敏度分析,2.灵敏度的计算,计算灵敏度的方法很多，有导数法、差商法、符号网络函数法、增量网络法、伴随网络法以及转置法等下面将着重介绍增量网络法和伴随网络法，而且只涉及线性系统的灵敏度分析2.灵敏度的计算,(1) 增量网络法,① 增量网络的形成,对线性网络,N,，,在网络的激励和拓扑结构一定的条件下，支路阻抗(或导纳)的微小变化必然会使支路电压和电流也有相应的微小变化a,),阻抗支路,(,b,),导纳支路,(,c,),独立源,(,d,),受控源,(1) 增量网络法,(2) 伴随网络法,伴随网络法的基础是特勒根(,Te llegen),定理该定理叙述了在具有相同拓扑结构的两个网络中，对应支路电流与电压之间的特定关系① 特勒根定理,② 伴随网络的构成,(2) 伴随网络法,③ 用伴随网络法计算灵敏度,(,a,),形成端口网络,(,b,),灵敏度的计算,,③ 用伴随网络法计算灵敏度,10.7.2 容差分析,容差分析通常包括最坏情况分析和统计分析。
1.最坏情况分析,2.统计分析,容差的统计分析方法有两种：蒙特卡罗(,Monte-Carlo),分析法和矩量法(,Moment Method)10.7.2 容差分析,10.8 电路的最优化设计,电路的最优化设计利用数学中最优化理论与方法，以计算机为手段，对电路进行辅助设计目前电路的最优化设计过程一般是先根据所需的电路性能指标要求，设计者给出初始方案(包括电路的结构和元件参数)；然后由计算机进行分析，求出初始电路的各项响应，并与指标要求相比较；若不能满足要求，则自动修改电路参数和结构， 经过反复计算、修改，直至满足性能要求为止10.8 电路的最优化设计电路的最优化设计利用数学中最优化,1.目标函数,电路的最优化设计，实际上是在一定约束条件下求函数极值的问题在有约束条件的情况下，它的数学描述可以是,,minF,(,P,),q,i,(,P,)≥0,i,=1，2…，,m,,h,j,P,=0,j,=1，2…，,n,1.目标函数,其中,P,(,p,1,，,p,2,，…，,p,n,),为设计变量，在电路设计中一般为元件参数向量q,i,(,P,),和,h,j,(,P,),为设计变量的约束方程。
F,(,P,),称为目标函数2.最优化基本原理,当目标函数建立之后，最优化设计问题就转为求变化元件参数向量,P,，,使,F,(,P,),为最小的问题，即求函数的极值问题其中P(p1，p2，…，pn)为设计变量，在电路设计中一般,3.最速下降法,在电路设计中，目标函数中的可调电路参数往往不止一个，属于多变量函数使用计算机分析时，一般采用迭代法，即用迭代法一步一步地把目标函数引向最小的方法3.最速下降法,4.最佳步长确定,由于在每步迭代时需要确定△,P,，,以确保,F,(,P,),逐渐减小，最速下降法指出,F,(,P,),的负梯度方向作为△,P,的方向于是在每步迭代中，目标函数只是步长,λ,的函数，对第,j,＋1,步迭代有,F,(,P,(,j,＋1),)=,F,(,P,(,j,),＋△,P,(,j,),)=,F,(,P,(,j,),＋,λ,(,j,),S,(,j,),),4.最佳步长确定,(1) 用外推法确定,λ,*,所在的区间,(2) 用黄金分割法确定,λ,*,(1) 用外推法确定λ*所在的区间,。