牛顿拉夫逊法&PQ分解法计算潮流MATLAB程序计算模型。
本程序模型基于MATLAB进行潮流计算，建议先安装matpower插件（MATLAB中非常重要的潮流计算的插件）。
本程序可进行牛拉法和PQ分解法潮流计算的切换，对比潮流计算的结果。很适合新手学习使用。
部分程序实例：
.rtcContent { padding: 30px; } .lineNode {font-size: 10pt; font-family: Menlo, Monaco, Consolas, "Courier New", monospace; font-style: normal; font-weight: normal; }tic %开始计时clear %清除工作空间变量clc %清屏%% ********************数据输入（全局标幺值p.u.）*********************%% 基本参数设置SB=100; %基准容量，单位MVAmode=2; % 1_极坐标下， 2_PQ分解法Tmax=100; %最大迭代次数limit=1.0e-5; %要求精度%% 母线参数：母线编号 节点类型 初始电压 初始相角(度) 负荷有功 负荷无功 基准电压(有名值)% 注意：1.母线填写顺序遵循：1_PQ节点-->2_PV节点-->3_平衡节点% 2.注入功率为+，消耗功率为-Bus=[1 1 1 0 -0.3 -0.18 115;2 1 1 0 -0.55 -0.13 10.5;3 2 1.1 0 0.5 0 115;4 3 1.05 0 0 0 115;];%% 线路参数：I侧母线 J侧母线 阻抗 1/2接地导纳Line=[1 2 0.1+0.4i 0.01528i;1 4 0.12+0.5i 0.0192i;2 4 0.08+0.4i 0.01413i;];%% 变压器参数：I侧母线 J侧母线 阻抗(归算I侧复数) 变比(I->J为k:1)Trans=[3 1 0.3i 1.1];%% 加接地无功补偿装置: 节点 导纳% 静止电容器，设置单组容量和总组数，以投入组数为自变量进行无功优化求解Cap=[2 0];%% **************************参数运算*******************************%% 变压器π型等效导纳参数%Trans_pi：1_I侧母线 2_J侧母线Trans_pi(:,1:2)=Trans(:,1:2);Trans_pi(:,3)=1/(Trans(:,3)./Trans(:,4));Trans_pi(:,4)=1/(Trans(:,3)./(1-Trans(:,4)));Trans_pi(:,5)=1/(Trans(:,3)./(Trans(:,4).^2-Trans(:,4)));%% 计算PQ节点数量m(去除PV节点和平衡节点),其中包含浮游节点，其初始P&Q=0n=size(Bus,1); %总节点数m=n-1;% 去除1个平衡节点数for i=1:nif Bus(i,2)==2 %除去PV节点数m=m-1;endend%% 形成P,Q,U,θ向量,以节点顺序排列for i=1:nP(Bus(i,1))=Bus(i,5);Q(Bus(i,1))=Bus(i,6);U(Bus(i,1))=Bus(i,3);cita(Bus(i,1))=Bus(i,4);
潮流计算结果：
