2022年matlab电力系统潮流计算.pdf

《2022年matlab电力系统潮流计算.pdf》由会员分享，可在线阅读，更多相关《2022年matlab电力系统潮流计算.pdf（22页珍藏版）》请在淘文阁 - 分享文档赚钱的网站上搜索。

1、matlab 电力系统潮流计算华中科技大学信 息 工 程 学 院课 程 设 计 报 告 书题目: 电力系统潮流计算专业: 电气工程及其自动化班级: 学号: 学生姓名 : 指导教师 : 2015 年 11 月 10 日学 号学生姓名专业 (班级 ) 电气工程及其自动化 () 设计题目电力系统潮流计算设计技术参数1.变压器的阻抗,励磁损耗。2.线路的阻抗 ,无功损耗。3.变压器及线路的功率损耗4.变压器及线路的电压降落5.各节点的功率与电压设计要求1、计算过程包括理论计算与MATLAB仿真两大部分。2、报告里面的理论分析(手工计算 )过程必须要有详细推导过程,再把仿真过程放到相应的部分中去。精品资

2、料 - - - 欢迎下载 - - - - - - - - - - - 欢迎下载 名师归纳 - - - - - - - - - -第 1 页，共 22 页 - - - - - - - - - - matlab 电力系统潮流计算工作量共 24 页;4000 字工作计划第 15 周:对系统模型分析与资料收集,找到计算的思路,总结计算方法 ; 第 16 周:对整个模型进行详细计算,并完成 word 文档。参考资料1 何仰赞 , 温增银电力系统分析(第三版 )M 、 华中科技大学 ,2002 2 刘同娟、MATLAB在电路分析中的应用、电气电子教学学报、2002 3 王守相 ,刘玉田电力系统潮流计算研究

3、现状-山东电力技术1996 年 05 期指导教师签字2015年 11 月 12 日精品资料 - - - 欢迎下载 - - - - - - - - - - - 欢迎下载 名师归纳 - - - - - - - - - -第 2 页，共 22 页 - - - - - - - - - - matlab 电力系统潮流计算信息工程学院课程设计成绩评定表学生姓名 : 学号: 专业(班级): 班课程设计题目 : 电力系统潮流计算指导教师评语 : 成绩: 指导教师 : 年月日精品资料 - - - 欢迎下载 - - - - - - - - - - - 欢迎下载 名师归纳 - - - - - - - - - -第

4、3 页，共 22 页 - - - - - - - - - - matlab 电力系统潮流计算摘要电力系统稳态分析包括潮流计算与静态安全分析。本文主要运用的事潮流计算, 潮流计算就是电力网络设计与运行中最基本的运算, 对电力网络的各种设计方案及各种运行方式进行潮流计算, 可以得到各种电网各节点的电压, 并求得网络的潮流及网络中的各元件的电力损耗, 进而求得电能损耗。本位就就是运用潮流计算具体分析 , 并有 MATLAB 仿真。关键词 : 电力系统潮流计算MA TLAB 仿真精品资料 - - - 欢迎下载 - - - - - - - - - - - 欢迎下载 名师归纳 - - - - - - -

5、- - -第 4 页，共 22 页 - - - - - - - - - - matlab 电力系统潮流计算Abstract Electric power system steady flow calculation and analysis of the static safety analysis、 This paper, by means of the calculation, flow calculation is the trend of the power network design and operation of the most basic operations of ele

6、ctric power network, various design scheme and the operation ways to tide computation, can get all kinds of each node of the power grid voltage and seek the trend of the network and the network of the components of the power loss, and getting electric power 、 The standard is to use the power flow ca

7、lculation and analysis, the specific have MATLAB simulation、Key words: Power system; Flow calculation; MATLAB simulation精品资料 - - - 欢迎下载 - - - - - - - - - - - 欢迎下载 名师归纳 - - - - - - - - - -第 5 页，共 22 页 - - - - - - - - - - matlab 电力系统潮流计算目录1 任务提出与方案论证 . 2 2 总体设计 . 2 2、1 潮流计算等值电路 . 3 2、2 建立电力系统模型 . 3 2、

8、3 模型的调试与运行 . 3 3 详细设计 . 错误!未定义书签。3、1 计算前提 . 3 3、2 手工计算 . 6 4 设计图及源程序 . 9 4、1MATLAB 仿真 . 9 4、2 潮流计算源程序 . 10 5 总结 . 16 参考文献 . 17 精品资料 - - - 欢迎下载 - - - - - - - - - - - 欢迎下载 名师归纳 - - - - - - - - - -第 6 页，共 22 页 - - - - - - - - - - matlab 电力系统潮流计算1 任务提出与方案论证潮流计算就是在给定电力系统网络结构、参数与决定系统运行状态的边界条件的情况下确定系统稳态运行状

9、态的一种基本方法, 就是电力系统规划与运营中不可缺少的一个重要组成部分。 可以说, 它就是电力系统分析中最基本、最重要的计算 , 就是系统安全、经济分析与实时控制与调度的基础。常规潮流计算的任务就是根据给定的运行条件与网路结构确定整个系统的运行状态, 如各母线上的电压 ( 幅值及相角 ) 、网络中的功率分布以及功率损耗等。潮流计算的结果就是电力系统稳定计算与故障分析的基础。在电力系统运行方式与规划方案的研究中, 都需要进行潮流计算以比较运行方式或规划供电方案的可行性、可靠性与经济性。同时, 为了实时监控电力系统的运行状态, 也需要进行大量而快速的潮流计算。因此 , 潮流计算就是电力系统中应用最

10、广泛、最基本与最重要的一种电气运算。在系统规划设计与安排系统的运行方式时, 采用离线潮流计算 ; 在电力系统运行状态的实时监控中 , 则采用在线潮流计算。就是电力系统研究人员长期研究的一个课题。它既就是对电力系统规划设计与运行方式的合理性、可靠性及经济性进行定量分析的依据, 又就是电力系统静态与暂态稳定计算的基础。潮流计算经历了一个由手工到应用数字电子计算机的发展过程, 现在的潮流算法都以计算机的应用为前提用计算机进行潮流计算主要步骤在于编制计算机程序,这就是一项非常复杂的工作。对系统进行潮流分析,本文利用MATLAB 中的 Simpower Systems工具箱设计电力系统 ,在 simul

11、ink 环境下 ,不仅可以仿真系统的动态过程,还可以对系统进行稳态潮流分析。2 总体设计SimpowerSystems使用 Simulink 环境,可以将该系统中的发电机、变压器,线路等模型联结起来 ,形成电力系统仿真模拟图。在加人测量模块,并对各元件的参数进行设置后,用measurement与 sink 中的仪器可以观察各元件的电压、电流、功率的大小。精品资料 - - - 欢迎下载 - - - - - - - - - - - 欢迎下载 名师归纳 - - - - - - - - - -第 7 页，共 22 页 - - - - - - - - - - matlab 电力系统潮流计算2、1 潮流计

12、算等值电路GGGGGGGGGGGGGGMW20MW1548 .0cos073. 016. 0136. 002NdxxxkV10kV110MW2580km32kmMWVAdYN16211,kV10MWVAdYN1611,MW30MW12070kmMW63485. 0cos06.0161.0134. 002Ndxxx8. 0cos075.0161. 0136.002NdxxxMW123MW50185. 0cos054. 0154. 0128. 002NdxxxMW2528. 0cos0591. 0157. 0128. 002NdxxxMWVA63MW35MVAdYYN10211,kV110kV11

13、0kV10kV3530kmkV35kV10MW15MWVAdYN10211,MWVAdYN63411,kV10MW105.10%5. 0%737.1500ssVIkWpkWpMWVAdYN20211,5.10%5. 0%896.1800ssVIkWpkWp5.10%5. 0%737.1500ssVIkWpkWpMW255.10%35. 0%1214400ssVIkWpkWp5.10%35. 0%1214400ssVIkWpkWp5.10%55. 0%501100ssVIkWpkWpMW80MW55.17%5. 6%5 .10%55. 0%632.13)31()32()21(00ssssVVVI

14、kWpkWp130km 2 、2 建立电力系统模型在 Simulink 中按照电力系统原型选择元件进行建模。所建立的模型与建立的方法在详细设计中详述。在电力系统模型的建立工程中主要涉及到的就是:元器件的选择及其参数的设置;发电机选型 ;变压器选择 ;线路的选择 ;负荷模型的选择;母线选择。 2 、3 模型的调试与运行建立系统模型 ,并设置好参数以后,就可以在Simulink 环境下进行仿真运行。运行的具体结果与分析也在详细设计中详述。3 详细设计3、1 计算前提首先就是发电机的参数计算,先对 5 个发电厂简化为5 台发电机来计算。发电机 G1: MVarQMWP45)8. 0tan(arcco

15、s606015411发电机 G2: MVarQMWP156)85. 0tan(arccos25225263422发电机 G3: 精品资料 - - - 欢迎下载 - - - - - - - - - - - 欢迎下载 名师归纳 - - - - - - - - - -第 8 页，共 22 页 - - - - - - - - - - matlab 电力系统潮流计算333 123636tan(arccos0.8)27PMWQMVar发电机 G4: 441 505050tan(arccos0.85)31PMWQMVar发电机 G5: 552255050tan(arccos0.8)37.5PMWQMVar其

16、次就是变电站的参数计算,我们还就是对7 个变电站简化为7 台变压器来计算。变压器 T1: MVAjSIjpSSVVXSVpRNNNsTNNsT)0800.00157. 0(100%406.791010161105.1010%450. 310)1016(11073100001322132323221变压器 T2:(双并联 ) MVAjSIjpSSVVXSVpRNNNsTNNsT)2000.00372.0()100%(27625.311010201105.102110%21346. 110)1020(110892110210002322232323222变压器 T3:(四并联 ) 22333232

17、2233003011121 11010100.09244(63 10 )%1110.5 11010105.0424463 10%4 ()(0.17600.8820)100sNTNsNTNNpVRSVVXSISpjSjMVA变压器 T4:(双并联 ) MVAjSSXXRRTTTT)1600.00314.0(27030.39217250.12101041414变压器 T5: 精品资料 - - - 欢迎下载 - - - - - - - - - - - 欢迎下载 名师归纳 - - - - - - - - - -第 9 页，共 22 页 - - - - - - - - - - matlab 电力系统潮流

18、计算MVAjSSXXRRTTTT)2205.00440.0(41168.2043680.0403053535变压器 T6:(两个三绕组变压器并联) 386.010)1010(3563213232362616TTTRRR75.6%21%25.0%21%75.10%21%)21()32()31(3)31()32()21(2)32()31()21(1ssssssssssssVVVVVVVVVVVV134.410%21153.010%21584.610%21233622262116NNsTNNsTNNsTSVVXSVVXSVVXMVAjIjPS)1100.00264.0()10100%(200606变

19、压器 T7:(双并联 ) MVAjSIjpSSVVXSVpRNNNsTNNsT)1100.00220.0()100%(2431.6101010355.102110%21306.010)1010(35502110210007322732323227再次就是传输线参数计算,5 条传输线的具体计算如下。根据教材查得kmr/21.00kmx/4 .00kmSb/108. 260线路 L1: 101101641012110.21 408.40.440162.8 10401.12 1010.67762LLLLLNRrlXxlBblSQB VMVar线路 L2: 精品资料 - - - 欢迎下载 - - -

20、- - - - - - - - 欢迎下载 名师归纳 - - - - - - - - - -第 10 页，共 22 页 - - - - - - - - - - matlab 电力系统潮流计算202202642022220.21 13027.30.4 130522.8 101303.64 1012.20222LLLLLNRrlXxlBblSQB VMVar线路 L3:(双回路 ) 30330364303233110.21707.3522110.4701422222.8 10703.92 1012.37162LLLLLNRrlXxlBblSQB VMVar线路 L4: 404404644042440

21、.21 6012.60.460242.8 10601.68 1011.01642LLLLLNRrlXxlBblSQB VMVar线路 L5:(双回路 ) 50550564505253110.21202.122110.420422222.8 10201.121010.06862LLLLLNRrlXxlBblSQB VMVar3、2 手工计算FLR1: 221112210112222111222222210()(3.45074.406)(0.02850.6562)11010(10.04420.1142)10.04420.1142()(8.416)(0.0700.1334)110(TTTNaTLLL

22、LNTTTNPSRjXjjMVAVSMWSSjQjMVAPQSRjXjjMVAVPQSRjXV222212114045)(1.34631.7625)(0.40329.5156)110206045200.40329.5156(39.596835.4844)25(4.482635.9144)bGTcbaLLjjMVASSSjjjMVASSSjQSjMVAFLR2: 精品资料 - - - 欢迎下载 - - - - - - - - - - - 欢迎下载 名师归纳 - - - - - - - - - -第 11 页，共 22 页 - - - - - - - - - - matlab 电力系统潮流计算22

23、223332222222222223032252156()(0.0925.042)(0.667936.6024)110(4.493134.1048)4.493134.1048()(27.352)(2.675.0854)110120TTTNcLLLNdGcTLPQSRjXjjMVAVSjMVAPQSRjXjjMVAVSSSSSjQS2(132.9792149.229)LjMVAFLR3: 222244422d2222333223404627()(1.72539.703)(0.10912.5101)110(133.5955149.9956)133.5955149.9956()(7.3514)(24

24、.5146.682)1103025TTTNLLLNeGdTPQSRjXjjMVAVSjMVAPQSRjXjjMVAVSSSSSj33(89.945130.0151)LLQSjMVAFLR4: 222255522222244422450545031()(0.36820.168)(0.10525.7687)110(92.7481133.9937)92.7487133.9937()(12.624)(27.65452.674)11080TTTNeLLLNfGeTLPQSRjXjjMVAVSjMVAPQSRjXjjMVAVSSSSSjQ4(34.9449107.3469)LSjMVAFLR5: 2727

25、07522525065226 315(0.3066.431)(0.05621.1812)3515(15.07820.3422)15.07820.3422(2.14)(0.38990.743)355(20.49451.1266)1537.5ThTLLihLLTSjjMVASSSj QjMVASjjMVASSSSj QjMVAS2226 22226 126 156 26 3(0.3864.34)(0.5145.7793)3520.65050.5451(0.3860.153)(0.13450.0533)3526.33698.7369(0.3866.584)(3.290556.1256)3535(25

26、.TTgfTGTTijjMVASjjMVASjjMVASSSSSSS5114194.12)jMVA计算每一个FLR 的功率分布与电压分布计算如下: FLR1: 精品资料 - - - 欢迎下载 - - - - - - - - - - - 欢迎下载 名师归纳 - - - - - - - - - -第 12 页，共 22 页 - - - - - - - - - - matlab 电力系统潮流计算221140 1.34645 31.762512.8970115115102.103010.0442 8.40.1442 160.8489102.1030101.2541TNbTLbabLPRQXVkVVVV

27、kVPRQXVkVVVVVkVFLR2: 功率分布 : *32*2323*23*2323()(0.0925.042) (132.9792149.229) 1418.6727.39257.042(4.881213.8097)()(27.352)(132.9792149.229) 1418.6727.39TbNLdNLTLTbNLTdNLTLTVVjjSSVjjMVAVVjjSSVZZZZZZZZZZ257.042(108.687122.62)jjMVA电压分布 : 12222(4.881213.8097)(2.675.0854)(7.55128.7243)7.5512 27.38.7243 52

28、2.424102.1030102.103( 2.424)104.527cLLLdbLSSSjjjMVAVkVVVVkVFLR3: 功率分布 : *343*3434*34*3434()(1.72539.703) (89.945130.0151) 1037.9279.07553.73(59.44416.846)()(7.3514)(89.945130.0151) 1037.9279.075GdTLeNLTLTLbNTeNLTLTVVjjSSVjjMVAVVjjSSVZZZZZZZZZZ53.73(31.81160.1256)jjMVA电压分布 : 13333(59.44419.846)(24.51

29、46.682)(83.95426.836)83.954 7.3526.836 149.404105.564396.16eLLLedLSSSjjjMVAVkVVVVkVFLR4: 功率分布 : 精品资料 - - - 欢迎下载 - - - - - - - - - - - 欢迎下载 名师归纳 - - - - - - - - - -第 13 页，共 22 页 - - - - - - - - - - matlab 电力系统潮流计算*534*4545*344*4545()(0.36820.168) (34.9449107.3469) 1037.927=12.96844.168(20.84319.689)(

30、)(12.624)(34.9449107.3469) 1037.927=1TGdLfNLTLTGdLTfNLTLTVVjjSSVjjMVAVVjjSSVZZZZZZZZZZ2.96844.168(1.39844.389)jjMVA电压分布 : 13343(59.44416.846)(24.5146.682)(83.95463.528)83.954 12.663.528 2424.464105.564381.10eLLLedLSSSjjjMVAVkVVVVkVFLR5: 这里我们先将f 点与发电机G5 当做电源 ,经过61TZ与63TZ构成两端供电网络以g 点作为运算负荷进行计算。64(0.38

31、64.134) (20.265670.9293)(22.093837)35(3.90025.1175)0.77210.718(0.3866.584)(20.265670.9293)(22.093837)35(16.506191.7905)0.77210.718TTjjSjMVAjjjSjMVAj电压分布 : 6316363636362625(16.642197.5698)16.64210.38697.56984.13410.9186373726.081420.26560.38670.9293( 0.153)0.116226.081426.197620.49452.11.12664TTTTgTT

32、igTLSSSjMVAVkVVVVVkVVVVV51.81526.197624.3826hiLkVVVV4 设计图及源程序4、1MATLAB仿真相关的原始数据输入格式如下: 1、B1 就是支路参数矩阵,第一列与第二列就是节点编号。节点编号由小到大编写。精品资料 - - - 欢迎下载 - - - - - - - - - - - 欢迎下载 名师归纳 - - - - - - - - - -第 14 页，共 22 页 - - - - - - - - - - matlab 电力系统潮流计算2、 对于含有变压器的支路,第一列为低压侧节点编号,第二列为高压侧节点编号,将变压器的串联阻抗置于低压侧处理,第三列

33、为支路的串列阻抗参数,第四列为支路的对地导纳参数,第五烈为含变压器支路的变压器的变比,第六列为变压器就是否就是否含有变压器的参数,其中 “1”为含有变压器, “0”为不含有变压器。3、B2 为节点参数矩阵,其中第一列为节点注入发电功率参数;第二列为节点负荷功率参数;第三列为节点电压参数;第六列为节点类型参数,其中 “1”为平衡节点 , “2”为 PQ 节点 , “3”为 PV 节点参数。4、X 为节点号与对地参数矩阵。其中第一列为节点编号,第二列为节点对地参数。4、2 潮流计算源程序close all clear all n=input( 请输入节点数 :n=); n1=input( 请输入支

34、路数 :n1=); isb=input( 请输入平衡节点号:isb=); pr=input( 请输入误差精度 :pr=); B1=input( 请输入支路参数:B1=); B2=input( 请输入节点参数:B2=); X=input( 节点号与对地参数:X=); Y=zeros(n); Times=1; %置迭代次数为初始值%创建节点导纳矩阵for i=1:n1 if B1(i,6)=0 %不含变压器的支路p=B1(i,1); q=B1(i,2); Y(p,q)=Y(p,q)-1/B1(i,3); Y(q,p)=Y(p,q); Y(p,p)=Y(p,p)+1/B1(i,3)+0、5*B1(i

35、,4); Y(q,q)=Y(q,q)+1/B1(i,3)+0、5*B1(i,4); else %含有变压器的支路p=B1(i,1); q=B1(i,2); Y(p,q)=Y(p,q)-1/(B1(i,3)*B1(i,5); Y(q,p)=Y(p,q); Y(p,p)=Y(p,p)+1/B1(i,3); Y(q,q)=Y(q,q)+1/(B1(i,5)2*B1(i,3); end end Y OrgS=zeros(2*n-2,1); DetaS=zeros(2*n-2,1); %将 OrgS、DetaS 初始化%创建 OrgS,用于存储初始功率参数h=0; j=0; for i=1:n %对 P

36、Q节点的处理精品资料 - - - 欢迎下载 - - - - - - - - - - - 欢迎下载 名师归纳 - - - - - - - - - -第 15 页，共 22 页 - - - - - - - - - - matlab 电力系统潮流计算if i=isb&B2(i,6)=2 h=h+1; for j=1:n OrgS(2*h-1,1)=OrgS(2*h-1,1)+real(B2(i,3)*(real(Y(i,j)*real(B2(j,3)-imag(Y(i,j)*imag(B2(j,3)+imag(B2(i,3)*(real(Y(i,j)*imag(B2(j,3)+imag(Y(i,j)

37、*real(B2(j,3); OrgS(2*h,1)=OrgS(2*h,1)+imag(B2(i,3)*(real(Y(i,j)*real(B2(j,3)-imag(Y(i,j)*imag(B2(j,3)-real(B2(i,3)*(real(Y(i,j)*imag(B2(j,3)+imag(Y(i,j)*real(B2(j,3); end end end for i=1:n %对 PV 节点的处理 ,注意这时不可再将h 初始化为 0 if i=isb&B2(i,6)=3 h=h+1; for j=1:n OrgS(2*h-1,1)=OrgS(2*h-1,1)+real(B2(i,3)*(re

38、al(Y(i,j)*real(B2(j,3)-imag(Y(i,j)*imag(B2(j,3)+imag(B2(i,3)*(real(Y(i,j)*imag(B2(j,3)+imag(Y(i,j)*real(B2(j,3); OrgS(2*h,1)=OrgS(2*h,1)+imag(B2(i,3)*(real(Y(i,j)*real(B2(j,3)-imag(Y(i,j)*imag(B2(j,3)-real(B2(i,3)*(real(Y(i,j)*imag(B2(j,3)+imag(Y(i,j)*real(B2(j,3); end end end OrgS %创建 PVU 用于存储 PV 节

39、点的初始电压PVU=zeros(n-h-1,1); t=0; for i=1:n if B2(i,6)=3 t=t+1; PVU(t,1)=B2(i,3); end end PVU %创建 DetaS,用于存储有功功率、无功功率与电压幅值的不平衡量h=0; for i=1:n %对 PQ 节点的处理if i=isb&B2(i,6)=2 h=h+1; DetaS(2*h-1,1)=real(B2(i,2)-OrgS(2*h-1,1); DetaS(2*h,1)=imag(B2(i,2)-OrgS(2*h,1); end end t=0; for i=1:n %对 PV 节点的处理 ,注意这时不可

40、再将h 初始化为 0 精品资料 - - - 欢迎下载 - - - - - - - - - - - 欢迎下载 名师归纳 - - - - - - - - - -第 16 页，共 22 页 - - - - - - - - - - matlab 电力系统潮流计算if i=isb&B2(i,6)=3 h=h+1; t=t+1; DetaS(2*h-1,1)=real(B2(i,2)-OrgS(2*h-1,1); DetaS(2*h,1)=real(PVU(t,1)2+imag(PVU(t,1)2-real(B2(i,3)2-imag(B2(i,3)2; end end DetaS %创建 I,用于存储节

41、点电流参数i=zeros(n-1,1); h=0; for i=1:n if i=isb h=h+1; I(h,1)=(OrgS(2*h-1,1)-OrgS(2*h,1)*sqrt(-1)/conj(B2(i,3); end end I %创建 Jacbi(雅可比矩阵 ) Jacbi=zeros(2*n-2); h=0; k=0; for i=1:n %对 PQ 节点的处理if B2(i,6)=2 h=h+1; for j=1:n if j=isb k=k+1; if i=j %对角元素的处理Jacbi(2*h-1,2*k-1)=-imag(Y(i,j)*real(B2(i,3)+real(Y

42、(i,j)*imag(B2(i,3)+imag(I(h,1); Jacbi(2*h-1,2*k)=real(Y(i,j)*real(B2(i,3)+imag(Y(i,j)*imag(B2(i,3)+real(I(h,1); Jacbi(2*h,2*k-1)=-Jacbi(2*h-1,2*k)+2*real(I(h,1); Jacbi(2*h,2*k)=Jacbi(2*h-1,2*k-1)-2*imag(I(h,1); else %非对角元素的处理Jacbi(2*h-1,2*k-1)=-imag(Y(i,j)*real(B2(i,3)+real(Y(i,j)*imag(B2(i,3); Jacb

43、i(2*h-1,2*k)=real(Y(i,j)*real(B2(i,3)+imag(Y(i,j)*imag(B2(i,3); Jacbi(2*h,2*k-1)=-Jacbi(2*h-1,2*k); Jacbi(2*h,2*k)=Jacbi(2*h-1,2*k-1); end if k=(n-1) % 将用于内循环的指针置于初始值,以确保雅可比矩阵换行k=0; end end end end end k=0; 精品资料 - - - 欢迎下载 - - - - - - - - - - - 欢迎下载 名师归纳 - - - - - - - - - -第 17 页，共 22 页 - - - - - -

44、- - - - matlab 电力系统潮流计算for i=1:n %对 PV 节点的处理if B2(i,6)=3 h=h+1; for j=1:n if j=isb k=k+1; if i=j %对角元素的处理Jacbi(2*h-1,2*k-1)=-imag(Y(i,j)*real(B2(i,3)+real(Y(i,j)*imag(B2(i,3)+imag(I(h,1); Jacbi(2*h-1,2*k)=real(Y(i,j)*real(B2(i,3)+imag(Y(i,j)*imag(B2(i,3)+real(I(h,1); Jacbi(2*h,2*k-1)=2*imag(B2(i,3);

45、 Jacbi(2*h,2*k)=2*real(B2(i,3); else %非对角元素的处理Jacbi(2*h-1,2*k-1)=-imag(Y(i,j)*real(B2(i,3)+real(Y(i,j)*imag(B2(i,3); Jacbi(2*h-1,2*k)=real(Y(i,j)*real(B2(i,3)+imag(Y(i,j)*imag(B2(i,3); Jacbi(2*h,2*k-1)=0; Jacbi(2*h,2*k)=0; end if k=(n-1) %将用于内循环的指针置于初始值,以确保雅可比矩阵换行k=0; end end end end end Jacbi %求解修正

46、方程 ,获取节点电压的不平衡量DetaU=zeros(2*n-2,1); DetaU=inv(Jacbi)*DetaS; DetaU %修正节点电压j=0; for i=1:n %对 PQ 节点处理if B2(i,6)=2 j=j+1; B2(i,3)=B2(i,3)+DetaU(2*j,1)+DetaU(2*j-1,1)*sqrt(-1); end end for i=1:n %对 PV 节点的处理if B2(i,6)=3 j=j+1; B2(i,3)=B2(i,3)+DetaU(2*j,1)+DetaU(2*j-1,1)*sqrt(-1); end end B2 %开始循环 * while

47、 abs(max(DetaU)pr OrgS=zeros(2*n-2,1); %! 初始功率参数在迭代过程中就是不累加的,所以在这里必须将其初始化为零矩阵精品资料 - - - 欢迎下载 - - - - - - - - - - - 欢迎下载 名师归纳 - - - - - - - - - -第 18 页，共 22 页 - - - - - - - - - - matlab 电力系统潮流计算h=0; j=0; for i=1:n if i=isb&B2(i,6)=2 h=h+1; for j=1:n OrgS(2*h-1,1)=OrgS(2*h-1,1)+real(B2(i,3)*(real(Y(i,

48、j)*real(B2(j,3)-imag(Y(i,j)*imag(B2(j,3)+imag(B2(i,3)*(real(Y(i,j)*imag(B2(j,3)+imag(Y(i,j)*real(B2(j,3); OrgS(2*h,1)=OrgS(2*h,1)+imag(B2(i,3)*(real(Y(i,j)*real(B2(j,3)-imag(Y(i,j)*imag(B2(j,3)-real(B2(i,3)*(real(Y(i,j)*imag(B2(j,3)+imag(Y(i,j)*real(B2(j,3); end end end for i=1:n if i=isb&B2(i,6)=3

49、h=h+1; for j=1:n OrgS(2*h-1,1)=OrgS(2*h-1,1)+real(B2(i,3)*(real(Y(i,j)*real(B2(j,3)-imag(Y(i,j)*imag(B2(j,3)+imag(B2(i,3)*(real(Y(i,j)*imag(B2(j,3)+imag(Y(i,j)*real(B2(j,3); OrgS(2*h,1)=OrgS(2*h,1)+imag(B2(i,3)*(real(Y(i,j)*real(B2(j,3)-imag(Y(i,j)*imag(B2(j,3)-real(B2(i,3)*(real(Y(i,j)*imag(B2(j,3)

50、+imag(Y(i,j)*real(B2(j,3); end end end OrgS %创建 DetaS h=0; for i=1:n if i=isb&B2(i,6)=2 h=h+1; DetaS(2*h-1,1)=real(B2(i,2)-OrgS(2*h-1,1); DetaS(2*h,1)=imag(B2(i,2)-OrgS(2*h,1); end end t=0; for i=1:n if i=isb&B2(i,6)=3 h=h+1; t=t+1; DetaS(2*h-1,1)=real(B2(i,2)-OrgS(2*h-1,1); DetaS(2*h,1)=real(PVU(t,