电力系统潮流计算程序详细汇总Word文件下载.docx
- 文档编号:7234930
- 上传时间:2023-05-08
- 格式:DOCX
- 页数:10
- 大小:15.68KB
电力系统潮流计算程序详细汇总Word文件下载.docx
《电力系统潮流计算程序详细汇总Word文件下载.docx》由会员分享,可在线阅读,更多相关《电力系统潮流计算程序详细汇总Word文件下载.docx(10页珍藏版)》请在冰点文库上搜索。
功、电压幅值、电压相位。
%“1”为平衡节点,“2”为PQ节点,“3”为PV节点参数。
n=input('
请输入节点数:
n='
;
n1=input('
请输入支路数:
n1='
isb=input('
请输入平衡节点号:
isb='
pr=input('
请输入误差精度:
pr='
B1=input('
请输入支路参数:
B1='
B2=input('
请输入节点参数:
B2='
Y=zeros(n;
N=1;
%建立节点导纳矩阵
fori=1:
n1
ifB1(i,8==0
p=B1(i,1;
q=B1(i,2;
Y(p,q=Y(p,q-1/(B1(i,3+B1(i,4-B1(i,5;
Y(q,p=Y(p,q;
Y(p,p=Y(p,p+1/(B1(i,3+B1(i,4+0.5*B1(i,6+B1(i,5;
Y(q,q=Y(q,q+1/(B1(i,3+B1(i,4+0.5*B1(i,6+B1(i,5;
else
Y(p,q=Y(p,q-1/((B1(i,3+B1(i,4*B1(i,7-B1(i,5;
Y(p,p=Y(p,p+1/(B1(i,3+B1(i,4+B1(i,5;
Y(q,q=Y(q,q+1/(B1(i,7^2*(B1(i,3+B1(i,4+B1(i,5;
end
Y
G=real(Y;
B=imag(Y;
PriS=zeros(2*n-2,1;
ImbS=zeros(2*n-2,1;
%创建PriS,用于存储初始功率参数
h=0;
j=0;
n
ifi~=isb&
B2(i,2==1
h=h+1;
forj=1:
PriS(2*h-1,1=PriS(2*h-1,1+B2(i,5*(G(i,j*B2(j,5-B(i,j*B2(j,6+B2(i,6*(G(i,j*B2(j,6+B(i,
j*B2(j,5;
PriS(2*h,1=PriS(2*h,1+B2(i,6*(G(i,j*B2(j,5-B(i,j*B2(j,6-B2(i,5*(G(i,j*B2(j,6+B(i,j*B
2(j,5;
B2(i,2==2
PriS(2*h-1,1=PriS(2*h-1,1+
B2(i,5*(G(i,j*B2(j,5-B(i,j*B2(j,6+B2(i,6*(G(i,j*B2(j,6+B(i,j*B2(j,5;
PriS(2*h,1=PriS(2*h,1+
B2(i,6*(G(i,j*B2(j,5-B(i,j*B2(j,6-B2(i,5*(G(i,j*B2(j,6+B(i,j*B2(j,5
PriS
U3=zeros(n-h-1,2;
%U3存储PV节点的初始电压
t=0;
ifB2(i,2==2
t=t+1;
U3(t,1=B2(i,5;
U3(t,2=B2(i,6;
U3
%ImbS于存储有功功率、无功功率和电压幅值的不平衡量
ImbS(2*h-1,1=B2(i,3-PriS(2*h-1,1;
ImbS(2*h,1=B2(i,4-PriS(2*h,1;
ImbS(2*h,1=U3(t,1^2+U3(t,2^2-B2(i,5^2-B2(i,6^2;
ImbS
I=zeros(n-1,1;
%I,存储节点电流参数
ifi~=isb
I(h,1=(PriS(2*h-1,1-PriS(2*h,1*sqrt(-1/conj(B2(i,5+B2(i,6*sqrt(-1;
end
I
Jacbi=zeros(2*n-2;
%Jacbi(雅可比矩阵
k=0;
ifB2(i,2==1
ifj~=isb
k=k+1;
ifi==j
Jacbi(2*h-1,2*k-1=-B(i,j*B2(i,5+G(i,j*B2(i,6+imag(I(h,1;
Jacbi(2*h-1,2*k=G(i,j*B2(i,5+B(i,j*B2(i,6+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;
Jacbi(2*h-1,2*k-1=-B(i,j*B2(i,5+G(i,j*B2(i,6;
Jacbi(2*h-1,2*k=G(i,j*B2(i,5+B(i,j*B2(i,6;
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
Jacbi(2*h-1,2*k-1=-B(i,j*B2(i,5+G(i,j*B2(i,6+imag(I(h,1;
Jacbi(2*h-1,2*k=G(i,j*B2(i,5+B(i,j*B2(i,6+real(I(h,1;
Jacbi(2*h,2*k-1=2*B2(i,6;
Jacbi(2*h,2*k=2*B2(i,5;
Jacbi(2*h-1,2*k-1=-B(i,j*B2(i,5+G(i,j*B2(i,6;
Jacbi(2*h-1,2*k=G(i,j*B2(i,5+B(i,j*B2(i,6;
Jacbi(2*h,2*k-1=0;
Jacbi(2*h,2*k=0;
Jacbi
%求解修正方程,获取节点电压的不平衡量
ImbU=zeros(2*n-2,1;
ImbU=inv(Jacbi*ImbS;
ImbU
%修正节点电压
j=j+1;
B2(i,5=B2(i,5+ImbU(2*j,1;
B2(i,6=B2(i,6+ImbU(2*j-1,1;
B2
whileabs(max(ImbU>
pr
PriS(2*h-1,1=PriS(2*h-1,1+B2(i,5*(G(i,j*B2(j,5-B(i,j*B2(j,6+B2(i,6*(G(i,j*B2(j,6+B(i,j*B2(j,5;
PriS(2*h,1=PriS(2*h,1+B2(i,6*(G(i,j*B2(j,5-B(i,j*B2(j,6-B2(i,5*(G(i,j*B2(j,6+B(i,j*B2(j,5;
%创建ImbS
n%对PQ节点的处理
ImbS(2*h,1=B2(i,4-PriS(2*h,1;
endendt=0;
fori=1:
n%对PV节点的处理ifi~=isb&
B2(i,2==2h=h+1;
t=t+1;
ImbS(2*h-1,1=B2(i,3-PriS(2*h-1,1;
ImbS(2*h,1=U3(t,1^2+U3(t,2^2-B2(i,5^2-B2(i,6^2;
endendImbS%创建II=zeros(n-1,1;
h=0;
nifi~=isbh=h+1;
I(h,1=(PriS(2*h-1,1-PriS(2*h,1*sqrt(-1/conj(B2(i,5+B2(i,6*sqrt(-1;
endendI%创建JacbiJacbi=zeros(2*n-2;
k=0;
nifB2(i,2==1h=h+1;
forj=1:
nifj~=isbk=k+1;
ifi==jJacbi(2*h-1,2*k-1=-B(i,j*B2(i,5+G(i,j*B2(i,6+imag(I(h,1;
Jacbi(2*h-1,2*k=G(i,j*B2(i,5+B(i,j*B2(i,6+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;
elseJacbi(2*h-1,2*k-1=-B(i,j*B2(i,5+G(i,j*B2(i,6;
Jacbi(2*h-1,2*k=G(i,j*B2(i,5+B(i,j*B2(i,6;
Jacbi(2*h,2*k-1=-Jacbi(2*h-1,2*k;
Jacbi(2*h,2*k=Jacbi(2*h-1,2*k-1;
endifk==(n-1k=0;
endendendendendk=0;
nifB2(i,2==2h=h+1;
ifi==jJacbi(2*h-1,2*k-1=-B(i,j*B2(i,5+G(i,j*B2(i,6+imag(I(h,1;
Jacbi(2*h,2*k=2*B2(i,5;
elseJacbi(2*h-1,2*k-1=-B(i,j*B2(i,5+G(i,j*B2(i,6;
Jacbi(2*h-1,2*k=G(i,j*B2(i,5+B(i,j*B2(i,6;
Jacbi(2*h,2*k-1=0;
Jacbi(2*h,2*k=0;
endifk==(n-1k=0;
endendendendendJacbiImbU=zeros(2*n-2,1;
ImbU=inv(Jacbi*ImbS;
ImbU%修正节点电压j=0;
nifB2(i,2==1j=j+1;
B2(i,5=B2(i,5+ImbU(2*j,1;
B2(i,6=B2(i,6+ImbU(2*j-1,1;
endendfori=1:
nifB2(i,2==2j=j+1;
endendB2N=N+1;
%迭代次数加1endN
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 电力系统 潮流 计算 程序 详细 汇总