1、 1 2 0.0922 0.0470; 2 3 0.4930 0.2511; 3 4 0.3660 0.1864; 4 5 0.3811 0.1941; 5 6 0.8190 0.7070; 6 7 0.1872 0.6188; 7 8 0.7114 0.2351; 8 9 1.0300 0.7400; 9 10 1.0440 0.7400; 10 11 0.1966 0.0650; 11 12 0.3744 0.1238; 12 13 1.4680 1.1550; 13 14 0.5416 0.7129; 14 15 0.5910 0.5260; 15 16 0.7463 0.5450; 1
2、6 17 1.2890 1.7210; 17 18 0.7320 0.5740; 2 19 0.1640 0.1565; 19 20 1.5042 1.3554; 20 21 0.4095 0.4784; 21 22 0.7089 0.9373; 3 23 0.4512 0.3083; 23 24 0.8980 0.7091; 24 25 0.8960 0.7011; 6 26 0.2030 0.1034; 26 27 0.2842 0.1447; 27 28 1.0590 0.9337; 28 29 0.8042 0.7006; 29 30 0.5075 0.2585; 30 31 0.97
3、44 0.9630; 31 32 0.3105 0.3619; 32 33 0.3410 0.5302; % 支路，阻抗NodeData = 2 100.00 60.00; 3 90.00 40.00; 4 120.00 80.00; 5 60.00 30.00; 6 60.00 20.00; 7 200.00 100.00; 8 200.00 100.00; 9 60.00 20.00; 10 60.00 20.00; 11 45.00 30.00; 12 60.00 35.00; 13 60.00 35.00; 14 120.00 80.00; 15 60.00 10.00; 16 60.
4、00 20.00; 17 60.00 20.00; 18 90.00 40.00; 19 90.00 40.00; 20 90.00 40.00; 21 90.00 40.00; 22 90.00 40.00; 23 90.00 50.00; 24 420.00 200.00; 25 420.00 200.00; 26 60.00 25.00; 27 60.00 25.00; 28 60.00 20.00; 29 120.00 70.00; 30 200.00 600.00; 31 150.00 70.00; 32 210.00 100.00; 33 60.00 40.00; % 节点，负荷U
5、B = 12.66; % 电压基准 kVSB = 10; % 功率基准 MVAZB = UB2/SB; % 阻抗基准 ohmBranchData（:,3,4） = BranchData（:,3,4） / ZB; % 阻抗标幺化NodeData（:,2,3） = NodeData（:,2,3） / SB / 1000;% 功率标幺化NN = 33; % 节点数A0 = zeros（NN）;for n = 1:NN-1 A0（BranchData（n,1）,BranchData（n,2） = 1;end % 形成 A0AssociatedMatrix=0;for n=2: AssociatedMa
6、trix（n,n）=1; temp=BranchData（n-1,1）; AssociatedMatrix（n,1:n-1）=AssociatedMatrix（temp,1:n-1）;endA0T = A0; % 形成 A0 的转置S = 0;-NodeData（:,2） - i*NodeData（:,3）; % 形成 SZL = 0;,3） + i*BranchData（:,4）; % 形成 ZLV = ones（NN,1）;V（1） = 1; % 各个节点电压赋初值IL（NN,1） = -conj（S（NN） / V（NN）; % 最末支路电流赋初值Delta = 1; % 收敛判据赋初值
7、TempV = V; % 赋初值，用于记忆上次迭代结果while Delta 1e-8 countnum=countnum+1; IN = conj（S ./ V）; % 节点注入电流 for n = 1: IL（NN-n） = A0（NN-n,NN-n+1:end） * IL（NN-n+1:end） - IN（NN-n）; end % 电流回代过程 for n = 2:NN V（n） = A0T（n,1:n-1） * V（1:n-1） - ZL（n） * IL（n）; end % 电压前推过程 Delta = max（abs（V-TempV）; % 更新收敛判据 TempV = V; % 记
8、忆迭代结果Vangle（:,1）=abs（V）;,2）=angle（V）/3.1415*180;for i=1: st=BranchData（i,1）; en=BranchData（i,2）; Sij（i,1）=V（st）*conj（V（st）-V（en）/ZL（i+1）; Sji（i,1）=V（en）*conj（V（en）-V（st）/ZL（i+1）;End实验五 输电网潮流计算实验 本实验针对高压输电网络具有多环、多电源的结构特性，要求根据功率平衡方程的泰勒展开式近似的求取系统的状态变量，使系统潮流计算能够由计算机自行完成，即根据已知的电力网的数学模型（节点导纳矩阵）及各节点参数，由计算程
9、序运行完成该电力系统的潮流计算。通过实验教学加深学生对复杂电力系统潮流计算计算方法的理解，学会运用电力系统的数学模型，掌握潮流计算的过程及其特点，熟悉各种常用应用软件，熟悉硬件设备的使用方法，加强编制调试计算机程序的能力，提高工程计算的能力，学习如何将理论知识和实际工程问题结合起来。 编制调试电力系统潮流计算的计算机程序。程序要求根据已知的电网的数学模型（节点导纳矩阵）及各节点参数，完成该电力系统的潮流计算，要求计算出节点电压、功率等参数。1.采用牛顿拉夫逊法编制程序编制潮流程序。2.先编制好的电力系统潮流计算的计算程序原代码由自备移动存储设备导入计算机。3.应的编程环境下对程序进行组织调试。
10、4.应用计算例题验证程序的计算效果。5.对调试正确的计算程序进行存储、打印。6.完成本次实验的实验报告。四、 实验程序%牛顿拉夫逊法进行潮流计算% B1矩阵：1、支路首端号；2、末端号；3、支路阻抗；4、支路对地电纳% 5、支路的变比；6、支路首端处于K侧为1，1侧为0% B2矩阵：1、该节点发电机功率；2、该节点负荷功率；3、节点电压初始值% 4、PV节点电压V的给定值；5、节点所接的无功补偿设备的容量% 6、节点分类标号：1为平衡节点；2为PQ节点；% 3为PV节点；clear;clc;n=4;%input（请输入节点数：n=）;nl=4;请输入支路数：nl=isb=1;请输入平衡母线节点
11、号：isb=pr=0.00001;请输入误差精度：pr=B1=2 4 0.1+0.4i 0.01528*2i 1 0 ; 3 4 0.3i 0 1./1.1 0 ; 1 4 0.12+0.5i 0.01920*2i 1 0; 1 2 0.08+0.4i 0.01413*2i 1 0 ; 请输入各节点参数形成的矩阵：B2=B2=0 0 1.05 1.05 0 1; 0 0.3+0.18i 1 0 0 2; 0.5 0 1.1 1.1 0 3; 0 0.55+0.13i 1 0 0 2;Y=zeros（n）;e=zeros（1,n）;f=zeros（1,n）;V=zeros（1,n）;sida=z
12、eros（1,n）;S1=zeros（nl）;% % %-nl%支路数 if B1（i,6）=0%左节点处于1侧 p=B1（i,1）;q=B1（i,2）; else %左节点处于K侧 p=B1（i,2）;q=B1（i,1）; end Y（p,q）=Y（p,q）-1./（B1（i,3）*B1（i,5）;%非对角元 Y（q,p）=Y（p,q）; %非对角元 Y（q,q）=Y（q,q）+1./（B1（i,3）*B1（i,5）2）+B1（i,4）./2;%对角元K侧 Y（p,p）=Y（p,p）+1./B1（i,3）+B1（i,4）./2; %对角元1侧 %求导纳矩阵disp（导纳矩阵 Y=disp（Y
13、）%-G=real（Y）;B=imag（Y）;%分解出导纳阵的实部和虚部 n%给定各节点初始电压的实部和虚部 e（i）=real（B2（i,3）; f（i）=imag（B2（i,3）; V（i）=B2（i,4）;%PV节点电压给定模值 n%给定各节点注入功率 S（i）=B2（i,1）-B2（i,2）; %i节点注入功率SG-SL B（i,i）=B（i,i）+B2（i,5）;%i节点无功补偿量 %-P=real（S）;Q=imag（S）; %分解出各节点注入的有功和无功功率ICT1=0;IT2=1;N0=2*n;N=N0+1;a=0; %迭代次数ICT1、a；不满足收敛要求的节点数IT2whil
14、e IT2=0 % N0=2*n 雅可比矩阵的阶数；N=N0+1扩展列 IT2=0;a=a+1; for i=1:n if i=isb%非平衡节点 C（i）=0;D（i）=0; for j1=1: C（i）=C（i）+G（i,j1）*e（j1）-B（i,j1）*f（j1）;%（Gij*ej-Bij*fj） D（i）=D（i）+G（i,j1）*f（j1）+B（i,j1）*e（j1）;%（Gij*fj+Bij*ej） end P1=C（i）*e（i）+f（i）*D（i）;%节点功率P计算ei（Gij*ej-Bij*fj）+fi（Gij*fj+Bij*ej） Q1=C（i）*f（i）-e（i）*D（
15、i）;%节点功率Q计算fi（Gij*ej-Bij*fj）-ei（Gij*fj+Bij*ej）%求i节点有功和无功功率P,Q的计算值 V2=e（i）2+f（i）2;%电压模平方 %- 以下针对非PV节点来求取功率差及Jacobi矩阵元素 - if B2（i,6）=3%非PV节点 DP=P（i）-P1;%节点有功功率差 DQ=Q（i）-Q1; %节点无功功率差 %- 以上为除平衡节点外其它节点的功率计算 -%-求取Jacobi矩阵 - for j1=1: if j1=isb&j1=i%非平衡节点&非对角元 X1=-G（i,j1）*e（i）-B（i,j1）*f（i）;% dP/de=-dQ/df X
16、2=B（i,j1）*e（i）-G（i,j1）*f（i）;% dP/df=dQ/de X3=X2; % X2=dp/df X3=dQ/de X4=-X1; % X1=dP/de X4=dQ/df p=2*i-1;q=2*j1-1; J（p,q）=X3;J（p,N）=DQ;m=p+1; % X3=dQ/de J（p,N）=DQ节点无功功率差 J（m,q）=X1;J（m,N）=DP;q=q+1; % X1=dP/de J（m,N）=DP节点有功功率差 J（p,q）=X4;J（m,q）=X2; % X4=dQ/df X2=dp/df elseif j1=i&j1=isb%非平衡节点&对角元 X1=-C
17、（i）-G（i,i）*e（i）-B（i,i）*f（i）;% dP/de X2=-D（i）+B（i,i）*e（i）-G（i,i）*f（i）;% dP/df X3=D（i）+B（i,i）*e（i）-G（i,i）*f（i）; % dQ/de X4=-C（i）+G（i,i）*e（i）+B（i,i）*f（i）;% dQ/df p=2*i-1;J（p,q）=X3;%扩展列Q m=p+1; J（m,q）=X1;J（p,q）=X4;%扩展列P J（m,q）=X2; end end else%-下面是针对PV节点来求取Jacobi矩阵的元素 - DP=P（i）-P1;% PV节点有功误差 DV=V（i）2-V2;% PV节点电压误差 % dP/de % dP/df X5=0;X6=0;J（p,q）=X5;J（p,N）=DV; % PV节点电压误差J（p,q）=X6; % PV节点有功误差 J（m,q
