《潮流分析matlab》由会员分享，可在线阅读，更多相关《潮流分析matlab（13页珍藏版）》请在金锄头文库上搜索。
1、潮流分析潮流分析 matlabmatlab请输入节点数：n=5请输入支路数：n1=5请输入平衡母线节点号：isb=5请输入误差精度：pr=0.0001请输入由支路参数形成的矩阵：B1=1 2 0.025+0.08i 0.07i 1 0;1 3 0.03+0.1i 0.09 1 0;2 3 0.02+0.06i 0.05i 1 0; 4 2 0.1905i 0 1.0522 1;5 3 0.1905i 0 1.0522 1请输入各节点参数形成的矩阵：B2=0 -0.8055-0.5320i 1 0 0 2;0 -0.18-0.12i 1 0 0 2;0 0 1 0 0 2;0 0.5 1.052
2、2 1.0522 0 3;0 0 1.0522 1.0522 0 1请输入由节点参数形成的矩阵：X=1 0;2 0;3 0;4 0;5 0;导纳矩阵 Y=6.3560 -20.5272i -3.5587 +11.3879i -2.7523 + 9.1743i 0 0 -3.5587 +11.3879i 8.5587 -31.5772i -5.0000 +15.0000i 0 + 4.9889i 0 -2.7523 + 9.1743i -5.0000 +15.0000i 7.7973 -29.3987i 0 0 + 4.9889i0 0 + 4.9889i 0 0 - 4.7414i 0 0 0
3、 0 + 4.9889i 0 0 - 4.7414i雅可比矩阵为：-20.5622 11.3879 9.1743 0 -6.4010 3.5587 2.752311.3879 -31.6372 15.0000 5.2493 3.5587 -8.5587 5.00009.1743 15.0000 -29.4237 0 2.7523 5.0000 -7.84230 5.2493 0 -5.2493 0 0 06.3110 -3.5587 -2.7523 0 -20.4922 11.3879 9.1743-3.5587 8.5587 -5.0000 0 11.3879 -31.5172 15.000
4、0-2.7523 -5.0000 7.7523 0 9.1743 15.0000 -29.3737雅可比矩阵为：-24.7747 13.7501 11.0245 0 -8.7903 3.9345 2.983113.5438 -36.5472 17.3684 5.6351 4.5948 -10.1178 5.731310.8460 17.3335 -33.7956 0 3.5781 5.8359 -8.93540 5.6351 0 -5.6351 0 0.5384 06.9176 -3.9345 -2.9831 0 -25.9551 13.7501 11.0245-4.5948 9.7878 -
5、5.7313 0.5384 13.5438 -36.8940 17.3684-3.5781 -5.8359 9.0080 0 10.8460 17.3335 -33.8576雅可比矩阵为：-24.1898 13.4226 10.7672 0 -8.4629 3.8896 2.957513.2489 -35.8499 17.0208 5.5803 4.4453 -9.9289 5.625310.6171 16.9918 -33.1727 0 3.4578 5.7123 -8.79750 5.5803 0 -5.5803 0 0.5003 06.8471 -3.8896 -2.9575 0 -25
6、.2549 13.4226 10.7672-4.4453 9.5702 -5.6253 0.5003 13.2489 -36.0917 17.0208-3.4578 -5.7123 8.7994 0 10.6171 16.9918 -33.1741雅可比矩阵为：-24.1809 13.4176 10.7633 0 -8.4576 3.8892 2.957413.2446 -35.8395 17.0155 5.5794 4.4428 -9.9265 5.623710.6138 16.9866 -33.1635 0 3.4557 5.7103 -8.79580 5.5794 0 -5.5794 0
7、 0.5000 06.8466 -3.8892 -2.9574 0 -25.2449 13.4176 10.7633-4.4428 9.5665 -5.6237 0.5000 13.2446 -36.0795 17.0155-3.4557 -5.7103 8.7958 0 10.6138 16.9866 -33.1635-各节点电压 U 为（节点从小到大排列）:1.0928 + 0.0984i 1.0646 + 0.0736i 1.0598 + 0.0705i 1.0520 - 0.0214i 1.0522 -各节点电压相角为（节点从小到大排列）:5.1431 3.9536 3.8076 -1
8、.1673 0-按公式计算全部线路功率，结果如下:1.0000 2.0000 0.3952 + 0.2502i1.0000 3.0000 0.4103 + 0.2818i2.0000 3.0000 0.0702 + 0.0377i2.0000 4.0000 0.5000 + 0.3984i3.0000 5.0000 0.3702 + 0.3585i2.0000 1.0000 -0.3902 - 0.3162i3.0000 1.0000 -0.3002 - 0.2646i3.0000 2.0000 -0.0700 - 0.0939i4.0000 2.0000 -0.5000 - 0.3301i5
9、.0000 3.0000 -0.3702 - 0.3137i各条支路的功率损耗 DS 为(顺序同您输入 B1 时一样):0.0050 - 0.0659i0.1101 + 0.0171i0.0002 - 0.0562i0.0000 + 0.0684i-0.0000 + 0.0449i平衡节点的功率:0.8055 + 0.5320i 0.1800 + 0.1200i 0 -0.5000 0 %本程序的功能是用牛顿-拉夫逊法进行潮流计算 直角坐标n=input(请输入节点数：n=);n1=input(请输入支路数：n1=);isb=input(请输入平衡母线节点号：isb=);pr=input(请输
10、入误差精度：pr=);%B1 数据说明：第一列和第二列是节点编号，节点好由小到大编写%对于含有变压器的的支路 ，第一列为低压侧节点编号，第二列为高压侧节点编号%第三列为支路的串联阻抗参数%第四列为支路的对地导纳参数%第五列为含变压器支路的变压器变比%第六列为变压器是否含有变压器B1=input(请输入由支路参数形成的矩阵：B1=);%B2 数据说明：第一列为节点注入发电功率参数；第二列为节点负荷功率参数%第三列为节点电压参数，第六列为节点类型参数B2=input(请输入各节点参数形成的矩阵：B2=);X=input(请输入由节点参数形成的矩阵：X=);Y=zeros(n); e=zeros(1
11、,n);f=zeros(1,n);V=zeros (1,n);O=zeros(1,n);S1=zeros(n1);for i=1:nif X(i,2)=0p=X(i,1);Y(p,p)=1./X(i,2);endendfor i=1:n1if B1(i,6)=0p=B1(i,1);q=B1(i,2);elsep=B1(i,2);q=B1(i,1);endY(p,q)=Y(p,q)-1./(B1(i,3)*B1(i,5);Y(p,q)=Y(p,q);Y(p,q)=Y(q,q)+1./(B1(i,3)*B1(i,5)2)+B1(i,4)./2;Y(p,p)=Y(p,p)+1./B1(i,3)+B1(i,4)./2;end G=real(Y);B=imag(Y);for i=1:ne(i)=real(B2(i,3);f(i)=imag(B2(i,3);V(i)=B2(i,4);endfor
