自動(dòng)控制裝置課程設(shè)計(jì)--電力系統(tǒng)潮流計(jì)算的程序設(shè)計(jì)

上傳人：奔*** IP屬地：河北更新時(shí)間：2024-03-01 格式：doc 頁(yè)數(shù)：27 大小：471.00KB 人氣指數(shù)：12 舉報(bào) 版權(quán)申訴

自動(dòng)控制裝置課程設(shè)計(jì)--電力系統(tǒng)潮流計(jì)算的程序設(shè)計(jì)_第1頁(yè)

已閱讀1頁(yè)，還剩26頁(yè)未讀，繼續(xù)免費(fèi)閱讀

版權(quán)說(shuō)明：本文檔由用戶提供并上傳，收益歸屬內(nèi)容提供方，若內(nèi)容存在侵權(quán)，請(qǐng)進(jìn)行舉報(bào)或認(rèn)領(lǐng)

文檔簡(jiǎn)介

1、　　《自動(dòng)控制裝置語(yǔ)言程序設(shè)計(jì)》　　題目: 電力系統(tǒng)潮流計(jì)算的程序設(shè)計(jì)　　系別: 自動(dòng)化學(xué)院 　　專業(yè): 07電氣工程及其自動(dòng)化　　(低壓電力智能控制方向)

2、;　　前言　　潮流計(jì)算是研究電力系統(tǒng)穩(wěn)態(tài)運(yùn)行情況的一種基本電氣計(jì)算，常規(guī)潮流計(jì)算的任務(wù)是根據(jù)給定的運(yùn)行條件和網(wǎng)路結(jié)構(gòu)確定整個(gè)系統(tǒng)的運(yùn)行狀態(tài)，如各母線上的電壓（幅值及相角）、網(wǎng)絡(luò)中的功率分布以及功率損耗等。潮流計(jì)算的結(jié)果是電力系統(tǒng)穩(wěn)定計(jì)算和故障分析的基礎(chǔ)。 　　具體表現(xiàn)在以下方面： <p

3、>　　(1)在電網(wǎng)規(guī)劃階段,通過(guò)潮流計(jì)算,合理規(guī)劃電源容量及接入點(diǎn),合理規(guī)劃網(wǎng)架,選擇無(wú)功補(bǔ)償方案,滿足規(guī)劃水平的大、小方式下潮流交換控制、調(diào)峰、調(diào)相、調(diào)壓的要求。 　　(2)在編制年運(yùn)行方式時(shí),在預(yù)計(jì)負(fù)荷增長(zhǎng)及新設(shè)備投運(yùn)基礎(chǔ)上,選擇典型方式進(jìn)行潮流計(jì)算,發(fā)現(xiàn)電網(wǎng)中薄弱環(huán)節(jié),供調(diào)度員日常調(diào)度控制參考,并對(duì)規(guī)劃、基建部門提出改進(jìn)網(wǎng)架結(jié)構(gòu),加快基建進(jìn)度的建議。 <

4、p>　　(3)正常檢修及特殊運(yùn)行方式下的潮流計(jì)算,用于日運(yùn)行方式的編制,指導(dǎo)發(fā)電廠開(kāi)機(jī)方式,有功、無(wú)功調(diào)整方案及負(fù)荷調(diào)整方案,滿足線路、變壓器熱穩(wěn)定要求及電壓質(zhì)量要求。 　　(4)預(yù)想事故、設(shè)備退出運(yùn)行對(duì)靜態(tài)安全的影響分析及作出預(yù)想的運(yùn)行方式調(diào)整方案。 　　總結(jié)為在電力系統(tǒng)運(yùn)行方式和規(guī)劃方案的研究中，都需要進(jìn)行潮流計(jì)算以比較運(yùn)行方式或規(guī)劃供電方案的

5、可行性、可靠性和經(jīng)濟(jì)性。同時(shí)，為了實(shí)時(shí)監(jiān)控電力系統(tǒng)的運(yùn)行狀態(tài)，也需要進(jìn)行大量而快速的潮流計(jì)算。因此，潮流計(jì)算是電力系統(tǒng)中應(yīng)用最廣泛、最基本和最重要的一種電氣運(yùn)算。在系統(tǒng)規(guī)劃設(shè)計(jì)和安排系統(tǒng)的運(yùn)行方式時(shí)，采用離線潮流計(jì)算；在電力系統(tǒng)運(yùn)行狀態(tài)的實(shí)時(shí)監(jiān)控中，則采用在線潮流計(jì)算。　　隨著計(jì)算機(jī)技術(shù)的迅速發(fā)展和普及，電子計(jì)算機(jī)已成為分析計(jì)算復(fù)雜電力系統(tǒng)各種運(yùn)行情況的主要工具。<

6、p>　　目錄　　目錄1　　一、課程設(shè)計(jì)目的和任務(wù)2　　二、設(shè)計(jì)主要儀器設(shè)備和材料3　　三、牛頓法潮流計(jì)算的步驟及計(jì)算程序框圖4<b&g

7、t;　　四、程序代碼7　　五、計(jì)算結(jié)果24　　六、小結(jié)與展望24　　七、致謝25　　八、主要參考文獻(xiàn)25　　一、課程設(shè)計(jì)目的

8、和任務(wù)　　設(shè)計(jì)目的：通過(guò)對(duì)電力系統(tǒng)分析的學(xué)習(xí)，了解電力系統(tǒng)的運(yùn)行和能量管理，計(jì)算電力系統(tǒng)潮流分布以及有功和無(wú)功功率的分點(diǎn)，掌握電力系統(tǒng)環(huán)網(wǎng)（牛頓-拉夫遜法、P-Q分解法）的計(jì)算方法。并熟練使用計(jì)算機(jī)編程。　　設(shè)計(jì)任務(wù)：如圖所示電力系統(tǒng)節(jié)點(diǎn)潮流及阻抗、導(dǎo)納分布圖，用直角坐標(biāo)牛頓拉夫遜法（極坐標(biāo)牛頓-拉夫遜法、P-Q分解法）計(jì)算機(jī)編程計(jì)算各線路潮流分布。寫(xiě)出

9、編程思想及流程框圖。并附錄原程序。　　圖1 IEEE 9節(jié)點(diǎn)電力系統(tǒng)線路圖　　IEEE9節(jié)點(diǎn)電力系統(tǒng)原始數(shù)據(jù)：圖中GEN1節(jié)點(diǎn)為平衡節(jié)點(diǎn)，V1=1.04，GEN2和GEN3為PV節(jié)點(diǎn)，V2=1.025，P2=1.63，V3=1.025，P3=0.85，其他節(jié)點(diǎn)為PQ節(jié)點(diǎn)，S7=1.25+j0.5，S8=0.9+j0.3，S9=1.0+j0.35。<

10、;/p>　　表1 系統(tǒng)基本參數(shù)數(shù)據(jù)　　表2 系統(tǒng)線路原始參數(shù)(標(biāo)幺值)　　表3 節(jié)點(diǎn)功率、電壓數(shù)據(jù)(標(biāo)幺值)　　二、設(shè)計(jì)主要儀器設(shè)備和材料　　表4 設(shè)計(jì)主要設(shè)備　　三、牛頓法潮流計(jì)算的步驟及計(jì)算

11、程序框圖　　牛頓法潮流計(jì)算的求解過(guò)程大致可以分為以下幾步：　　（1）輸入系統(tǒng)的原始數(shù)據(jù)，包括系統(tǒng)的信息，各線路和變壓器所在節(jié)點(diǎn)的編號(hào)和等值電路中的參數(shù)；各負(fù)荷所在節(jié)點(diǎn)的編號(hào)及其所取用的有功功率和無(wú)功功率；作為PQ節(jié)點(diǎn)的發(fā)電機(jī)所在節(jié)點(diǎn)和編號(hào)及所給定的有功功率和無(wú)功功率，PV節(jié)點(diǎn)的編號(hào)和所接發(fā)電機(jī)的有功功率和電壓有效值的給定值，平衡節(jié)點(diǎn)的編號(hào)和給定的電壓有效

12、值。　?。?）形成節(jié)點(diǎn)導(dǎo)納矩陣。　?。?）給定各PQ節(jié)點(diǎn)的電壓初值和除平衡節(jié)點(diǎn)外各節(jié)點(diǎn)電壓相位的初值并組成待求量的初始向量和。　?。?）置迭代次數(shù)k=0。　?。?）應(yīng)用、及PV節(jié)點(diǎn)和平衡節(jié)點(diǎn)所給定的電壓，按式1：　　計(jì)算各PQ節(jié)點(diǎn)

13、的有功功率誤差和無(wú)功功率誤差以及各PV節(jié)點(diǎn)的有功功率誤差，并組成功率誤差向量和。　?。?）按式2：　　中的收斂判據(jù)判斷最大的功率誤差是否小于容許值，如果滿足則轉(zhuǎn)到第（11）步，否則進(jìn)行下一步。　?。?）應(yīng)用、按式3和式4：　　計(jì)算雅可

14、比矩陣元素，并形成式5中的雅可比矩陣。　?。?）解修正方程式5，得出和。　?。?）應(yīng)用式6：　　計(jì)算各節(jié)點(diǎn)電壓和相位的修正值，即新的初值。　?。?0）置k=k+1返回第（5）步繼續(xù)進(jìn)行下一輪迭代。<p&g

15、t;　?。?1）按式7：　　計(jì)算平衡節(jié)點(diǎn)的發(fā)電機(jī)（或注入）有功功率和無(wú)功功率，并按式8：　　計(jì)算各PV節(jié)點(diǎn)發(fā)電機(jī)（或注入）的無(wú)功功率；按式9和式10：　　計(jì)算各元件兩端的功率和損耗，最后輸出計(jì)算結(jié)果。　　與上述步驟相對(duì)應(yīng)的計(jì)算程序框圖

16、如下圖所示：　　四、程序代碼　　// 文件一，NEquation.h: interface for the NEquation class.//　　////////////////////////////////////////////////////////////////

17、//////　　#if !defined(AFX_NEQUATION_H__7D6DCBC7_B370_4F2E_A720_4FB8BA7B55D9__INCLUDED_)　　#define AFX_NEQUATION_H__7D6DCBC7_B370_4F2E_A720_4FB8BA7B55D9__INCLUDED_<p&

18、gt;　　#if _MSC_VER > 1000　　#pragma once　　#endif// _MSC_VER > 1000　　#define CHECKERROR　　class NEquation <

19、;b>　　{　　public:　　NEquation();　　virtual ~NEquation();　　public:　　v

20、oidSetSize(int size);　　double&Data(int lhs, int rhs);　　double&Value(int lhs);　　intRun();　　protected:<p&g

21、t;　　intm_nNumber;　　double*m_nDataBuffer;　　double*m_nValue;　　voidFactorial();　　voidForward();　　voidBack

22、ward();　　virtual doubleZEROD();　　virtual doubleZEROV();　　};　　NEquation::NEquation()

23、　　{　　m_nDataBuffer=NULL;　　m_nValue=NULL;　　m_nNumber=0;　　}　　NEquation::~NEquatio

24、n()　　{　　if(m_nDataBuffer != NULL) delete []m_nDataBuffer;　　if(m_nValue != NULL) delete []m_nValue;　　}</b&

25、gt;　　double NEquation::ZEROD(){return 0.0;}　　double NEquation::ZEROV(){return 0.0;}　　void NEquation::SetSize(int size)　　{</b

26、>　　if(size < 1)return; 　　if(m_nDataBuffer != NULL) delete []m_nDataBuffer;　　if(m_nValue != NULL) delete []m_nValue;　　m_nDataBuffer

27、 = new double[size * size];　　m_nValue = new double[size];　　for(int i=0; i< size*size; i++)m_nDataBuffer[i] = ZEROD();　　for(i=0; i< size; i++)m_nValue[i

28、] = ZEROV();　　m_nNumber = size;　　}　　double & NEquation::Data(int lhs, int rhs)　　{&

29、lt;p>　　#if defined(CHECKERROR)　　if((lhs<m_nNumber)&&(lhs>=0)&&(rhs<m_nNumber)&&(rhs>=0))　　#endif　　r

30、eturn m_nDataBuffer[lhs * m_nNumber + rhs];　　#if defined(CHECKERROR)　　else　　printf("Error");　　return m_nDat

31、aBuffer[0];　　#endif　　}　　double & NEquation::Value(int lhs)　　{&l

32、t;p>　　#if defined(CHECKERROR)　　if((lhs<m_nNumber)&&(lhs>=0))　　#endif　　return m_nValue[lhs];　　#if def

33、ined(CHECKERROR)　　else　　printf("Error");　　return m_nValue[0];　　#endif<p&g

34、t;　　}　　int NEquation::Run()　　{　　Factorial();　　Forward();　　Backward();</p&g

35、t;　　return 1;　　}　　double reverse( double ff){return (double)1.0/ff;};　　void NEquation::Factorial()

36、　　{　　int i, j, k;　　for(i=0; i<m_nNumber; i++)　　{　　//規(guī)格化format line;<

37、p>　　Data(i, i) = reverse(Data(i, i));　　for(j = i+1; j< m_nNumber; j++) Data(i, j) = Data(i, i) * Data(i, j);　　//消去　　for(j = i+1; j&

38、lt;m_nNumber; j++)　　{　　for(k=i+1; k<m_nNumber; k++)　　{　　Data(j,k) = Data(j, k) - Data(j,i) * Da

39、ta(i, k);　　}　　}　　}　　}　　void NEquation::

40、Forward()　　{　　int i, j;　　for(i=0; i< m_nNumber; i++)　　{

41、　　Value(i) = Data(i, i) * Value(i);　　for(j = i+1; j < m_nNumber; j++)　　Value(j) = Value(j) - Data(j, i) * Value(i);　　}<p

42、>　　}　　void NEquation::Backward()　　{　　int i, j;　　for(i = 1; i < m_nNumber;

43、 i++)　　{　　for(j = 1; j < i+1; j++)　　Value(m_nNumber - i -1) -= Data(m_nNumber -i-1, m_nNumber -j) * Value(m_nNumber -j); &l

44、t;p>　　}　　}　　#endif// !defined(AFX_NEQUATION_H__7D6DCBC7_B370_4F2E_A720_4FB8BA7B55D9__INCLUDED_)　　//////////文件二，LF.h

45、//////////////　　#ifndef _LF_LCR_　　#define _LF_LCR_　　#include "DataDesign.h"　　int MakeY( int nB, int nL, Line* sL, double** YG,

46、 double** YB );　　int CalDeltaPQ( int npv, int npq, Bus* bus, double** YG, double** YB, int* p_Jtobus, double* deltaf );　　int FormJacobian( int npv, int npq, Bus* bus, double** YG

47、, double** YB, int* p_Jtobus, double** Jac );　　#endif//_LF_LCR_　　////////文件三，數(shù)據(jù)設(shè)計(jì)，DataDesign.h //////////　　#ifndef _DATA_LCR_INCLUDE_

48、#define _DATA_LCR_INCLUDE_　　struct Line　　{　　int Num,NumI,NumJ,Type;　　double R,X,B,K;　　}

49、;　　struct Bus　　{　　int Num ;　　double Volt,Phase,GenP,GenQ,LoadP,LoadQ;<p&g

50、t;　　int Type;　　};　　struct Shunt　　{　　int Num,NumI;　　double

51、G,B;　　};　　#endif//_DATA_LCR_INCLUDE_　　/////文件四，F(xiàn)ormF.cpp //////////　　#include <math.h>　　#

52、include <fstream.h>　　#include "LF.h"　　#include "DataDesign.h"　　#include <stdio.h>　　//form delta p and de

53、lta q　　int CalDeltaPQ( int npv, int npq, Bus* bus, double** YG, double** YB, int* p_Jtobus, double* deltaf )　　{　　int i,j;//循環(huán)變量

54、;　　int k;//節(jié)點(diǎn)修正變量　　double i1,i2;　　double dphase;　　for(i=0;i<npv+npq+npq;i++)　　{

55、　　if(i<npv+npq) //form deltaP　　{　　i1=0;　　k=p_Jtobus[i];　　for(j=0;j<npv+npq+1;j++)<

56、;/p>　　{　　dphase=bus[k].Phase-bus[j].Phase;　　i1+=bus[j].Volt*(YG[k][j]*cos(dphase)+YB[k][j]*sin(dphase));　　}</b&g

57、t;　　deltaf[i]=bus[k].GenP-bus[k].LoadP-bus[k].Volt*i1;　　}　　else　　{ //form delta

58、Q　　i2=0;　　k=p_Jtobus[i];　　for(j=0;j<npv+npq+1;j++)　　{　　dphase=bus[k].Phase

59、-bus[j].Phase;　　i2+=bus[j].Volt*(YG[k][j]*sin(dphase)-YB[k][j]*cos(dphase));　　}　　deltaf[i]=bus[k].GenQ-bus[k].LoadQ-bus[k].Volt*i2;</p&

60、gt;　　}　　}　　for(i=0;i<npv+npq*2;i++)　　deltaf[i]*=(-1);　　return 1;&

61、lt;/p>　　}　　/////文件五，計(jì)算雅可比矩陣，F(xiàn)ormJ.cpp ///////////　　#include <math.h>　　#include <fstream.h>

62、　#include "LF.h"　　#include "DataDesign.h"　　#include <malloc.h>　　#include <string.h>　　//////////////生成雅可比

63、矩陣///////////////////　　int FormJacobian( int npv, int npq, Bus* bus, double** YG, double** YB, int* p_Jtobus, double** Jac )　　{

64、int i,j;　　int np = npv+npq;　　double dphase;　　////////生成 H 矩陣////////////　　double** H;　　H = (double**)new

65、char[(np)*sizeof(void*)];　　for( i=0;i<np;i++ )　　{　　H[i]=new double[np];　　memset(H[i],0,sizeof(double)*(np) );//清零

66、　　}　　for(i=0;i<np;i++)　　{　　double temp=0;　　for(j=0;j<np;j++)

67、　　{　　if(i!=j)　　{　　dphase=bus[i+1].Phase-bus[j+1].Phase;　　H[i][j]=bus[

68、i+1].Volt*bus[j+1].Volt*(YG[i+1][j+1]*sin(dphase)-YB[i+1][j+1]*cos(dphase));　　temp+=H[i][j];　　}　　}

69、　　double dphase1=bus[i].Phase-bus[0].Phase;　　H[i][i]=temp*(-1)-bus[i+1].Volt*bus[0].Volt*(YG[i+1][0]*sin(dphase1)-YB[i+1][0]*cos(dphase));　　}<p&

70、gt;　　ofstream fout1("out.txt",ios::app);　　/////////生成 L 矩陣////////////　　double** L;　　L = (double**)new char[(np)*sizeof(void*)];<p&g

71、t;　　for( i=0;i<np;i++ )　　{　　L[i]=new double[np];　　memset(L[i],0,sizeof(double)*(np) );　　}</p

72、>　　for(i=npv;i<np;i++)　　{　　for(j=npv;j<np;j++)　　{　　if(i!=j)

73、　　{　　L[i][j]=H[i][j];　　}　　}　　L[i][i]=-1*H[i][i]-2*bus[i+1]

74、.Volt*bus[i+1].Volt*YB[i+1][i+1];　　}　　/////////////生成 J 矩陣///////////////////　　double** J;　　J = (double**)new char[(n

75、p)*sizeof(void*)];　　for( i=0;i<np;i++ )　　{　　J[i]=new double[np];　　memset(J[i],0,sizeof(double)*(np) );

76、　　}　　for(i=npv;i<np;i++)　　{　　double temp=0;　　for(j=0;j<np+1;j++)<p&g

77、t;　　{　　if(i!=j)　　{　　dphase=bus[i+1].Phase-bus[j+1].Phase;　　J[i][j]=-1*bus[i+1].

78、Volt*bus[j+1].Volt*(YG[i+1][j+1]*cos(dphase)+YB[i+1][j+1]*sin(dphase));　　temp+=-1*J[i][j];　　}　　}

79、J[i][i]=temp;　　}　　/////////////////生成 N 矩陣 /////////////////　　double** N;　　N = (double**)new char[(np)*sizeof(void*)]

80、;　　for( i=0;i<np;i++ )　　{　　N[i]=new double[np];　　memset(N[i],0,sizeof(double)*(np) );

81、;　　}　　for(i=0;i<np;i++)　　{　　for(j=npv;j<np;j++)　　{

82、if(i!=j)　　{　　N[i][j]=-1*J[i][j];　　}　　}　　N[i][i]=

83、J[i][i]+2*bus[i+1].Volt*bus[i+1].Volt*YG[i+1][i+1];　　}　　for(i=0;i<np;i++)　　for(j=0;j<np;j++)　　Jac[i][j]=H[i][j]

84、;　　for(i=np;i<np+npq;i++)　　for(j=0;j<np;j++)　　Jac[i][j]=J[i-np+2][j];　　for(i=0;i<np;i++)　　for(j=np;j<np

85、+npq;j++)　　Jac[i][j]=N[i][j-np+2];　　for(i=np;i<np+npq;i++)　　for(j=np;j<np+npq;j++)　　Jac[i][j]=L[i-np+2][j-np+2];<p&

86、gt;　　//////////////////輸出雅可比矩陣///////////////////////　　ofstream fout("out.txt",ios::app);　　static int dd=0; dd++;　　fout << "----------

87、----"<<"第"<<dd<<"次迭代后的"<<"雅可比矩陣--------------------" << endl;　　for(i=0; i<np+npq; i++)　　{&l

88、t;/p>　　for(j=0; j<np+npq; j++ )　　fout << Jac[i][j] << "\t";　　fout << endl;　　}<

89、p>　　fout.close();　　return 1;　　}　　/////////////文件六，節(jié)點(diǎn)導(dǎo)納矩陣的計(jì)算，MakeY.cpp /////////////　　#include &

90、lt;math.h>　　#include <fstream.h>　　#include "LF.h"　　#include <stdio.h>　　#include<iomanip.h><

91、;p>　　int MakeY( int nB, int nL, Line* sL, double** YG, double** YB )　　{　　int i,j;　　for(i=0;i<nB;i++)&

92、lt;p>　　for(j=0;j<nB;j++)　　{　　YG[i][j]=0.0;　　YB[i][j]=0.0;　　};　　for(i=0;

93、i<nL;i++)　　{　　int start,end;　　double r,x,b,k,g,b1;　　start=sL[i].NumI;　　end=sL[i].NumJ;</p&

94、gt;　　r=sL[i].R;　　x=sL[i].X;　　b=sL[i].B;　　k=sL[i].K;　　g=r/(r*r+x*x);　　b1=-x/(r*r+x*x);<p&g

95、t;　　if(sL[i].Type==2)//線路中含有變壓器　　{　　YG[start][end]=-1*g/k;　　YG[end][start]=-1*g/k;　　YB[start][end]=-1*b1/k;

96、;　　YB[end][start]=-1*b1/k;　　YG[start][start]+=g/k/k;　　YG[end][end]+=g;　　YB[start][start]+=0.5*b+b1/k/k;　　YB[end][end]+=0.5*b+

97、b1;　　continue;　　}　　if(sL[i].Type==1)//線路中不含變壓器　　{　　YG[sta

98、rt][end]=-1*g;　　YG[end][start]=-1*g;　　YB[start][end]=-1*b1;　　YB[end][start]=-1*b1;　　YG[start][start]+=g;　　YG[end][e

99、nd]+=g;　　YB[start][start]+=0.5*b+b1;　　YB[end][end]+=0.5*b+b1;　　continue;　　}

100、printf("Error Input !\n");　　return 0;　　}　　/////////////輸出節(jié)點(diǎn)導(dǎo)納矩陣/////////////////　　ofstre

101、am fout("out.txt");　　fout << "--------------節(jié)點(diǎn)導(dǎo)納矩陣--------------------" << endl;　　for(i=0;i<nB;i++)　　{&

102、lt;/p>　　for(j=0;j<nB;j++)　　fout <<setw(2)<< YG[i][j] << "+j" << YB[i][j] <<" ";　　fout << endl;<

103、/p>　　}　　fout.close();　　return 1;　　}　　//////////////文件七，基于牛頓法的潮流計(jì)算

104、程序，LF.cpp ///////////////　　#include <stdio.h>　　#include <malloc.h>　　#include <string.h>　　#include <math.h></p

105、>　　#include <fstream.h>　　#include "NEquation.h"　　#include "LF.h"　　int main( int argc, char* argv[] )<p

106、>　　{　　///////////////////讀入數(shù)據(jù)///////////////////////　　int nB, nL, nSH;　　char filename[64];　　printf("請(qǐng)輸入原始數(shù)據(jù)的文

107、件名:(0 for default value)\n");　　scanf("%s", filename);　　if( strcmp(filename,"0")==0 )　　{　　str

108、cpy( filename, "in.txt" );　　}　　ifstream fin;　　fin.open(filename);　　char buffer[128];

109、fin.getline(buffer,128);　　fin >> nB >> nL >> nSH;//讀入系統(tǒng)的節(jié)點(diǎn)數(shù)、線路數(shù)和分支數(shù)　　double** YG, **YB;　　YG = (double**)new char[(nB)*sizeof(void*)];<

110、/p>　　for( int i=0;i<nB;i++ )　　{　　YG[i]=new double[nB];　　memset(YG[i],0,sizeof(double)*(nB) );<b&g

111、t;　　}　　YB = (double**)new char[(nB)*sizeof(void*)];　　for( i=0;i<nB;i++ )　　{　　YB[i]=new double[nB];</p

112、>　　memset(YB[i],0,sizeof(double)*(nB) );　　}　　Line* line;　　line = new Line[nL];　　Bus* bus;</

113、b>　　bus = new Bus[nB];　　Shunt* shunt;　　shunt = new Shunt[nSH];　　fin.getline(buffer,128);　　fin.getline(buffer,128

114、);　　int npv=0, npq=0;　　for( i=0; i<nB; i++ )　　{　　///讀入系統(tǒng)的節(jié)點(diǎn)號(hào)、節(jié)點(diǎn)類型、節(jié)點(diǎn)電壓、相角、節(jié)點(diǎn)負(fù)載有功功率、節(jié)點(diǎn)負(fù)載無(wú)功率、節(jié)點(diǎn)輸入有功功率、節(jié)點(diǎn)輸入無(wú)功功率///<

115、;/p>　　fin >> bus[i].Num >> bus[i].Type >> bus[i].Volt >> bus[i].Phase >> bus[i].LoadP >> bus[i].LoadQ >> bus[i].GenP >> bus[i].GenQ;　　bus

116、[i].Phase =bus[i].Phase*3.141592/180.0;　　bus[i].Num -=1;　　if( bus[i].Type==2)　　npv++;　　if( bus[i].Type==1)

117、　　npq++;　　};　　fin.getline(buffer,128);　　fin.getline(buffer,128);　　for( i=0; i<nL; i++

118、 )　　{　　/////////////讀入系統(tǒng)線路標(biāo)號(hào)、線路起始節(jié)點(diǎn)、線路終止節(jié)點(diǎn)、線路類型、線路電阻、線路電抗、線路對(duì)地導(dǎo)納///////////////　　fin >> line[i].Num >> line[i].NumI>> li

119、ne[i].NumJ >> line[i].Type >> line[i].R >> line[i].X >> line[i].B>> line[i].K;　　line[i].NumI--;　　line[i].NumJ--;　　}<

120、;/b>　　fin.getline(buffer,128);　　fin.getline(buffer,128);　　for( i=0; i<nSH; i++ )　　{<b&g

121、t;　　}　　fin.close();　　////////////組成節(jié)點(diǎn)導(dǎo)納矩陣////////////////　　MakeY( nB, nL, line, YG, YB );　　/////////////////內(nèi)存分配//////////////

122、/////　　double* deltaf = new double[npv+npq*2];　　int* p_Jtobus = new int[npv+npq+npq];　　int nJ=0;　　for( i=0; i<npv+

123、npq+1; i++ )　　{　　if( bus[i].Type==2 )//PV 節(jié)點(diǎn)　　{　　p_Jtobus[nJ]=i;　　n

124、J++;　　}　　}　　for( i=0; i<npv+npq+1; i++ )　　{

125、if( bus[i].Type==1 )//PQ 節(jié)點(diǎn)　　{　　p_Jtobus[nJ]=i;　　nJ++;　　}

126、　　}　　for( i=0; i<npv+npq+1; i++ )　　{　　if( bus[i].Type==1 )//PQ 節(jié)點(diǎn)　　{</p&

127、gt;　　p_Jtobus[nJ]=i;　　nJ++;　　}　　}　　double** Jac;

128、;　　Jac = (double**)new char[(nJ)*sizeof(void*)];　　for( i=0;i<nJ;i++ )　　{　　Jac[i]=new double[nJ];　　memset(Jac[i],0,

129、sizeof(double)*(nJ) );　　}　　double* deltaX;　　deltaX = new double[nJ];　　double mdelta=0;

130、　int it=0;　　do　　{　　it++;　　mdelta = 0.0;　　//form

131、delta p and delta q　　CalDeltaPQ( npv, npq, bus, YG, YB, p_Jtobus, deltaf );　　/////////////形成雅可比矩陣////////////////　　FormJacobian( npv, npq, bus, YG, YB, p_Jto

132、bus, Jac );　　//////////////計(jì)算負(fù)載潮流/////////////////　　NEquation equ;　　equ.SetSize(nJ);　　for( i=0; i<nJ; i++ )<b

133、>　　{　　equ.Value(i) = deltaf[i];　　for( int j=0; j<nJ; j++ )　　equ.Data(i,j) = Jac[i][j];　　}

134、　　equ.Run();　　/////////////修正電壓幅值和電壓相角///////////　　for( i=0; i<nJ; i++ )　　{　　deltaX[i] = equ.Value(i);<

135、;/p>　　int busi = p_Jtobus[i];　　if( i<npv+npq )　　bus[busi].Phase -= equ.Value(i);　　else　　bus[busi].Vo

136、lt -= equ.Value(i)*bus[busi].Volt;　　if( fabs(equ.Value(i))>fabs(mdelta) )　　mdelta = equ.Value(i);　　}　　}while( fabs

137、(mdelta)>1e-5 );//精度控制　　//////////////////系統(tǒng)潮流計(jì)算///////////////// 　　/////////////////各線路潮流計(jì)算////////////////　　double* p_Pij, *p_Qij, *p_Pji, *p_Qji;<

138、/p>　　p_Pij = new double[nL];　　p_Qij = new double[nL];　　p_Pji = new double[nL];　　p_Qji = new double[nL];　　for( i=0; i<nL

眾賞文庫(kù)> 全部分類> 畢業(yè)設(shè)計(jì)

溫馨提示

1. 本站所有資源如無(wú)特殊說(shuō)明，都需要本地電腦安裝OFFICE2007和PDF閱讀器。圖紙軟件為CAD,CAXA,PROE,UG,SolidWorks等.壓縮文件請(qǐng)下載最新的WinRAR軟件解壓。
2. 本站的文檔不包含任何第三方提供的附件圖紙等，如果需要附件，請(qǐng)聯(lián)系上傳者。文件的所有權(quán)益歸上傳用戶所有。
3. 本站RAR壓縮包中若帶圖紙，網(wǎng)頁(yè)內(nèi)容里面會(huì)有圖紙預(yù)覽，若沒(méi)有圖紙預(yù)覽就沒(méi)有圖紙。
4. 未經(jīng)權(quán)益所有人同意不得將文件中的內(nèi)容挪作商業(yè)或盈利用途。
5. 眾賞文庫(kù)僅提供信息存儲(chǔ)空間，僅對(duì)用戶上傳內(nèi)容的表現(xiàn)方式做保護(hù)處理，對(duì)用戶上傳分享的文檔內(nèi)容本身不做任何修改或編輯，并不能對(duì)任何下載內(nèi)容負(fù)責(zé)。
6. 下載文件中如有侵權(quán)或不適當(dāng)內(nèi)容，請(qǐng)與我們聯(lián)系，我們立即糾正。
7. 本站不保證下載資源的準(zhǔn)確性、安全性和完整性, 同時(shí)也不承擔(dān)用戶因使用這些下載資源對(duì)自己和他人造成任何形式的傷害或損失。

自動(dòng)控制裝置課程設(shè)計(jì)--電力系統(tǒng)潮流計(jì)算的程序設(shè)計(jì)

文檔簡(jiǎn)介

溫馨提示

最新文檔

評(píng)論

自動(dòng)控制裝置課程設(shè)計(jì)--電力系統(tǒng)潮流計(jì)算的程序設(shè)計(jì)

文檔簡(jiǎn)介

溫馨提示

最新文檔

評(píng)論

免費(fèi)下載