MATLAB牛顿拉夫逊法算潮流分析（国外英文资料）.docVIP

MATLAB牛顿拉夫逊法算潮流分析 The function of this program is to use Newton Raphson method for power flow calculation % B1 matrix: 1, branch head number, 2, end number; 3, branch impedance; 4, line to earth electricity (or transformer admittance); 5, the ratio of the branch; 6, the first end of the branch is at the K side, and the 1 side is 0; the 1 side is the second; % 7, line / transformer identification (0/1) transformer parameter, when the branch head end is on the K side, the mark is 1, return to the end side, 0 count to the first end side % B2 matrix: 1, the node generator power; 2, the node load power; % 3, PQ node voltage initial value, or the given value of balance node and PV node voltage 4, the node is connected with reactive compensation, parallel capacitance (sense) of the electrical charge % 5, node classification label: 1 is balanced node (should be node 1); 2 is PQ node; 3 is PV node; Clear; Isb=1;%input (enter the balanced bus node number: isb=); Pr=1e-5;%input (\ input error precision: pr=); %--------------------------------------------------- N=10;%input (enter node number: n=); Nl=10;%input (enter the number of branches: nl=); B1=[1 2 3.4+12.8i 1.4e-4i 100; 14 5.1+19.2i 2.1e-4i 100; 23 4.25+16i 1.75e-4i 100; 45 4.25+16i 7e-4i 100; 13 5.1+19.2i 2.1e-4i 100; 64 5.95+22.4i 9.8e-4i 100; 27 1.78+53.89i 0 38.5/231 01; 38 1.49+48.02i 0 11/231 01; 49 1.49+48.02i 0 11/231 01; 510 2.46+70.17i 0 38.5/231 0 1]%input (` enter the matrix formed by the branch parameter]: B1=); B2=[0 225.501; 0022002; 0022002; 0022002; 0022002; 120023103; 0 61.11+37.87i 3502; 0 47.53+29.46i 1002; 0 54.32+33.66i 1002; 0, 40.74+25.25i, 350, 2],%input (enter the matrix of each node parameter: B2=); %------------------------------------------------------------- %n=4;%input (enter node number: n=); %nl=4;%input (enter the number of branches: nl=); %B1=[1 2 4+16i 0100; %1 3 4+16i 0100; %2 3 2+8i 0100; %2 4, 1.49+48.02i 0, 11/110 0, 1],%input ([enter the matrix formed by the branch parameter]: B1