电力系统分析潮流计算最终完整版文档格式.docx
《电力系统分析潮流计算最终完整版文档格式.docx》由会员分享，可在线阅读，更多相关《电力系统分析潮流计算最终完整版文档格式.docx（21页珍藏版）》请在冰豆网上搜索。
%第一列与第二列就是节点编号。
节点编号由小到大编写
%对于含有变压器得支路，第一列为低压侧节点编号，第二列为高压侧节点编号
%第三列为支路得串列阻抗参数，含变压器支路此值为变压器短路电抗
%第四列为支路得对地导纳参数，含变压器支路此值不代入计算
%第五烈为含变压器支路得变压器得变比，变压器非标准电压比
%第六列为变压器就是否就是否含有变压器得参数，其中“1”为含有变压器，“0”为不含有变压器
%B2为节点参数矩阵
%第一列为节点注入发电功率参数
%第二列为节点负荷功率参数
%第三列为节点电压参数
%第四列
%第五列
%第六列为节点类型参数，“1”为平衡节点，“2”为PQ节点，“3”为PV节点参数
%X为节点号与对地参数矩阵
%第一列为节点编号
%第二列为节点对地参数
%默认算例
%n=4;
%n1=4;
%isb=4;
%pr=0、00001;
%B1=[120、1667i00、88641;
130、1302+0、2479i0、0258i10;
140、1736+0、3306i0、0344i10;
340、2603+0、4959i0、0518i10];
%B2=[001002;
0-0、5-0、3i1002;
0、201、05003;
0-0、15-0、1i1、05001];
%X=[10;
20、05i;
30;
40];
clear;
clc;
num=input（'
就是否采用默认数据？
（1-默认数据;
2-手动输入）'
）;
ifnum==1
n=4;
n1=4;
isb=4;
pr=0、00001;
B1=[120、1667i00、88641;
B2=[001002;
X=[10;
else
n=input（'
请输入节点数:
n='
n1=input（'
请输入支路数:
n1='
isb=input（'
请输入平衡节点号:
isb='
pr=input（'
请输入误差精度:
pr='
B1=input（'
请输入支路参数:
B1='
B2=input（'
请输入节点参数:
B2='
X=input（'
节点号与对地参数:
X='
end
Times=1;
%迭代次数
%创建节点导纳矩阵
Y=zeros（n）;
fori=1:
n1
ifB1（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（i,5）*（B1（i,5）-1））/B1（i,3）;
end
Y（i,i）=Y（i,i）+X（i,2）;
%计及补偿电容电纳
disp（'
导纳矩阵为：
'
disp（Y）;
%显示导纳矩阵
%初始化OrgS、DetaS
OrgS=zeros（2*n-2,1）;
DetaS=zeros（2*n-2,1）;
%创建OrgS，用于存储初始功率参数
h=0;
j=0;
n%对PQ节点得处理
ifi~=isb&
B2（i,6）==2%不就是平衡点&
就是PQ点
h=h+1;
forj=1:
n
%公式8-74
%Pi=ei*（Gij*ej-Bij*fj）+fi*（Gij*fj+Bij*ej）
%Qi=fi*（Gij*ej-Bij*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（B2（j,3））+imag（Y（i,j））*real（B2（j,3）））;
end
n%对PV节点得处理，注意这时不可再将h初始化为0
B2（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;
ifB2（i,6）==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）;
%delQi
B2（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;
%delUi
%DetaS
%创建I，用于存储节点电流参数
i=zeros（n-1,1）;
ifi~=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;
ifB2（i,6）==2
forj=1:
ifj~=isb
k=k+1;
ifi==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%非对角元素得处理
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）=-Jacbi（2*h-1,2*k）;
Jacbi（2*h,2*k）=Jacbi（2*h-1,2*k-1）;
ifk==（n-1）%将用于内循环得指针置于初始值，以确保雅可比矩阵换行
k=0;
n%对PV节点得处理
升级会员 