牛顿拉夫逊法PQ分解法计算潮流MATLAB程序计算模型。

牛顿拉夫逊法&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:n    if Bus(i,2)==2 %除去PV节点数        m=m-1;    endend%% 形成P,Q,U,θ向量,以节点顺序排列for i=1:n    P(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);

潮流计算结果：