電分實(shí)驗(yàn)報(bào)告

實(shí)驗(yàn)報(bào)告實(shí)驗(yàn)課程：電力系統(tǒng)分析（下）學(xué)生姓名：學(xué)號(hào)：專業(yè)班級(jí)：電氣152班2018年6月17日南昌大學(xué)實(shí)驗(yàn)報(bào)告學(xué)生姓名：學(xué)號(hào)：_專業(yè)班級(jí)：電氣152實(shí)驗(yàn)類型：口驗(yàn)證口綜合■設(shè)計(jì)口創(chuàng)新實(shí)驗(yàn)日期：實(shí)驗(yàn)成績(jī)：一、實(shí)驗(yàn)名稱配電網(wǎng)潮流計(jì)算實(shí)驗(yàn)二、實(shí)驗(yàn)?zāi)康模罕緦?shí)驗(yàn)針對(duì)配電網(wǎng)具有輻射式開環(huán)結(jié)構(gòu)的特性，基于前推回代的方法編制程序使系統(tǒng)潮流計(jì)算能夠由計(jì)算機(jī)自行完成，即已知配電網(wǎng)首端節(jié)點(diǎn)電壓和各末端節(jié)點(diǎn)的功率，由末端向首端計(jì)算網(wǎng)絡(luò)功率分布，再由首端向末端計(jì)算節(jié)點(diǎn)電壓分布，迭代直至收斂條件。通過(guò)實(shí)驗(yàn)教學(xué)加深學(xué)生對(duì)配電網(wǎng)潮流計(jì)算方法的理解，掌握樹狀輻射式配電網(wǎng)潮流計(jì)算的基本算法。熟悉各種常用應(yīng)用軟件，熟悉硬件設(shè)備的使用方法，加強(qiáng)編制調(diào)試計(jì)算機(jī)程序的能力，提高工程計(jì)算的能力，學(xué)習(xí)如何將理論知識(shí)和實(shí)際工程問(wèn)題結(jié)合起來(lái)。三、實(shí)驗(yàn)內(nèi)容：編制調(diào)試配電網(wǎng)潮流計(jì)算的計(jì)算機(jī)程序。程序要求根據(jù)已知的配電網(wǎng)參數(shù)，完成該配電系統(tǒng)的潮流計(jì)算，要求計(jì)算出節(jié)點(diǎn)電壓、功率等參數(shù)。先編制好的電力系統(tǒng)潮流計(jì)算的計(jì)算程序原代碼由自備移動(dòng)存儲(chǔ)設(shè)備導(dǎo)入計(jì)算機(jī)。應(yīng)的編程環(huán)境下對(duì)程序進(jìn)行組織調(diào)試。應(yīng)用計(jì)算例題驗(yàn)證程序的計(jì)算效果。對(duì)調(diào)試正確的計(jì)算程序進(jìn)行存儲(chǔ)、打印。完成本次實(shí)驗(yàn)的實(shí)驗(yàn)報(bào)告。S4四、實(shí)驗(yàn)數(shù)據(jù)

如圖所示一個(gè)5節(jié)點(diǎn)的配電網(wǎng)系統(tǒng)，S1=S2=2+j2(MV?A),S3=S4=S5=4+j4(MV?A),假定所有負(fù)荷均為恒功率負(fù)荷，節(jié)點(diǎn)1為參考節(jié)點(diǎn)相角為0。計(jì)算電網(wǎng)個(gè)節(jié)點(diǎn)電壓及支路功率分布，收斂條件為e<10-6。五、程序框圖[0.54+0.65j;0.62+0.5j;0.6+0.35j;0.72+0.75j;1.0+0.55j;23456六、實(shí)驗(yàn)程序functionDSPowerFlowCalc(linedata,BusPower,U)clcifnargin<3linedata=1232557[0.54+0.65j;0.62+0.5j;0.6+0.35j;0.72+0.75j;1.0+0.55j;23456BusPower=[00.6+0.45j0.4+0.3j0.4+0.28j0.6+0.4j0.4+0.3j0.5+0.35j0.5+0.4j];%節(jié)點(diǎn)負(fù)荷功率U=[10.510101010101010];%節(jié)點(diǎn)電壓end[S,U,Zm,tailbus,headbus,branchnum]=DataProcess(linedata,BusPower,U);forn=1:100fprintfd第%d輪迭代結(jié)果\n',n);fprintf(各支路的首段功率S/(kV?A)\n');S=forward(S,U,Zm,tailbus,headbus,branchnum);U0=U;fprintf('各節(jié)點(diǎn)電壓U/kV\n');U=backward(S,U,Zm,tailbus,headbus,branchnum);%收斂判據(jù)ifmax(abs(U-U0)/U(1))<1e-4break;endendendfunction[S,U,Zm,tailbus,headbus,branchnum]=DataProcess(linedata,BusPower,U)%數(shù)據(jù)預(yù)處理branchnum=length(linedata(:,1));busnum=max([linedata(:,1)',linedata(:,2)']);S=zeros(busnum);forn=1:busnumS(n,n)=BusPower(n);endiflength(U)<2%當(dāng)U只給出了電源點(diǎn)的電壓時(shí)，%根據(jù)電源點(diǎn)電壓高于其他節(jié)點(diǎn)電壓5%,計(jì)算其他節(jié)點(diǎn)的初始電壓U(2:busnum)=ones(1,busnum-1)*(U(1)/1.05);endZm=zeros(busnum);%Zm(i,j)表示節(jié)點(diǎn)ij之間的阻抗linknum=zeros(1,busnum);forn=1:branchnumtemp1=linedata(n,1);temp2=linedata(n,2);%將各支路參數(shù)保存到一個(gè)矩陣中Zm(temp1,temp2)=linedata(n,3);Zm(temp2,temp1)=Zm(temp1,temp2);%記錄各個(gè)節(jié)點(diǎn)所連接的支路數(shù)linknum(temp1)=linknum(temp1)+1;linknum(temp2)=linknum(temp2)+1;end%[tailbus,headbus]=PathPlan(Zm,linknum,branchnum);[tailbus,headbus]=PathPlan1(linedata,busnum,branchnum);endfunction[tailbus,headbus]=PathPlan(Zm,linknum,branchnum)%通過(guò)逐條拆除支路的方法確定支路計(jì)算的順序%tailbus(i)表示前推過(guò)程中，計(jì)算的第i條支路的尾節(jié)點(diǎn)%headbus(i)表示前推過(guò)程中，計(jì)算的第i條支路的首節(jié)點(diǎn)%回代時(shí)路徑剛好相反tailbus=zeros(1,branchnum);headbus=tailbus;form=1:branchnumtailbus(m)=find(linknum(2:end)==1,1)+1;headbus(m)=find(Zm(tailbus(m),:)~=0,1);linknum(tailbus(m))=linknum(tailbus(m))-1;linknum(headbus(m))=linknum(headbus(m))-1;endendfunction[tailbus,headbus]=PathPlan1(linedata,busnum,branchnum)tailbus=zeros(1,branchnum);headbus=tailbus;linknum=zeros(1,busnum);%首先追加從根節(jié)點(diǎn)headbus(branchnum)=1;linknum(1)=linknum(1)+1;temp1=find(linedata(:,1)==1);temp2=find(linedata(:,2)==1);if~isempty(temp1)tailbus(branchnum)=linedata(temp1(1),2);linknum(linedata(temp1(1),2))=linknum(linedata(temp1(1),2))+1;linedata(temp1(1),:)=[];elseif~isempty(temp2)tailbus(branchnum)=linedata(temp2(1),1);linknum(linedata(temp2(1),1))=linknum(linedata(temp2(1),1))+1;linedata(temp2(1),:)=[];endforn=2:branchnumnowbranchnum=length(linedata(:,1));form=1:nowbranchnum%linknum(linedata(m,1))~=0,說(shuō)明該支路的一個(gè)節(jié)點(diǎn)已經(jīng)存在路徑中iflinknum(linedata(m,1))~=0headbus(branchnum-n+1)=linedata(m,1);tailbus(branchnum-n+1)=linedata(m,2);linknum(linedata(m,1))=linknum(linedata(m,1))+1;linknum(linedata(m,2))=linknum(linedata(m,2))+1;linedata(m,:)=[];break%跳出本次遍尋elseiflinknum(linedata(m,2))~=0headbus(branchnum-n+1)=linedata(m,2);tailbus(branchnum-n+1)=linedata(m,1);linknum(linedata(m,1))=linknum(linedata(m,1))+1;linknum(linedata(m,2))=linknum(linedata(m,2))+1;linedata(m,:)=[];breakendendendendfunctionS=forward(S,U,Zm,tailbus,headbus,branchnum)%前推過(guò)程，從葉節(jié)點(diǎn)向電源點(diǎn)計(jì)算form=1:branchnump=tailbus(m);q=headbus(m);dS=(abs(sum(S(p,:)))/U(p))A2*Zm(q,p);S(q,p)=sum(S(p,:))+dS;fprintf('%s%s','a'+q-1,'a'+p-1)disp(S(q,p))endendfunctionU=backward(S,U,Zm,tailbus,headbus,branchnum)%回代過(guò)程，從電源點(diǎn)向葉節(jié)點(diǎn)計(jì)算form=branchnum:-1:1q=headbus(m);%首節(jié)點(diǎn)p=tailbus(m);%尾節(jié)點(diǎn)DU=(real(S(q,p))*real(Zm(q,p))+imag(S(q,p))*imag(Zm(q,p)))/U(q);dU=(real(S(q,p))*imag(Zm(q,p))-imag(S(q,p))*real(Zm(q,p)))/U(q);U(p)=sqrt((U(q)-DU)A2+dUA2);fprintf('%s*,*a*+p-1)disp(U(p))endend七、實(shí)驗(yàn)結(jié)果

第1輪迭代結(jié)果各支路的首段功率S/(kV?A)hh(J.5037?0一4021ieg0.5024+0.3513ie10.1025+0.3011ib(?15292?L()780icd040L40.2fiOSibc0.80750.5857icib3.5385+2.6338i各節(jié)點(diǎn)電壓U/kVb10.1553c10.0772d10.0435c9674f9.9103g9.9223h10.0909第2掄迭代姑果各支路的首段功率bhTOC\o"1-5"\h\zerbecdbc己b各節(jié)點(diǎn)電壓U/kVb10.1557c10.0776d10.0439TOC\o"1-5"\h\ze9.9677I9.9JO7r土92對(duì)h10.0913fx?S/(kV?A)0.5036+S/(kV?A)0.5036+0.402010.5025+0.351310一，1。2|5+0.3014I1.5295+I.078210.4014+0,280810.8074+0.585713.5356+2.6304i針對(duì)輻射狀配電網(wǎng)的的特點(diǎn)，應(yīng)用前推回代法進(jìn)行潮流計(jì)算。同時(shí)利用葉節(jié)點(diǎn)只與一條支路聯(lián)接的特點(diǎn)來(lái)確定前推回代的支路順序。該方法可以處理多分支輻射網(wǎng)，無(wú)需復(fù)雜的的網(wǎng)絡(luò)編號(hào)、不用形成導(dǎo)納矩陣，迭代次數(shù)少，結(jié)果精確。九、調(diào)試記錄利用fprintf輸出結(jié)果時(shí)不能輸出電壓的虛部，改成disp可以完整輸出電壓。南昌大學(xué)實(shí)驗(yàn)報(bào)告學(xué)生姓名：學(xué)號(hào)：—專業(yè)班級(jí)：電氣152實(shí)驗(yàn)類型：口驗(yàn)證口綜合■設(shè)計(jì)口創(chuàng)新實(shí)驗(yàn)日期：實(shí)驗(yàn)成績(jī)：一、實(shí)驗(yàn)題目：輸電網(wǎng)潮流計(jì)算實(shí)驗(yàn)二、實(shí)驗(yàn)?zāi)康模罕緦?shí)驗(yàn)針對(duì)高壓輸電網(wǎng)絡(luò)具有多環(huán)、多電源的結(jié)構(gòu)特性，要求根據(jù)功率平衡方程的泰勒展開式近似的求取系統(tǒng)的狀態(tài)變量，使系統(tǒng)潮流計(jì)算能夠由計(jì)算機(jī)自行完成，即根據(jù)已知的電力網(wǎng)的數(shù)學(xué)模型（節(jié)點(diǎn)導(dǎo)納矩陣）及各節(jié)點(diǎn)參數(shù)，由計(jì)算程序運(yùn)行完成該電力系統(tǒng)的潮流計(jì)算。通過(guò)實(shí)驗(yàn)教學(xué)加深學(xué)生對(duì)復(fù)雜電力系統(tǒng)潮流計(jì)算計(jì)算方法的理解，學(xué)會(huì)運(yùn)用電力系統(tǒng)的數(shù)學(xué)模型，掌握潮流計(jì)算的過(guò)程及其特點(diǎn)，熟悉各種常用應(yīng)用軟件，熟悉硬件設(shè)備的使用方法，加強(qiáng)編制調(diào)試計(jì)算機(jī)程序的能力，提高工程計(jì)算的能力，學(xué)習(xí)如何將理論知識(shí)和實(shí)際工程問(wèn)題結(jié)合起來(lái)。三、實(shí)驗(yàn)內(nèi)容：編制調(diào)試電力系統(tǒng)潮流計(jì)算的計(jì)算機(jī)程序。程序要求根據(jù)已知的電網(wǎng)的數(shù)學(xué)模型（節(jié)點(diǎn)導(dǎo)納矩陣）及各節(jié)點(diǎn)參數(shù)，完成該電力系統(tǒng)的潮流計(jì)算，要求計(jì)算出節(jié)點(diǎn)電壓、功率等參數(shù)。米用牛頓拉夫遜法編制程序編制潮流程序。先編制好的電力系統(tǒng)潮流計(jì)算的計(jì)算程序原代碼由自備移動(dòng)存儲(chǔ)設(shè)備導(dǎo)入計(jì)算機(jī)。應(yīng)的編程環(huán)境下對(duì)程序進(jìn)行組織調(diào)試。應(yīng)用計(jì)算例題驗(yàn)證程序的計(jì)算效果。對(duì)調(diào)試正確的計(jì)算程序進(jìn)行存儲(chǔ)、打印。完成本次實(shí)驗(yàn)的實(shí)驗(yàn)報(bào)告。四、實(shí)驗(yàn)數(shù)據(jù)4^在上圖所示的簡(jiǎn)單電力系統(tǒng)中，系統(tǒng)中節(jié)點(diǎn)1、2為PQ節(jié)點(diǎn)，節(jié)點(diǎn)3為PV節(jié)點(diǎn)，節(jié)點(diǎn)4為平衡節(jié)點(diǎn)，已給定P1s+jQ1s=-0.30-j0.18P2s+jQ2s=-0.55-j0.13P3s=0.5V3s=1.10V4s=1.05Z0°容許誤差e=10-5。五、程序框圖六、實(shí)驗(yàn)程序functionPowerFlowNRRec(linedata,busdata,error)%節(jié)點(diǎn)標(biāo)號(hào)無(wú)特殊要求clcifnargin<1linedata=[140.10.40.015281;1300.301.1;120.120.50.019201;420.080.400.014131];busdata=[10-0.30-0.1810;40-0.55-0.1310;320.501.100;21001.050;];endifnargin<3error=10A-5;endclc%優(yōu)化節(jié)點(diǎn)排序pqbus=find(busdata(:,2)==0);pvbus=find(busdata(:,2)==2);vcbus=find(busdata(:,2)==1);bussort=[pqbus'pvbus'vcbus'];Num=busdata(bussort,1);%對(duì)節(jié)點(diǎn)進(jìn)行重新排列，使其排列順序?yàn)镻Q節(jié)點(diǎn)->PV節(jié)點(diǎn)->平衡節(jié)點(diǎn)linedata=NumResort(linedata,Num);P=busdata(:,3);Q=busdata(:,4);U=busdata(:,5);Us=U;deta=busdata(:,6);e=U.*cos(deta);f=U.*sin(deta);busnum=length(busdata(bussort,1));pqnum=length(pqbus);%生成節(jié)點(diǎn)導(dǎo)納矩陣Y=BuildY(linedata);fork=1:100[dP,dQ,dU2]=CalcCollection(P,Q,Us,e,f,Y,busnum,pqnum);ifmax(abs([dPdQ]))<errorbreak;endJ=CalcJ(e,f,Y,busnum,pqnum);PQU=zeros(busnum-1,1);forii=1:pqnumPQU(2*ii-1)=dP(ii);PQU(2*ii)=dQ(ii);endforii=pqnum+1:busnum-1PQU(2*ii-1)=dP(ii);PQU(2*ii)=dU2(ii);enddef=Jordan(J,PQU);forii=1:busnum-1f(ii)=f(ii)+def(2*ii-1);e(ii)=e(ii)+def(2*ii);endendifk>=100disp('不收斂！')；return;end[Sbus,Sij,Svc]=CalcS(linedata,e,f,Y,busnum);disp('編號(hào):');disp(Num');disp('節(jié)點(diǎn)電壓:');disp(e'+1j*f');disp('平衡節(jié)點(diǎn)功率：')；disp(Svc);disp('線路功率');forn=1:length(Sij(:,1))fprintf('%o%o',Num(Sij(n,1)),Num(Sij(n,2)))disp(Sij(n,3))enddisp('節(jié)點(diǎn)功率');forn=1:busnumfprintf('%o',Num(n))disp(Sbus(n))endendfunction[Sbus,Sij,Svc]=CalcS(linedata,e,f,Y,busnum)Sbus=zeros(1,busnum);G=real(Y);B=imag(Y);Qi=zeros(1,busnum);Ud=e+1j*f;%公式11-57forii=1:busnumforjj=1:busnumSbus(ii)=Sbus(ii)+Ud(ii)*conj(Y(ii,jj))*conj(Ud(jj));endendQiSvc=Sbus(busnum);Sij=zeros(length(linedata(:,1)),3);forn=1:length(linedata(:,1))ii=linedata(n,1);jj=linedata(n,2);Sij(2*n-1,1)=ii;Sij(2*n-1,2)=jj;Sij(2*n-1,3)=Ud(ii)*(conj(Ud(ii))...*conj(1j*linedata(n,5)+linedata(n,6)*(linedata(n,6)-1)/(linedata(n,3)+1j*linedata(n,4)))...+(conj(Ud(ii))-conj(Ud(jj)))*conj(-Y(ii,jj)));ii=linedata(n,2);jj=linedata(n,1);Sij(2*n,1)=ii;

Sij(2*n,2)=jj;Sij(2*n,3)=Ud(ii)*(conj(Ud(ii))...*conj(1j*linedata(n,5)+(1-linedata(n,6))/(linedata(n,3)+1j*linedata(n,4)))...+(conj(Ud(ii))-conj(Ud(jj)))*conj(-Y(ii,jj)));endendfunctionlinedata1=NumResort(linedata,Num)linedata1=linedata;forn=1:length(Num)n1=find(linedata(:,1)==Num(n));n2=find(linedata(:,2)==Num(n));linedata1(n1,1)=n*ones(length(n1),1);linedata1(n2,2)=n*ones(length(n2),1);endendfunctionJ=CalcJ(e,f,Y,busnum,pqnum)G=real(Y);B=imag(Y);J=zeros(2*(busnum-1));H=zeros(pqnum,pqnum);N=M=L=forH;H;H;iifor=1:pqnumjj=1:pqnumifii?=jjH(ii,jj)=-B(ii,jj)*e(ii)+G(ii,jj)*f(ii);N(ii,jj)=G(ii,jj)*e(ii)+B(ii,jj)*f(ii);M(ii,jj)=-N(ii,jj);L(ii,jj)=H(ii,jj);J(2*ii-1,2*jj-1)=H(ii,jj);J(2*ii-1,2*jj)N=M=L=forH;H;H;iifor=1:pqnumjj=1:pqnum+1j*(G(ii,ii)*f(ii)+B(ii,ii)*e(ii));forjj=1:busnumifii?=jjlii=Iii+G(ii,jj)*e(jj)-B(ii,jj)*f(jj)...+1j*(G(ii,jj)*f(jj)+B(ii,jj)*e(jj));endendH(ii,ii)=-B(ii,ii)*e(ii)+G(ii,ii)*f(ii)+imag(Iii);N(ii,ii)=G(ii,ii)*e(ii)+B(ii,ii)*f(ii)+real(Iii);M(ii,ii)=-G(ii,ii)*e(ii)-B(ii,ii)*f(ii)+real(Iii);L(ii,ii)=-B(ii,ii)*e(ii)+G(ii,ii)*f(ii)-imag(Iii);J(2*ii-1,2*ii-1)=H(ii,ii);J(2*ii-1,2*ii)=N(ii,ii);J(2*ii,2*ii-1)=M(ii,ii);J(2*ii,2*ii)=L(ii,ii);endH=zeros(pqnum,busnum-pqnum-2);N=H;M=H;L=H;forii=1:pqnumforjj=pqnum+1:busnum-1ifii?=jjH(ii,jj)=-B(ii,jj)*e(ii)+G(ii,jj)*f(ii);N(ii,jj)=G(ii,jj)*e(ii)+B(ii,jj)*f(ii);M(ii,jj)=-N(ii,jj);L(ii,jj)=H(ii,jj);J(2*ii-1,2*jj-1)=H(ii,jj);J(2*ii-1,2*jj)=N(ii,jj);J(2*ii,2*jj-1)=M(ii,jj);J(2*ii,2*jj)=L(ii,jj);endendendH=zeros(busnum-pqnum-2,pqnum);N=H;R=H;S=H;forii=pqnum+1:busnum-1forjj=1:pqnumifii?=jjH(ii,jj)=-B(ii,jj)*e(ii)+G(ii,jj)*f(ii);N(ii,jj)=G(ii,jj)*e(ii)+B(ii,jj)*f(ii);R(ii,jj)=0;S(ii,jj)=0;J(2*ii-1,2*jj-1)=H(ii,jj);J(2*ii-1,2*jj)=N(ii,jj);J(2*ii,2*jj-1)=R(ii,jj);

J(2*ii,2*jj)=S(ii,jj);endendendH=N=R=S=forzeros(busnum-pqnum-2);H;H;H;iifor=pqnum+1:busnum-1jj=pqnum+1:busnum-1ifii?=jjH(ii,jj)=-B(ii,jj)*e(ii)+G(ii,jj)*f(ii);N(ii,jj)=G(ii,jj)*e(ii)+B(ii,jj)*f(ii);R(ii,jj)=0;S(ii,jj)=0;J(2*ii-1,2*jj-1)endH=N=R=S=forzeros(busnum-pqnum-2);H;H;H;iifor=pqnum+1:busnum-1jj=pqnum+1:busnum-1+1j*(G(ii,ii)*f(ii)+B(ii,ii)*e(ii));forjj=1:busnumifii?=jjlii=Iii+G(ii,jj)*e(jj)-B(ii,jj)*f(jj)...+1j*(G(ii,jj)*f(jj)+B(ii,jj)*e(jj));endendH(ii,ii)=-B(ii,ii)*e(ii)+G(ii,ii)*f(ii)+imag(Iii);N(ii,ii)=G(ii,ii)*e(ii)+B(ii,ii)*f(ii)+real(Iii);R(ii,ii)=2*f(ii);S(ii,ii)=2*e(ii);J(2*ii-1,2*ii-1)=H(ii,ii);J(2*ii-1,2*ii)=N(ii,ii);J(2*ii,2*ii-1)=R(ii,ii);J(2*ii,2*ii)=S(ii,ii);endendfunction[dP,dQ,dU2]=CalcCollection(P,Q,Us,e,f,Y,busnum,pqnum)G=real(Y);B=imag(Y);Pt=P;Qt=Q;dP=zeros(1,busnum);dQ=zeros(1,busnum);dU2=zeros(1,busnum);forii=1:pqnumPi=0;Qi=0;forjj=1:busnumPi=Pi+e(ii)*(G(ii,jj)*e(jj)-B(ii,jj)*f(jj))...+f(ii)*(G(ii,jj)*f(jj)+B(ii,jj)*e(jj));Qi=Qi+f(ii)*(G(ii,jj)*e(jj)-B(ii,jj)*f(jj))...-e(ii)*(G(ii,jj)*f(jj)+B(ii,jj)*e(jj));enddP(ii)=P(ii)-Pi;dQ(ii)=Q(ii)-Qi;Pt(ii)=Pi;Qt(ii)=Qi;endforii=pqnum+1:busnum-1Pi=0;forjj=1:busnumPi=Pi+e(ii)*(G(ii,jj)*e(jj)-B(ii,jj)*f(jj))...+f(ii)*(G(ii,jj)*f(jj)+B(ii,jj)*e(jj));enddP(ii)=P(ii)-Pi;Pt(ii)=Pi;dU2(ii)=Us(ii)A2-(e(ii)A2+f(ii)A2);endendfunctionY=BuildY(linedata)nf=linedata(:,1);nt=linedata(:,2);r=linedata(:,3);x=linedata(:,4);b=linedata(:,5);k=linedata(:,6);branchnum=length(nf);busnum=max([nf',nt']);y=ones(branchnum,1)./(r+1j*x);Y=zeros(busnum);forn=1:branchnumY(nf(n),nt(n))=Y(nf(n),nt(n))-k(n)*y(n);Y(nt(n),nf(n))=Y(nf(n),nt(n));Y(nf(n),nf(n))=Y(nf(n),nf(n))+k(n)A2*y(n)+1j*b(n);Y(nt(n),nt(n))=Y(nt(n),nt(n))+y(n)+1j*b(n);endendfunctions=Jordan(A,b)%約當(dāng)消元ifnargin<1A=[123;275;149];b=[16-3]';endifrank(A)~=rank([A,b])error('A矩陣的秩和增廣矩陣的秩不相同，方程不存在唯一解');end[~,n]=size(A);A(:,n+1)=b;fork=1:nA(k,k:end)=A(k,k:end)/A(k,k);%規(guī)格化forq=1:k-1%消上三角ifA(q,k)~=0A(q,k:end)=A(q,k:end)-A(k,k:end).*A(q,k);endendforp=k+1:n%消下三角ifA(p,k)~=0A(p,k:end)=A(p,k:end)-A(k,k:end).*A(p,k);end

endends=A(:,end);end七、實(shí)驗(yàn)結(jié)果編號(hào):節(jié)點(diǎn)電壓:0.98160.OO8E5i0,958了0.l084i1.0924+0.l290i1.0500平衡節(jié)點(diǎn)功率：0.3679+■0.264710.2462。，。1I6i0.010610.029310.0931i0.136110.104510.140610_1602j線路功率-0.24000.2462。，。1I6i0.010610.029310.0931i0.136110.104510.140610_1602j節(jié)點(diǎn)功率-0.3000-0.ISOOi4-0.5500-0.1300i30.5000+0.［）934i0.3679卜0M64了iA?八、實(shí)驗(yàn)分析通過(guò)在程序中對(duì)節(jié)點(diǎn)排列順序進(jìn)行優(yōu)化，這樣無(wú)需對(duì)網(wǎng)絡(luò)進(jìn)行特定編號(hào)。同時(shí)為了提高速度，沒(méi)有通過(guò)求逆的方式來(lái)求解線性方程組，而是利用約當(dāng)消元的方法來(lái)求解線性方程組。八、調(diào)試記錄在計(jì)算公式時(shí)，沒(méi)有理解以下符號(hào)的意義導(dǎo)致代碼出錯(cuò)：是節(jié)點(diǎn)i的電壓幅值（誤：節(jié)點(diǎn)i的電壓）；為（誤：節(jié)點(diǎn)i的所有對(duì)地導(dǎo)納之和）；為（誤：）。在編寫節(jié)點(diǎn)優(yōu)化部分代碼時(shí)，沒(méi)有分清行號(hào)與列號(hào)的字母導(dǎo)致編程出錯(cuò)。南昌大學(xué)實(shí)驗(yàn)報(bào)告學(xué)生姓名：學(xué)號(hào)：6專業(yè)班級(jí)：電氣152實(shí)驗(yàn)類型：口驗(yàn)證口綜合■設(shè)計(jì)口創(chuàng)新實(shí)驗(yàn)日期：實(shí)驗(yàn)成績(jī)：一、實(shí)驗(yàn)題目水火電廠間有功功率負(fù)荷的經(jīng)濟(jì)分配二、實(shí)驗(yàn)?zāi)康模罕緦?shí)驗(yàn)通過(guò)對(duì)水火電廠間有功功率負(fù)荷的經(jīng)濟(jì)分配的計(jì)算機(jī)程序編制與調(diào)試，獲得能計(jì)算有功負(fù)荷經(jīng)濟(jì)分配的計(jì)算機(jī)程序，即根據(jù)已知的水火電廠的耗量特性、發(fā)電用水量及各時(shí)間段的有功負(fù)荷功率，由計(jì)算機(jī)程序自動(dòng)完成各時(shí)段各電廠的有功出力的計(jì)算。通過(guò)實(shí)驗(yàn)教學(xué)加深學(xué)生對(duì)電力系統(tǒng)經(jīng)濟(jì)運(yùn)行概念的理解，學(xué)會(huì)運(yùn)用數(shù)學(xué)模型進(jìn)行水火電廠間的有功功率負(fù)荷的經(jīng)濟(jì)分配計(jì)算，熟悉各種常用應(yīng)用軟件，熟悉硬件設(shè)備的使用方法，加強(qiáng)編制調(diào)試計(jì)算機(jī)程序的能力，提高工程計(jì)算的能力，學(xué)習(xí)如何將理論知識(shí)和實(shí)際工程問(wèn)題結(jié)合起來(lái)。三、實(shí)驗(yàn)內(nèi)容：編制調(diào)試水、火電廠間有功功率負(fù)荷經(jīng)濟(jì)分配的計(jì)算機(jī)程序。程序要求根據(jù)已知參數(shù)，完成水火電廠間有功功率負(fù)荷的經(jīng)濟(jì)分配計(jì)算，要求計(jì)算出各電廠滿足經(jīng)濟(jì)運(yùn)行要求所分配的有功負(fù)荷等參數(shù)。熟悉電力系統(tǒng)經(jīng)濟(jì)運(yùn)行的計(jì)算方法，按照計(jì)算方法編制程序。將事先編制好的計(jì)算程序原代碼由自備移動(dòng)存儲(chǔ)設(shè)備導(dǎo)入計(jì)算機(jī)。在相應(yīng)的編程環(huán)境下對(duì)程序進(jìn)行組織調(diào)試。應(yīng)用計(jì)算例題驗(yàn)證程序的計(jì)算效果。對(duì)調(diào)試正確的計(jì)算程序進(jìn)行存儲(chǔ)、打印。完成本次實(shí)驗(yàn)的實(shí)驗(yàn)報(bào)告。四、實(shí)驗(yàn)數(shù)據(jù)一個(gè)火電廠和一個(gè)水電廠并聯(lián)運(yùn)行?；痣姀S的燃料消耗特性為：水電廠的耗水量特性為：水電廠的給定日用水量為。系統(tǒng)的日負(fù)荷變化如下：時(shí)負(fù)荷為350MW，時(shí)，負(fù)荷為700MW，時(shí)負(fù)荷為500MW?；痣姀S容量為600MW，水電廠容量為450MW。試確定水、火電廠間的功率經(jīng)濟(jì)分配。五、程序框圖六、實(shí)驗(yàn)程序functionOptimizeDispatckclc,clear,closeallsymsgamaxPLDxPHxsW=2+0.8*PHx+1.5e-3*PHxA2;PHs=(0.4-0.8*gamax+0.0007*PLDx)/(0.003*gamax+0.0007);PTs=(0.8*gamax-0.4+0.003*gamax*PLDx)/(0.003*gamax+0.0007);gama=0;PLD=[350,700,500];T=[8,10,6];WH=1.5e7;%d=0.2;kp=0.15;ki=0.01;kd=0.01;W=0;W0=0;W1=0;plot(1:0.01:100,WH,'r');gridonholdonforn=1:100W0=W1;W1=W;W=0;form=1:length(PLD)PH(m)=subs(PHs,{gamax,PLDx},{gama,PLD(m)});PT(m)=subs(PTs,{gamax