下载文档原格式
计算机辅助电路仿真技术实验报告篇一：电路计算机仿真实验报告电路计算机仿真分析实验报告实验一直流电路工作点分析和直流扫描分析一、实验目的1、学习使用Pspice软件，熟悉它的工作流程，即绘制电路图、元件类别的选择及其参数的赋值、分析类型的建立及其参数的设置、Probe窗口的设置和分析的运行过程等。
2、学习使用Pspice进行直流工作点分析和直流扫描分析的操作步骤。
二、原理与说明对于电阻电路，可以用直观法（支路电流法、节点电压法、回路电流法）列写电路方程，求解电路中各个电压和电流。
PSPICE软件是采用节点电压法对电路进行分析的。
使用PSPICE软件进行电路的计算机辅助分析时，首先在capture环境下编辑电路，用PSPICE的元件符号库绘制电路图并进行编辑、存盘。
然后调用分析模块、选择分析类型，就可以“自动”进行电路分析了。
需要强调的是，PSPICE软件是采用节点电压法“自动”列写节点电压方程的，因此，在绘制电路图时，一定要有参考节点（即接地点）。
此外，一个元件为一条“支路”（branch）,要注意支路（也就是元件）的参考方向。
对于二端元件的参考方向定义为正端子指向负端子。
三、示例实验应用PSPICE求解图1-1所示电路个节点电压和各支路电流。
图1-1 直流电路分析电路图R22图1-2 仿真结果四、选做实验1、实验电路图（1）直流工作点分析，即求各节点电压和各元件电压和电流。
（2）直流扫描分析，即当电压源Us1的电压在0-12V之间变化时，求负载电阻RL中电流IRL随电压源Us1的变化曲线。
I图1-3 选做实验电路图2、仿真结果图1-4 选做实验仿真结果3、直流扫描分析的输出波形图1-5 选做实验直流扫描分析的输出波形4、数据输出V_Vs1 I++00 ++00 ++00 ++00 ++00 ++00 ++00 ++00 ++00++00 ++00 ++00 ++00从图1-3可以得到IRL与USI的函数关系为：IRL=+US1=+五、思考题与讨论：1、根据图1-1、1-3及所得仿真结果验证基尔霍夫定律。
四川⼤学电⼒系统计算机辅助分析实验报告⽬录1.潮流计算基础实验 (2)⼀：实验名称(Title) (2)⼆：实验⽬的(Purpose) (2)三：实验⼯具/材料(Tools/Materials) (2)四：实验内容（Content） (2)五：实验模型（Model） (3)六：实验结果(Results) (4)七：实验分析(Analysis) (6)1.母线电压幅值及相⾓ (7)2.交流线结果 (7)3.功率及损耗 (7)2.潮流计算分析实验 (8)⼀：实验名称(Title) (8)⼆：实验⽬的(Purpose) (8)三：实验内容（Content） (8)四：实验数据/结果(Data/Results) (9)1. 不同潮流迭代算法的仿真效率 (9)2. 调整给GEN2的电压幅值 (10)3. 调整GEN2的有功输出 (11)4. 调整STNC-230负荷的有功消耗 (13)五：实验分析/讨论(Analysis/Discussion) (14)1. 分析不同潮流迭代算法的仿真效率 (14)2. 调整GEN2节点电压幅值，分析其对潮流的影响 (14)3. 调整发电机GEN2的有功输出，分析其对潮流的影响 (17)4. 调整STNC-230负荷的有功消耗，分析其对潮流的影响 (20)六：实验⼼得（Experience ） (21)1.潮流计算基础实验⼀：实验名称(Title)WSCC 9节点潮流计算基础实验⼆：实验⽬的(Purpose)1. 掌握潮流计算的概念、原理和计算数据要求；2. 熟练使⽤ PSASP 建⽴电⼒系统的潮流计算模型，并完成潮流计算；3. 掌握潮流计算结果数据整理和分析。
三：实验⼯具/材料(Tools/Materials)电⼒系统分析综合程序（Power System Analysis Software Package）电⼒系统分析综合程序（Power System Analysis Software Package）简称PSASP。
电力系统计算机辅助分析综合训练实验报告姓名：学科、专业：电气工程及其自动化学号：完成日期：2012.7.15大连理工大学Dalian UniversityofTechnology目录实验一PSASP的基本使用 (1)实验二5节点系统的潮流计算与PSASP验证 (4)实验三IEEE14节点系统的潮流计算与潮流调整 (16)实验四潮流调整 (26)实验五对称短路计算和不对称短路计算 (39)实验一PSASP的基本使用一、实验目的学习掌握PSASP仿真软件图模一体化的操作界面，掌握电力系统各元件的数学模型，为进行电力系统分析计算打下基础。
二、实验内容利用一个简单5节点系统首先进行潮流估算，再以PSASP仿真软件进行仿真计算，通过图形和文本两种方式进行仿真，验证计算结果的正确性。
三、实验原理《电力系统分析综合程序》(PowerSystemAnalysisSoftwarePackage,PSASP)是一套历史长久、功能强大、使用方便的电力系统分析程序，是高度集成和开放具有我国自主知识产权的大型软件包。
PSASP基于电网基础数据库、固定模型库以及用户自定义模型库的支持，可进行电力系统(输电、供电和配电系统)的各种计算分析。
全国各省市、香港地区电力规划设计、生产调度运行、科研教育等超过400家用户，应用于多项大型电力系统工程计算分析，应用于多所大学作为科研和教学的有力工具。
功能：潮流计算、短路计算、暂稳计算、小扰动稳定、电压稳定、静态安全分析、最优潮流等。
PASASP7.0图模一体化平台，“图”有单线图、地理位置接线图、厂站主接线图；“模”指模型及数据。
可边绘图边建数据，也可根据已有数据进行图形的快速绘制。
PASASP将数据分为四类：基础数据、计算数据、结果数据、用户自定义数据等。
基础数据：基于物理特性的电力系统元件库和元件公用参数库计算数据：根据不同计算的需要，与计算密切相关的数据结果数据：各种计算结果的数据，与相关计算数据相对应实验步骤：1、建立网络结构，输入元件参数与发电和负荷数据（基础数据可分数据组）。
四川大学电气信息学院实验报告书课程名称：电力系统分析的计算机算法实验项目：暂态稳定实验专业班级：电气工程及其自动化专业103班级实验时间：2017年5月12日气信息学院专业中心实验室目录一：实验目的(Purpose)1.掌握暂稳计算的概念、原理和计算数据要求；2.熟练使用PSASP 建立电力系统的暂稳计算模型，并完成暂稳计算；3.掌握暂稳计算结果的数据图形化整理和分析；4.利用PSASP，分析电力系统暂态稳定受故障条件、受不同运行方式的影响及其规律；5.通过计算机仿真，巩固《电力系统分析理论》所学，对“电力系统暂态稳定”及其影响因素加深理论和实践认识；6.学习技术图形的绘制二：实验内容（Content）1.不同故障地点对该系统暂稳的影响2.不同故障类型对该系统暂稳的影响3.不同故障持续时间对该系统暂稳的影响4.调整系统的初始运行方式对该系统暂稳的影响三：实验数据/分析(Data/Analysis)Ⅰ.系统的初始运行方式不变1.不同故障地点对该系统暂稳的影响故障类型：三相短路故障线路：GEN2-230~STNC-230故障开始时间：1s切除时间分别选取不同故障地点30%，50%，70%，观察功角差及有功、频率变化情况。
最后每一种情况再确定极限切除时间。
a)功角差变化情况时的功角差摇摆曲线；图2是在线路50%处故障时的功角差摇摆曲线；图3是在线路70%处故障时的功角差摇摆曲线。
从仿真曲线可以看出s切除故障后，故障在线路50%处和在70%处，系统均能保持暂态稳定，但线路50%切除故障的相对功角变化幅度更大；而在线路30%处切除故障后，功角差无限增大，系统将失去暂态稳定。
b)有功，频率变化情况不同切除故障时间的时域仿真结果如图4-图12所示，图4-图6是在线路30%处故障时的各发电机的有功、频率变化曲线；图7-图9是在线路50%处故障时的各发电机的有功、频率变化曲线；图10-图12是在线路70%处故障时的各发电机的有功、频率变化曲线。
昆明理工大学《电力系统计算机辅助分析》上机实验（指导书）主编唐岚电力工程学院二〇〇七年一月目录前言 ................................................................................... I I 实验一MATLAB软件的基本操作 (1)实验二MATLAB 程序的基本结构 (2)实验三电力系统计算中常用的数值算法 (4)实验四电力网络的数学模型 (5)实验五潮流计算 (6)实验六短路电流计算 (7)实验七静态稳定和暂态稳定计算 (8)附录：学生实验报告表头格式 (9)前言《电力系统计算机辅助分析》是电气工程及自动化专业的专业必修核心课程，是《电力系统分析基础》的后续课程。
其主要任务是：使学生深入学习电力系统潮流、短路、稳定计算的计算机算法，提高学生应用计算机对电力系统进行分析和计算的能力。
《电力系统计算机辅助分析》课程教学的难点在于如何让学生掌握将《电力系统分析基础》中所学的各种电力系统元件的数学模型用网络方程联系起来并转化为程序代码的方法，进而用相应的数值算法求解之。
这要求学生学好先修课程《电力系统分析基础》、《计算方法》和《程序设计基础》。
但由于后两门课程开设较早，学生掌握情况差异也很大，所以极不利于教学。
因此，我们结合选用的《电力系统分析》教材，以MATLAB作为软件工具来实现各种电力系统计算程序。
这样可以利用MATLAB软件强大的数值计算能力和相对简单的编程语言，以减少教学过程中的不利因素，在促进学生掌握《电力系统计算机辅助分析》核心内容的同时，也让学生初步学习了MATLAB这个在科学和工程领域应用十分广泛的软件。
作为《电力系统计算机辅助分析》课程的辅助教学材料，其内容从MATLAB 软件的基本操作，到电力系统三大基本计算，共由七个上机实验组成。
所有七个实验均为课程教学过程中必做的基本实验。
下一步将考虑加入初步制订电力系统运行方式的综合实验，可供课程设计时选做。
昆明理工大学《电力系统计算机辅助分析》上机实验（指导书）电力工程学院二〇一四年十月目录实验一MATLAB软件的基本操作及程序的基本结构.2实验二电力系统计算中常用的数值算法及电力网络的数学模型 (7)实验三短路电流计算 (10)实验四潮流计算 (9)实验五静态稳定和暂态稳定计算 (12)附录：学生实验报告表头格式 (13)实验一 MATLAB 软件的基本操作及程序的基本结构 上机实验目的1、了解MATLAB 软件的基本构成并熟悉集成界面菜单命令与各工作窗口；2、掌握MATLAB 的两种运行方式；3、学习程序编辑器的使用，掌握常用数据类型与运算符；4、掌握m 文件的两种格式：脚本与函数；5、掌握顺序结构、循环结构、分枝结构和程序流程控制语句；6、掌握用M 文件编辑器进行简单调试的方法。
7、掌握用MATLAB 绘制二维图形的基本方法。
预习要求1、利用图书馆和互联网资源，了解MATLAB 软件的基本信息。
列举出MATLAB 软件的主要特点和应用领域。
2、列出实现三种结构的和程序流程控制的关键字和语法规则；3、学习MATLAB 函数的编写方法。
上机实验内容1、启动MATLAB 软件，识别出常用的四个窗口――命令窗口（Command Window ）、历史命令窗口（Command History ）、工作空间浏览器（Workspace Browser ）。
熟悉Matlab 工作环境。
2、上机求解下面的电路问题，脚本文件已经给出。
已知某个ＲＣ电路的端电压的表达式为：26,0~10t u e t -==区间时，试绘制电压的波形。
提示：在MATLAB 的M 文件编辑器中键入以下命令语句，并保存为ex.m ：％％％％％％％％％％％％％％％％％％％％％％％％％％％％％％％％％％％％％Clear,clc,close t=0:0.5:10; u=6*exp(-2*t);plot(t,u);title(‘RC 电路的电压响应曲线’);xlabel(‘时间/s’);ylabel(‘电压/v’) grid%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%5、 对于如下电路图，已知()200cos(90)os u t t ω=+V ，122R =Ω，50L mH =，233R =Ω，100C F μ=，自己编写程序，画出()c u t 和()s u t 的波形（四个周期）。
南京理工大学电力系统计算机辅助分析实验报告作者: 赵朋学号：0608190255 学院(系):动力工程学院专业: 电气工程及其自动化指导教师：郭新红2009 年5 月1 日实验一七节点电力系统数据库建立一实验目的(1)熟悉PSASP系统的操作环境(2)掌握在文本方式下数据建立和编辑(3) 熟悉系统单线图的编辑二实验内容根据实验指导书内容创建数据库，并完成系统单线图的绘制。
三实验原理七节点系统图如下所示：四实验步骤1 创建数据库目录进入PSASP主画面，点击“创建”，输入数据目录“D:/ZJS”，系统基准容量100兆伏安。
2 录入基础数据进入文本方式数据编辑环境，根据所给参数建立基础数据和公用参数。
根据实验指导书提示完成数据组选择，母线、交流线、两绕组变压器、直流线的相关参数的录入，发电机及其调节器数据的建立，并在PSS栏中编辑参数，然后建立区域数据。
3 系统单线图编辑进入图形方式数据编辑环境，在单线图编辑窗口中，绘制各种元件。
依照单线图的绘制方法，按以下顺序完成该系统单线图的绘制：①7条母线②4条交流线和2个电抗器③4台变压器④1条支流线⑤3台发电机⑥1个负荷实验二潮流计算一实验目的(1)理解PSASP潮流计算的流程；(2)掌握在文本和图形两种方式下潮流计算的方法；(3)掌握潮流计算结果输出方法。
二实验原理潮流计算实验是利用PSASP的潮流计算功能实现的，它的流程和结构如下：三实验步骤第一步：建立七节点系统的数据库由于在实验一中已经建立了七节点系统数据库，实验二是在实验一的基础上进行的，因此这里可以直接利用实验一的数据库建立潮流计算作业。
第二步：建立潮流计算作业1.定义以下2个方案：2.定义潮流作业作业1：一般方式的潮流计算作业2：具有控制功能的潮流计算3.在文本和图形方式下进行潮流计算①进行一般方式的潮流计算1②进行具有控制功能的潮流计算24.潮流计算结果输出将作业1和作业2的结果先以报表形式输出至文件保存,具体要求见实验指导书。
电力系统分析的计算机算法实验报告实验项目 : PSD-BPA潮流计算实验学院 : 电气信息学院专业 :任课老师:班级: 2010级学号:姓名:图 9-20 三机系统网络结构图表9-2 网络参数表9-3 节点参数表9-4 发电机参数一、根据表格中的数据对系统的IEE90文件进行修改程序清单：.=========================================================================.一、文件头部分一般为控制语句，用于指定：.1。
程标识（如"IEEE9"）以及工程名称（如"IEEE_9BUS_TEST_SYSTEM"）；(POWERFLOW,CASEID=IEEE9,PROJECT=IEEE_9BUS_TEST_SYSTEM).2。
法选项，这里将采用两次PQ分解法，然后转入牛顿_拉夫逊法，最大迭代次数为15次/ SOL_ITER,DECOUPLED=2,CURRENT=0,NEWTON=20,OPITM=0 \.3。
定基准功率/ MVA_BASE=100. \.4。
成的文件名称，如用于稳定计算的BSE文件名"IEEE90.BSE"以及用于画图的MAP文件名为"IEEE90.MAP"/ NEW_BASE,FILE=IEEE90.BSE \/ PF_MAP,FILE=IEEE90.MAP \./ NETWORK_DATA, FILE=, RXCHECK=ON\.5。
输出选项，这里将全部输出潮流的原始数据、计算结果/ P_INPUT_LIST,FULL \/ P_OUTPUT_LIST,FULL \.6。
析模式选择，这里按第四级来选择，具体的含义参见手册/ P_ANALYSIS_REPORT,LEVEL=4 \.==========================================================================.二、网络数据部分，包括节点数据、支路数据、变压器数据等等.==========================================================================..1.节点数据，包括普通PQ节点、PV节点以及平衡节点等等BS GEN1 16.501 999. 999. 1.01BE GEN2 18.001 163. 999. 1.01BE GEN3 13.801 85. 999. 1.01B GEN1-230230.01B STNB-230230.02 90. 30.B GEN3-230230.01B STNC-230230.01 100. 35.0 20.B GEN2-230230.01 10.B STNA-230230.01 125. 50. 20..2.变压器数据,包括普通的变压器、移相器、带调节的变压器等等T GEN1 16.5 GEN1-230230. .0576 16.5 242.T GEN2 18.0 GEN2-230230. .06250 18.0 242.T GEN3 13.8 GEN3-230230. .0586 13.8 242..3.支路数据，同一路径上的平行线路需要用回路标志来区分L GEN1-230230. STNB-230230. .0 170 .0920 .1580 L GEN1-230230. STNA-230230. .0100 .0850 .1760 L STNB-230230. GEN3-230230. .0390 .1700 .3580 L STNA-230230. GEN2-230230. .0320 .1610 .1530 L GEN3-230230. STNC-230230. .0119 .1008 .2090 L GEN2-230230. STNC-230230. .0085 .0720 .1490 .==========================================================================.三、结束.==========================================================================/CHANGES\PZ 01 1.0 1.0(END)二、执行计算，输出结果文件最后的输出结果：详细输出列表从输出结果可以看出最后结果收敛，具体数据已经详细算出。
附录：学生实验报告表头格式电力工程学院学生实验报告实验课程名称：《电力系统计算机辅助分析》开课实验室：计算中心408 2014 年 12 月 11 日legend('us','uc')%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%一、讨论1、MATLAB软件是什么？有什么特点？答：Matlab是一个高级的矩阵/阵列语言，它包含控制语句、函数、数据结构、输入和输出和面向对象编程特点。
用户可以在命令窗口中将输入语句与执行命令同步，也可以先编写好一个较大的复杂的应用程（M文件）后再一起运行。
特点：1. 编程效率高，2. 用户使用方便，3. 扩充能力强，交互性好，4. 移植性和开放性很好，5. 语句简单，内涵丰富，6. 高效方便的矩阵和数组运算，7. 方便的绘图功能。
3. 如何获得命令或函数的帮助文档？答：MATLAB 的各个函数，不管是内建函数、M文件函数、还是MEX文件函数等，一般它们都有M文件的使用帮助和函数功能说明，各个工具箱通常情况下也具有一个与工具箱名相间的M文件用来说明工具箱的构成内容等，在MATLAB命令窗口中，可以通过指令来获取这些纯文本的帮助信息。
通常能够起到帮助作用，获取帮助信息的指令有help、lookfor、which、doc、get、type等。
实验二电力系统计算中常用的数值算法及电力网络的数学模型1、用因子表法求解电路方程中的节点电压；A= [1/10+1/15 -1/10 -1/15 -1-1/10 1/10+1/50+1/30 -1/50 0-1/15 -1/50 1/15+1/50+1/30 01 0 0 0 ];B=[0 0 0 10];[n,m]=size(A);for i=1:nA(i,i)=1/A(i,i);for j=i+1:nA(i,j)=A(i,j)*A(i,i);endfor k=i+1:nfor j=i+1:nA(k,j)=A(k,j)-A(k,i)*A(i,j);endendenddisp('矩阵A的因子表为：');disp(A)for i=1:nB(i)=B(i)*A(i,i);for j=i+1:nB(j)=B(j)-A(j,i)*B(i);endendfor i=n-1:-1:1for j=i+1:-1:2B(j-1)=B(j-1)-A(j-1,i+1)*B(i+1);endenddisp('在因子表的基础上求解线性方程组的解为：x=');disp(B)2、将图2.1中的电压源分别改为2V、4V、6V、8V时，用因子表进行求解，体会该算法带来的益处；①电压源改为2V②电压源改为4V③电压源改为6V④电压源改为8V3、上机验证课本例17.6－17.7的计算结果；①例17.6②例17.74、设计系统变更时修改导纳矩阵的函数（参考课本P.323），并调试通过，给出导纳矩阵修正函数的调用规则及其主要的实现代码。
% n=input('请输入节点数:n=');% nl=input('请输入支路数:nl=');% B=input('请输入由支路参数形成的矩阵:B=');% X=input('请输入由节点号及其对地阻抗形成的矩阵:X=');n=4;nl=4;B=[1 2 0.08+0.4i 0 1 0;2 3 0.1+0.4i 0 1 0;3 4 0.3i 0 1 0;1 3 0.12+0.5i 0 1 0];X=[1 0;2 0;3 0;4 0]Y=zeros(n);for i=1:nif X(i,2)~=0;p=X(i,1);Y(p,p)=1./X(i,2);endendfor i=1:nlif B(i,6)==0p=B(i,1); q=B(i,2);elsep=B(i,2); q=B(i,1);endY(p,q)=Y(p,q)-1./(B(i,3)*B(i,5));Y(q,p)=Y(p,q);Y(q,q)=Y(q,q)+1./(B(i,3)*B(i,5)^2)+B(i,4)./2;Y(p,p)=Y(p,p)+1./B(i,3)+B(i,4)./2;enddisp('导纳矩阵Y=:');disp(Y)B1=input('请输入追加支路参数形成的矩阵:B1=');% B1=[3 4 0.3i 0 1.2 0];% B1=[1 2 0.3i 0 1 0];% B1=[1 4 0.3i 0 1 0];% B1=[3 5 0.3i 0 1 0];z=B1(3);% if B1(6)==0% p=B1(1); q=B1(2);% else% p=B1(2); q=B1(1);% endp=B1(1); q=B1(2);if B1(6)==1B1(5)=1/(B1(5));endif B1(2)>nY(p,p)=Y(p,p)+1/z; Y(p,q)=-1/z; Y(q,p)=-1/z; Y(q,q)=(1/z)+B1(4);end思考题1.计算短路电流的目的是什么？答：选择电气设备,计算各种保护的定值,制定运行方案等2.在计算短路电流时，是否用到了节点阻抗阵中的所有元素？据此，可以如何简化短路计算程序？答：不需要，可以利用高斯消元法消去部分中间节点和联络节点，化简成等效网络再进行计算。
实验四潮流计算1.用牛顿拉夫逊法求解如下非线性方程组，注意保留迭代过程中的中间值。
clcclear allsyms x1 x2f1=x1.^2+x2.^2-5*x1+1;f2=x1.^2+x2.^2-3*x2-3;J(1,1)=diff(f1,x1);J(1,2)=diff(f1,x2);J(2,1)=diff(f2,x1);J(2,2)=diff(f2,x2);% Xn=input('X0=')Xn=[3;3];a=0.01;b=1;t=0; %迭代次数while (b)fn1=subs(f1,{x1 x2},{Xn(1,1) Xn(2,1)}); fn2=subs(f2,{x1 x2},{Xn(1,1) Xn(2,1)});fn(1,1)=fn1;fn(2,1)=fn2;Js(1,1)=subs(J(1,1),{x1 x2},{Xn(1,1) Xn(2,1)});Js(1,2)=subs(J(1,2),{x1 x2},{Xn(1,1) Xn(2,1)}); Js(2,1)=subs(J(2,1),{x1 x2},{Xn(1,1) Xn(2,1)});Js(2,2)=subs(J(2,2),{x1 x2},{Xn(1,1) Xn(2,1)}); Ja=-1*Js;Xns=double(Ja\fn)Xn=Xns+XnXnsa(1)=abs(Xns(1))Xnsa(2)=abs(Xns(2))Xnsm=max(Xnsa)t=t+1;if Xnsm<ab=0;else b=1;enddisp(Xn);enddisp('迭代次数=')disp(t);disp('最终结果为')disp(Xn);2.上机验证课本例19.1和例19.2的计算结果；3.将两种算法的程序计算程序进行合并clc实验五静态稳定和暂态稳定计算上机验证例20.1与例20.2的计算结果2、用MATLAB数值工具箱中的函数（解微分方程、求特征值）替代例程中的相应部分，重复验证过程。
clcS0=0.583+0.361i;V0=1;Xd=2.489;w0=100*pi;Tj=10;zn=0;Eq=sqrt((V0+imag(S0)*Xd./V0)^2+(real(S0)*Xd./V0)^2);dtj0=atan(real(S0)*Xd./(V0*(V0+imag(S0)*Xd./V0)));Psl=Eq*V0./Xd;C=w0;d=-1./Tj*Eq*cos(dtj0)*V0./Xd;E=-zn./Tj;Kp=(Psl-real(S0))./real(S0);A=[0 C;d E][D]=eig(A)%用此函数求矩阵的特征值%disp('该系统的静态稳定极限：Psl=');disp(Psl);disp('该系统的静态储备系数为：Kp=');disp(Kp);disp('该系统的线性微分方程的特征根为：');disp(D);if zn>=0if(D(1)==conj(D(2)))||((real(D(1))<0)&&(real(D(2))<0))disp('该系统是静态稳定的。
');elsedisp('该系统是静态不稳定的。
');endelseif zn<0disp('该系统是静态不稳定的。
');endA =。
