參考Matpower 在python上復(fù)現(xiàn)了 N-R 潮流計算 （原創(chuàng)程序：全網(wǎng)獨一份）

摘要：電力系統(tǒng)潮流計算是進行電力系統(tǒng)穩(wěn)態(tài)研究的基礎(chǔ)，一般是通過輸入的 P，Q，V，求解所有節(jié)點的 (V.real and V.imag) or (V.max and Vangle)，目前潮流計算工具或軟件主要有matlab的Matpower工具箱，以及一眾圖形化軟件。這些工具箱or軟件靈活性一般，難以隨意定制化功能。而目前市面上的python版本潮流計算代碼，主要還是基于for循環(huán)實現(xiàn)，沒有充分發(fā)揮numpy的向量化計算優(yōu)勢。本文主要基于Matpower工具箱中的算法，充分運用高等電力網(wǎng)絡(luò)分析理論，盡可能的將 For Loop -> Vector，在簡化程序的同時，提高效率（起碼編程效率提高了，計算效率提沒提不好說）。

一、需要掌握什么理論知識？

電網(wǎng)絡(luò)分析：節(jié)點-支路關(guān)聯(lián)矩陣、節(jié)點參數(shù)、支路參數(shù)等等。

高數(shù)，數(shù)值分析（牛拉法），矩陣論。等等。

二、程序如何使用？

1、requirement：numpy 、python3、（pandas、scipy：僅用于數(shù)據(jù)IO），沒啥版本要求，裝新版的總沒錯。

2、算例數(shù)據(jù)來源，matpower的/data/case...m

3、算例數(shù)據(jù)轉(zhuǎn)換： .m? -> .mat

4、運行：PF_NR.py

可以快樂的進行潮流計算了?。?/p>

5、算例分析

python版程序與matlab版程序運行結(jié)果幾乎完全相同。

三、結(jié)果和展望

程序問題：

1、新增、刪除節(jié)點數(shù)據(jù)只有接口，沒功能。

2、未考慮：發(fā)電機功率越限后的節(jié)點類型轉(zhuǎn)換。

3、僅編寫了直角坐標(biāo)下的雅克比矩陣

4、無償分享程序，emmm，僅有一些必要的注釋。

5、請大家多多點贊呦！！

6.程序鏈接

鏈接：https://pan.baidu.com/s/13CC6dP36Seo3PmRS1FU80w?

提取碼：ryp1