潮流計(jì)算上機(jī)報(bào)告華電

潮流計(jì)算上機(jī)報(bào)告華電院系：電氣與電子工程學(xué)院劉寶柱課程設(shè)計(jì)的目的與要求培養(yǎng)學(xué)生的電力系統(tǒng)潮流計(jì)算機(jī)編程能力，掌握計(jì)算機(jī)潮流計(jì)算的相關(guān)知識二、設(shè)計(jì)正文(詳細(xì)內(nèi)容見附錄)地支路30。y地支路30。Y(1,1)=(,Y(1,2)=(,Y(1,3)=(,Y(1,4)=(,Y(2,1)=(,Y(2,2)=(,Y(2,4)=(,Y(3,1)=(,Y(3,3)=(,Y(4,1)=(,Y(4,2)=(,Y(4,4)=(,112233P(0)=PP(0)=P[e(0)4(Ge(0)Bf(0))+f(0)4(Gf(0)11s11s11jj1j111jjj=1j=111s11s11jj1j111jjj=1j=11j11j1P(0)=PP(0)=0.52596Q(0)=QQ(0)=0.019622s2,22s2P(0)=PP(0)=0.5U(0)2=U2U(0)2=0.033s3,33s34．根據(jù)求的第一次迭代時(shí)雅可比矩陣各元素的公式計(jì)算雅可比矩陣各個(gè)元素的具體「8.32104-1.03803-2.302080.56148-4.03334L00.992658.06298-0.56148-2.3020800-2.302080.561484.82612-1.0662900-0.56148-3.66667-2.3020801.0182104.52691003.66667000]-3.666670002.20000」33：e(1)=e(0)+Ae(0)=0.993491e(1)=e(0)+Ae(0)=0.9764322e(1)=e(0)+Ae(0)=1.1000033111222333返回第三步重新迭代，并校驗(yàn)收斂與否，令:=10-4。經(jīng)過四次迭代后，收斂條件滿U=0.98459-j0.007641U=0.95859-j0.109102U=1.09229-j0.130023U=1.05000-j0.0000044S~=-0.24437-j0.00094S~=0.04437-j0.11804S~=0.23821-j0.0243221S~=0.31179+j0.1269424S~=-0.04634-j0.1262441S~=-0.32153-j0.17564429．網(wǎng)損為:S~=0.01787-j0.04487何特點(diǎn)多慢少高少高多軟件計(jì)算北京城市電網(wǎng)的潮流，你應(yīng)該做哪些工作(收集哪些數(shù)據(jù)，如何整理，計(jì)算結(jié)果如何分析)b.對于平衡節(jié)點(diǎn)要了解節(jié)點(diǎn)的電壓大小相位、及節(jié)點(diǎn)所能提供的最大最小有功無功功率c.PV節(jié)點(diǎn)要知道節(jié)點(diǎn)電壓大小注入有功功率及節(jié)點(diǎn)所能提供的最大和最小無功功.率d.PQ節(jié)點(diǎn)要知道節(jié)點(diǎn)的注入有功和無功功率B所有支路的數(shù)據(jù)：b.各支路的電阻、電感、電納c.各變壓器的變比。整理的結(jié)果寫成本實(shí)驗(yàn)中所要求的格式(原始數(shù)據(jù)的txt文檔)，再用本實(shí)驗(yàn)所編制的程序進(jìn)行求解，得到各節(jié)點(diǎn)電壓、相位，各線路傳輸功率、損耗，平衡節(jié)點(diǎn)注入功率等數(shù)值。3.計(jì)算結(jié)果分析：可以大致了解整個(gè)電網(wǎng)的運(yùn)行狀況并按照電力系統(tǒng)知識依次進(jìn)行所需的調(diào)整。如調(diào)節(jié)各節(jié)點(diǎn)電壓的大小以防止超出規(guī)定限度，調(diào)整各線路的有功無功分布狀況分析整個(gè)系統(tǒng)的網(wǎng)損是否達(dá)到標(biāo)準(zhǔn)設(shè)計(jì)中遇到的問題和解決的辦法。1.在設(shè)計(jì)最開始遇到的問題就是由于以前C語言所學(xué)習(xí)的知識遺忘較多，對于程序的使用不太熟練，通過老師講解、看書、詢問同學(xué)逐漸掃除了這個(gè)2.設(shè)計(jì)的時(shí)候?qū)τ陔娏ο到y(tǒng)整體的知識脈絡(luò)不是很清晰，編寫程序時(shí)所定義的變量容易弄混淆，編寫的程序產(chǎn)生了較多的錯(cuò)誤，經(jīng)過仔細(xì)學(xué)習(xí)和對每時(shí)也對潮流計(jì)算的過程有了課程設(shè)計(jì)總結(jié)或結(jié)論該門課程的主要目的是培養(yǎng)學(xué)生的電力系統(tǒng)潮流計(jì)算機(jī)編程能力，掌握計(jì)算機(jī)潮流計(jì)算的相關(guān)知識。潮流計(jì)算是電力系統(tǒng)非常重要的分析計(jì)算，是電力系統(tǒng)最基本的計(jì)算，主要用以研究系統(tǒng)規(guī)劃和運(yùn)行中提出的各種問題。現(xiàn)代電力系統(tǒng)潮流就算過程比較復(fù)雜，手算有很大的困難，而計(jì)算機(jī)作為一種處理數(shù)據(jù)的工具，具有計(jì)算速度快，準(zhǔn)確率高的優(yōu)點(diǎn)，因此現(xiàn)在電力系統(tǒng)的潮流大多也都采用計(jì)算機(jī)的算法求解。對于寫程序而言，C語言是一種簡單但很實(shí)用的語言，利用C語言編程實(shí)現(xiàn)潮流計(jì)算不僅簡單易行，而且還能加深對電力系統(tǒng)潮流這一物理感念的理解。通過這次潮流對c++的面向?qū)ο缶幊碳捌溥^程有了更深入的了解。對于潮流計(jì)算，我更深刻的掌握了牛拉法解潮流方程的方法以及對數(shù)學(xué)高次多元方程線性化在迭代求解的方法有了更深入的認(rèn)識。通過這次課程設(shè)計(jì)，自己動(dòng)手編寫潮流計(jì)算的程序，仔細(xì)分析了潮流計(jì)算各個(gè)步驟，理論與實(shí)踐相結(jié)合,對潮流計(jì)算的思想和方法以及過程有了更深入的理解。在這兩周的學(xué)習(xí)里，我不僅對C語言的學(xué)習(xí)有了更深層次的理解，還對自己在《電力系統(tǒng)分析》這門課程里所學(xué)到的知識進(jìn)行了查缺補(bǔ)漏，對電力系統(tǒng)的運(yùn)附錄(設(shè)計(jì)流程圖、程序、表格、數(shù)據(jù)等)設(shè)非平衡節(jié)點(diǎn)電壓初值e(0)iii客比矩陣各元素客比矩陣各元素ijijijijijijiiiiimaxmax迭代是否收斂maxmax的新電壓iii計(jì)算平衡節(jié)點(diǎn)的功率及線路2.程序設(shè)計(jì):includeincludeincludeinclude#defineNBUS5#defineNLINE7/*Globalvariables*/VAbledt{R{intNum;doubleVolt,Phase,GenP,GenQ,LoadP,LoadQ;{GBvoidmain){FILE*fp;FILE*fpout;doubledddddddtrxgbttLLe,ps,qs,shsh,sum,w;USNBUSYBNBUSNBUSV2],U[100],dPQ[100],PQ[100][2],JJ[100][100]={0},H[100][100],J[100][100],N[LsPNBUSNBUSsQNBUSNBUSdspdsqsumgen,dp;psqsdtheinputdata}fscanffpdd%d",&nBUS,&nL,&nSH);iinBUSiBusiNumsBusiTypesBusiVoltsBusiPhasesBus[i].GenP=sBus[i].GenQ=sBus[i].LoadP=sBus[i].LoadQ=;fscanf(fp,"%d,%lf,%lf,%lf,%lf,%lf,%lf,%d",&i1,&d1,&d2,&d3,&d4,&d5,&d6,&i2);sBusi.Num=i1;sBus[i].Volt=d1;sBus[i].Phase=d2;sBus[i].GenP=d3;sBus[i].GenQ=d4;sBusiLoadPdsBusiLoadQd;sBus[i].Type=i2;iinLiumsLiNumIsLiNumJsLiRsLiXsLiBsLiKfscanf(fp,"%2d%3d%3d%lf%lf%lf%lf",&i1,&i2,&i3,&d1,&d2,&d3,&d4);sLi.Num=i1;sL[i].NumI=i2;sL[i].NumJ=i3;sL[i].R=d1;sL[i].X=d2;sL[i].B=d3;sL[i].}iinSHisSHiNumsSHiNumIsSHiGsSHi].B=;fscanffpd3d%lf",&i1,&i2,&d1);sSHiNumisSHiNumIisSHiBd1;}trixforiinBUSiforjjnBUSj+){ijijforllnLl{sLlRxdYGiiYGii]+g;YGjjYGjj]+g;YBiiYBii+b+sL[l].B;YBjjYBjj+b+sL[l].B;YG[i][j]=YG[i][j]-g;YG[j][i]=YG[j][i]-g;YBijYBij-b;YBjiYBji-b;}{YG[i][i]=YG[i][i]+g/sL[l].B/sL[l].B;YB[i][i]=YB[i][i]+b/sL[l].B/sL[l].B;YGjj=YG[j][j]+g;YBjj=YB[j][j]+b;YGi[j]=YG[i][j]-g/sL[l].B;YGj[i]=YG[j][i]-g/sL[l].B;YBi[j]=YB[i][j]-b/sL[l].B;YBj[i]=YB[j][i]-b/sL[l].B;}}}{Vi]=;Vi]=;}inBUSi{Vi][1]=sBus[i-1].Volt;Vi]=;}inBUSi{Vi][1]=sBus[i-1].Volt;Vi][0]=sBus[i-1].Phase;}ndlxtxtwNULL{printf"Cannotopenthefilenamed'電壓初值.txt'\n");}fprintffp-\n");inBUSifprintffpYd(%,%\n",i,V[i][0],V[i][1]);{量*/Si{{PQiPQiSj{PQiVi[1]*V[j][1]*(YG[i][j]*cos(V[i][0]-Vj+YB[i][j]*sin(V[i][0]-V[j][0]));PQiVi[1]*V[j][1]*(YG[i][j]*sin(V[i][0]-Vj-YB[i][j]*cos(V[i][0]-V[j][0]));}}}Si{{PQhPQisBusiGenPsBusiLoadPPQhPQisBusiGenQsBusiLoadQ}{PQhPQisBusiGenPsBusiLoadP}}txtwNULL{printfCannotopenthefilenamed偏移量.txt'\n");}fprintffp\n");fprintffpdPQdnidPQi克比矩陣*/Si{Sj{{HijPQi][1]-V[i][1]*V[i][1]*YB[i][j];Jij]=-PQ[i][0]-V[i][1]*V[i][1]*YG[i][j];Nij]=-PQ[i][0]+V[i][1]*V[i][1]*YG[i][j];Lij]=-PQ[i][1]-V[i][1]*V[i][1]*YB[i][j];}{HijVi][1]*V[j][1]*(YG[i][j]*sin(V[i][0]-V[j][0])-YBijcosVi0]-V[j][0]));JijViVj][1]*(YG[i][j]*cos(V[i][0]-Vj+YB[i][j]*sin(V[i][0]-V[j][0]));Nij]=V[i][1]*V[j][1]*(YG[i][j]*cos(V[i][0]-Vj+YB[i][j]*sin(V[i][0]-V[j][0]));LijViVj][1]*(YG[i][j]*sin(V[i][0]-V[j][0])-YBijcosVi0]-V[j][0]));}}}Si{{Sj{pe{hlHijhlNij}pe{hlHij}}Sj{pe{hlJijhlLij}pe{hlJij}}}{Sj{pe{hlHijhlNij}pe{hlHij}}}}lxtxtwNULL{printfCannotopenthefilenamed'雅克比矩陣.txt'\n");}fprintffp--\n");fprintffp");{fpdtti}{fprintffpd,i);{fptJJij}}{}{{d{}}{}{{JikJJiLHkHkt}}{}{JJkjJJkjJJkk];}kk{{JJijJJijJJi[k]*JJ[k][j];/*使下三角陣的元}iJJikdPQk}}U[h-1]=dPQ[h-1];forihi>=0;i--){{}U[i]=(dPQ[i]-t);}}twNULL{rintferrexit}{fUi}Si{{Vi=U[h];Vi=U[h];}{Vi=U[h];}}xtwNULL{printfCannotopenthefilenamed電壓值.txt'\n");}fprintffpn");SifprintffpYdniVi0],V[i][1]);w;{{wUi];}}{}}xtxtwNULL{printfCannotopenthefilenamed'電壓終值.txt'\n");}fprintffp\n");SifprintffpYdniVi0]*180/,V[i][1]);Si{{Sj{psViVj[1]*(YG[i][j]*cos(V[i][0]-Vj+YB[i][j]*sin(V[i][0]-V[j][0]));qsViVjYGij*sin(V[i][0]-V[j][0])-YBijcosVi0]-V[j][0]));}}}Si{Sj{}}forllnLl{sLlRxdsPijViVigVi][1]*V[j][1]*(g*cos(V[i][0]-Vj+b*sin(V[i][0]-V[j][0]));jgsinViVj-b*cos(V[i][0]-V[j][0])));sPjiVjVjgVi][1]*V[j][