潮流计算程序(极坐标）——python
本人电气工程研一萌新，初开学学习python，拿潮流计算练手写了个潮流计算程序，与大家一起交流分享进步。）
程序算例是IEEE30节点，网上可以找IEEE标准测试系统数据
我就不摆各种公式了（相信来找潮流计算程序的都是带任务来的，公式啥的都懂，教材也有）
直接进入正题
用到的库是numpy和pandas，一个是矩阵运算用，一个主要用来读取数据
首先是写入数据，整理节点原始数据
利用pandas读excel数据，再利用numpy形成矩阵，真的超快（我刚开始写的时候傻傻的在pycharm里一个个数据敲打...）ps：原始数据格式就是代码里备注的
然后就是整理线路数据
然后就是节点导纳矩阵了，这里还有一个关键点是，python不同matlab，它的索引是从0开始，这是特别需要注意的！
接下来就是要求有功和无功不平衡量。
因为我懒得重新编号，在python中对照新旧节点数据纠错属实不方便，所以在后面的不平衡量矩阵这里我放置里两列，然后筛选出需要的节点，比如有功不平衡量需要把平衡节点拿掉，第一行放置除平衡节点的节点号，第二行是它们的不平衡量。
然后雅克比矩阵也是需要做出相应的改变，需要先找出相应节点（笔者认为这里需要仔细理解雅克比矩阵中各元素的意义，不到单纯的理解对角和非对角哦，那只是特殊情况，因为书上的节点是按照PQPU平衡节点排的）
（篇幅原因，K矩阵和L矩阵放不下了，隐藏起来了，其实都是模仿着H和N写的）
得到雅克比矩阵后，就是要求除幅值相角修正量了，就是简单地雅克比矩阵求逆然后乘不平衡量矩阵，然后把各个数据分配。
最后就是迭代输出结果了
最后贴出结果与标准测试系统结果对比，节点电压都一致，就是相角还是会有0.000x的误差，应该是差不多的了。
结果应该是目前没办法做到完全一致，应该是处于合理误差内（自我安慰），笔者做潮流计算的目的本来也是为了锻炼自己的编程能力嘛，编程就是在一次次报错中进步的嘛。
这是我第一次发文分享，如果对各位读者有帮助的希望各位看官点点赞，后续写好了其他的电力系统编程项目（在写机组组合优化调度ing）再跟大家交流分享。
同时，如果哪里说的不好不对希望大家多多指正，谢谢。
最后最后，其实这适用于所有节点的，只需输入4个excel表格原始数据
原创作品哦~
潮流计算程序(极坐标）——python相关推荐
- 潮流课设matlab编程,10kV配电网络简单潮流计算程序——电力系统分析课设2018（matlab...
10kV配电网络简单潮流计算程序--电力系统分析课设2018(matlab 10kV配电网络简单潮流计算程序--电力系统分析课设2018(matlab) 10kV配电网络简单潮流计算程序--电力系统分 ...
- 电力系统潮流C语言,电力系统中的潮流计算程序
<电力系统中的潮流计算程序>由会员分享,可在线阅读,更多相关<电力系统中的潮流计算程序(7页珍藏版)>请在人人文库网上搜索. 1.B2=input(请输入各节点参数形成的矩阵: ...
- 电力系统潮流计算程序实现C语言版（动态节点+稀疏技术）
概述 电力系统潮流计算程序采用的是牛顿-拉夫逊法(直角坐标) 潮流计算程序的系统节点数是可以动态变化的,根据节点大小分配储存空间 为了减少储存空间和加快计算速度,根据电力系统节点导纳矩阵稀疏的特点,计 ...
- 环形网络潮流计算matlab，利用matlab编程计算任意环形网络牛拉法潮流计算程序
环形网络潮流计算matlab,利用matlab编程计算任意环形网络牛拉法潮流计算程序,程序通用性强,通过修改参数可以得到任意节点和网络的环形网络牛拉法潮流计算. YID:696064261479453 ...
- 环形网络潮流计算matlab 利用matlab编程计算任意环形网络牛拉法潮流计算程序
环形网络潮流计算matlab 利用matlab编程计算任意环形网络牛拉法潮流计算程序,程序通用性强,通过修改参数可以得到任意节点和网络的环形网络牛拉法潮流计算. YID:856064261479453 ...
- 环形网络潮流计算matlab 利用matlab编程计算任意环形网络牛拉法潮流计算程序，程序通用性强
环形网络潮流计算matlab 利用matlab编程计算任意环形网络牛拉法潮流计算程序,程序通用性强,通过修改参数可以得到任意节点和网络的环形网络牛拉法潮流计算. 现有:6960642614794538 ...
- 环形网络潮流计算matlab 利用matlab编程计算任意环形网络牛拉法潮流计算程序，程序通用性强，通过修改参数
环形网络潮流计算matlab 利用matlab编程计算任意环形网络牛拉法潮流计算程序,程序通用性强,通过修改参数可以得到任意节点和网络的环形网络牛拉法潮流计算. 现有:6960642614794538 ...
- 用matlab计算潮流牛拉法,Matlab牛拉法潮流计算程序
Matlab牛拉法潮流计算程序 V(k)=sqrt(e(k)^2+f(k)^2); %计算各节点电压的模值 sida(k)=atan(f(k)./e(k))*180./pi; %计算各节点电压的角度 ...
- 潮流计算程序————节点导纳矩阵及其算法(Python语言描述)
目录 代码 例子 代码 import numpy as npdef jiedianmatrix(B, X, n):n1 =B.shape[0]Y = np.mat(np.zeros
最新文章
- 通过webbrowser控件获取验证码
- Hybris Commerce Cloud backoffice的一些使用截图 - OAuth工作中心
- Linux 命令之 declare -- 声明或显示 shell 变量
- java写一个99到0_Java中一个普通的循环为何从10开始到99连续相乘会得到0？
- 互联网晚报 | 3月1日 星期二 | 2022年苹果可能推出三款Apple Watch；“国家中小学智慧教育平台”投入试运行...
- 数据分析应用统计学之基本统计量【平均数、众数、中位数、四分位数】
- java分库校验商户流水号是否重复,asp中用数据库生成不重复的流水号
- nginx普通配置/负载均衡配置/ssl/https配置
- java 8 map reduce_java8之Lambda表达式 4：MapReduce开发案例
- 台达变频器485通讯接线图_台达变频器怎么接线 台达变频器接线图详解
- python 计算置信区间_python 计算数据中值的置信区间
- java导出繁体字word_利用简繁体字转换功能实现简体字文档转换为繁体字文档的方法...
- 自己尝试使用简单数据集实现决策树 代码——《机器学习实战》
- 【k8s】 k8s | 一 | 介绍
- 怎樣制作线段动画_教您怎样才能演示图形平移过程？
- 欧拉筛法（线性筛）的学习理解
- A*算法项目实践之一：栅格法的使用与障碍物栅格的生成
- 实现一个函数输入123456789，输出123，456，789”
- php段落首行缩进2字符,WordPress文章段落自动缩进两个字符方法
- valueOf()方法详解