潮流计算(matlab)实例计算.docx

精品文档
精品文档
潮流例题：根据给定的参数或工程具体要求（如图），收集和查阅资料；学****相关软件（软件自选:
本设计选择
Matlab进行设计）
O4-j(IO5
().1
精品文档
精品文档
(H+);y(1,4)=1/(+);y(1,5)=1/(+);
y(2,3)=1/(+);y(2,5)=1/(+);
y(3,4)=1/(+);
y(4,5)=0;fori=1:5forj=i:5y(j,i)=y(i,j);endend
Y=0;
%求互导纳fori=1:5forj=1:5ifi~=jY(i,j)=-y(i,j);endendend%求自导纳fori=1:5
Y(i,i)=sum(y(i,:));endY%Y为导纳矩阵G=real(Y);B=imag(Y);%原始节点功率S(1)=+;
S(2)=--;
S(3)=--;
S(4)=--;
S(5)=0;
P=real(S);Q=imag(S);%赋初值
U=ones(1,5);U(5)=;e=zeros(1,5);
ox=ones(8,1);fx=ones(8,1);
count=0%计算迭代次数
whilemax(fx)>1e-5fori=1:4
forj=1:4
H(i,j)=0;N(i,j)=0;M(i,j)=0;L(i,j)=0;oP(i)=0;oQ(i)=0;endendfori=1:4
精品文档
精品文档
forj=1:5
oP(i)=oP(i)-U(i)*U(j)*(G(i,j)*cos(e(i)-e(j))+B(i,j)*sin(e(i)-e(j)));oQ(i)=oQ(i)-U(i)*U(j)*(G(i,j)*sin(e(i)-e(j))-B(i,j)*cos(e(i)-e(j)));end
oP(i)=oP(i)+P(i);oQ(i)=oQ(i)+Q(i);
end
fx=[oP,oQ]';
%求雅克比矩阵
%当i~=j时候求H,N,M,L如下：
fori=1:4
forj=1:4
ifi~=jH(i,j)=-U(i)*U(j)*(G(i,j)*sin(e(i)-e(j))-B(i,j)*cos(e(i)-e(j)));N(i,j)=-U(i)*U(j)*(G(i,j)*cos(e(i)-e(j))+B(i,j)*sin(e(i)-e(j)));
L(i,j)=H(i,j);
M(i,j)=-N(i,j);
end
end
end
H,N,M,L
%当i=j时H,N,M,L如下：
fori=1:4
forj=1:5
ifi~=j
H(i,i)=H(i,i)+U(i)*U(j)*(G(i,j)*sin(e(i)-e(j))-B(i,j)*cos(e(i)-e(j)));N(i,i)=N(i,i)-U(i)*U(j)*(G(i,j)*cos(e(i)-e(j))+B(i,j)*sin(e(i)-e(j)));M(i,i)=M(i,i)-U(i)*U(j)*(G(i,j)*cos(e(i)-e(j))+B(i,j)*sin(e(i)-e(j)));L(i,i)=L(i,i)-U(i)*U(j)*(G(i,j)*sin(e(i)-e(j))-B(i,j)*cos(e(i)-e(j)));endend
N(i,i尸N(i,i)-2*(U(i)『2*G(i,i);
L(i,i尸L(i,i)+2*(U(i)『2*B(i,i);
end
J=[H,N;M,L]%J为雅克比矩阵ox=-((inv(J))*fx);
fori=1:4
oe(i)=ox(i);oU(i)=ox(i+4)*U(i);
end
fori=1:4
e(i)=e(i)+oe(i);U(i)=U(i)+oU(i);
end
count=count+1;
end
ox,U,e,count
%求节点注入的净功率
i=5;
forj=1:5
精品文档
精品文档
P(i)=U(i)*UO)*(G(i,j)*cos(e(i)-eQ))+B(i,j)*sin(e(i)-e(j)))+P(i);
Q(i)=U(i)*UQ)*(G(i,j)*sin(e(i)-eQ))-B(i,j)*cos(e(i)-eO)))+Q(i);end
S(5)=P(5)+Q(5)*sqrt(-1);
%求节点注入电流
I=Y*U'
.运行结果
Y值:
Window
■TNew to