q分解法潮流计算matlab小程序

《q分解法潮流计算matlab小程序》由会员上传分享，免费在线阅读，更多相关内容在学术论文-天天文库。

1、-Q分解法潮流计算matlab小程序1(2010-06-2314:43:08)转载▼标签：杂谈题B是东南人孕陈怡编的电力系统分析里，五节点的算例，本程序的通用性釘待提高。各元件的有关数据如下：(1)发电机G1的端电压为lp.u，发出的有功、无功可调节；(2)发电机G2的端电压为lp.u,按指定的有功P=0.5发电，取8=0.0001;(3)其中，G1为平衡节点，编号为5;G2为PV节点，编号为4;其余为PQ节点，编号为1、2、3。其它数据见上图所示。(4)木程序均釆用标么值进行计算，角度为度()，但程序运行屮均是用弧度进行计在编程开始吋，先创建了MATLAB数掘，如下所示，仅在M

2、ATLAB命令窗□巾有效:systemdata.dat表＞5“‘系统数1R”：552311branch.dat表示“支路数据”:表4-1支路数据支路状态节点号节点号RXB(K)1120.0250.080.141130.030.10.181230.020.060.1032400.19051.052233500.19051.0522generator.dat农示“发电机数据”:表4-2发电机数据140.500.01500.01500.01500.0150140.82151.000.01500.01500.01500.0150140.82loads.dat丧不•“负荷数掘”:表4-3负

3、荷数裾节点状态节点号PQ11-0.8055-0.53212-0.1800-0.12013001.00.5pvnode.dat表示“节点数据”：44.600000.00001balance.dat表示“平衡节点数据”：51.0000下面对木程序作简单的说明：其巾，node_counts表;‘节点、数”，branch_counts表不•“女路数（含并联文路、接地女路）”，generatorcounts表不“发电机数”，loadcounts表不“货荷数”，pvnodecounts表不“PV'i*/点数”，balcounts表示“平衡节点数”，branch表示“文路数則•”，genera

4、tor表示“发电机数损”，loads表示“负荷数拋”，pvnode表示“Vf点数据”，balance表术“Y•衡Vf点数据”，Y表示“Vf点导纳矩阵”，B1表示“B”’，B2表示“B"”，U表示“节点电压”，dPQ表示“冇功无功的不平衡量”，K表示“迭代次数”，eps表示“收敛精度”，Ski、S』I表示“支路潮流”，p表示“注入有功”，q表示“注入无功”。主程序：%==========^算例为5节点5支路网（东南大学陈怡编）P-Q法潮流汁算土程序/，丄，*1«*1««9/%M函数在Workspace看不到炎M的位function[UPQ，SIJ，SJI，BalanceNodeP

5、Q，LossIJPQ，LossP，LossQ]=flowmain(systemdata，branch，generator，loads,pvnode,balance)clearclc%足义全局变fi2010-1globalnode_counts;%节点数globalbranch_counts;%i路数(念并联jJd路、接地支路)globalgenerator_counts;%发电机数globalload_counts;%负倚数globalpvnode_counts;%PV节点数globalbal_counts;%平衡节点globalbranch;%文路数掘globalgenerat

6、or;%发电机数据globalloads;%负荷数掘globalpvnode;%PV节点数掘globalbalance;%平衡节点数拋globalY;%节点导纳矩阵globalB1;%globalB2;globalU;%P点电压globaldPQ；globalKKK;%迭代次数globaleps;%收敛精度globalSIJ;%支路潮流globalp;%注入有功globalq;%注入无功%下述A法仅在MATLAB命令窗II屮有效loadsystemdata.dat;%系统数据loadbranch.dat;%支路数据loadgenerator.dat;%发啦机数据loadloads

7、.dat;%负荷数据loadpvnode.dat;%PV节点数据loadbalance.dat;%平衡节点数据%赋初值node_counts=systemdata(1);branch_counts=systemdata(2);generator_counts=systemdata(3);load_counts=systenidata(4);pvnode_counts=systemdata(5);bal_counts=systemdata(6);%%%Vf点异纳矩阵Y的形成Y=pr