电力系统直流潮流计算研究【IEEE9节点】（Matlab代码实现）

💥💥💞💞欢迎来到本博客❤️❤️💥💥

🏆博主优势：🌞🌞🌞博客内容尽量做到思维缜密，逻辑清晰，为了方便读者。

⛳️座右铭：行百里者，半于九十。

📋📋📋本文目录如下：🎁🎁🎁

目录

💥1 概述

📚2 运行结果

2.1 IEEE9节点

🎉3 参考文献

🌈4 Matlab代码实现

💥1 概述

直流潮流发的特点是用电力系统的交流潮流（有功功率和无功功率）用等值的直流电流来代替。甚至只用直流电路的解析法来分析电力系统的有功潮流，而根不如你不考虑无功分布对有功的影响。这样一来计算速度加快，但计算的准确度有所降低，本方法适用于对潮流计算准确度要求不高的计算场景。
下面先对直流潮流法的原理进行简单介绍：

上图为直流法的等值图，在上图所示的输电线路中，有功潮流为：

在这里插入图片描述

为了快速计算的需要，将上式进行了三项简化：
（1）考虑一般高压电网中线路的电阻远小于电抗，对地电导也可以忽略即 Gii=0 Gij =0
（2）按照标幺值计算时，节点电压与其额定电压相差不大，故有：Ui≈Uj≈1.0；
（3）线路两端的电压相角差（θi-θj）较小，所以有：

在这里插入图片描述

这样，上式前两项均为零，只剩第三项

在这里插入图片描述

这就相当于线路两端的直流电位分别为θi和θj。线路的直流电阻是Xij。则用矩阵表示为如下式所示。

式中：B0为正常运行时网络的节点电纳矩阵；
θ为网络中各节点的电压相位角的向量；
P为节点注入的有功功率向量；

📚2 运行结果

2.1 IEEE9节点

clc
close all
clear 
%% 算例
mpc = case9;

%% 潮流计算
[theta1,P_branch,M,Z,slackbus] = DCpowerflow(mpc);

%% 输出结果
disp('=============================');
disp('支路潮流矩阵')
disp('=============================');
disp('')
disp([num2str(P_branch)]);



disp('=============================');
disp('节点相位矩阵')
disp('=============================');
disp('')
disp([num2str(theta1)]);

clc
close all
clear
%% 算例
mpc = case9;

%% 潮流计算
[theta1,P_branch,M,Z,slackbus] = DCpowerflow(mpc);

%% 输出结果
disp('=============================');
disp('支路潮流矩阵')
disp('=============================');
disp('')
disp([num2str(P_branch)]);

disp('=============================');
disp('节点相位矩阵')
disp('=============================');
disp('')
disp([num2str(theta1)]);