电力系统分析潮流计算最终完整版文档格式.doc

1、%第一列和第二列是节点编号。节点编号由小到大编写%对于含有变压器的支路，第一列为低压侧节点编号，第二列为高压侧节点编号%第三列为支路的串列阻抗参数，含变压器支路此值为变压器短路电抗%第四列为支路的对地导纳参数，含变压器支路此值不代入计算%第五烈为含变压器支路的变压器的变比，变压器非标准电压比%第六列为变压器是否是否含有变压器的参数，其中“1”为含有变压器，“0”为不含有变压器%B2为节点参数矩阵%第一列为节点注入发电功率参数%第二列为节点负荷功率参数%第三列为节点电压参数%第四列%第五列%第六列为节点类型参数，“1”为平衡节点，“2”为PQ节点，“3”为PV节点参数%X为节点号和对地参数矩阵%

2、第一列为节点编号%第二列为节点对地参数%默认算例% n=4;% n1=4;% isb=4;% pr=0.00001;% B1=1 2 0.1667i 0 0.8864 1;1 3 0.1302+0.2479i 0.0258i 1 0;1 4 0.1736+0.3306i 0.0344i 1 0;3 4 0.2603+0.4959i 0.0518i 1 0;% B2=0 0 1 0 0 2;0 -0.5-0.3i 1 0 0 2;0.2 0 1.05 0 0 3;0 -0.15-0.1i 1.05 0 0 1;% X=1 0;2 0.05i;3 0;4 0;clear;clc;num=input

3、（是否采用默认数据？（1-默认数据;2-手动输入）;if num=1 n=4; n1=4; isb=4; pr=0.00001; B1=1 2 0.1667i 0 0.8864 1; B2=0 0 1 0 0 2; X=1 0;else n=input（请输入节点数:n= n1=input（请输入支路数:n1= isb=input（请输入平衡节点号:isb= pr=input（请输入误差精度:pr= B1=input（请输入支路参数:B1= B2=input（请输入节点参数:B2= X=input（节点号和对地参数:X=endTimes=1; %迭代次数%创建节点导纳矩阵Y=zeros（n）;

4、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,4）; Y（q,q）=Y（q,q）+1/B1（i,3）+0.5*B1（i,4）; else %含有变压器的支路 Y（p,q）=Y（p,q）-B1（i,5）/B1（i,3）; Y（p,p）=Y（p,p）+B1（i,5）/B1（i,3）+（1-B1（i,5）/B1（i,3）; Y（q,q）=Y（q,q）+B1（i,5）/B1（i,3）+（B1（

5、i,5）*（B1（i,5）-1）/B1（i,3）; end Y（i,i）=Y（i,i）+X（i,2）; %计及补偿电容电纳disp（导纳矩阵为：disp（Y）; %显示导纳矩阵 %初始化OrgS、DetaSOrgS=zeros（2*n-2,1）;DetaS=zeros（2*n-2,1）; %创建OrgS，用于存储初始功率参数h=0;j=0;n %对PQ节点的处理 if i=isb&B2（i,6）=2 %不是平衡点&是PQ点 h=h+1; for j=1:n %公式8-74 %Pi=ei*（Gij*ej-Bij*fj）+fi*（Gij*fj+Bij*ej） %Qi=fi*（Gij*ej-Bij*

6、fj）-ei*（Gij*fj+Bij*ej） 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（

7、B2（j,3）+imag（Y（i,j）*real（B2（j,3）; endn %对PV节点的处理，注意这时不可再将h初始化为0B2（i,6）=3 %不是平衡点&是PV点 %公式8-75-a 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）;%创建PVU 用于存储PV节点的初始电压PVU=zeros（n-h-1,1）;t=0; if B2（i,6

8、）=3 t=t+1; PVU（t,1）=B2（i,3）;%创建DetaS，用于存储有功功率、无功功率和电压幅值的不平衡量B2（i,6）=2 DetaS（2*h-1,1）=real（B2（i,2）-OrgS（2*h-1,1）; %delPi DetaS（2*h,1）=imag（B2（i,2）-OrgS（2*h,1）; %delQiB2（i,6）=3 DetaS（2*h-1,1）=real（B2（i,1）-OrgS（2*h-1,1）; %delPi DetaS（2*h,1）=real（PVU（t,1）2+imag（PVU（t,1）2-real（B2（i,3）2-imag（B2（i,3）2; %de

9、lUi% DetaS%创建I，用于存储节点电流参数i=zeros（n-1,1）; if i=isb I（h,1）=（OrgS（2*h-1,1）-OrgS（2*h,1）*sqrt（-1）/conj（B2（i,3）;%conj求共轭%创建Jacbi（雅可比矩阵）Jacbi=zeros（2*n-2）;k=0; if B2（i,6）=2 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）=-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 %非对角