%?function:?計(jì)算雅可比矩陣
function?[?J?]?=?Jacobi（?nmUcitaBGPiQi?）
%雅可比矩陣的計(jì)算
%分塊?H?N?K?L
%i!=j時(shí)
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時(shí)由于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

?屬性????????????大小?????日期????時(shí)間???名稱
-----------?---------??----------?-----??----

?????文件????????832??2018-12-20?10:39??牛頓拉夫遜法&PQ分解法計(jì)算潮流\Jacobi.m

?????文件????????265??2018-12-21?10:30??牛頓拉夫遜法&PQ分解法計(jì)算潮流\Polar_Correct.m

?????文件???????4576??2018-12-21?15:42??牛頓拉夫遜法&PQ分解法計(jì)算潮流\PowerFlow.m

?????文件????????329??2018-12-21?10:30??牛頓拉夫遜法&PQ分解法計(jì)算潮流\PQ_Correct.m

?????文件????????558??2018-12-20?10:54??牛頓拉夫遜法&PQ分解法計(jì)算潮流\Unbalanced.m

?????目錄??????????0??2018-12-26?14:15??牛頓拉夫遜法&PQ分解法計(jì)算潮流

-----------?---------??----------?-----??----

?????????????????6560????????????????????6