Matlab程序_牛顿拉夫逊法&PQ;分解法计算潮流

1.Matlab程序，附带超详细批注 2.分别用牛顿拉夫逊法和PQ分解法求解潮流 3.算例采用华科大电分下册例11-5作为验证，计算结果完全一致

% function: 计算雅可比矩阵
function [ J ] = Jacobi（ nmUcitaBGPiQi ）
%雅可比矩阵的计算
%分块 H N K L
%i!=j时
for i=1:n-1
	for j=1:n-1
		H（ij）=-U（i）*U（j）*（G（ij）*sin（cita（i）-cita（j））-B（ij）*cos（cita（i）-cita（j）））;
	end
end
for i=1:n-1
	for j=1:m
		N（ij）=-U（i）*U（j）*（G（ij）*cos（cita（i）-cita（j））+B（ij）*sin（cita（i）-cita（j）））;
	end
end
for i=1:m
	for j=1:n-1
		K（ij）=U（i）*U（j）*（G（ij）*cos（cita（i）-cita（j））+B（ij）*sin（cita（i）-cita（j）））;
	end
end
for i=1:m
	for j=1:m
		L（ij）=-U（i）*U（j）*（G（ij）*sin（cita（i）-cita（j））-B（ij）*cos（cita（i）-cita（j）））;
	end
end
%i==j时由于m<=n-1
for i=1:n-1
	H（ii）=U（i）^2*B（ii）+Qi（i）;
end
for i=1:m
	N（ii）=-U（i）^2*G（ii）-Pi（i）;
end
for i=1:m
	K（ii）=U（i）^2*G（ii）-Pi（i）;
end
for i=1:m
	L（ii）=U（i）^2*B（ii）-Qi（i）;
end

%合成雅可比矩阵
J=[H N;K L];

end

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----

     文件        832  2018-12-20 10:39  牛顿拉夫逊法&PQ分解法计算潮流\Jacobi.m

     文件        265  2018-12-21 10:30  牛顿拉夫逊法&PQ分解法计算潮流\Polar_Correct.m

     文件       4576  2018-12-21 15:42  牛顿拉夫逊法&PQ分解法计算潮流\PowerFlow.m

     文件        329  2018-12-21 10:30  牛顿拉夫逊法&PQ分解法计算潮流\PQ_Correct.m

     文件        558  2018-12-20 10:54  牛顿拉夫逊法&PQ分解法计算潮流\Unbalanced.m

     目录          0  2018-12-26 14:15  牛顿拉夫逊法&PQ分解法计算潮流

----------- ---------  ---------- -----  ----

                 6560                    6