MATLAB潮流计算设计

第第PAGE1NUMPAGES6060各种修改.调整以外，还要注意输入和输出的方便性和灵活性，加强人机联系，以便使计算人员能及时监视计算过程并适当地控制计算的进行。用以检查系统各元件是否过负荷.各点电压是否满足要求，功率的分布和分配是否合PVPPQ0。一个独立的电力网中只设一个平衡点。3计算支路功率分布，PVLagrangePROCESSINGSYSTEMSNETWORKSnf1(x1, xn)f2(x1,

)fn(x1, xn)

x1x

x2x

xnx

f(x(0)x(0),x(0)x(0),x(0)x(0))

f(x(0)x(0),x(0)x(0),x(0)x(0))

f(x(0)x(0),x(0)x(0),x(0)x(0))

1nx(0)122

nf(x(0),x(0) (0))n

x(0)

x(0)

x(0)0n n

01020303 01020303n0102f(x(0),x(0) (0)) n0102

x(0)

x(0) 2

x(0)

0202

03 0301f(x(0),x(0) x(0)) 01

x(0)

x(0) n

x(0)n n

0

00f(x(0),x(0) x(0))00

x(0)

1f(x(0),x(0),x(0))

*x

0

0 2

(x(0),x(0) x(0)

nn

0

02n 2n11

12nx(1)x(0)x12n

k+1第第PAGE8NUMPAGES60f(x(k),x(k) x(k))

kk

kk kk

x(k)

1f(x(k),x(k) x(k)

*x(k

k

k 2

(x(k),x(k) x(k)

nn

k

kx(k

x(k x(k

x(k1)x(k)x(k得到修正

,……

,对各量进行修正 （i=1,2,……，n）迭代过程一直进行到满足收敛判据mx|f(x(k),x(k),，x(k))|} 2.1.mPVm+1，m+2，…，n-1PVn.

YijgijjBij，本文中节点电压以直角坐标形式表示，即Vieijfi第第PAGE9NUMPAGES60 Piei(GijejBijfj)fi(GijfjBijej)

Qi

fi(GijejBijfj)ei(GijfjBijej

12m

nnein

nnfj)fin

(Gif

Bije

))QiQisQiQisfi(GijejBifj)ei(GifjBijej

nneij

(Gije

fj)

nnj

(Gif

Bije

(im1,m2,,nV2V2V2V2(e2f2

.nVnenjfn修正方程WJVW1W1V1WJJJV 2

W3

*V3

MATLAB1Y计算有功功率的不平衡量ΔPi,从而求出Pi/Ui(i1,2,3...,n,isPAGE10NUMPAGES60J解修正方程式,求各节点的电压的变化量ei=ei+Δei(i=1,2,3...计算无功功率的不平衡量ΔQi,从而求出Qi/Ui解修正方程式,求各节点的电压大小的变化量

(i=1,2,3...,n,i≠s)求各节点的电压大小的新值UiUi

(i=1,2,3...,n,i≠s) in SSUS

UiPi SijUiIijUi[Uiyi0(UiUj)yij]Pij SjiUjIjiUj[Uiyj0(UjUi)yji]Pji SijSijS

Pijn节点导纳矩阵的非对角元Yiji,j

NVhV1

i,j特殊说明所有的参数都用标幺值)为Zij,则线路导纳为YijZijN纳为Y

1N/(Zij*N*N),Y0

Y0ijJY2.2

其中计算△P;△Q;△VwhileP=P(k)-e(k)*(dotG(k,:)*e-dotB(k,:)*f)-f(k)dotG(k,:)*f+dotQ=Q(k)-f(k)*(dotG(k,:)*e-dotB(k,:)*f)+e(k)dotG(k,:)*f+dotb(k,:)*e);W(k)=[PQ]';k=k+1;whileL=P(k)-e(k)*(dotG(k,:)*e-dotB(k,:)*f)-f(k)dotG(k,:)*f+dotb(k,:)*e);W(k)=[L(k)v(k)]';V(k)=[Li=1j=1whilei=^m+1ifj=nelseifi=jJ(i,j)=[-(dotG(i,:)*e-dotB(i,:)*f)-G(i,i)*e(i)-B(i,i)*f(i)-(dotB(i,:)*e)+B(i,i)*e(i)-G(i,i)*f(i);dotG(i,:)*f+dotB(i,:)*e+B(i,i)*e(i)-G(i,i)*f(i)-(dotG(i,:)*e-dotB(i,:)*f)+G(i,i)*e(i)+B(i,i)*f(i)];J(i,j)=[-(G(i,j)*e(i)+B(i,j)f(i))B(i,j)*e(i)-G(i,j)*f(i);B(i,j)*e(i)-G(i,j)*f(i)whilei=^nifj=nelseifi=jJ(i,j)=[-(dotG(i,:)*e-dotB(i,:)*f)-G(i,i)*e(i)-B(i,i)*f(i)-(dotB(i,:)*e)+B(i,i)*e(i)-G(i,i)*f(i);-2*e(i)-2*f(i)];J(i,j)=[-(G(i,j)*e(i)+B(i,j)f(i))B(i,j)*e(i)-G(i,j)*f(i);00];MATLABM文件来解决；MATLABMATLAB提供MATLABIEEE-6BUS标准试验系统的潮流计算进行仿真,具体计算中系统节点分三PVQδ未PV;n1.12.1RX112214316425535643746|V(0)||V(0)||V(0)||V(0)| |V(0)|1.100,|V(0)| (0)(0)(0)(0)(0)(0) P150.0,Q15.0,P230.0,Q218.0,P355.0,Q3P450.0,Q45.0,P5 2.3IEE－62.20+000+3.2520i0.5765-00-0.5765+0000.4449-0+6.8353i-0.4449+0+00+6.8353i1.1124-0-0.5583+0-0.5765+1.3085i-0.4449+01.0214-0-0.4339+00-0.5583+00.9922-

nnein

nnfj)fin

(Gif

Bije

))QiQisQiQisfi(GijejBifj)ei(GifjBijej计算Pi

nneij

(Gije

fj)

nnj

(Gif

Bije

(im1,m2,,nV2V2V2V2(e2f2 计算Pi

W1W1V1WJJJV 2

W3

*V3

解之可得V1V2由V1e1V2f1Vn2en1Vn1e1e1e1,f1

f1f1en1en1en1,fn1

fn1fn1whilewhilewhileholdonholdonYlabel('P(*)Q(+)')whilei~=m+1ify==nelseifwhilei~=nify==nelseifwhileh~=na节点1中P(*),Q(+)随迭代次数k的变 b节点2中P(*),Q(+)随迭代次数k的变c节点3中P(*),Q(+)随迭代次数k的变 d节点4中P(*),Q(+)随迭代次数k的变e5P(*),Vk2.4既反映在建模方法的技术层面，同样更高地反映到对建模过程和模型配置管理的要求。（Interoperability（FDMS仿真模型互操作性的实现需要遵循一定标准，HLA标准定义了对象模型模板（OMT（FOMCORBA、DCOM的自动化程度可以大大提高。仿真模型互操作性是近十年逐渐发展起来研究课题，VV&AVV&AVV&AVV&A调管理问题是其发展应用中必须解决的一些问题。仿真系统的规模越来越大，复杂VV&AVV&AVV&AVV&AVV&AVV&AVV&AVV&AVV&AVV&AVV&AVV&AVV&AVV&AVV&AVV&AVV&AVV&AVV&AVV&AVV&A80SIMNET布交互仿真技术的快速发展一方面是得益于分布式计算技术和通信网络技术快速发90IEEE1278真应用交互协议数据单元（PDU）和编码、网络要求、仿真工程管理、校核验证与验（IEEE1278同时发展的聚合级仿真（ALSP定义了聚合级仿真应用之间互联进行大规模作战仿真的所应遵循的体系1995转变的计划。HLAHLA信息接口的规定，IEEE20009HLA[1]IEEE1516其成为工业标准。分布交互仿真标准化工作在新世纪的开始将进一步深入地开展，作，同时大大提高了仿真系统开放性、扩展性和可管理性。一般地，仿真系统体系HLARTI；仿真模型层包括完成仿真应用所需的各种仿真模型设计及其实现，仿真模型的OMT（FOM）和仿真对象模型（SOM；（UNICORE、MOLNimrod/G、EcoGridNinf、Bricks国也有中科院、清华大学等多个单位在进行网格方面研究。基于网格的仿真可以充DARPA(DefenseAdvancedResearchProjects1,386100,298（4）WebVRMLXML设计开发中得到广泛的应用，人们可以利用具有友好的人机界面仿真应用来进行工技术学习的难度大大降低，将有力地促进科学技术水平提高。全精确地代表真实对象,存在一个可信性(Credibility1962,BiggsandCawthorne(1962(Bloodhound13,16(SCS)型可信性技术委员会(TCMC)这是一个重要的里程碑。七十年代中期以来,模型可信性成为许多重要学术会议的主1981(SCSC)和冬季仿真会议(WSC)关于模型可信性的专题讨论;1984(IFAC)于模型可信性的专题讨论;过去几年的“高性能计算机会议”中都有关于仿真模型校核(Verification)和验证(Validation)的讨论;美国军事运筹学会(MORS)自1989期以来召开了多次有关模型校核、验证和确认(Accreditation)的小型讨论会年11VVA;对基于知识的系统的VVA,9193年召开过两次重要会议(EUROVAV)。九十年代以来,许多政府、民间部门和学术机构都成立了相应的组织,以VVA(DoD)1991真办公室(DMSO)VVADMSO1993VVABMDOSDIO),都先后制VVA(NIST1992VVAVVA;IEEE前些年发表了一个软件VV(ANSI)准指导委员会采用;计算机仿真学会(SCS)1993比较粗浅和分散,尚属于起步阶段,甚至人们对有关仿真可信性的概念和一般方法论TCMC术迅速发展和需求日益迫切的形势极不相称,需大力加强。许多学者和工程技术人员在系统仿真可信性研究方面做了大量工作,BalciandSargent(1984a)汇集了19841967年,FishmanandKiviat(1967)就指出了评估仿真模型时牵涉到两个方面:模型校核和模型验证,其中,“校核”指检查计算机程序的逻辑和代码,而“验证”指证实仿真模型和真实对象有相同的行为。NaylorandFinger(1967)认识到难以建立一个广泛适用判断模型和真实对象之间的一致性。基于经济模型的验证,他们提出了三种历史方法论,即理性主义(Rationalism)(Empiricism)Mckenny(1967)ShrankandHolt(1967)NaylorandFinger时指出,模型是否有效与仿真结果的用途有关,而且认为仿真验证要回答的不是模型表达真实对象的真假,而是仿真结果中的误差是否大到淹没其有用性的程度。不管他进一步认为验证的目的是比较真实对象和仿真模型所具有的一些特性,而不一定是针对仿真的建模机理。Mihram(1972)将模型开发过程分为五个步骤,即系统分析、系统合成、模型校核、模型验证和模型分析,实质上已经把仿真可信性评估作为仿真工作Garratt(1974)Schellenberger(1974)模型有效性的准则,包括技术有效性、运行有效性和动态有效性。表面看来这是一个Tytula(1978)NaylorandfingerVanHorn对应、十分接近。KheirandHolmes(1978)认为应尽可能将系统分成多个子系统,这样证,最后验证全系统模型。而且,针对导弹系统模型的验证,提出了一个多层次、多步,TCMC1979准术语,为人们所普遍接受和使用。他们将仿真环境分为三个基本要素,即:真实对象(Reality)、概念(Conceptual)模型和计算(Computerized)模型,联系这些要素的是仿真的三个过程,即分析、编程和计算机仿真,而评估这些过程的可信性的步骤为:模型认可(Qualification)义。Schruben(1980)认为为一个非理论的过程建立一个严格有效的模型是不可能的,可信性反映在人们对基于仿真结果进行决策的信心上,而且强调仿真建模者必须建立清晰的文档,供有关人员确认仿真模型时参考。Grider(1981)提出了半实物仿真的校核和验证问题,指出仿真模型校核不仅仅是要校核计算机仿真程序,还应包括仿真系统中的硬件设备。Kheir(1981)TCMCSargent(1984,1988)实际含义及评估方法由问题属性、系统属性和模型属性决定,并给出了这些属性的内的数据是充足和正确的。kheir(1985)实对象和仿真要求,因此用户应参与可信性评估过程,而且他们的意见是主要的。中,不能等仿真工程结束后才进行。另外,他还强调了用户参与的重要性。“如何表示和度量仿真的可信性?用什么方法对仿真模型进行VVA?”Hermann(1967)指出可以从不同角度验证仿真模型和真实对象之间的一致性,如二者的典型事件类。一类是主观判断比较。如Hermann(1967)认为可根据直觉对仿真模型进行外观Schruben(1980)Turing和仿真输出进行比较,如果他们不能区别两类输出,则认为仿真结果是有效的;否则,需要研究他们是如何能区分的,这可以为修正仿真模型提供有价值的反馈信息,或者是确认使他们进行区分的因素对仿真用途而言没有影响,从而认为模型对该用途而言VanHorn(1971)的结果可能在模型中使用的参数值、假设和结构的某个范围内都是有效的,通过灵敏度(Sensitivity)分析可确定这些模型特征的一个范围,使得在其中由仿真获得的结果都是有效的。如果没有理由怀疑模型的特征超出这个范围,则对模型有效性的信心大大增强。Miller(1974)提出并演示的灵敏度分析方法稍微不同,他计算了与模型的一些参考条件相关的一组灵敏度系数,并估计模型输入参数的概率分布,将此误差分布转换系统的输出,包括静态输出(随机变量)和动态输出(随机过程或时间序列)。X2检验、KolmogorovSmirnovCramer-VonMisesGarratt(1974)(MANOVA)Teorey(1975)TZucker(1981X2检验方法验证复杂系统的仿真模型,即随机向量的比较,认X2检验,也可以先通过所谓标准化消除变量之间的相关性,然后使用X2检验BalciandSargent(1982a,1982bWright(1972)方法,如方差分析法。SchatzoffandTillman(1975)就在检验不同的调度算法时,用假设检验方法验证从一个真实计算机系统和从该系统的仿真模型得到的结果之间的一致性。KheirandHolmes(1978Theil(TIC)作为性能指标比较导弹系统仿真试验和飞行试验的数据,后来,RowlandandHolmes(1978)提出用适合于稀疏随机数TICGriner(1978)Baird(1980)Bayes(RF)FishmanandKiviat(1967)了检验两个频谱等价性的步骤,用于比较一个排队模型和其代表的真实过程。Gallantetal.(1974)描述了一个方法,结合使用Schuster重复试验的结果。MontgomeryandConard(1980)使用频谱方法比较导弹系统仿真和飞行试验数据,周宪民(1988)和杨军、刘藻珍(1994)采用这种方法验证导弹系统仿真模abcabc可信性与仿真模型的用途有关,又可理解为可用性,关心的是偏差是否大到淹灭其有用性的程度。应将VVA只有在少数情况下,才有可能对可信性的某些方面作客观定量表示和评估,而大多VVAVVA。努力提高VVA以上诸方面都依赖于研究的具体问题,难以建立一个统一的规范,每一个仿真工程都是一个挑战。由以上方法性研究工作可以得出如下结论:电力系统原型是一个由各种电气设备和控制装置构成的整体，其数学模型包括描述设备动态特征的微分方程和描述设备之间电气联系的代数方程;设备电气联接关系或离散的时变参数;工干预和自动控制;一般可将电力系统数学模型用高维非线性、非自治的分时段微分一代数方程组(ntiilgbaiqutinos)来描述。电力系统安全运20世纪50目的由物理原型抽象出数学模型。建模过程应考虑以下几个问题:①依据研究问题的侧重点不同，对原型进行不同程度的简化;②电力系统由若干子系统组成，应保证整dqo坐标的接口;③在保证精度前提下，数学模型应有利于数值求解，例如忽略发电机定子暂态对电力系统机电暂态影响很小，解是采用一定的数值方法和软件技术来设计仿真程序。仿真程序的编写应注意:①选择数值稳定性高、误差可以控制、计算量少和尽可能节省内存空间的算法;②充分利用计算机硬件和软件的新技术，提高仿真系统的性价比;③注意提高仿真程序的可靠性、可扩展性和调试维修方便性④仿真程序应具有直观友好的人机界面。数字仿真的独特优势也是促使其快速发展的重要因素:①不受系统规模和复杂性的限制;组建几十台机、上百条支路的动模实验室己经非常困难，而优秀的电力系统数字远强于前者。②保证被研究系统的安全性;许多试验，例如稳定性破坏试验，直接在济性;动模试验需要专职人员进行大量准备工作，而且存在设备折旧费用，数字仿真只需少数计算人员参与，过程短，花费少。④可用于电力系统规划和前瞻性研究;在电力系统仿真就是通过建立适当的数学模型来模拟实际电路的一种研究方法,随着电力系统的不断扩大和网络化,实际电力网络拓扑系统变得越来越复杂,而这时候高效的模拟仿真也变得越来越重要,随着计算机技术的不断发展,电力系统仿真软件已成为电力系统工作者进行电力系统规划、保护、调度及故障研究的重要工具。FACTSNETOMAC国电力公司(PTI)PSS/E(PowerSystemSimulatorforEngineering),Math-worksMATLAB(5.2PowerSystemBlockset(PSB)工具箱,以及中国电力科学院给出的基于DOSPSASP。上述几种软件各有各的特点:PSS/E多种输出设备相连,输入输出可根据用户要求进行设计,它要求使用者有一定的编程EMTPPSASPPSASP（1）目前，MATLABMATLABSimulinkSimulinkSimulinkMATLAB（2）MTALABPSP(电力系统软件包)1）PSP简介PSP库内模型分为电源组、元件组、电力电子元件(PowerGUD1998年，MathworkSMATLABVersion5.2PowerSysterBloek(PSB)6.5Sin1POwersystems(SPS)，仿真功能得到了很大的提高。SPS提供了一个现代的分析工Siulink环境下通过鼠标选取相应的模型拖放到设计窗口，然后进行连接并设置各模块的参数。以变压器为例:首先将变压器拖放到设计窗可以通过PoI进行系统的稳态潮流计算和初始设置，这些可以通过设置“odFlondMehinnitiliztionPV节点或平衡节点(或电压相角用电力系统工具箱提供的powerZsys()可以提取从给定电源到输出端子的状态方程模Bode图进行电力系统的频域分析，寻找系统的振荡模式。PSB是一个图形工具，利用它可以在仿真环境中对电力系统进行设计与仿真Simulink(PowerSystemsBlockset)功能非常强Connectors(ElectricalElements(3.1powerlib（3）及库中元件很多电阻(电容、电感)的仿真电阻元件在Elements库中,双击Elements10/11/2016Library:Powerlib/Elements(SeriesorParallelRLCBranch)。仿真时用鼠标把RLC串联支路拖放至标题为untitledRLCRR、L、CLC0 (3-PhaseProgrammableVoltageSource),受控电压源(ControlledVoltageSource)及电源参数设定窗口,需设定的参数有:电压（peakamplitude，相数(phase),频率(frequency)，瞬时时间(sampletime)。Powergui改变初始状态,开始从起始条件模拟.状态参数的文件名就是建立电容电压和电感UcLi.显示负荷的波动情况,并能设置包括三相电机在内的三相电力系统的初始条件,模拟过程在稳定状态下进行.这一选项适于含以下几种电机类型的电路:简化的同步发电机、同步电机和异步电机(鼠笼)模块.自动生成系统的状态-空间模型(SS)(若安装了控制系统工具箱),自动打开LTI该报表保存在扩展名为.rep用Powergui当把母线设置为相对相(phase-to-phase)Powergui模块的Steady-StateVoltagesandCurrents及功率损耗.这样在仿真系统动态运行过程的基础上,就能快捷的完成电力系统的简PowerguiLoadFlowandMachineInitializationUpdateFlow系统稳态计算中,负荷以恒定的功率表示,而在某些计算中有用固定阻抗或负荷电压静态特性来表示负荷的.SimPowerSystems点类型的不同,已知参数的情况也就不同.SimPowerSystems适用任何情况,应根据不同的情况选择正确的负荷仿真方式.PQ在仿真模拟图中,母线B1—B4PQPV然而在母线为PVRLCPVSynchronousMachinepu模块可以做为同步发电机也可以做为模拟电动机系统.于是利用这一特性,通过对参3.53.5与三相RLCPowerguiUpdateLoadFlow仿真系统中负荷模块的特殊处理.然而上述电动机模型在仿真系统中依然存在一些问题,该负载只能与非电流源元件(如线路)相连.这是因为在SimPowerSystems统中,电动机模块是由等效电流源进行模拟的.而系统中还有其它使用电流源模拟的元件,比如变压器.那么当需要模拟负载通过降压变压器从系统吸收功率时,将二者直接连接.仿真系统会提示出错信息:两个电流源不能串联运行.也就是说在SimPowerSystems元件串联运行时,可以在电动机的出口处并上一个小负载,改变纯粹的串联结构,从而建立一个新的负载替代固定阻抗的负载模块,如图3.6作了.RLC0.002%),所以对系统的潮流分析不会有太大的影响.用以上方法搭建的电力系统模型,可以使负载输出恒定在所设定的值,在这时观察系统的潮流变化发现,母线上的电压降落在±5%以内,满足系统对电压损耗的要求.这就使得潮流分析达到更理想的值.3.7IEE6MATLABpowergui3.83.8powerguimachine5/synchronousmachinepustandard,bustypepv6/synchronousmachinepustandard,在bustype3.93.9updateloadflow563.10

3.10windowsbutton；单选按钮(radiobutton)；检查框(cheekbos)；列表框(listbox)；下拉菜单(popupmenu)；滑块(slider)；编辑框(editbox)；静态文字(statictext)；框架(frame)matlabfigure命令中完成的。例如figwin=figure('position',[100100300],'name','uicontrol:pushbutton','numbertitle','off');语句生成如图4.14.1figure 20],'string','tuichu','CallBack',['close(figwin)'])其中'position'是生成窗'string'callback‘style’是其控件属性。语句是在以上窗口中添加按钮控件生成窗口如下；

4.2figurefigwin=figure('position',[100100600330],'name','uicontrol:AxesHandle=axes('position',[0.40.150.50.7],'box','on');axis([010010])axisequalholdonaxisequalholdonaxisequalholdonaxisequalholdonaxisequalholdonaxisequalholdplot([12],[8plot([34],[8plot([89],[8plot([12],[2plot([67],[2plot([89],[2plot([1.21.2],[2plot([1.83.2],[7.7plot([3.23.2],[7.7plot([3.83.8],[7.7plot([3.85.5],[7.7plot([8.28.2],[7.7plot([8.88.8],[2plot([1.81.8],[2.3plot([4.56.2],[2.3plot([6.26.2],[2.3plot([6.86.8],[2.3plot([6.88.2],[2.3plot([8.28.2],[2.3plot([6.58.2],[7.7plot([1.81.8],[7.7plot([1.83.5],[2.3plot([1.51.5],[2plot([1.53.5],[2.3plot([3.53.5],[8holdplot([1.51.5],[2plot([8.58.5],[2plot([6.56.5],[2plot([1.51.5],[8plot([3.53.5],[8plot([8.58.5],[8axis([0100framatbox=uicontrol(figwin,'style','frame','position',[15205305]);framatbox1=uicontrol(figwin,'style','frame','position',[2105250framatbox2=uicontrol(figwin,'style','frame','position',[4905110text=uicontrol(figwin,'style','text','position',[55313515text=uicontrol(figwin,'style','text','position',[55311515text=uicontrol(figwin,'style','text','position',[5539515text=uicontrol(figwin,'style','text','position',[5537515text=uicontrol(figwin,'style','text','position',[5535515text=uicontrol(figwin,'style','text','position',[5533515text=uicontrol(figwin,'style','text','position',[2805540text=uicontrol(figwin,'style','text','position',[2809010text=uicontrol(figwin,'style','text','position',[28023010text=uicontrol(figwin,'style','text','position',[4009010text=uicontrol(figwin,'style','text','position',[4509010text=uicontrol(figwin,'style','text','position',[33523010text=uicontrol(figwin,'style','text','position',[45023010text=uicontrol(figwin,'style','text','position',[41020025text=uicontrol(figwin,'style','text','position',[3458535text=uicontrol(figwin,'style','text','position',[3805540text=uicontrol(figwin,'style','text','position',[42025040radio1=uicontrol(figwin,'style','pushbutton','position',[190290100.8],[8.410],'string','','callback',['line([35],[2.82.8]);line([4.85],[2.95],[2.710],'string','','callback',['line([2.73.7],[4.36.9]);line([3.66.9]);line([3.83.7],[6.710],'string','','callback',['line([11],[46]);line([0.91],[5.86]);line([1.110],'string','','callback',['line([78],[2.52.5]);line([7.88],[2.68],[2.410],'string','','callback',['line([8.58.5],[46]);line([8.58.6],[48.4],[410],'string','','callback',['line([57],[8.28.2]);line([6.87],[8.17],[8.310],'string','','callback',['line([23],[7.57.5]);line([22.2],[7.52.2],[7.510],'string','','callback',['line([35],[1.31.3]);line([33.2],[1.33.2],[1.310],'string','','callback',['line([1.41.4],[46]);line([1.41.3],[41.5],[410],'string','','callback',['line([78],[1.41.4]);line([77.2],[1.47.2],[1.410],'string','','callback',['line([9.19.1],[46]);line([9.09.1],[5.89.1],[5.810],'string','','callback',['line([57],[7.37.3]);line([55.2],[7.35.2],[7.310],'string','','callback',['line([23],[8.48.4]);line([2.83],[8.33],[8.5push1=uicontrol(figwin,'style','pushbutton','position',[300107020],'string','清除push2=uicontrol(figwin,'style','pushbutton','position',[380107020],'string','退出push3=uicontrol(figwin,'style','pushbutton','position',[220107020],'string','计算text=uicontrol(figwin,'style','text','position',[7531045text=uicontrol(figwin,'style','text','position',[14031045text=uicontrol(figwin,'style','text','position',[529055text=uicontrol(figwin,'style','text','position',[12529015text=uicontrol(figwin,'style','text','position',[12527015text=uicontrol(figwin,'style','text','position',[12525015text=uicontrol(figwin,'style','text','position',[12523015text=uicontrol(figwin,'style','text','position',[12521015text=uicontrol(figwin,'style','text','position',[12519015text=uicontrol(figwin,'style','text','position',[12517015text=uicontrol(figwin,'style','text','position',[12515015text=uicontrol(figwin,'style','text','position',[12513015text=uicontrol(figwin,'style','text','position',[12511015text=uicontrol(figwin,'style','text','position',[1259015text=uicontrol(figwin,'style','text','position',[1257015text=uicontrol(figwin,'style','text','position',[1255015text=uicontrol(figwin,'style','text','position',[1253015text=uicontrol(figwin,'style','text','position',[1251015e1=uicontrol(figwin,'style','edit','position',[52028030e2=uicontrol(figwin,'style','edit','position',[52026030e3=uicontrol(figwin,'style','edit','position',[52024030e4=uicontrol(figwin,'style','edit','position',[52022030e5=uicontrol(figwin,'style','edit','position',[52020030e6=uicontrol(figwin,'style','edit','position',[52018030e7=uicontrol(figwin,'style','edit','position',[55028030e8=uicontrol(figwin,'style','edit','position',[55026030e9=uicontrol(figwin,'style','edit','position',[55024030e10=uicontrol(figwin,'style','edit','position',[55022030e11=uicontrol(figwin,'style','edit','position',[55020030e12=uicontrol(figwin,'style','edit','position',[55018030e13=uicontrol(figwin,'style','edit','position',[52016060text=uicontrol(figwin,'style','text','position',[52531050text=uicontrol(figwin,'style','text','position',[50028520text=uicontrol(figwin,'style','text','position',[50026520text=uicontrol(figwin,'style','text','position',[50024520text=uicontrol(figwin,'style','text','position',[50022520text=uicontrol(figwin,'style','text','position',[49020525text=uicontrol(figwin,'style','text','position',[50018520text=uicontrol(figwin,'style','text','position',[50016520text=uicontrol(figwin,'style','text','position',[50013520text=uicontrol(figwin,'style','text','position',[50011520text=uicontrol(figwin,'style','text','position',[5009520text=uicontrol(figwin,'style','text','position',[5007520text=uicontrol(figwin,'style','text','position',[5005520text=uicontrol(figwin,'style','text','position',[5003520text=uicontrol(figwin,'style','text','position',[5001520text=uicontrol(figwin,'style','text','position',[5201520text=uicontrol(figwin,'style','text','position',[7511045text=uicontrol(figwin,'style','text','position',[7529045text=uicontrol(figwin,'style','text','position',[14029045text=uicontrol(figwin,'style','text','position',[7527045text=uicontrol(figwin,'style','text','position',[14027045text=uicontrol(figwin,'style','text','position',[7525045text=uicontrol(figwin,'style','text','position',[14025045text=uicontrol(figwin,'style','text','position',[7523045text=uicontrol(figwin,'style','text','position',[14023045text=uicontrol(figwin,'style','text','position',[7521045text=uicontrol(figwin,'style','text','position',[14021045text=uicontrol(figwin,'style','text','position',[7519045text=uicontrol(figwin,'style','text','position',[14019045text=uicontrol(figwin,'style','text','pos