pq分解法在潮流分析中的应用

《pq分解法在潮流分析中的应用》由会员上传分享，免费在线阅读，更多相关内容在行业资料-天天文库。

1、PQ分解法在潮流分析中的应用一、实验目的1、了解典型节点系统的潮流和分析计算过程。2、学会运用matlab对电力系统潮流进行计算的方法及原理，我们主要介绍PQ分解法。二、实验要求分析在负荷侧改变无功补偿容量时对潮流分布的影响和对网损的影响，也可进行其他的例如改变回路数，断路器状态等操作观察分析系统潮流变化及网损变化。三、实验过程对于典型的节点系统，分别有IEEE14、30、57、118等，在实验中我们选取了典型30节点的系统。1、计算程序：在matlab中建立m文件程序，输入30个节点以及各支路的数据，其格式如下：bas

2、emva=100;accuracy=0.001;accel=1.8;maxiter=100;%IEEE30-BUSTESTSYSTEM(AmericanElectricPower)%BusBusVoltageAngle---Load-----------Generator-----StaticMvar%NocodeMag.DegreeMWMvarMWMvarQminQmax+Qc/-Qlbusdata=[111.060.00.00.00.00.0000221.0430.021.7012.740.00.0-40500301

3、.00.02.41.20.00.0000401.060.07.61.60.00.0000521.010.094.219.00.00.0-40400601.00.00.00.00.00.0000701.00.022.810.90.00.0000821.010.030.030.00.00.0-30400901.00.00.00.00.00.00001001.00.05.82.00.00.0-62419以上是30个节点对应的数据，其中基准容量取为100MVA，精度值为0.001，加速因子是1.8，算法的最大迭代次数为100次。

4、%Linecode%BusbusRX1/2B=1forlines%nlnrp.u.p.u.p.u.>1or<1tr.tapatbusnllinedata=[120.01920.05750.026401130.04520.18520.020401240.05700.17370.018401340.01320.03790.004201250.04720.19830.020901260.05810.17630.018701460.01190.04140.004501570.04600.11600.010201670.02670

5、.08200.008501680.01200.04200.004501690.00.20800.00.9786100.556000.9699110.2080019100.1100014120.256000.93212130.140001lfybus%formthebusadmittancematrixlfDecouple%Loadflowsolutionbydecouplemethodbusout%Printsthepowerflowsolutiononthescreenlineflow%Computesanddispl

6、aysthelineflowandlosses以上是网络支路中的数据，通过首、末端节点来表示，包括了支路的电阻、电抗、0.5的电纳和变比标幺值，其中标幺值为1时表示线路，大于1时表示降压变压器，小于1时表示升压变压器。其后的一系列是对函数的声明，lfDecouple是形成节点导纳矩阵程序，lfDecouple是PQ分解法程序，busout是输出母线的潮流结果程序，算出各个节点的电压的幅值和相位，lineflow是计算各条支路潮流分布和网损的程序。节点导纳矩阵程序lfbus：%ThisprogramobtainsthBus

7、AdmittanceMatrixforpowerflowsolution%Copyright(c)1998byH.Saadatj=sqrt(-1);i=sqrt(-1);nl=linedata(:,1);nr=linedata(:,2);R=linedata(:,3);X=linedata(:,4);Bc=j*linedata(:,5);a=linedata(:,6);nbr=length(linedata(:,1));nbus=max(max(nl),max(nr));Z=R+j*X;y=ones(nbr,1)./Z;

8、%branchadmittanceforn=1:nbrifa(n)<=0a(n)=1;elseendYbus=zeros(nbus,nbus);%initializeYbustozero%formationoftheoffdiagonalelementsfork=1:nbr;Ybus(nl(k),nr(k))=Ybu