離線潮流計(jì)算課程設(shè)計(jì)

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

2023年全國(guó)碩士研究生考試考研英語(yǔ)一試題真題（含答案詳解+作文范文）_第1頁(yè)

已閱讀1頁(yè)，還剩22頁(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、　　離線潮流計(jì)算課程設(shè)計(jì)　　班級(jí) 學(xué)號(hào) 姓名 　　一、課程設(shè)計(jì)的目的　　離線潮流計(jì)算課程設(shè)計(jì)是“電力系統(tǒng)穩(wěn)態(tài)分析”課程教學(xué)的重要實(shí)踐環(huán)節(jié)，目的是使學(xué)生加深理解所學(xué)知識(shí)，鍛煉應(yīng)用理論解決實(shí)際問(wèn)題的能力。

2、　　課程設(shè)計(jì)內(nèi)容　　1、手算算例1（精度取10 -2），求潮流計(jì)算結(jié)果。　　2、編寫(xiě)潮流計(jì)算程序，計(jì)算算例2 。　　三、課程設(shè)計(jì)要求　　計(jì)算機(jī)語(yǔ)言：編程語(yǔ)

3、言可任選，建議使用C語(yǔ)言，使用VC6.0平臺(tái) 。　　潮流計(jì)算方法：牛頓—拉夫遜法和P-Q分解法。　　編寫(xiě)潮流計(jì)算程序，并在計(jì)算機(jī)上調(diào)試通過(guò)，運(yùn)行結(jié)果正確。　　要求潮流計(jì)算程序具有一定的通用性，基礎(chǔ)數(shù)據(jù)以數(shù)據(jù)文件的形式輸入。　　撰寫(xiě)課程設(shè)計(jì)報(bào)告，內(nèi)容包括：

4、;　　潮流計(jì)算的意義　　潮流計(jì)算幾種主要方法的比較　　手算算例1的計(jì)算過(guò)程及結(jié)果　　程序框圖、源程序及符號(hào)說(shuō)明　　算例及計(jì)算結(jié)果&l

5、t;b>　　、驗(yàn)收方式　　以答辯的形式進(jìn)行驗(yàn)收，演示程序并回答有關(guān)問(wèn)題。　　、算例　　節(jié)點(diǎn)①為平衡節(jié)點(diǎn)，U1=1∠0°，節(jié)點(diǎn)②為PV節(jié)點(diǎn)，P2=0.4，U2=1.1，節(jié)點(diǎn)③為PQ節(jié)點(diǎn)，S3 = 0.8+ j 0.6，網(wǎng)絡(luò)結(jié)構(gòu)和參數(shù)如圖1所示

6、。　　GN1節(jié)點(diǎn)為平衡節(jié)點(diǎn)，V9=1.04，δ9=0.0，GEN2和GEN3為PV節(jié)點(diǎn)，V7=1.025，P7=1.63，V8=1.025，P8=0.85，其它節(jié)點(diǎn)為PQ節(jié)點(diǎn)，S1=1.25+j0.5，S3=0.9+j0.3，S5=1.0+j0.35，ε=0.0001，網(wǎng)絡(luò)結(jié)構(gòu)如圖2所示。　　圖1 算例1</p

7、>　　圖2 算例2　　Txt 文檔原始數(shù)據(jù)　　Line　　1 1 1 1 2 0.01 0.085000001 0.088　　2 1 1 1 6 0.032000

8、002 0.160999998 0.152999997　　3 1 1 2 3 0.017000001 0.092 0.079000004　　4 2 1 2 9 0 0.057999998 1　　5 1 1 3 4 0.039000001 0.170000002 0.179000005<

9、;p>　　6 1 1 4 5 0.012 0.101000004 0.104999997　　7 2 1 4 8 0 0.059 1　　8 1 1 5 6 0.018999999 0.071999997 0.075000003　　9 2 1 6 7 0 0.063000001 1<

10、;p>　　Node　　1 1 -1.25 -0.5　　2 1 0.000000 0.000000　　3 1 -0.9 -0.3　　4 1 0.000000 0.000000　　5 1 -

11、1.000000 -0.35　　6 1 0.000000 0.000000　　7 2 1.63 1.025　　8 2 0.85 1.025　　9 0 1.04 0.000000　　上機(jī)程序：編程

12、;　　#include <stdafx.h>　　#include <string.h>　　#include <stdio.h>　　#include <stdlib.h>　　#incl

13、ude <conio.h>　　#include <float.h>　　#include <math.h>　　using namespace System;　　#define LINE_TOTAL_NUM 9<

14、;p>　　#define NODE_TOTAL_NUM 9　　#define NODE_DATA_FILENAME "node.txt"　　#define LINE_DATA_FILENAME "line.txt"　　#define B "B.txt&quo

15、t;　　#define G "G.txt"　　#define deltaPQ "deltapq.txt"　　#define JACCOBI "jaccobi.txt"　　#define PQ "pq.t

16、xt"　　#define line_S "linedeltaS.txt"　　#define U "u.txt"　　#define line_pq "lineS.txt"　　struct NodeData

17、　　{　　unsigned int Index; //node index　　unsigned int Type; //node type: PQ:1,PV:2,balance point:0　　float FirstInp

18、ut; //PQ or PV:active power,balance point:V　　float SecondInput; //PQ:reactive power,PV:V,balance point:angle 　　};　　struct NodeDat

19、a gNodeData[NODE_TOTAL_NUM]; // define a global variable of node 確定節(jié)點(diǎn)讀入結(jié)構(gòu)體　　struct LineData　　{　　unsigned int Index; //line index</

20、p>　　unsigned int Type; //line type: Line:1,tranformer:2,line connected to ground:3　　unsigned int Status; //line status: on:1,off:0　　unsigned int SrcNode; //sour

21、ce node index of this line　　unsigned int DestNode;//destination node index of this line　　float Resistance; //resistance of the line　　float Reactance;

22、 //reactance of the line　　float Conductance; //Line: conductance B/2,transformer:change ratio　　};　　struct LineData gLineData[LINE_TOTA

23、L_NUM]; //確定支路節(jié)點(diǎn)導(dǎo)納結(jié)鉤體　　struct NodeVoltage　　{　　float Voltage;　　float Angle;　　};</b

24、>　　struct NodeVoltage gNodeVoltage[NODE_TOTAL_NUM];　　struct DeltaPQ　　{ //定義不平衡功率的結(jié)構(gòu)體　　float Del

25、taP; //有功　　float DeltaQ; //無(wú)功　　};　　struct DeltaPQ gDeltaPQ[NODE_

26、TOTAL_NUM]; //定義不平衡功率數(shù)組gUnbPower　　void SolveEquation(unsigned int Dimension,float** FactorMatrix,float* ConstVector);　　//void SolveEquation(unsigned int Dimension,

27、float** FactorMatrix,float* ConstVector);　　float gY_G[NODE_TOTAL_NUM][NODE_TOTAL_NUM];　　float gY_B[NODE_TOTAL_NUM][NODE_TOTAL_NUM];　　int main(array<System

28、::String ^> ^args)　　{　　Console::WriteLine(L"Hello World");　　float node_P[9],node_Q[9],line_deltap[9],line_deltaq[9];

29、;　　unsigned int i,j,k,n,m; //cycle variable;　　FILE* fp; //file structure　　unsigned int PVNum=0; //the number of PV node　　fp=fopen(NODE_DATA_FILE

30、NAME,"r"); //open and read the node data　　for(i=0;i<NODE_TOTAL_NUM;i++)　　{　　fscanf(fp," %d %d %f %f ",&gNodeDat

31、a[i].Index,　　&gNodeData[i].Type,　　&gNodeData[i].FirstInput, &l

32、t;/p>　　&gNodeData[i].SecondInput);　　//initiate the NodeVoltage according to node type 讀入節(jié)點(diǎn)數(shù)據(jù)PQ或PV　　switch(gNodeData[i].Type)　　{&

33、lt;/b>　　case 1: //PQ node　　gNodeVoltage[i].Voltage = 1.0;　　gNodeVoltage[i].Angle = 0.0;　　break;&

34、lt;p>　　case 2: //PV node　　gNodeVoltage[i].Voltage = gNodeData[i].SecondInput;　　gNodeVoltage[i].Angle = 0.0;　　PVNum++;&

35、lt;p>　　break;　　case 0: //balance node　　gNodeVoltage[i].Voltage = gNodeData[i].FirstInput;　　gNodeVoltage[i].Angle = gNodeData[i].Sec

36、ondInput;　　break;　　default:　　printf("there is no such a node type, data has mistake");&l

37、t;b>　　return 0;　　}　　}　　fclose(fp);　　fp=fopen(LINE_DATA_FILENAME,"r"); //ope

38、n and read the line data 讀節(jié)點(diǎn)導(dǎo)納數(shù)據(jù)　　for(i=0;i<LINE_TOTAL_NUM;i++)　　{　　fscanf(fp," %u %u %u %u %u %f %f %f",&gLineData[i].

39、Index,　　&gLineData[i].Type,　　&gLineData[i].Status,　　&gLineData[i].SrcNode,　　&gLineData[i].DestNode,

40、　　&gLineData[i].Resistance,　　&gLineData[i].Reactance,　　&gLineData[i].Conductance);　　}　　fclose(fp);</p&

41、gt;　　float line_g[9], line_b[9],line_g0[9], line_b0[9];　　for(i=0;i<NODE_TOTAL_NUM;i++)　　{ line_g[i]=0;　　line_b[i]=0;

42、　line_g0[i]=0;　　line_b0[i]=0;　　for(j=0;j<NODE_TOTAL_NUM;j++)　　{　　gY_G[i][j] = 0.0;　　gY_B[i][j

43、] = 0.0;　　}　　}　　for(i=0;i<LINE_TOTAL_NUM;i++)　　{　　float

44、 tAbs;　　float LineG,LineB;　　unsigned int tSrcNode,tDestNode;　　tSrcNode = gLineData[i].SrcNode - 1;//counting from zero　　tDestNode = gLin

45、eData[i].DestNode - 1;　　tAbs = pow(gLineData[i].Resistance,2)+pow(gLineData[i].Reactance,2);//由阻抗求導(dǎo)納　　LineG = gLineData[i].Resistance/tAbs;　　LineB = (-1)*gLi

46、neData[i].Reactance/tAbs;　　if(gLineData[i].Type)　　{　　if(gLineData[i].Type == 1)//it means Line 線路　　{</b&

47、gt;　　gY_G[tSrcNode][tSrcNode]+= LineG;　　gY_B[tSrcNode][tSrcNode] += LineB + gLineData[i].Conductance;　　gY_G[tDestNode][tDestNode] += LineG;<p&

48、gt;　　gY_B[tDestNode][tDestNode] += LineB + gLineData[i].Conductance;　　gY_G[tSrcNode][tDestNode] += (-1)*LineG;　　gY_B[tSrcNode][tDestNode] += (-1)*LineB;　　gY_

49、G[tDestNode][tSrcNode] += (-1)*LineG;　　gY_B[tDestNode][tSrcNode] += (-1)*LineB;　　line_g[i]=LineG;　　line_b[i]=LineB;　　line_b0[i]=gLineDat

50、a[i].Conductance;　　}　　else if(gLineData[i].Type == 2) //it means Transformer 變壓器　　{　　float VolRat

51、io = gLineData[i].Conductance;　　gY_G[tSrcNode][tSrcNode] += LineG/pow(VolRatio,2);　　gY_B[tSrcNode][tSrcNode] += LineB/pow(VolRatio,2);　　gY_G[tDestNode][tDest

52、Node] += LineG;　　gY_B[tDestNode][tDestNode] += LineB;　　gY_G[tSrcNode][tDestNode] += (-1)*LineG/VolRatio;　　gY_B[tSrcNode][tDestNode] += (-1)*LineB/VolRatio;&l

53、t;/p>　　gY_G[tDestNode][tSrcNode] += (-1)*LineG/VolRatio;　　gY_B[tDestNode][tSrcNode] += (-1)*LineB/VolRatio;　　line_g[i]=LineG/VolRatio;　　line_b[

54、i]=LineB/VolRatio;　　}　　else if(gLineData[i].Type == 3)//it means line connected to ground　　{　　unsi

55、gned int tNotZeroNode;　　if(gLineData[i].SrcNode != 0)　　tNotZeroNode = gLineData[i].SrcNode-1;　　else　　tNotZeroNode = g

56、LineData[i].DestNode-1;;　　gY_G[tNotZeroNode][tNotZeroNode] += LineG;　　gY_B[tNotZeroNode][tNotZeroNode] += LineB;　　}else<p&g

57、t;　　return 0;　　}　　}　　//form deltaPQ　　//form Jacobin Matrix

58、　　//N-R method　　fp=fopen(G,"w"); //write the ymatrix data to file　　for(i=0;i<9;i++)　　{　　for(j=0;j<9;

59、j++)　　fprintf(fp," %f ",gY_G[i][j]);　　fprintf(fp,"\n");　　}　　fprintf(fp,"\n\n");</p&g

60、t;　　fclose(fp);　　fp=fopen(B,"w"); //write the ymatrix data to file　　for(i=0;i<9;i++)　　{　　for

61、(j=0;j<9;j++)　　fprintf(fp," %f ",gY_B[i][j]);　　fprintf(fp,"\n");　　}　　fprintf(fp,"\n\n&quo

62、t;);　　fclose(fp);　　for(m=0;m<5;m++)　　{ 　　unsigned int tSrcNode,tDestNode;　　float p1,p2,q1,q2

63、,line_p[9][9],line_q[9][9];　　for(i=0;i<9;i++)　　for(j=0;j<9;j++)　　{　　line_p[i][j]=line_q[i][j]=0;<

64、;p>　　}　　fp=fopen(line_S,"w");　　fprintf(fp,"------ 線上的網(wǎng)損功率----- ");　　fprintf(fp,"\n\n");<p&

65、gt;　　for(i=0;i<LINE_TOTAL_NUM;i++)　　{ 　　tSrcNode = gLineData[i].SrcNode - 1;//counting from zero　　tDestNode = gLineData[i].DestNod

66、e - 1;　　p1= line_g[i]*(gNodeVoltage[tSrcNode].Voltage*gNodeVoltage[tSrcNode].Voltage-gNodeVoltage[tSrcNode].Voltage*gNodeVoltage[tDestNode].Voltage*cos(gNodeVoltage[tSrcNode].Angle-gNodeVoltage[tDest

67、Node].Angle))-gNodeVoltage[tSrcNode].Voltage*gNodeVoltage[tDestNode].Voltage* line_b[i]*sin(gNodeVoltage[tSrcNode].Angle-gNodeVoltage[tDestNode].Angle);　　q1= -(line_b[i]*(gNodeVoltage[tSrcNode].Volta

68、ge*gNodeVoltage[tSrcNode].Voltage-gNodeVoltage[tSrcNode].Voltage*gNodeVoltage[tDestNode].Voltage*cos(gNodeVoltage[tSrcNode].Angle-gNodeVoltage[tDestNode].Angle))+gNodeVoltage[tSrcNode].Voltage*gNodeVoltage[tDestNode].Vol

69、tage* line_g[i]*sin(gNodeVoltage[tSrcNode].Angle-gNodeVoltage[tDestNode].Angle)+gNodeVoltage[tSrcNode].Voltage*gNodeVoltage[tSrcNode].Voltage* line_b0[i]);　　p2= line_g[i]*(gNodeVoltage[tDestNode].Vol

70、tage*gNodeVoltage[tDestNode].Voltage-gNodeVoltage[tSrcNode].Voltage*gNodeVoltage[tDestNode].Voltage*cos(gNodeVoltage[tDestNode].Angle-gNodeVoltage[tSrcNode].Angle))-gNodeVoltage[tSrcNode].Voltage*gNodeVoltage[tDestNode].

71、Voltage* line_b[i]*sin(gNodeVoltage[tDestNode].Angle-gNodeVoltage[tSrcNode].Angle);　　q2= -(line_b[i]*(gNodeVoltage[tDestNode].Voltage*gNodeVoltage[tDestNode].Voltage-gNodeVoltage[tSrcNode].Voltage*gN

72、odeVoltage[tDestNode].Voltage*cos(gNodeVoltage[tDestNode].Angle-gNodeVoltage[tSrcNode].Angle))+gNodeVoltage[tSrcNode].Voltage*gNodeVoltage[tDestNode].Voltage* line_g[i]*sin(gNodeVoltage[tDestNode].Angle-gNodeVoltage[tSrc

73、Node].Angle)+gNodeVoltage[tDestNode].Voltage*gNodeVoltage[tDestNode].Voltage* line_b0[i]);　　line_deltap[i]=p1+p2;　　line_deltaq[i]=q1+q2;　　line_p[tSrcNode][tD

74、estNode]=p1;　　line_p[tDestNode][tSrcNode]=p2;　　line_q[tSrcNode][tDestNode]=q1;　　line_q[tDestNode][tSrcNode]=q2;　　fprintf(fp," 第%d 條

75、 %6.3f + j%6.3f ",i+1,line_deltap[i],line_deltaq[i]);　　fprintf(fp,"\n");　　}　　fprintf(fp,"\n\n");<

76、;p>　　fclose(fp);　　fp=fopen(line_pq,"w");　　fprintf(fp,"------ 線路潮流----- ");　　fprintf(fp,"\n\n");　　fprint

77、f(fp,"有功 1 2 3 4 5 6 7 8 9\n");　　for(i=0;i<9;i++)　　{　　fprintf(fp,&

78、quot; %d ", i+1 );　　for(j=0;j<9;j++)　　{　　fprintf(fp," %6.3f ",line_p[i][j]);　　}</b&g

79、t;　　fprintf(fp,"\n");　　}　　fprintf(fp,"\n\n");　　fprintf(fp," 無(wú)功 1 2 3

80、4 5 6 7 8 9\n");　　for(i=0;i<9;i++)　　{　　fprintf(fp," %d ", i+1 );

81、;　　for(j=0;j<9;j++)　　{　　fprintf(fp," %6.3f ",line_q[i][j]);　　}　　fprintf(fp,"\n&

82、quot;);　　}　　fprintf(fp,"\n\n");　　fclose(fp);　　float tempP[NODE_TOTAL_NUM]={0.0},tempQ[NODE_TOTAL_NUM]={0.0}

83、; //設(shè)兩個(gè)變量數(shù)組用于存儲(chǔ)兩個(gè)累加和　　for(i=0;i<9;i++)　　{　　for(j=0;j<9;j++)　　{<

84、;/p>　　tempP[i]+=gNodeVoltage[j].Voltage*(gY_G[i][j]*cos(gNodeVoltage[i].Angle-gNodeVoltage[j].Angle)+ 　　gY_B[i][j]*sin(gNodeVoltage[i].Angle-gNodeVoltage[j].Angle));</

85、p>　　tempQ[i]+=gNodeVoltage[j].Voltage*(gY_G[i][j]*sin(gNodeVoltage[i].Angle-gNodeVoltage[j].Angle)-　　gY_B[i][j]*cos(gNodeVoltage[i].Angle-gNodeVoltage[j].Angle));&l

86、t;b>　　}　　}　　for(i=0;i<9;i++) //利用公式求解每個(gè)節(jié)點(diǎn)的不平衡功率　　{</p

87、>　　if(gNodeData[i].Type == 1) //如果是PQ節(jié)點(diǎn)則△P和△Q都要求　　{　　gDeltaPQ[i].DeltaP=gNodeData[i].FirstInput-gNodeVoltage

88、[i].Voltage*tempP[i];　　gDeltaPQ[i].DeltaQ=gNodeData[i].SecondInput-gNodeVoltage[i].Voltage*tempQ[i];　　node_P[i]=gNodeData[i].FirstInput;　　node_Q[i]=gNodeData

89、[i].SecondInput;　　}　　if(gNodeData[i].Type==2) //如果是PV節(jié)點(diǎn)只需求△P　　{<p&g

90、t;　　gDeltaPQ[i].DeltaP=gNodeData[i].FirstInput-gNodeVoltage[i].Voltage*tempP[i];　　node_P[i]=gNodeData[i].FirstInput;　　node_Q[i]=gNodeVoltage[i].Voltage*tempQ[i];<p

91、>　　}　　if(gNodeData[i].Type==0) 　　{　　node_P[i]=gNodeVoltage[i].Voltage*tempP[i];　　node_Q[i]=gNo

92、deVoltage[i].Voltage*tempQ[i];　　}　　}　　float B_Array[14];　　for(i=0;i<NODE_TOTAL_NUM;i++)

93、 //向數(shù)組中放入不平衡功率的值　　{　　if(gNodeData[i].Type==1) //PQ節(jié)點(diǎn)兩個(gè)都要輸入　　{

94、　　B_Array[i]=gDeltaPQ[i].DeltaP; 　　B_Array[i+NODE_TOTAL_NUM-1]=gDeltaPQ[i].DeltaQ; 　　}　　if(gNodeData[i].Type==2)

95、 //PV節(jié)點(diǎn)只輸入△P　　B_Array[i]=gDeltaPQ[i].DeltaP; 　　if(gNodeData[i].Type==0) //平衡節(jié)點(diǎn)不輸入數(shù)據(jù)　　{　　}<

96、;/b>　　}　　fp=fopen(deltaPQ,"w"); 　　for(i=0;i<14;i++)　　{　　fpri

97、ntf(fp," %f ",B_Array[i]);　　fprintf(fp,"\n");　　}　　fprintf(fp,"\n\n");　　fclose(fp);</

98、p>　　float gJaccobi[14][14];　　float H[8][8],N[8][6],J[6][8],L[6][6];　　float s;　　for(i=0;i<8;i++)　　for

99、(j=0;j<8;j++)　　{　　{ s=0;　　for(k=0;k<9;k++)　　s=s+gNodeVoltage[k].Voltage*(gY_G[i][k]*sin(gNo

100、deVoltage[i].Angle-gNodeVoltage[k].Angle)-gY_B[i][k]*cos(gNodeVoltage[i].Angle-gNodeVoltage[k].Angle));　　}　　if(i==j)<

101、;b>　　{ 　　H[i][i]=gNodeVoltage[i].Voltage*(s-gNodeVoltage[i].Voltage*(gY_G[i][i]*sin(gNodeVoltage[i].Angle-gNodeVoltage[i].Angle)-gY_B[i][i]*cos(gNodeVoltage[i].Angle-gNodeVoltage[i].Angl

102、e)));　　}　　else H[i][j]=-gNodeVoltage[i].Voltage*gNodeVoltage[j].Voltage*(gY_G[i][j]*sin(gNodeVoltage[i].Angle-gNodeVoltage[j].Angle)-gY_B[i][j]*cos(gNodeVolt

103、age[i].Angle-gNodeVoltage[j].Angle));　　}　　for(i=0;i<8;i++)　　for(j=0;j<6;j++)　　{

104、　　s=0;　　for(k=0;k<9;k++)　　{　　s=s+gNodeVoltage[k].Voltage*(gY_G[i][k]*cos(gNodeVoltage[i].Angle-gNodeVoltage

105、[k].Angle)+gY_B[i][k]*sin(gNodeVoltage[i].Angle-gNodeVoltage[k].Angle));　　}　　if(i==j)　　{

106、　　N[i][i]=-gNodeVoltage[i].Voltage*(s-gNodeVoltage[i].Voltage*(gY_G[i][i]*cos(gNodeVoltage[i].Angle-gNodeVoltage[i].Angle)+gY_B[i][i]*sin(gNodeVoltage[i].Angle-gNodeVoltage[i].Angle)))-2*gNodeVoltage[i].Voltag

107、e*gNodeVoltage[i].Voltage*gY_G[i][i];　　}　　else N[i][j]=-gNodeVoltage[i].Voltage*gNodeVoltage[j].Voltage*(gY_G[i][j]*cos(gNodeVoltage[i].Angle-gNodeVoltage[j]

108、.Angle)+gY_B[i][j]*sin(gNodeVoltage[i].Angle-gNodeVoltage[j].Angle));　　}　　for(i=0;i<6;i++)　　for(j=0;j<8;j++)<

109、;b>　　{　　s=0;　　for(k=0;k<9;k++)　　{ 　　s=s+gNodeVoltage[k].Voltage*(gY_G[i][k]*cos(gN

110、odeVoltage[i].Angle-gNodeVoltage[k].Angle)+gY_B[i][k]*sin(gNodeVoltage[i].Angle-gNodeVoltage[k].Angle));　　}　　if(i==j)&l

111、t;b>　　{ 　　J[i][i]=-gNodeVoltage[i].Voltage*(s-gNodeVoltage[i].Voltage*(gY_G[i][i]*cos(gNodeVoltage[i].Angle-gNodeVoltage[i].Angle)+gY_B[i][i]*sin(gNodeVoltage[i].Angle-gNodeVoltage[i].An

112、gle)));　　}　　else J[i][j]=gNodeVoltage[i].Voltage*gNodeVoltage[j].Voltage*(gY_G[i][j]*cos(gNodeVoltage[i].Angle-gNodeVoltage[j].Angle)+gY_B[i][j]*sin(gNodeVol

113、tage[i].Angle-gNodeVoltage[j].Angle));　　}　　for(i=0;i<6;i++)　　for(j=0;j<6;j++)　　{

114、　　s=0;　　for(k=0;k<9;k++)　　{ 　　s=s+gNodeVoltage[k].Voltage*(gY_G[i][k]*sin(gNodeVoltage[i].Angle-gNodeVoltag

115、e[k].Angle)-gY_B[i][k]*cos(gNodeVoltage[i].Angle-gNodeVoltage[k].Angle));　　}　　if(i==j)　　{

116、　　L[i][i]=-gNodeVoltage[i].Voltage*(s-gNodeVoltage[i].Voltage*(gY_G[i][i]*sin(gNodeVoltage[i].Angle-gNodeVoltage[i].Angle)-gY_B[i][i]*cos(gNodeVoltage[i].Angle-gNodeVoltage[i].Angle)))+2*gNodeVoltage[i].Voltag

117、e*gNodeVoltage[i].Voltage*gY_B[i][i];　　}　　else L[i][j]=-gNodeVoltage[i].Voltage*gNodeVoltage[j].Voltage*(gY_G[i][j]*sin(gNodeVoltage[i].Angle-gNodeVoltage[j]

118、.Angle)-gY_B[i][j]*cos(gNodeVoltage[i].Angle-gNodeVoltage[j].Angle));　　}　　for(i=0;i<8;i++)　　{

119、for(j=0;j<8;j++)　　gJaccobi[i][j]=H[i][j];　　}　　for(i=0;i<8;i++)　　{　　for(

120、j=8;j<14;j++)　　gJaccobi[i][j]=N[i][j-8];　　}　　for(i=8;i<14;i++)　　{　　for(

121、j=0;j<8;j++)　　gJaccobi[i][j]=J[i-8][j];　　}　　for(i=8;i<14;i++)　　{　　for(j

122、=8;j<14;j++)　　gJaccobi[i][j]=L[i-8][j-8];　　}　　fp=fopen(JACCOBI,"w"); //write the Jmatrix data to file

123、fprintf(fp," ----------%d 次迭代之后的雅各比---------- ", m+1);　　fprintf(fp,"\n");　　for(i=0;i<14;i++)　　{

124、;　　for(j=0;j<14;j++)　　fprintf(fp," %f ", gJaccobi[i][j]);　　fprintf(fp,"\n");　　}　　fclose(fp);&l

125、t;/p>　　fp=fopen(U,"w"); 　　fprintf(fp,"----- %d次迭代之后的節(jié)點(diǎn)電壓和相角-----" ,m+1);　　fprintf(fp,"\n\n");　　for(i=0;i<9;i

126、++)　　{　　fprintf(fp," %6.3f ,%6.3f ",gNodeVoltage[i].Voltage,gNodeVoltage[i].Angle*180/3.14 );　　fprintf(fp,"\n");</

127、p>　　}　　fclose(fp);　　float *Jgb[14]={0};　　for(i=0;i<14;i++)　　Jgb[i]=gJaccobi[i];

128、　　SolveEquation(14,Jgb,B_Array);　　for(i=0;i<8;i++)　　gNodeVoltage[i].Angle=gNodeVoltage[i].Angle-B_Array[i];　　for(i=0;i<6;i++)　　gNo

眾賞文庫(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ì)自己和他人造成任何形式的傷害或損失。

離線潮流計(jì)算課程設(shè)計(jì)

文檔簡(jiǎn)介

溫馨提示

最新文檔

評(píng)論

離線潮流計(jì)算課程設(shè)計(jì)

文檔簡(jiǎn)介

溫馨提示

最新文檔

評(píng)論

免費(fèi)下載