電力系統(tǒng)潮流計(jì)算機(jī)算法的開發(fā)
1.7萬字 61頁
包括程序代碼



摘 要

電力系統(tǒng)潮流計(jì)算是電力系統(tǒng)穩(wěn)態(tài)分析中最重要的基本電氣計(jì)算，它不僅是動(dòng)態(tài)穩(wěn)定分析的基礎(chǔ)，而且也是電力系統(tǒng)安全經(jīng)濟(jì)運(yùn)行和規(guī)劃設(shè)計(jì)的基礎(chǔ)。本文是在全面介紹牛頓-拉夫遜法和P-Q分解法計(jì)算潮流原理的基礎(chǔ)上，利用Matlab語言程序設(shè)計(jì)編制用計(jì)算機(jī)計(jì)算電力系統(tǒng)潮流的程序。P-Q分解法利用了電力系統(tǒng)的一些特有的運(yùn)行特性,對(duì)牛頓-拉夫遜法做了簡(jiǎn)化，同時(shí)采用稀疏技術(shù)，從而大大降低了對(duì)計(jì)算機(jī)內(nèi)存的要求，提高了計(jì)算速度，不僅可用于離線計(jì)算，而且可用于在線計(jì)算。該設(shè)計(jì)具有計(jì)算方法可靠、滿足收斂性、基本上能達(dá)到設(shè)計(jì)要求。

關(guān)鍵詞: 電力系統(tǒng)潮流計(jì)算,牛頓-拉夫遜法,P-Q分解法,Matlab
ABSTRACT

Power flow calculation is the most important and fundamental algorithm in power system, which is the basis of calculating system stability and safe economy running. This article is was introducing comprehensively Newton - Rough abdicates the law and in the P-Q resolution computation tidal current principle foundation, calculates the electrical power system tidal current using the Matlab language programming establishment with the computer the procedure. The P-Q resolution has used the electrical power system some unique movement characteristic, abdicated the law to Newton - Rough to make the simplification, simultaneously used the sparse technology, thus reduced greatly to the computer memory request, enhanced the computation speed, not only available in off-line computation, moreover available in online computation. This design has the computational method reliably, the satisfied astringency, basically can meet the design requirements.
Key word: Electrical power system tidal current computation, Newton - Rough abdicates the law, P-Q resolution, Matlab
目 錄

摘 要 I
ABSTRACT II
第一章 前 言 1
1.1 電力系統(tǒng)敘述 1
1.2 潮流計(jì)算簡(jiǎn)介 1
1.3 潮流計(jì)算的意義及其發(fā)展 2
1.4 本文的主要工作 4
第二章 電力網(wǎng)絡(luò)的數(shù)學(xué)模型 5
2.1 導(dǎo)納矩陣的原理及計(jì)算方法 5
2.1.1 自導(dǎo)納和互導(dǎo)納的確定方法 5
2.1.2 節(jié)點(diǎn)導(dǎo)納矩陣的性質(zhì)及意義 6
2.1.3 非標(biāo)準(zhǔn)變比變壓器等值電路 8
2.2 網(wǎng)絡(luò)方程的解法 8
第三章 簡(jiǎn)單電力系統(tǒng)的潮流計(jì)算 12
3.1 開式網(wǎng)絡(luò)的電壓和功率分布計(jì)算 12
3.2 簡(jiǎn)單閉式網(wǎng)絡(luò)的功率分布計(jì)算 15
第四章 電力系統(tǒng)潮流計(jì)算的計(jì)算機(jī)算法 24
4.1 概述 24
4.2 潮流計(jì)算的基本方程 24
4.2.1 節(jié)點(diǎn)的分類 24
4.2.2 潮流計(jì)算的約束條件 25
4.2.3 基本方程式 26
4.3 牛頓-拉夫遜法潮流計(jì)算 28
4.3.1 牛頓-拉夫遜法概要 28
4.3.2 牛頓-拉夫遜法的基本方程式 32
4.3.3 計(jì)算步驟和程序框圖 36
4.4 PQ分解法潮流計(jì)算 38
4.4.1 PQ分解法的基本方程式 38
4.4.2 計(jì)算步驟和程序框圖 40
4.4.3 實(shí)例 43
第五章 總 結(jié) 46
致 謝 47
參考文獻(xiàn) 48
附 錄 49



參考文獻(xiàn)
[1] 于永利. 基于廣義Tellegen定理的電力系統(tǒng)潮流算法研究[C]. 遼寧工學(xué)院, 2007: 1 ~ 4.
[2] 武曉朦, 張飛廷. 電力系統(tǒng)的P-Q分解法潮流計(jì)算[J]. 現(xiàn)代電子技術(shù), 2002. 8.(20): 105.
[3] 陳珩. 電力系統(tǒng)穩(wěn)態(tài)分析(第三版) [M]. 北京: 中國(guó)電力出版社, 2007: 111 ~ 123.
[4] 陳立新, 楊光宇. 電力系統(tǒng)分析[M]. 北京: 中國(guó)電力出版社, 2005: 100~135.
[5] Fernando L A.Computational complexity in power systems[J]. IEEE Trans on PAS, 1976, 95(4) : 1028 - 1037.
[6] 劉天琪, 邱曉燕. 電力系統(tǒng)分析理論[M]. 北京: 科學(xué)出版社, 2005: 70 ~ 98
[7] B Stott, O Alsac. Fast Decoupled Load Flow. IEEE Transactions on Power Apparatus and Systems, 1974, 93(3): 859-869.
[8] Borkowska B. Probabilistic Load Flow. IEEE Trans on Power Apparatus and Systems, 1974, PAS-93:752 - 759.
[9] Carpentier J. Optimal Power Flows. Int J Elec Power and Energy Syst, April



附 錄
%本程序的功能是用PQ分解法進(jìn)行潮流計(jì)算
function y=DDpq(n,m,pr,a,t0,v,P,Q)%
n=input('請(qǐng)輸入節(jié)點(diǎn)數(shù)：n=');% n為節(jié)點(diǎn)數(shù)
m=input('請(qǐng)輸入PQ節(jié)點(diǎn)數(shù)：m=');% m為PQ節(jié)點(diǎn)數(shù)
pr=input('請(qǐng)輸入誤差精度：pr=');% pr為誤差精度
a=input('請(qǐng)輸入矩陣：a=');% a為支路矩陣
t0=input('請(qǐng)輸入網(wǎng)絡(luò)結(jié)點(diǎn)的電壓相角迭代初值：t0=');% t0為網(wǎng)絡(luò)節(jié)點(diǎn)電壓相角迭代初值
v=input('請(qǐng)輸入節(jié)點(diǎn)的電壓迭代初值和pv結(jié)點(diǎn)的給定電壓值：v=');% v為網(wǎng)絡(luò)結(jié)點(diǎn)的電壓迭代初值和pv結(jié)點(diǎn)的給定電壓值
P=input('請(qǐng)輸入網(wǎng)絡(luò)結(jié)點(diǎn)的給定有功功率值：P=');% P為網(wǎng)絡(luò)結(jié)點(diǎn)的給定有功功率值
Q=input('請(qǐng)輸入網(wǎng)絡(luò)結(jié)點(diǎn)的給定無功功率值：Q=');% Q為網(wǎng)絡(luò)結(jié)點(diǎn)的給定無功功率值
Y=zeros(n,n);
for d=1:length(a(:))/4
if a(d,1)==0
Y(a(d,2),a(d,2))=Y(a(d,2),a(d,2))+1/a(d,3);
else
if a(d,4)>0
Y(a(d,1),a(d,1))=Y(a(d,1),a(d,1))+1/a(d,3);
Y(a(d,1),a(d,2))=-1/(a(d,3)*a(d,4));
Y(a(d,2),a(d,1))=-1/(a(d,3)*a(d,4));
......