基于MATLAB-GUI的线性动态电路分析

河南科技学院机电学院本科学位论文I摘 要在电路分析中,对线性动态电路的分析方法一般是经典法和拉普拉斯变换法。不过对于较高阶的动态电路,用这些方法就会比较繁琐。而本文设计了一种新的方法-伴随网 络法。 这种方法是把动态电路的过渡过程时间t 0T 划分成若干时间间隔t，把动态元件电感 L和电容C用相应的离散模型来取代，经过代换后的电路称为原电路的伴随网络。因而对该时间间隔内相应伴随网络的分析可视为稳态电路分析。不过这种方法需要多次重复计算,若由人来计算就会相当的复杂。所以我们利用MATLAB软件来进行编程，把计算工作交 给计算机。这样就可以方便快捷地算出结果。为了使软 件的界面更加友好，我 们 利用MATLAB自带的GUI工具做了一个图形用户界面。 用户只需要按照界面的要求进行元件参数和支路编号的输入，运行后即可在屏幕上显示GUI界面下的输出电压波形和电流波形。关键词：伴随网络法，图形用户界面，离散模型IIAbstractGenerally, the method of the linear dynamic circuit analysis is the Classic and Laplace transform method in circuit analysis. But these methods will be complex to the high-level linear dynamic circuit. In this article we design a new method-accompanying network method. This method divided the time (t0T) of transitional process of dynamic circuit into a number of intervals t, The dynamic component inductance L and capacitance C with the corresponding discrete model to replace, after substitution of the circuit , the original circuit is called accompanying network. So the analysis of accompanying network could be deemed as steady-state analysis. However, the times of calculation are so many. It will be too complex to calculate it by people. Therefore, we use MATLAB to program. By doing this, it will be convenient to obtain the results. At the same time, in order to get a more friendly interface, I also use the tools (GUI) of MATLAB made a graphical user interface. The user only needs to input the parameter of the components and the number of branches according to the requirement of the interface. After the running, it could output the voltage and current waveforms on the GUI interface.Key words：accompanying network , graphic user interface, discrete model河南科技学院机电学院本科学位论文III目 录1 绪论 .12 线性动态电路分析的常用方法 .22.1 经典法 .22.2 拉普拉斯变换法 .23 电路分析中常用的仿真软件 .23.1 PSPICE 和 SABER.23.1.1 PSPICE 简介 .23.1.2 SABER 简介 .33.2 MATLAB.43.2.1 MATLAB 简介 .43.2.2 MATLAB-GUI 简介 .54 伴随网络法 .74.1 动态元件的离散模型 .74.1.1 电容 C 的伴随网络模型 .74.1.2 电感 L 的伴随网络模型 .84.2 伴随网络法分析瞬态电路的步骤 .85 伴随网络法的 MATLAB 编程实现 .85.1 程序流程图 .95.1.1 主程序流程图 .95.1.2 GUI 程序流程图 .105.2 算例分析 .105.2.1 输入参数 .115.2.2 GUI 输出参数 .126 结语 .13致谢 .14参考文献 .15附录 .16IV河南科技学院机电学院本科学位论文11 绪论通常情况下，我们对一个线性动态电路进行瞬态分析往往采用的是经典法，即列写出电路的时域微分方程并求解,得出电感电流或者电容电压,然后再根据各支路间的约束关系进一步求得其它各支路的电压或电流。众所周知在经典法中，如果电路中只含有一个动态元件(电感或电容) ,那么列出的方程是一阶微分方程,其求解过程相对简单; 但如果电路中含有两个或两个以上的 动态元件，则所列出的方程将会是二阶或高阶微分方程，其求解过程会很麻烦，工作量很大。有鉴于此，人们又考虑采用积分变换的方法，利用拉普拉斯变换把已知的时域函数变换为频域函数，从而把时域的微分方程化为频域的代数方程。求出频域函数后，再作拉普拉斯反变换 ，返回 时域，从而求得满足电路初始条件的原微分方程的解，且不需要确定积 分常数 1。但 该方法同样 也存在运算量大的问题，尤其是对于含有多个动态元件的高阶复杂动态电路，并且计算得出的结果很不直观，无法清晰地看出在过渡 过程中各电量随时间变化的规律。本文就是从这个意义上出发，探索出一种线性电路瞬态分析的新方法伴随网络法。伴随网络法从建立电路方程开始，就 设 法避开微分方程。它把 动态电路的过渡过程时间 t0T 划分成若干时间间隔 t，把动态元件电感 L 和电容C 用相 应的离散模型来取代，经过代换后的电路称 为原电路的伴随网络。对于每一个时间间隔 t 而言，在伴随网络中，不再含有 动态 元件，取而代之的是动态元件 L、C 的离散模型，因而对该时间间隔下相应伴随网络的分析可视为稳态电路分析。也就是说，通过伴随网络法，可以将瞬 态电路分析 归结为一系列不同离散时刻下电阻网络的稳态分析。在对每一个离散时刻下的电阻网络进行方程列写时采用的是改进结点法，所列出的线性方程组采用的是电路方程的矩阵型式。其中结点电压方程矩阵的建立是采用的直接添加法，方程系数矩阵及右端项的元素并非整体形成，而是采用逐次扫描电路中的每一个支路元件，分别添加它们对方程的贡献的方法。建立方程伊始，方程系数矩阵及右端项全部为零元素所填充，每扫描一个元件，就将它对方程的贡献填到合适的位置，这样逐次扫描，逐次添加，直至网 络中每一条支路均被扫描，网络方程便被建立起来。当然整个网络方程的建立和求解都是借助于计算机并由MATLAB 软件编程实现的，充分利用了该软件 强大的数值计算功能和图形输出能力，力求更加方便直观。整个 电路矩阵的建立也采用的是框架式结构，即电路规模、 结点数、支路数以及各支路元件参数均由使用者自行输入，待分析支路的编号也由使用者输入，运行后即可直接显示出待求支路的电压和电流波形。同时为了使输入输出界面更加友好，我还利用MATLAB自带的GUI工具作2了一个图形用户界面，用户只需要按照界面的要求进行网络拓扑参数和待分析支路编号的输入，运行后即可在屏幕上显示GUI界面下的输出电压和电流波形。简单方便，清晰直观。2 线性动态电路分析的常用方法我们知道，对于动态电 路的分析常用的方法有经典法和拉普拉斯变换法。对于经典法，我们常用来求解一些低阶的微分方程。运用拉普拉斯变换法可以求解一些较高阶的微分方程。2.1 经典法对于仅含一个动态元件的简单电路,结合所给电路模型,根据KCL (基尔霍夫电流定律) 、KVL (基尔霍夫 电压定律)和支路的VCR (电压电流关系)建立描述电路的方程,建立的方程是以时间为自变量的线性常微分方程,然后求解常微分方程,从而得到电路所求变量(电压或电流)。2.2 拉普拉斯变换法对于具有多个动态元件的高阶线性动态电路,用直接求解微分方程的方法比较困难。例如对于一个n阶 方程,直接求解时需要知道变量及其各阶导数在t=0 +时刻的值,而电路中给定的初始状态是各电容电压和电感电流在t=0 +时刻的值,从这些值求得所需初始条件的工作量很大。于是人们通过拉氏变换,把已知的时域函数变换为频域函数,从而把时域的微分方程组转化为频域的代数方程组,求解代数方程组得到频域解后，再经 拉普拉斯反变换返回时域，可以得到满足电路初始条件的原微分方程组的解,而不需要确定积分常数。因此，拉普拉斯变换是求解高阶复杂动态电路的有效而且重要的方法之一。3 电路分析中常用的仿真软件目前国内外常用的电路仿真软件有 PSPICE、SABER 和 MATLAB 等一些软件。这些 软件都有各自强大的功能，在不同领域有着各自的特点,现在将其介绍如下:3.1 PSPICE 和 SABER3.1.1 PSPICE 简介用于模拟电路仿真的SPICE （Simulation Program with Integrated Circuit Emphasis）软件是于1972年由美国加州大学伯克利分校的计算机辅助设计小组利用FORTRAN 语言开发而成的，主要用于大规模集成电路的计算机辅助设计。河南科技学院机电学院本科学位论文3SPICE的正式版SPICE.2G在1975年正式推出，但是该程序的运行环境至少为小型机。1985年，加州大学伯克利分校用C 语言对SPICE软件进行了改写，并由MICROSIM公司推出。PSPICE是该公司在SPICE.2G版本的基础上升级并应用到PC机上的仿真软件。1998年著名的EDA商业软件开发商ORCAD公司与MICROSIM公司正式合并，自此MICROSIM公司的PSPICE产品正式并入ORCAD 公司的商业EDA 系统中。不久之后，ORCAD 公司正式推出了ORCAD PSPICE Release 10.5，与传统的SPICE软件相比，PSPICE 10.5在三大方面实现了重大 变革：首先，在对模拟电路进行直流、交流和瞬态等基本 电路特性分析的基础上，实现了蒙特卡罗分析、最坏情况分析以及优化设计等较为复杂的电路特性分析；第二，不但能够对模拟电路进行仿真，而且能够对数字 电路、数 /模混合电路进行仿真；第三，集成度大大提高，电 路图绘制完成后可直接进行电路仿真，并且可以随时分析观察仿真结果。现在国内使用较多的是PSPICE8.0版本，它能 够工作在 WINDOWS XP环境下，占用硬盘空间较小，整个软件有原理图编辑，电 路仿真，激励编辑，元器件编辑，波形 图等几部分组成，使用时是一个整体。 虽然目前 PSPICE软件在高校中已经非常流行，但是它也存在着明显的缺点。比如：仿真模型都是针对小功率电子器件的；软件采用变步长算法，对于周期性开关状态变化的电力电子电路而言，将耗费大量时间在步长上面；输出数据后的格式兼容性也不理想。3.1.2 SABER 简介SABER是美国Analogy公司开发、 现由Synopsys公司经营的系统仿真软件，是当今世界上功能强大的电力电子仿真软件之一。可用于电子、电力电子、机 电一体化、机械、光电、光学、控制等不同类型系统构成的混合系统仿真，这也是SABER的最大特点。SABER 作为混合仿真系统，可以兼容模拟、数字、控制量的混合仿真，便于在不同层面上分析和解决问题，我 们 从以下几个方面对SABER仿真软件进行介绍：(1).原理图输 入和仿真。SABER Sketch是SABER 的原理 图输入工具，通过它可以直接进入SABER仿真引擎。在SABER Sketch中，用户能够创建自己的原理图，启动 SABER完成各种仿真（偏置点分析、DC分析、 AC分析、瞬 态分析、温度分析、参数分析、傅立叶分析、蒙特卡诺分析、噪声分析、应力分析、失真分析等），可以直接在原理图上查看仿真结果，SABER Sketch及其仿真功能可以帮助用户完成混合信号、混合技术（ 电气、液 压等）系统的仿真分析。(2).数据可视 化和分析。Cosmos Scope是SABER的波形查看和仿真结果分析4工具，它的测量工具有50多种标准的测量功能，可以对波形进行准确的定量分析。它的专利工具波形计算器，可以对波形进行多种数学操作。(3).模型库。SABER拥有市场上最大的电气、混合信号、混合技 术模型库，它具有很大的通用模型库和较为精确的具体型号的器件模型，其元件模型库中有4700多种器件模型，500多种通用模型，能够满足航空、汽车和电源设计的需求。(4).建模。不同类型的设计需要不同类型的模型，SABER提供了完整的建模功能，可以满足各种仿真与分析的需求。其建模 语言主要有MAST、VHDLAMS、 Fortran，建模工具包括 StateAMS、5维的图表建模工具TLU，SABER可以对SPICE、 SIMULINK模型进行模型转换，同时SABER还拥有强大的参数提取工具，可以通过协同仿真实现模型复用。与PSPICE相比，SABER是功能更为强大的仿真软件，它可以仿真电力电子元件、电 路和系统，不仅具有PSPICE 的功能，而且具有更丰富的元件 库和更精致的仿真描述能力，还能结合数学控制方程模块工作。SABER还可以仿真电力传动、机械、热 力、流体等其他运动过程。SABER 的仿真真 实性很好，从仿真的电路到实际的电路实现，期间参数基本不用修改。但是SABER应用的主要困 难是软件操作较为复杂，市场价格高昂，比较适合于大型企业和科研院所应用，中小企业一般都是委托研究。3.2 MATLAB3.2.1 MATLAB 简介MATLAB 是矩阵实验室（Matrix Laboratory）的简称，它是由 MathWorks 公司开发出来的一套功能强大的工程应用软件，它和 Mathematica、Maple 并称为三大数学软件。它在 数值计算 方面的功能可以说是首屈一指。MATLAB 可以进行矩阵运算、 绘制函数和数据、实现算法、创建用户界面、 连接其它编程语言的程序等，主要应用于工程计算、控制设计、信号处理与通 讯、 图像处理、 信号检测、金融建模设计与分析等领域。MATLAB 软件主要包括以下 5 部分：MATLAB 开 发环境。 这是一套可以帮助你使用 MATLAB 函数和文件的工具，这些工具中的大多数都采用图形用户界面，它包括 MATLAB 桌面和命令窗口、命令历史、编辑器和调试器、帮助 浏览器、工作空间、文件和查询路径。MATLAB 数学函数库。这是一个涵盖了从求和、正弦、余弦这样一些基本函数，到矩阵求逆、矩阵特征值、贝赛尔函数和快速傅里叶变换等这样一些复杂函数的运算法则的大集合。MATLAB 语言。这是一种高级矩阵语言，它具有可以控制（数据）流状态、河南科技学院机电学院本科学位论文5函数、数据结构、输入输出和面向对象编程的特征。它允许快速创建小程序，也可以创建大的复杂的应用程序。图形处理。MATLAB 可方便地将矢量和矩 阵用图形来显示，也可注释和打印图形。它包括一些高级（ 绘图）函数， 这些高级函数可用于二 维和三维数据可视化、图像 处理、动画和图像再现。它也包括一些低 级 函数，可使用户定制图形界面和在你的 MATLAB 应用软件上创建完整的图形用 户界面。MATLAB 应用程序接口。这是一个允许你将用 C 语言和 Fortran 语言编写好的程序与 MATLAB 程序链接的库，它可以从 MATLAB 动态链接库中调用程序、将 MATLAB 作为一个 计算引擎、读或写 MAT 文件。MATLAB 主要有以下优点；（1）友好的工作平台和 编程环境；（2）简单易用的程序语言；（3）强大的科学计算及数据处理能力；（4）出色的图形处理功能；（5）应用广泛的模块集合工具箱；（6）实用的程序接口和发布平台；（7）应用软件开发（包括用户界面）。3.2.2 MATLAB-GUI 简 介GUI 的全称是 Graphical User Interfaces，其汉语意思是 图形用户界面。是指由窗口、光标、按键、菜 单、文字 说明等对象（Objects）构成的一个用户界面。用户通过一定的方法（如鼠标或键盘）选择、激活这些图形对象，使计算机产生某种动作或变化，比如：进行计算、绘图等。假如 读者所从事的数据分析、解方程、计算结果可视工作比较单一，那么一般不会考 虑 GUI 的制作。但是如果读者想向别人提供应用程序，想进行某种技术、方法的演示，想制作一个能 够让用户反复使用且操作简单的专用工具，那么图形用户界面也许是最好的选择之一。这些功能是 MATLAB 的一大亮点。本文就是将 电路分析的 结果做一个 GUI。可以使用户很方便直观的去观看运算的结果。6图 1. 伴随网络法的 GUI 界面1.支路选择对话框的功能如图1所示,左边的一个对话框上面写着请选择支路。我们把它称作支路选择对话框,它是用MATLAB-GUI 中的listbox制作而成的。GUI中的listbox是一个列表框，其显 示是由srting属性定 义的一系列列表选项,用 户能够选择其中的一个或者多个选项。在本设计中其使用方法如下：在整个程序运行时,会提醒用户将电路的相关规模和数据输入程序。然后计算机将会自动地运算,并且把运算的结果存储起来。当程序运行到调用 GUI 界面时。支路 选择对话框会自行回 调已经存储好的电路中各个支路的编号。即支路对话框中将会列出各个支路的编号。在本设计中我们定义一次只能选择列表项中的一个选项，来避免选项间 的冲突。所以用 户在查看 电路的波形时只能选择其中的一个支路。如果选择多个的 话,系统会提示你选择方式错误。2.坐标图的功能在GUI界面中,右边的是一个坐标平面,我们称它为坐标图。坐标图是用GUI中的axes 制作而成的。axes 可以调用每一个支路的电压或电流数据。并将其画成波形。坐标轴的横轴为时间轴 ,其纵轴的单位是电压v或是电流A。画出的波形一般比较的平滑。用户可以清楚的看到各个时刻的电压或电流值,并且总体的电压或电流走向也一目了然。当然具体是电压波形还是电流的波形,这就要看你选择了。Axes是一个二维的坐标图形,其调用的数据是经过运算保存后的所选择支路的数据,这些数据可以通过一个命令plot画成一个图形。并且通过选择GUI下面的电压或电流按钮,程序将自动改变坐标轴的单位。河南科技学院机电学院本科学位论文73.电压电流选择对话框电压电流选择对话框的功能是选择所要查看的电压波形或电流波形。该对话框是由一个button group和两个radio button所组 成。button group是一个按钮组合框。组 合框是一个图形封 闭区域,它可以把相关的控件(单选按钮或多选按钮等)组合在一起,使得界面的按钮功能和相互间的关系更容易理解。Button group 中包含了两个 radio button。radio button 是单选按钮。单选按钮一般是以组出现的。一组单选 按钮之间的关系是相互排斥的,也就是说任何时候只有一个按钮有效。在本设计 中只用到了两个单选按钮。上面的一个是电流波形选择按钮,下面一个是电压波形选择按钮。在用户没有做出选择的时候,系统默认的是电压波形。让 axes 回调所选支路的电压数据,并 给坐标轴辅以电压单位。当用户点击选择电流波形时,将触动 axes 回调所选支路的 电流数据,并给坐标轴以电流单位。4 伴随网络法伴随网络法的数学模型是把动态电路的过渡过程时间t 0 T 划分成若干时间间隔t, 把动态 元件电感L 和电容C用相应的离散模型来代替，经过代换后的电路称为原电路的伴随网络电路。 对于每一个时间间隔t而言,在伴随网络中,我们认为不再含有动态元件,取而代之的是动态元件L, C的离散模型,因而对该时间间隔下相应伴随网络的分析可视为稳态电路分析。也就是说,通过伴随网络法,可以将瞬态电路分析归结为一系列不同离散时刻下电阻网络的稳态分析。当然这种思路的背景主要是基于计算机辅助分析的基础上，因为在这种过程中要重复计算很多次，单靠人工是无法实现 的，不 过通过计算机就可以既方便又快捷的算出结果。4.1 动态元件的离散模型4.1.1 电容 C 的伴随网络模型图2.电容C及其伴随网络模型如图 2 所示，图中左边电容 C 的电压 u（t）与电流 i(t)为关联参考方向，则：+-Ci(t)u(t) +-Un+1In+1G=2C/h(2C/h)Un+I8i(t)=Cdu(t)/dt （1）在某一离散时刻 tn+1,in+1=Cun+1 （2）我们采用梯形积分公式：un+1=un+(un+un+1)h/2 （3）故而由式（2）、 （3）可得：in+1=un+12C/h-(un2C/h+in) （4）由此得出图 2 中右边所示电容的伴随网络模型，该离散模型是由一个电导和一个电流源并联而成。电导 G 的取值为 2C/h，电 流源的取值为(2C/h)u n+in，其中参数 C 为电 容值（法拉）， h 为划分后的时间间隔 t（秒）。4.1.2 电感 L 的伴随网络模型图3.电感L及其伴随网络模型如图 3 所示，图中左边电容 L 的电压 u（t）与电流 i(t)为关联参考方向，则：u(t)=Ldi(t)/dt （5）在某一离散时刻 tn+1，un+1=Lin+1 （6）我们再次采用梯形积分公式：in+1=in+(in+in+1)h/2 （7）故而由式（6）、 （7）可得：in+1=un+1h/2L+(unh/2L+in) （8）由此得出图 3 中右边所示电感的伴随网络模型，该离散模型同样是由一个电导和一个电流源并联而成。 电导 G 的取值为 h/2L，电流源的取值为(h/2L)un+in，其中参数 L 为电感值（亨利），h 为划分后的时间间隔 t（秒）。4.2 伴随网络法分析瞬态电路的步骤(1)根据所给电 路的网络数据,对初始状态进行分析;(2)根据网络结 构预估过渡过程时间T,并将t 0T 划分(通常可等分)为若干时间间隔t = h,取时间步长为h,形成t 1 = h时的伴随网络,计 算出t 1时刻的网络响应;+-L+-Un+1In+1G=h/2L(h/2L)Un+Ii(t)u(t)河南科技学院机电学院本科学位论文9(3)根据上一时 刻计算出的结果,修正伴随网络的参数,从而获得下一时刻的伴随网络;(4)利用步骤(3)形成的伴随网 络,建立电路方程并求解,得出该时刻的响应;(5)如果 t2 T,分析结束,输 出结果;否则增加时间步长 h,转向步骤(3) 。5 伴随网络法的 MATLAB 编程实现本文编程采用的是目前电路分析领域使用较多的软件MATLAB, 该软件具有强大的数学运算功能和图形处理能力,其数据处理十分有效、精细,运行速度较快,且数据格式兼容性较好,便于数据的后处理与分析。105.1 程序流程图5.1.1 主程序流程图图 4.伴随网络法主程序流程图伴随网络法主程序流程图如图 4 所示，另外还有高斯消元子程序这里就不河南科技学院机电学院本科学位论文11再赘述了，所对应的 MATLAB 源程序在文后的附录中给出。5.1.2 GUI 程序流程图图 5.GUI 程序流程图5.2 算例分析有一线性网络如图 4 所示,如果电容元件的初始值为 UC(0+)= 0V,电感元件的初始值 iL(0+)=0A, R1=R2=R3=2, gm=1S, L=1H, C=2F。开关 K 在 t=0 时刻闭合,且 K闭合前电路已处于稳态。求: t0 时各支路的电压和 电流。12R1R2UR3CL+_ E+ _ gmUK0图6.算例电路图5.2.1 输入参数从网络结构我们可以看出,该电路共有5个结点, 7个支路(每个元件均算作一个支路) 。其中包含3个电阻, 1个独立电压源, 1个受控电流源,另外还有两个动态元件电容C 和 电感L。运行程序后 ,屏幕提示信息如下 :请输入电路规模请输入待求网络的节点数: 5n = 5请输入待求网络的支路数: 7b = 7请输入待求网络的动态元件数: 2q = 2请输入过渡过程的时间(单位:秒) : 15T = 15请输入第1支路的始节点,终节点,元件类型,参数1,参数2,控制支路始节点,终节点(没有则 置零) :- 1, 2, R, 2, 0, 0, 0请输入第2支路的始节点,终节点,元件类型,参数1,参数2,控制支路始节点,终节点(没有则 置零) :- 2, 3, R, 2, 0, 0, 0请输入第3支路的始节点,终节点,元件类型,参数1,参数2,控制支路始节点,终节点(没有则 置零) :- 3, 0, R, 2, 0, 0, 0请输入第4支路的始节点,终节点,元件类型,参数1,参数2,控制支路始节点,终节点(没有则 置零) :- 1, 0, E, 5, 0, 0, 0请输入第5支路的始节点,终节点,元件类型,参数1,参数2,控制支路始节点,终节点河南科技学院机电学院本科学位论文13(没有则 置零) :- 2, 0, VCCS, 1, 0, 2, 3请输入第6支路的始节点,终节点,元件类型,参数1,参数2,控制支路始节点,终节点(没有则 置零) :- 4, 0, C, 2, 0, 0, 0请输入第7支路的始节点,终节点,元件类型,参数1,参数2,控制支路始节点,终节点(没有则 置零) :- 3, 4, L, 1, 0, 0, 0当所有支路元件信息均输入后,程序开始进入循环,并给出各时间间隔内所对应的结点电压和困难支路电流。最后,程序会调用GUI,并显示各个支路的编号，以便供用户查看各个支路。5.2.2 GUI 输出参数本题以电容C所在第六支路为例进行仿真，其 结果图 7和图8所示。左边为所有支路，我们选中了第六支路，右边为电压电流的波形 图。用 户可以根据需要进行选择所要查看的电压或电路波形，图下方有一个返回健。用户在查看完之后就可以返回。根据实际的计算 ,我们发现仿真的结果和实际计算的结果完全一致。当然我们还可以把该方法应用到实际的教学过程中，相信教学的效率和效果将会大大的提高，这就是多媒体教学的特色。14图7. GUI的电压波形仿真图图8.GUI的电流波形仿真图6 结语本文是以MATLAB GUI为电路分析平台， 应用微积 分的数学思想，把高阶线性动态电路的瞬态分析转化为若干个连续的微小时间段内的线性电阻电路的稳态分析。 单就方法而言,伴随网络法更方便、直 观。而且将MATLAB可视化的功能应用于程序中，可以使用户 直观体验到电网络分析的最终结果，程序独到的编写使得网络拓扑可以由使用者自行输入,因此也具有较好的灵活性。电路元件涵盖了电路分析中常见的电阻(导) ,独立电压源,独立电流源,四种受控源,电感和电容等元件,过渡过程的时间以及采样间隔亦由使用者自行定义,并且可以由使用者选择需要分析的支路进而输出该支路的电压、电流波形,波形显示清晰直观。在程序河南科技学院机电学院本科学位论文15调试过程中,可以根据所选电路的网络拓扑结构预先估计一下过渡过程的时间和对应的采样间隔,选取几组不同的值试运行一下,从中挑出最优解。当然设计中也存在着一些不足和需要进一步完善的地方，如：支路元件中没有考虑到耦合电感元件，所采用的激励源均为直流电源等。今后可考虑将交流电源也添加进来，再加入耦合 电感， 这样的话该设计就将更加完美。致谢首先感谢我的指导老师。从论文的选题，研究思路，知 识的创新，论文的开题，撰文和修改等过程，无不浸透着导师的心血和汗水，也蕴含着导师对我的深切关怀。 论文所取的成果渗透着导师的辛劳，也凝聚着导师的智慧。 导师渊博的学识，敏 锐的科学洞察力，严谨的治学态度，强烈的 创新思想和对科研教育事业的执著追求都给我留下难忘的印象，并将激励我在今后的工作中勇敢地面对困难和挑战。此外本论文的顺利完成也离不开各位同学和朋友的关心和帮助，在此感谢他们一直以来对我的支持和帮助。16参考文献1 邹玲.电路理 论M.武 汉：华中科技大学出版社， 20062 楼顺天,陈生潭,雷虎民. MATLAB 5. x程序设计语言M. 西安: 西安电子科技大学出版社, 2000:1163 程少庚,崔杜武,刘小河.电网络分析M.西安:西安理工大学出版社, 2003: 65-68.4 童梅.电路的计算机辅助分析Pspice和MATLABM.北京:机械工业出版社, 2005:7-8.5 李培芳, 孙盾,李玉玲,李军. 电路学习辅导与题解M.北京:清华大学出版社, 2006:77-78.6 游磊,梁颖 .Matlab与Pspice在电路分析中的应用.成都航空职业技术学院学报J.2004(3) : 7 黄忠霖.控制系 统 MATLAB 计算及仿真M.北京：国防工业出版社，20068 张永瑞.电 路基础教程M.北京：科学出版社， 20059 邱关源.电 路(第四版)M.北京:高等教育出版社 ,2003:289-290.10Charles K.Alexander. Fundamentals of Electric CircuitsM.北京:清华大学出版社,2000.11周金萍.Matlab 6实践与提高M.北京:中国电力出版社,2002.12陈怀琛.MATLAB 6.1实 用指南M.北京:电子工业出版社,2000.13施阳,李俊 .MATLAB语 言工具箱TOOLBOX 实用指南M. 北京: 清华大学出版社,1998. 14薛定宇,陈 阳泉.基于MATLAB/Simulink的系统仿真技术与应用.清华大学出版社,200215James W.Nilsson, Susan A.Riedel. Introductory Circuit for Electrical and Computer Engineering. Publishing House of Electronic Iindustry,2007.1617河南科技学院机电学院本科学位论文17181920附录1. 伴随网络法 MATLAB 主程序：%*用伴随网络法（梯形公式）对线性电路进行瞬态分析*%clc %清除命令窗口，并使光标返回clear %清除工作空间的变量%*原始数据的输入和数据库的建立（各支路电压和电流均取关联参考方向*%输入格式： 支路的始结点， 终结点，元件类型，参数，控制支路始 结点， 终结点（没有则置零）回车disp(请输 入电路规模);n=input(请输入待求网络的结点数:);b=input(请输入待求网络的支路数:);q=input(请输入待求网络的动态元件数:);T=input(请输 入过渡过程的 预估时间（单位:秒）:);e_net=cell(b,7); %建立单元阵列%输入 时先输入非动态元件所在支路参数，然后再输入电容和电感所在支路参数%输入 顺序为：电阻( 导)，电压源，电流源， CCCS，VCCS，VCVS，CCVS，然后是电容和电感c=请输入第;d=支路的始节点，终节点，元件 类型，参数1,参数2,控制支路始节点,终节点(没有则置零):;for i=1:bdisp(c,num2str(i),d);e_net(i,:)=input(-);endFs=out;in;st;vl;vl0;c_out;c_in;database=cell2struct(e_net,Fs,2); %将单元阵列转换为结构阵列并建立支路元件数据库%*用改进节点法建立待求网络方程*%Y=zeros(n-1,n-1);18V=zeros(n-1,1);Y1=zeros(n-1,n-1);V1=zeros(n-1,1);Q=zeros(n-1,1);%*首先添加非动态元件的贡献进入节点电压方程*%for i=1:b-qk=database(i).out;l=database(i).in;v=database(i).vl;g=database(i).c_out;h=database(i).c_in;sdim=size(Y,1);switch database(i).stcase R %电阻(导)元件所在支路对伴随网络方程的贡献if l=0Y(k,k)=Y(k,k)+1/v;elseY(k,k)=Y(k,k)+1/v;Y(l,l)=Y(l,l)+1/v;Y(k,l)=Y(k,l)-1/v;Y(l,k)=Y(l,k)-1/v;endcase E %独立电压源所在支路对伴随网络方程的贡献Y(k,sdim+1)=1;Y(sdim+1,k)=1;V(sdim+1,1)=v;if l=0Y(l,sdim+1)=-1;Y(sdim+1,l)=-1;endcase I %独立电流源所在支路对伴随网络方程的贡献V(k,1)=V(k,1)-v;if l=0V(l,1)=V(l,1)+v;end河南科技学院机电学院本科学位论文19case CCCS %流控流源所在支路对伴随网络方程的贡献Y(k,sdim+1)=v;Y(g,sdim+1)=1;Y(sdim+1,g)=1;if h=0Y(sdim+1,h)=-1;Y(h,sdim+1)=-1;endif l=0Y(l,sdim+1)=-v;endcase VCCS %压控流源所在支路 对伴随网络方程的贡献Y(k,g)=Y(k,g)+v;if h=0Y(k,h)=Y(k,h)-v;endif l=0Y(l,g)=Y(l,g)-v;endif (l=0 & h=0)Y(l,h)=Y(l,h)+v;endcase VCVS %压控压源所在支路对伴随网络方程的贡献Y(k,sdim+1)=1;Y(sdim+1,k)=-1;Y(sdim+1,g)=v;if h=0Y(sdim+1,h)=-v;endif l=0Y(l,sdim+1)=-1;Y(sdim+1,l)=1;endcase CCVS %流控压源所在支路 对伴随网络方程的贡献Y(k,sdim+2)=1;20Y(sdim+2,k)=1;Y(g,sdim+1)=1;Y(sdim+1,g)=1;Y(sdim+2,sdim+1)=-v;if h=0Y(h,sdim+1)=-1;Y(sdim+1,h)=-1;endif l=0Y(l,sdim+2)=-1;Y(sdim+2,l)=-1; endendendsdim1=size(Y,1);sdim2=size(V,1);if sdim2pV1=Q;endfor j=b-q+1:b %将动态元件的贡献添入伴随网络方程，并开始循环k=database(j).out;l=database(j).in;v1=database(j).vl;v2=database(j).vl0;河南科技学院机电学院本科学位论文21g=database(j).c_out;h=database(j).c_in;sdimy=size(Y,1);switch database(j).stcase C %电容所在支路对伴随网络方程的贡献if (t=0 & l=0)Y(k,sdimy+1)=1;Y(l,sdimy+1)=-1;Y(sdimy+1,k)=1;Y(sdimy+1,l)=-1;V(sdimy+1,1)=v2;elseif (t=0 & l=0)Y(k,sdimy+1)=1;Y(sdimy+1,k)=1;V(sdimy+1,1)