程序名称##

极坐标下牛顿—拉夫逊潮流计算(matlab版+python版)

程序功能(对象)

适用于任意大小的纯交流电网,支持节点和支路的增删；
适用于接入多个风电、光伏等分布式电源；
子函数包含：节点导纳矩阵计算，雅克比矩阵计算(n-1+m维)；
误差分析：输出形式对标matpower内runpf()函数，误差<10^-3；

程序函数说明

主程序：Test_NL_PowerFlow
输入算例子程序：所有case文件
节点导纳子程序：createYbus文件
雅克比矩阵子程序：Jacobi文件

程序计算步骤及流程图

参数初始化，读取网络参数，在此标幺化

ac_data = case33bw_60;%IEEE33节点为例
ac_bus = ac_data.bus;
ac_branch = ac_data.branch;
ac_gen = ac_data.gen;
ac_dg = ac_data.dg;

形成节点导纳矩阵n*n维

Ybus = createYbus(ac_baseMVA, ac_bus, ac_branch);

%循环列写n-1+m个功率不平衡方程，节点注入功率计算

           for j=1:busNumdetal_ij=(Bus_V(i,3)-Bus_V(j,3))*pi/180;  Pii(i, 2)=Pii(i, 2)+Bus_V(j, 2)*(real(Ybus(i, j))*cos(detal_ij)+imag(Ybus(i, j))*sin(detal_ij));    Qii(i, 2)=Qii(i, 2)+Bus_V(j, 2)*(real(Ybus(i, j))*sin(detal_ij)-imag(Ybus(i, j))*cos(detal_ij));            endPii(i, 2) = Bus_V(i, 2)*Pii(i, 2);Qii(i, 2) = Bus_V(i, 2)*Qii(i, 2);Pi = Pii([ac_pq;ac_pv], :);                               %除去平衡节点且重置节点位置Qi = Qii([ac_pq;ac_pv], :);

计算功率偏差，判断是否收敛

         dP = Pacs - Pi(:,2);dQ = Qacs - Qi(:,2);

形成雅克比矩阵

         [J,H,N,K,L] = Jacobi(V,  Y, ac_pq ,ac_pv,Pi,Qi);

输出节点电压和支路功率
(Bus_V第二列为电压幅值,第三列为电压相角； S_branch第三列为Pij，S_branch第四列为Pji，S_branch第五列为支路ij损耗)

    Ui = Bus_V(:, 2) .* exp(1j * (Bus_V(:, 3)*pi/180));S_branch(i , 3) = Ui(from)*conj(Ui(from))*conj(yi0(from, to))+Ui(from)*(conj(Ui(from))-conj(Ui(to)))*conj(-Ybus(from, to));S_branch(i , 4) =Ui(to)*conj(Ui(to))*conj(yi0(to, from))+Ui(to)*(conj(Ui(to))-conj(Ui(from)))*conj(-Ybus(to, from));S_branch(i , 5) = real(S_branch(i , 3)+ S_branch(i , 4));

计算流程图如下图所示

基于牛顿拉夫逊法电能损耗计算流程图

程序误差分析

接入DG之前，matlab版计算结果与matpower计算结果对比(结果一致)，误差主要存在于收敛精度判断

matlab版18节点接入DG后的计算结果：(Bus_V为标幺值)
python版接入DG之后计算结果：(与MATLAB版结果一致)

程序连接

CSND借鉴版: https://download.csdn.net/download/WConstelltion/12311925
matlab版：https://download.csdn.net/download/WConstelltion/85045068
通用性程序请私信“深入”探讨

【极坐标下牛顿—拉夫逊潮流计算(matlab版+python版)】相关推荐

牛顿-拉夫逊法原理讲解以及python算例实现
文章目录一.基本原理二.原理深入讲解 1 .切线是曲线的线性逼近 2. 牛顿-拉夫逊法 3 .牛顿-拉夫逊法是否总是收敛(总是可以求得足够近似的根)? 三.pythons实例总结一.基本原理 ...
2021-06-16 节点电压为极坐标下的牛顿-拉夫逊法潮流计算学习
以节点电压为极坐标下的牛顿-拉夫逊法潮流计算节点电压用极坐标表示的牛顿-拉夫逊法的潮流计算节点电压用极坐标表示的牛顿-拉夫逊法的潮流计算当节点电压使用极坐标时,节点电压表示为: 节点功率方式 ...
极坐标 c语言,牛顿拉夫逊迭代法极坐标潮流计算C语言程序
<牛顿拉夫逊迭代法极坐标潮流计算C语言程序>由会员分享,可在线阅读,更多相关<牛顿拉夫逊迭代法极坐标潮流计算C语言程序(31页珍藏版)>请在人人文库网上搜索. 1.*利用牛 ...
牛顿法（牛顿拉夫逊）配电网潮流计算matlab程序
牛顿法配电网潮流计算matlab程序传统牛顿-拉夫逊算法,简称牛顿法,是将潮流计算方程组F(X)=0,进行泰勒展开.因泰勒展开有许多高阶项,而高阶项级数部分对计算结果影响很小,当忽略一阶以上部分时, ...
牛顿拉夫逊基波潮流计算通用型程序，runpf函数的替换
牛顿拉夫逊基波潮流计算通用型程序,runpf函数的替换,可提供matlab版和python版 ID:4480672886448715
电力系统潮流【牛顿-拉夫逊法】（4节点、5节点、6节点、9节点）（Matlab代码实现）
目录 1 概述 2 电力系统潮流计算概述 2.1 电力潮流发展进程 2.2牛顿拉夫逊法潮流计算 3 仿真结果 4 Matlab代码及文章讲解
牛顿拉夫逊PQ节点法c语言编程,牛顿拉夫逊介绍(原理、计算方法、程序)..doc
4 3 牛顿-拉夫逊法概述 3.1 牛顿-拉夫逊法基本原理电力系统潮流计算是电力系统分析中的一种最基本的计算,是对复杂电力系统正常和故障条件下稳态运行状态的计算.潮流计算的目标是求取电力系统在给定运 ...
牛顿-拉夫逊简单解释
核心: 牛顿-拉夫逊法在数学上是求解非线性代数方程式的有效方法. 其要点是把非线性方程式的求解过程变成反复地对相应的线性方程式进行求解的过程,即通常所称的逐次线性化过程.是数值计算普遍使用的重要方法, ...
C语言实现牛顿-拉夫逊newton raphson求近似根roots(附完整源码)
实现牛顿-拉夫逊newton raphson求近似根实现以下几个相关的接口实现牛顿-拉夫逊newton raphson求近似根的完整源码(定义,实现,main函数测试) 实现以下几个相关的接口 d ...

【极坐标下牛顿—拉夫逊潮流计算(matlab版+python版)】

程序名称##

程序功能(对象)

程序函数说明

程序计算步骤及流程图

程序误差分析

程序连接

【极坐标下牛顿—拉夫逊潮流计算(matlab版+python版)】相关推荐

最新文章

热门文章