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

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

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

版權說明：本文檔由用戶提供并上傳，收益歸屬內容提供方，若內容存在侵權，請進行舉報或認領

文檔簡介

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

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

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

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

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

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

7、t;　　四、程序代碼7　　五、計算結果24　　六、小結與展望24　　七、致謝25　　八、主要參考文獻25　　一、課程設計目的

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

9、編程思想及流程框圖。并附錄原程序。　　圖1 IEEE 9節(jié)點電力系統(tǒng)線路圖　　IEEE9節(jié)點電力系統(tǒng)原始數(shù)據(jù)：圖中GEN1節(jié)點為平衡節(jié)點，V1=1.04，GEN2和GEN3為PV節(jié)點，V2=1.025，P2=1.63，V3=1.025，P3=0.85，其他節(jié)點為PQ節(jié)點，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ù)(標幺值)　　表3 節(jié)點功率、電壓數(shù)據(jù)(標幺值)　　二、設計主要儀器設備和材料　　表4 設計主要設備　　三、牛頓法潮流計算的步驟及計算

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

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

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

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

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

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ù)設計，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é)點修正變量　　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>　　}　　/////文件五，計算雅可比矩陣，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é)點導納矩陣的計算，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é)點導納矩陣/////////////////　　ofstre

101、am fout("out.txt");　　fout << "--------------節(jié)點導納矩陣--------------------" << 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;　　}　　//////////////文件七，基于牛頓法的潮流計算

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("請輸入原始數(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é)點數(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é)點號、節(jié)點類型、節(jié)點電壓、相角、節(jié)點負載有功功率、節(jié)點負載無功率、節(jié)點輸入有功功率、節(jié)點輸入無功功率///<

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)線路標號、線路起始節(jié)點、線路終止節(jié)點、線路類型、線路電阻、線路電抗、線路對地導納///////////////　　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é)點導納矩陣////////////////　　MakeY( nB, nL, line, YG, YB );　　/////////////////內存分配//////////////

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é)點　　{　　p_Jtobus[nJ]=i;　　n

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

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

126、　　}　　for( i=0; i<npv+npq+1; i++ )　　{　　if( bus[i].Type==1 )//PQ 節(jié)點　　{</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 );　　//////////////計算負載潮流/////////////////　　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)潮流計算///////////////// 　　/////////////////各線路潮流計算////////////////　　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

眾賞文庫> 全部分類> 畢業(yè)設計

溫馨提示

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

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

文檔簡介

溫馨提示

最新文檔

評論

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

文檔簡介

溫馨提示

最新文檔

評論

免費下載