MATLAB牛顿Raphson能量流分析方法.docVIP

MATLAB牛顿Raphson能量流分析方法 %本程序的功能是用牛顿——拉夫逊法进行潮流计算 % B1矩阵：1、支路首端号； 2、末端号； 3、支路阻抗； 4、线路对地电纳 (或变压器导纳)； % 5、支路的变比； 6、支路首端处于K侧为1，1侧为0； % 7、线路/变压器标识（0/1）变压器参数当支路首端处于K侧标识为1时归算至末端侧，0归算至首端侧 % B2矩阵：1、该节点发电机功率； 2、该节点负荷功率； % 3、PQ节点电压初始值，或平衡节点及PV节点电压的给定值 % 4、节点所接无功补偿并联电容（感）的电纳 % 5、节点分类标号：1为平衡节点（应为1号节点）；2为PQ节点；3为PV节点； clear; isb=1; %input(请输入平衡母线节点号：isb=); pr=1e-5; %input(请输入误差精度：pr=); %--------------------------------------------------- n=10;%input(请输入节点数：n=); nl=10;%input(请输入支路数：nl=); B1=[1 2 3.4+12.8i 1.4e-4i 1 0 0; 1 4 5.1+19.2i 2.1e-4i 1 0 0; 2 3 4.25+16