From 7bde55ceb86007ed8efaa4de39ba1f7447ac84f1 Mon Sep 17 00:00:00 2001 From: dmy Date: Wed, 14 May 2014 10:00:41 +0800 Subject: [PATCH] =?UTF-8?q?=E5=AE=8C=E6=88=90=E4=BA=86=E7=94=B5=E5=8E=8B?= =?UTF-8?q?=E7=9A=84=E9=83=A8=E5=88=86?= MIME-Version: 1.0 Content-Type: text/plain; charset=UTF-8 Content-Transfer-Encoding: 8bit Signed-off-by: dmy --- DrawGap.asv | 6 ----- FormAA.asv | 12 --------- FormAlphaP.asv | 7 ----- FormG.asv | 11 -------- FormG.m | 6 ++--- FormH.asv | 31 ---------------------- FormLw.m | 11 ++------ FormLz.m | 10 ++----- FormYY.asv | 5 ---- Modification.asv | 18 ------------- OPF.m | 20 +++++++++----- admmatrix.asv | 31 ---------------------- func_deltG.asv | 34 ------------------------ func_deltG.m | 4 +-- modifyadmmatrix.asv | 18 ------------- openfile.asv | 64 --------------------------------------------- 16 files changed, 22 insertions(+), 266 deletions(-) delete mode 100644 DrawGap.asv delete mode 100644 FormAA.asv delete mode 100644 FormAlphaP.asv delete mode 100644 FormG.asv delete mode 100644 FormH.asv delete mode 100644 FormYY.asv delete mode 100644 Modification.asv delete mode 100644 admmatrix.asv delete mode 100644 func_deltG.asv delete mode 100644 modifyadmmatrix.asv delete mode 100644 openfile.asv diff --git a/DrawGap.asv b/DrawGap.asv deleted file mode 100644 index c9d0645..0000000 --- a/DrawGap.asv +++ /dev/null @@ -1,6 +0,0 @@ -function DrawGap(plotGap) -x=find(plotGap); -ts=size(x,2); - - -end \ No newline at end of file diff --git a/FormAA.asv b/FormAA.asv deleted file mode 100644 index 0d82338..0000000 --- a/FormAA.asv +++ /dev/null @@ -1,12 +0,0 @@ -function AA=FormAA(L_1Z,deltG,U_1W,Hcoma,deltH) -tOnes=eye(14); -tZeros=zeros(14); -AA=[ - tOnes,L_1Z,tZeros,tZeros,tZeros,zeros(14,10); - tZeros,tOnes,tZeros,tZeros,-deltG',zeros(14,10); - tZeros,tZeros,tOnes,U_1W,tZeros,zeros(14,10); - tZeros,tZeros,tZeros,tOnes,deltG,zeros(14,10); - tZeros,tZeros,tZeros,tZeros,Hcoma,deltH; - zeros(10,14),zeros(10,14),zeros(10,14),zeros(10,14),deltH',zeros(10,10); - ]; -end \ No newline at end of file diff --git a/FormAlphaP.asv b/FormAlphaP.asv deleted file mode 100644 index 3a8f5d9..0000000 --- a/FormAlphaP.asv +++ /dev/null @@ -1,7 +0,0 @@ -function AlphaP=FormAlphaP(Init_L,deltL,Init_U,deltU) -ti=deltL(delt) -t1=-Init_L./deltL'; -t2=-Init_U./deltU'; -t3=[t1,t2]; -t4=t3() -end \ No newline at end of file diff --git a/FormG.asv b/FormG.asv deleted file mode 100644 index fba473d..0000000 --- a/FormG.asv +++ /dev/null @@ -1,11 +0,0 @@ -function Mat_G=FormG(Volt,PVi,PGi,PG,QG,PD,QD,Loadi) - -Mat_G=[ - PG(PGi); - QG(PVi); - sparse(PD(Loadi)); - sparse(QD(Loadi)); - Volt'; - - ]; -end \ No newline at end of file diff --git a/FormG.m b/FormG.m index 7cf7e72..3c8ab70 100644 --- a/FormG.m +++ b/FormG.m @@ -1,10 +1,10 @@ -function Mat_G=FormG(Volt,PVi,PGi,PG,QG,PD,QD,Loadi,Vbi) +function Mat_G=FormG(Volt,PD,QD,Loadi,Vbi,mVolt,RealVolt,bigM,sigma) Mat_G=[ sparse(PD(Loadi)); sparse(QD(Loadi)); - Volt'; - Volt'; + Volt'-bigM*Vbi-mVolt'-sigma*RealVolt'; + Volt'+bigM*Vbi+mVolt'+sigma*RealVolt'; Vbi; ]; end \ No newline at end of file diff --git a/FormH.asv b/FormH.asv deleted file mode 100644 index 5983fcf..0000000 --- a/FormH.asv +++ /dev/null @@ -1,31 +0,0 @@ -function Mat_H=FormH(Busnum,GB,AngleIJMat,Volt,PG,PD,QG,QD,Y) -t1=real(GB).*cos(AngleIJMat)+imag(GB).*sin(AngleIJMat); -t2=Volt'*Volt; -t3=t1.*t2; -t4=sum(-t3,2);%P -t5=real(GB).*sin(AngleIJMat)-imag(GB).*cos(AngleIJMat); -t6=t2.*t5; -t7=sum(-t6,2);%Q -t8=PG-PD; -t9=QG-QD; -%Mat_H=([(PG-PD)',(QG-QD)'])'+([t4',t7'])'; -Mat_H(1:2:2*Busnum)=t8(1:Busnum)+t4(1:Busnum); -Mat_H(2:2:2*Busnum)=t9(1:Busnum)+t7(1:Busnum); -Mat_H=Mat_H'; -%% -QDcos=textread('300glys.txt'); -t=QD(PD==0 &&); -aa=QD; -QD(QD~=0)=PD(QD~=0)./tan(QDcos); -QD(PD==0)=t; -%% -%%%%一下是学姐给的公式 -AngleIJ=AngleIJMat-angle(GB); -%dP=PG-PD-diag(Volt)*Y*cos(AngleIJ)*Volt'; -dP=PG-PD-diag(Volt)*Y.*cos(AngleIJ)*Volt'; -dQ=QG-QD-diag(Volt)*Y.*sin(AngleIJ)*Volt'; -%Mat_H(1:2:2*Busnum)=dP(1:Busnum);暂时改一下 20111227 -%Mat_H(2:2:2*Busnum)=dQ(1:Busnum);暂时改一下 20111227 -Mat_H=[dP;dQ;]; - -end \ No newline at end of file diff --git a/FormLw.m b/FormLw.m index 01dfdbd..3a08974 100644 --- a/FormLw.m +++ b/FormLw.m @@ -1,11 +1,7 @@ -function Lw=FormLw(Mat_G,Init_U,GenU,Busnum,PVQU,PD,PD0,QD0,Loadi,KK,PF,noDataTransCapacity) +function Lw=FormLw(Mat_G,Init_U,Busnum,PD0,QD0,Loadi) KK=999; -%PU=GenU(:,2);%发电机有功上界 - VoltU=(1.1)*ones(1,Busnum); -%VoltU=10*ones(1,Busnum); PDU=PD0(Loadi); -% PDU=noDataTransCapacity; PDU(PDU>0)=1.200*PDU(PDU>0); PDU(PDU<0)=0.800*PDU(PDU<0); PDU(PDU==0)=0.400; @@ -13,7 +9,6 @@ realPD=PD0(Loadi); indPD=find(realPD>0); PDU(indPD(3:12:end))=1.55*realPD(indPD(3:12:end)); PDU(indPD(9:12:end))=1.05*realPD(indPD(9:12:end)); -%PDU=10*ones(length(Loadi),1); QDU=QD0(Loadi); QDU(QDU>0)=1.200*QDU(QDU>0); QDU(QDU<0)=0.800*QDU(QDU<0); @@ -22,9 +17,7 @@ realQD=QD0(Loadi); indQD=find(realQD>0); QDU(indQD(3:12:end))=1.55*realQD(indQD(3:12:end)); QDU(indQD(9:12:end))=1.05*realQD(indQD(9:12:end)); -% PF=0.85; -% QDU=1.0*PD(Loadi).*sqrt(1 -PF.^2)./PF; -t1=([PDU',QDU',VoltU,VoltU,1*ones(1,Busnum)])'; +t1=([PDU',QDU',0*VoltU,50*ones(1,Busnum),1*ones(1,Busnum)])'; t2=Mat_G+Init_U'-t1; Lw=t2; diff --git a/FormLz.m b/FormLz.m index b2b8984..d6c46f8 100644 --- a/FormLz.m +++ b/FormLz.m @@ -1,8 +1,5 @@ -function Lz=FormLz(Mat_G,Init_L,GenL,Busnum,PVQL,PD,PD0,QD0,Loadi,KK,PF) -KK=999; - +function Lz=FormLz(Mat_G,Init_L,Busnum,PD0,QD0,Loadi) VoltL=(0.9)*ones(1,Busnum); -%VoltL=-10*ones(1,Busnum); PDL=PD0(Loadi); PDL(PDL>0)=0.800*PDL(PDL>0); PDL(PDL<0)=1.200*PDL(PDL<0); @@ -11,7 +8,6 @@ realPD=PD0(Loadi); indPD=find(realPD>0); PDL(indPD(3:12:end))=0.95*realPD(indPD(3:12:end)); PDL(indPD(9:12:end))=0.45*realPD(indPD(9:12:end)); -%PDL=-10*ones(length(Loadi),1); QDL=QD0(Loadi); QDL(QDL>0)=0.800*QDL(QDL>0); QDL(QDL<0)=1.200*QDL(QDL<0); @@ -20,9 +16,7 @@ realQD=QD0(Loadi); indQD=find(realQD>0); QDL(indQD(3:12:end))=0.95*realQD(indQD(3:12:end)); QDL(indQD(9:12:end))=0.95*realQD(indQD(9:12:end)); -% QDL=0*PD(Loadi).*sqrt((1-PF.^2))./PF; -t1=([PDL',QDL',VoltL,VoltL,0*ones(1,Busnum)])'; +t1=([PDL',QDL',-5000*ones(1,length(VoltL)),0*VoltL,0*ones(1,Busnum)])'; t2=Mat_G-Init_L'-t1; Lz=t2; - end \ No newline at end of file diff --git a/FormYY.asv b/FormYY.asv deleted file mode 100644 index 85d2a13..0000000 --- a/FormYY.asv +++ /dev/null @@ -1,5 +0,0 @@ -function FormYY(Init_L,Lul.Lz,Init_U,Luu,Lw,Lz,LxComa) -t=[ - - -] -end \ No newline at end of file diff --git a/Modification.asv b/Modification.asv deleted file mode 100644 index ec3486f..0000000 --- a/Modification.asv +++ /dev/null @@ -1,18 +0,0 @@ -function [Init_Z,Init_L,Init_W,Init_U,Init_Y,deltZ,deltL,deltW,deltU,deltX,deltY,PG,QG,Volt,UAngel]=Modification(Init_Z,Init_L,Init_W,Init_U,Init_Y,deltZ,deltL,deltW,deltU,deltX,deltY,PG,QG,Volt,UAngel,PVi,ContrlCount,Balance,Busnum,PGi) -AlphaP=FormAlphaP(Init_L,deltL,Init_U,deltU); -fprintf('AlphaP %f\n',AlphaP); -AlphaD=FormAlphaD(Init_Z,deltZ,Init_W,deltW); -fprintf('AlphaD %f\n',AlphaD); - -Init_Z=Init_Z+AlphaD*deltZ'; -Init_L=Init_L+AlphaP*deltL'; -Init_W=Init_W+AlphaD*deltW'; -Init_U=Init_U+AlphaP*deltU'; -Init_Y=Init_Y+AlphaD*deltY'; -PG(PVi)=PG(PGi)+deltX(size(PGi,1)); -QG(PVi)=QG(PVi)+deltX(size(PGi,1)+1,1:(size(PVi,1)+size(PGi,1))); -t=deltX(size(PVi,1)+size(PGi,1)+1,ContrlCount)'; -t(2*Balance-1)=0; -Volt=Volt+AlphaP*t(2:2:2*Busnum); -UAngel=UAngel+AlphaP*t(1:2:2*Busnum); -end \ No newline at end of file diff --git a/OPF.m b/OPF.m index e84a22d..bf94e32 100644 --- a/OPF.m +++ b/OPF.m @@ -6,7 +6,7 @@ clear %% thesis=ForThesis(1,62); [kmax,Precision,UAngel,Volt,Busnum,PVi,PVu,Balance,Y,Angle,P0,Q0,r,c,GB,Linei,Linej,Transfori,Transforj,GenU,GenL,GenC,PG,QG,PD,QD,CenterA,PGi,PVQU,PVQL,Liner,Linex,Lineb,Transforr,Transforx,Branchi,Branchg,Branchb,Transfork0]= ... - pf('C:\lz\结题后重新计算\重构结果\78\木材厂堤Ⅰ线_2229880_2013-10-29\木材厂堤Ⅰ线_2229880_2013-10-29_iPso_newFile.txt'); + pf('E:\算例\东际911_2751267_2012-09-05\newFIle9-2.txt'); % pf('E:\算例\柳金Ⅰ926_21671693_2012-09-06\newFIle16.txt'); %pf('D:\Project\青秀降损项目\最小化潮流\最小潮流算例\原始\津头站津视922（3-1）_0.5_120%.txt'); %pf('D:\Project\最小化潮流\最小潮流算例\仙海919.txt'); @@ -47,10 +47,16 @@ plotGap=zeros(1,60); ContrlCount=size(Loadi,1)*2+Busnum*2+Busnum; kmax=60; Precision=Precision/1; +sigma=0.03; %% 加误差 -PD0(Loadi)=PD0(Loadi).*(1+normrnd(0,0.05,length(Loadi),1)); -QD0(Loadi)=QD0(Loadi).*(1+normrnd(0,0.05,length(Loadi),1)); +PD0(Loadi)=PD0(Loadi).*(1+normrnd(0,sigma,length(Loadi),1)); +QD0(Loadi)=QD0(Loadi).*(1+normrnd(0,sigma,length(Loadi),1)); +mPD=PD0; +mQD=QD0; +RealVolt=Volt0; +mVolt=Volt0.*(1+normrnd(0,sigma,1,length(Volt0))); Vbi=sparse(0.1*ones(Busnum,1)); +bigM=50; while(abs(Gap)>Precision) if KK>kmax break; @@ -77,11 +83,11 @@ while(abs(Gap)>Precision) %% 形成方程矩阵 Luu=Init_U'.*Init_W'+Init_u*ones(RestraintCount,1); Lul=Init_L'.*Init_Z'-Init_u*ones(RestraintCount,1); - Mat_G=FormG(Volt,PVi,PGi,PG,QG,PD,QD,Loadi,Vbi); + Mat_G=FormG(Volt,PD,QD,Loadi,Vbi,mVolt,RealVolt,bigM,sigma); Mat_H=FormH(Busnum,Volt,PG,PD,QG,QD,Y,UAngel,r,c,Angle,Loadi); Ly=Mat_H; - Lz=FormLz(Mat_G,Init_L,GenL,Busnum,PVQL,PD,RealPD,RealQD,Loadi,KK,PF); - Lw=FormLw(Mat_G,Init_U,GenU,Busnum,PVQU,PD,RealPD,RealQD,Loadi,KK,PF); + Lz=FormLz(Mat_G,Init_L,Busnum,PD0,QD0,Loadi); + Lw=FormLw(Mat_G,Init_U,Busnum,PD0,QD0,Loadi); Lx=FormLx(deltF,deltH,Init_Y,deltG,Init_Z,Init_W); YY=FormYY(Lul,Lz,Ly,Luu,Lw,Lx); %% 开始解方程 @@ -93,6 +99,6 @@ while(abs(Gap)>Precision) Gap=(Init_L*Init_Z'-Init_U*Init_W'); KK=KK+1; end -fprintf('%f\n',sum(full(Vbi))); +fprintf('目标函数: %f\n',sum(full(Vbi))); toc diff --git a/admmatrix.asv b/admmatrix.asv deleted file mode 100644 index 1896345..0000000 --- a/admmatrix.asv +++ /dev/null @@ -1,31 +0,0 @@ -function [G,B,GB,Y,r,c,Angle] = admmatrix(Busnum,Linei,Linej,Liner,Linex,Lineb,Transfori... - ,Transforj,Transforr,Transforx,Transfork0,Branchi,Branchb) -%************************************************************************** -% 程序功能 : 子函数——形成节点导纳矩阵Y -% 编 者： -% 编制时间：2010.12 -%************************************************************************** -%% 支路导纳计算 -G = -sparse(Linei,Linej,Liner./(Liner.^2+Linex.^2),Busnum,Busnum) - sparse(Linej,Linei,Liner./(Liner.^2+Linex.^2),Busnum,Busnum); -G = G - sparse(1:Busnum,1:Busnum,sum(G,2)'); % 计算各线路支路电导 -B = sparse(Linei,Linej,Linex./(Liner.^2+Linex.^2),Busnum,Busnum)+sparse(Linej,Linei,Linex./(Liner.^2+Linex.^2),Busnum,Busnum); -B = B - sparse(1:Busnum,1:Busnum,sum(B,2)')+sparse(Linei,Linei,Lineb,Busnum,Busnum)+sparse(Linej,Linej,Lineb,Busnum,Busnum); -%% 变压器支路计算 -if Transfori>0 - mr = Transforr./(Transforr.^2+Transforx.^2); % 计算变压器支路电导 - mx = -Transforx./(Transforr.^2+Transforx.^2); % 计算变压器支路电纳 - G = G-sparse(Transfori,Transforj,mr./Transfork0,Busnum,Busnum)-sparse(Transforj,Transfori,mr./Transfork0,Busnum,Busnum)... - +sparse(Transfori,Transfori,mr./Transfork0./Transfork0,Busnum,Busnum)+sparse(Transforj,Transforj,mr,Busnum,Busnum); - B = B-sparse(Transfori,Transforj,mx./Transfork0,Busnum,Busnum)-sparse(Transforj,Transfori,mx./Transfork0,Busnum,Busnum)... - +sparse(Transfori,Transfori,mx./Transfork0./Transfork0,Busnum,Busnum)+sparse(Transforj,Transforj,mx,Busnum,Busnum); -end -%% 接地支路计算 -if Branchi>0 % 判断有无接地支路 - B = B+sparse(Branchi,Branchi,Branchb,Busnum,Busnum); -end -%% 化作极坐标形式 -GB = G+B.*1i; %将电导，电纳合并，写成复数形式 -Y = abs(GB); %求节点导纳幅值 -[r,c] = find(Y); -Angle = angle(GB(GB~=0)); %求节点导纳角度 -%Angle=angle(GB); \ No newline at end of file diff --git a/func_deltG.asv b/func_deltG.asv deleted file mode 100644 index 5355449..0000000 --- a/func_deltG.asv +++ /dev/null @@ -1,34 +0,0 @@ -function deltG=func_deltG(Busnum,PVi,PGi,Loadi,PD,QD) -sizeLoadi=size(Loadi,1); -dg3_dPD=sparse(1:size(Loadi,1),1:size(Loadi,1),ones(size(Loadi,1),1),size(Loadi,1),size(Loadi,1)); -dg4_dPD=sparse(size(Loadi,1),length(Loadi)); -dg5_dPD=sparse(size(Loadi,1),Busnum); -dg6_dPD=dg5_dPD; -dg7_dPD=sparse(sizeLoadi,Busnum); -%% -dg3_dQD=sparse(length(Loadi),length(Loadi)); -dg4_dQD=sparse(1:size(Loadi,1),1:size(Loadi,1),ones(size(Loadi,1),1),size(Loadi,1),size(Loadi,1)); -dg5_dQD=sparse(size(Loadi,1),Busnum); -dg6_dQD=dg5_dQD; -dg7_dQD=sparse(sizeLoadi,Busnum); -%% -dg3_dx=sparse(2*Busnum,sizeLoadi); -dg4_dx=sparse(2*Busnum,length(Loadi)); -dg5_dx=[sparse(1:Busnum,1:Busnum,ones(Busnum,1),Busnum,Busnum); - sparse(Busnum,Busnum); - ]; -dg6_dx=dg5_dx; -dg7_dx=sparse(2*Busnum,Busnum); -%% -dg3_dvbi=sparse(Busnum,sizeLoadi); -dg4_dvbi=sparse(Busnum,length(Loadi)); -dg5_dvbi=sparse(Busnum,Busnum); -dg6_dvbi=sparse(Busnum,Busnum); -dg6_dvbi=sparse(Busnum,Busnum); -%% -deltG=[dg3_dPD,dg4_dPD,dg5_dPD,dg6_dPD; - dg3_dQD,dg4_dQD,dg5_dQD,dg6_dQD; - dg3_dx,dg4_dx,dg5_dx,dg6_dx; - dg3_dvbi,dg4_dvbi,dg5_dvbi,dg6_dvbi; -]; -end \ No newline at end of file diff --git a/func_deltG.m b/func_deltG.m index ee01ea9..d439ceb 100644 --- a/func_deltG.m +++ b/func_deltG.m @@ -22,8 +22,8 @@ dg7_dx=sparse(2*Busnum,Busnum); %% dg3_dvbi=sparse(Busnum,sizeLoadi); dg4_dvbi=sparse(Busnum,length(Loadi)); -dg5_dvbi=sparse(Busnum,Busnum); -dg6_dvbi=sparse(Busnum,Busnum); +dg5_dvbi=sparse(-eye(Busnum,Busnum)); +dg6_dvbi=sparse(eye(Busnum,Busnum)); dg7_dvbi=sparse(eye(Busnum,Busnum)); %% deltG=[dg3_dPD,dg4_dPD,dg5_dPD,dg6_dPD,dg7_dPD; diff --git a/modifyadmmatrix.asv b/modifyadmmatrix.asv deleted file mode 100644 index a917e8d..0000000 --- a/modifyadmmatrix.asv +++ /dev/null @@ -1,18 +0,0 @@ -function [new_G,new_B,GB,Y,r,c,Angle] = modifyadmmatrix(ii,jj,G,GG,B,BB) -%************************************************************************** -% 程序功能 : 子函数——形成节点导纳矩阵Y -% 编 者： -% 编制时间：2010.12 -%************************************************************************** -%% 支路导纳计算 -new_G=G; -new_G(ii,jj)=new_G(ii,jj)-G(ii,jj); -new_G(jj,ii)=new_G(jj,ii)-G(jj,ii); -new_G(ii,ii)=new_G(ii,ii)+G(ii,jj); -new_G(jj,jj)=new_G(jj,jj)+G(ii,jj); - -%% 化作极坐标形式 -GB = G+B.*1i; %将电导，电纳合并，写成复数形式 -Y = abs(GB); %求节点导纳幅值 -[r,c] = find(Y); -Angle = angle(GB(GB~=0)); %求节点导纳角度 \ No newline at end of file diff --git a/openfile.asv b/openfile.asv deleted file mode 100644 index 5095547..0000000 --- a/openfile.asv +++ /dev/null @@ -1,64 +0,0 @@ -function [Busnum,Balance,PQstandard,Precision,Linei,Linej,Liner,Linex,Lineb,kmax,Transfori ,... - Transforj,Transforr,Transforx,Transfork0,Branchi,Branchb,Pointpoweri,PG,QG,PD,QD,PVi,PVu,Gen,GenU,GenL,GenC,CenterA,PGi,PVQU,PVQL] = openfile(FileName) -%************************************************************************** -% 程序简介 : 子函数——读取潮流计算所需数据 -% 编 者： -% 编制时间 ：2010.12 -%************************************************************************** -data = dlmread(FileName); % 一次读入全部数据 -Busnum= data(1,1); % 节点数 -PQstandard = data(1,3); % 基准容量 -kmax = data(1,4); %最大迭代次数 -Precision = data(2,1); % 精度 -Balance = data(3,2); % 生成1到节点号的列向量 -CenterA=data(1,5); %中心参数 -LineNum=data(1,2); %支路数 -%% 各参数矩阵分块 -zeroRow = find(data(:,1)==0); %查找第一列元素为零的行号 -line = data(zeroRow(1)+1:zeroRow(2)-1,:); % 形成线路参数矩阵 -ground = data(zeroRow(2)+1:zeroRow(3)-1,:); % 形成对地支路参数矩阵 -tran = data(zeroRow(3)+1:zeroRow(4)-1,:); % 形成变压器参数矩阵 -buspq = data(zeroRow(4)+1:zeroRow(5)-1,:); % 形成节点功率参数矩阵 -PV = data(zeroRow(5)+1:zeroRow(6)-1,:); % 形成pv节点功率参数矩阵 -Gen=data(zeroRow(6)+1:zeroRow(7)-1,:); -%% 线路参数矩阵分块 -Linei = line(:,2); % 节点i -Linej= line(:,3); % 节点j -Liner = line(:,4); % 线路电阻 -Linex = line(:,5); % 线路电抗 -Lineb = line(:,6); % b/2 -%% 对地支路参数矩阵 -Branchi = ground(:,1); % 对地支路节点号 -Branchb = ground(:,2); % 对地支路的导纳 -%% 变压器参数矩阵 -Transfori = tran(:,2); % 节点i -Transforj= tran(:,3); % 节点j -Transforr = tran(:,4); % 变压器电阻 -Transforx= tran(:,5); % 变压器电抗 -Transfork0 = tran(:,6); % 变压器变比 -%% 节点功率参数矩阵 -Pointpoweri = buspq(:,1); -PG=buspq(:,2); % 发电机有功 -QG=buspq(:,3); % 发电机无功 -PD=buspq(:,4); % 负荷有功 -QD=buspq(:,5); % 负荷无功 -%% pv节点功率参数矩阵 -PVi = PV(:,1); % PV节点的节点号 -PVu = PV(:,2); % PV节点电压 -PVQL=PV(:,3);%PV节点无功下限 -PVQU=PV(:,4); %PV节点无功上限 -%% 发电机参数 -%GenU=Gen(:,[1 5 6]); -%GenL=Gen(:,[1 7 8]); -GenC=Gen(:,[1 2:4]); -t=GenC(:,2); -GenC(:,2)=GenC(:,4); -GenC(:,4)=t; -t=Gen(:,[1 5]); -%GenL=[t,PVQL(PVi)]; -GenL=t;%有功下界 -t=Gen(:,[1 6]); -%GenU=[t,PVQU(PVi)]; -GenU=t;%有功上届 -PGi=Gen(:,1);%发电机节点号 -end \ No newline at end of file