基于MATLAB的潮流计算Word格式文档下载.docx

1、a=zeros（n,1）;Ps=zeros（n,1）;Qs=zeros（n,1）;P=zeros（n,1）;Q=zeros（n,1）;p=zeros（n-1,1）;q=zeros（m,1）;aa=zeros（n-1,1）;u=zeros（m,1）;k=0;I=sqrt（-1）;%由于后面用到 ij 作为循环变量故重定义虚数单位Y=6.311-I*20.3722 -2.7523+I*9.1743 -3.5587+I*11.3879 0 0;-3.5587+I*11.3879 8.5587-I*30.9993 -5+I*15 I*4.9889 0;-2.7523+I*9.1743 -5+I*15

2、7.7523-I*28.7557 0 I*4.9889;0 I*4.9889 0 -I*5.2493 0;0 0 I*4.9889 0 -I*5.2493;%初值的设定%U=1;1;1;a=0;0;0;Ps=-0.8055;-0.18;0.5;Qs=-0.532;-0.12;3%求失配功率P（k）和Q（k）%for i=1:n-1s=0;for j=1:ns=s+U（j）*（real（Y（i,j）*cos（a（i）-a（j）+imag（Y（i,j）*sin（a（i）-a（j）;endP（i）=U（i）*s;s=s+U（j）*（real（Y（i,j）*sin（a（i）-a（j）-imag（Y（i

3、,j）*cos（a（i）-a（j）;Q（i）=U（i）*s;detp（i）=Ps（i）-P（i）;mdetq（i）=Qs（i）-Q（i）;% 求解雅可比矩阵 %while （max（detp）=0.0001 | max（detq）=0.0001） %设定精确度H=zeros（n-1,n-1）;N=zeros（n-1,m）;M=zeros（m,n-1）;L=zeros（m,m）;if i=j4H（i,j）=U（i）2*imag（Y（i,j）+Q（i）;elseH（i,j）=-U（i）*U（j）*（real（Y（i,j）*sin（a（i）-a（j）-imag（Y（i,j）*cos（a（i）-a（j

4、）;N（i,j）=-U（i）2*real（Y（i,j）-P（i）;N（i,j）=-U（i）*U（j）*（real（Y（i,j）*cos（a（i）-a（j）+imag（Y（i,j）*sin（a（i）-a（j）;M（i,j）=U（i）2*real（Y（i,j）-P（i）;M（i,j）=U（i）*U（j）*（real（Y（i,j）*cos（a（i）-a（j）+imag（Y（i,j）*sin（a（i）-a（j）;L（i,j）=U（i）2*imag（Y（i,j）-Q（i）;L（i,j）=-U（i）*U（j）*（real（Y（i,j）*sin（a（i）-a（j）-imag（Y（i,j）*cos（a（i）-

5、a（j）;5% 求解修正方程得到修正量 aa,u %J=-inv（H,N;M,L）*detp;detq;aa=J（1:n-1,:u=J（n:n-1+m,:n-1 %计算 a（节点电压相角）,U（节点电压大小）a（i）=a（i）+aa（i）;U（i）=U（i）+U（i）*u（i）;detp %迭代过程中失配功率变化情况k=k+1;n-1 %计算功率误差 p,q6A=a*180/pi; %对相角进行转换由弧度制装换为角度display（最终的计算结果如下: AU各节点的电压为:fprintf（节点%d: %d（%d）nn,i,U（i）,A（i）;失配功率为:detpdetq迭代的次数为: %dn,

6、k） %输出迭代次数2运行结果 fengjian请输入 PQ 节点个数:ans =-0.8055-0.18000.5000-0.34200.25750.40757ans =0.03470.0022-0.0094-0.0178-0.0138-0.0407-0.04471.0e-003 *0.5087-0.1291-0.2188-0.3127-0.4229-0.4416A =-7.2049-5.6096-5.3697-0.0056U = 1.00271.02631.03171.00008节点 1: 1.002729e+000（-7.204950e+000）节点 2: 1.026321e+000（-

7、5.609609e+000）节点 3: 1.031663e+000（-5.369712e+000）节点 4: 1（-5.622073e-003）节点 5: 1（0） detp =1.0e-006 *0.1069 -0.0529 -0.0548 -0.0083detq =1.0e-007 *-0.7838 -0.5117 -0.5255 3二 PQ 解耦迭代1 程序代码n=5;m=3;% 各参数初始化 % %导纳矩阵 %电压矢量 %相角矢量9%节点导纳矩阵的生成%初值的设定%求失配功率P（k）和Q（k）%10%生成 B和 B矩阵%Bp=zeros（n-1,n-1）;Bpp=zeros（m,m）;

8、n-1;Bp（i,j）=imag（Y（i,j）;end;m;Bpp（i,j）=imag（Y（i,j）;%aa=（-1）*inv（Bp）*detp（1）/U（1）;detp（2）/U（2）;detp（3）/U（3）;detp（4）/U（4）;u=（-1）*inv（Bpp）*detq（1）/U（1）;detq（2）/U（2）;detq（3）/U（3）;% 求 解 修 正 方 程 得 到 修 正 量aa,u %U（i）=U（i）+u（i）;11% detp%（记录迭代次数）生成的节点导纳矩阵为Y12算得矩阵 B为Bp算得矩阵 B”为Bpp2运行结果 fengjian1A =-7.2049-5.609

9、7-0.0057U =1.0027生成的节点导纳矩阵为Y =6.3110 -20.3722i -2.7523 + 9.1743i -3.5587 +11.3879i 0 013-3.5587 +11.3879i 8.5587 -30.9993i -5.0000 +15.0000i 0 + 4.9889i 0-2.7523 + 9.1743i -5.0000 +15.0000i 7.7523 -28.7557i 0 0 + 4.9889i0 0 + 4.9889i 0 0 - 5.2493i 00 0 0 + 4.9889i 0 0 - 5.2493i算得矩阵 B为Bp =-20.3722 9.

10、1743 11.3879 011.3879 -30.9993 15.0000 4.98899.1743 15.0000 -28.7557 00 4.9889 0 -5.2493算得矩阵 B”为Bpp =-20.3722 9.1743 11.387911.3879 -30.9993 15.00009.1743 15.0000 -28.7557 1.002732e+000（-7.204948e+000） 1.026322e+000（-5.609674e+000） 1.031663e+000（-5.369742e+000） 1（-5.676340e-003）detp =1.0e-004 *-0.2747 0.2761 0.0562 -0.0119-0.3272 0.0311 0.2931 814