一、直流潮流计算原理
直流潮流发的特点是用电力系统的交流潮流（有功功率和无功功率）等值的直流电流来代替。甚至只用直流电路的解析法来分析电力系统的有功潮流，而不考虑无功分布对有功的影响。这样一来计算速度加快，但计算的准确度有所降低，本方法适用于对潮流计算准确度要求不高的计算场景。
下面先对直流潮流法的原理进行简单介绍：
上图为直流法的等值图，在上图所示的输电线路中，有功潮流为：
为了快速计算的需要，将上式进行了三项简化：
（1）考虑一般高压电网中线路的电阻远小于电抗，对地电导也可以忽略即 Gii=0 Gij =0
（2）按照标幺值计算时，节点电压与其额定电压相差不大，故有：Ui≈Uj≈1.0；
（3）线路两端的电压相角差（θi-θj）较小，所以有：
这样，上式前两项均为零，只剩第三项
这就相当于线路两端的直流电位分别为θi和θj。线路的直流电阻是Xij。则用矩阵表示为如下式所示。
式中：B0为正常运行时网络的节点电纳矩阵；θ为网络中各节点的电压相位角的向量；P为节点注入的有功功率向量。
二、算例
以IEEE9节点系统为算例，系统参数如下。
系统结构如下
节点参数如下
支路参数如下
三、matlab程序
`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)]);`
。。。。。。。略
详细程序转链接：https://mp.weixin.qq.com/s?__biz=Mzg4MTc1MjE2Mg==&mid=2247483821&idx=1&sn=ff63c6c9d71ee7f8829daa708ad7c3aa&chksm=cf6060a0f817e9b6c3e2c488a4ea5a5c6c8d40be05c517e6657d2d0ec89aac1fce3b49e0f082&token=1687758221&lang=zh_CN#rd
感谢阅读，欢迎关注“电磁matlab”