基于直角坐标系的牛顿法潮流计算

基于直角坐标系的牛顿法潮流计算 linedata = [ 1 4 0.1 0.4 0.01528 1; 1 3 0 0.3 0 1.1; 1 2 0.12 0.5 0.01920 1; 4 2 0.08 0.40 0.01413 1]; busdata = [ 1 0 -0.30 -0.18 1 0; 4 0 -0.55 -0.13 1 0; 3 2 0.5 0 1.10 0; 2 1 0 0 1.05 0;];

function PowerFlowNRRec（linedatabusdataerror）
	%节点标号无特殊要求
	clc
	if nargin < 1
		% linedata = [
		% 1	2	0.1		0.4		0.01528 1;
		% 1	3	0		0.3		0 		1.1;
		% 1	4	0.12	0.5		0.01920 1;
		% 2	4	0.08	0.40	0.01413 1];
		% busdata = [
		% 1 0	-0.30 	-0.18 	1    	0;
		% 2 0	-0.55 	-0.13 	1 		0;
		% 3 2	0.5   	0		1.10	0;
		% 4 1	0		0		1.05	0;];
		
		linedata = [
			1	4	0.1		0.4		0.01528 1;
			1	3	0		0.3		0 		1.1;
			1	2	0.12	0.5		0.01920 1;
			4	2	0.08	0.40	0.01413 1];
		busdata = [
			1 0	-0.30 	-0.18 	1    	0;
			4 0	-0.55 	-0.13 	1 		0;
			3 2	0.5   	0		1.10	0;
			2 1	0		0		1.05	0;];
	end
	if nargin < 3
		error = 10^-5;
	end
	clc


	%优化节点排序
	pqbus = find（busdata（:2）==0）;
	pvbus = find（busdata（:2）==2）;
	basebus = find（busdata（:2）==1）;

	bussort = [pqbus‘ pvbu