【正文】 e power grid voltage and seek the trend of the network and the network of the ponents of the power loss, and getting electric standard is to use the power flow calculation and analysis, the specific have MATLAB words: Power system。本文主要運(yùn)用的事潮流計(jì)算，潮流計(jì)算是電力網(wǎng)絡(luò)設(shè)計(jì)與運(yùn)行中最基本的運(yùn)算，對(duì)電力網(wǎng)絡(luò)的各種設(shè)計(jì)方案及各種運(yùn)行方式進(jìn)行潮流計(jì)算，可以得到各種電網(wǎng)各節(jié)點(diǎn)的電壓，并求得網(wǎng)絡(luò)的潮流及網(wǎng)絡(luò)中的各元件的電力損耗，進(jìn)而求得電能損耗。for(i=0。fp1=fopen(”“,”w+“)。for(i=0。printf(”n=====ping heng jie dian gong lv =====n“)。int i,j。float ds[N][N]。float dp[N][N]={0}。float q[N][N]={0}。float qs=0,qv1=0,qv2=0。else goto top。ielse { kp=0。iguass(NM,M,y2,B,x2)。imax)max=fabs(dqu[i])。iguass(1,N1,y1,B,x1)。imax)max=fabs(dpu[i])。float x1[N1],x2。for(i=0。} getch()。b[i][j]=lb。g[i][j]=lg。b[j][j]+=lb+lb1。g[j][j]+=lg。lb=lx/(lr*lr+lx*lx)。lb1=tb。lr=tR。jHeadnode1。i=0。int i,j,k。4 0。2 0。0 5+0i 0 3。021i 1 0 0 2。4 2 0 1。1 3 + 0 1 0。end end B2 Times=Times+1。end end for i=1:n if B2(i,6)==3 j=j+1。for i=1:n if B2(i,6)==2 j=j+1。DetaU=inv(Jacbi)*DetaS。end if k==(n1)k=0。Jacbi(2*h,2*k1)=0。else Jacbi(2*h1,2*k1)=imag(Y(i,j))*real(B2(i,3))+real(Y(i,j))*imag(B2(i,3))。Jacbi(2*h,2*k1)=2*imag(B2(i,3))。if i==j Jacbi(2*h1,2*k1)=imag(Y(i,j))*real(B2(i,3))+real(Y(i,j))*imag(B2(i,3))+imag(I(h,1))。for i=1:n if B2(i,6)==3 h=h+1。end if k==(n1)k=0。Jacbi(2*h,2*k1)=Jacbi(2*h1,2*k)。else Jacbi(2*h1,2*k1)=imag(Y(i,j))*real(B2(i,3))+real(Y(i,j))*imag(B2(i,3))。Jacbi(2*h,2*k1)=Jacbi(2*h1,2*k)+2*real(I(h,1))。if i==j Jacbi(2*h1,2*k1)=imag(Y(i,j))*real(B2(i,3))+real(Y(i,j))*imag(B2(i,3))+imag(I(h,1))。for i=1:n if B2(i,6)==2 h=h+1。h=0。I(h,1)=(OrgS(2*h1,1)OrgS(2*h,1)*sqrt(1))/conj(B2(i,3))。h=0。DetaS(2*h,1)=real(PVU(t,1))^2+imag(PVU(t,1))^2real(B2(i,3))^2imag(B2(i,3))^2。t=t+1。for i=1:n if i~=isbamp。DetaS(2*h,1)=imag(B2(i,2))OrgS(2*h,1)。B2(i,6)==2 h=h+1。end end end OrgS %創(chuàng)建DetaS h=0。for j=1:n OrgS(2*h1,1)=OrgS(2*h1,1)+real(B2(i,3))*(real(Y(i,j))*real(B2(j,3))imag(Y(i,j))*imag(B2(j,3)))+imag(B2(i,3))*(real(Y(i,j))*imag(B2(j,3))+imag(Y(i,j))*real(B2(j,3)))。end end end for i=1:n if i~=isbamp。for j=1:n OrgS(2*h1,1)=OrgS(2*h1,1)+real(B2(i,3))*(real(Y(i,j))*real(B2(j,3))imag(Y(i,j))*imag(B2(j,3)))+imag(B2(i,3))*(real(Y(i,j))*imag(B2(j,3))+imag(Y(i,j))*real(B2(j,3)))。for i=1:n if i~=isbamp。%！!初始功率參數(shù)在迭代過程中是不累加的，所以在這里必須將其初始化為零矩陣 h=0。B2(i,3)=B2(i,3)+DetaU(2*j,1)+DetaU(2*j1,1)*sqrt(1)。B2(i,3)=B2(i,3)+DetaU(2*j,1)+DetaU(2*j1,1)*sqrt(1)。DetaU %修正節(jié)點(diǎn)電壓 j=0。end end end end end Jacbi %求解修正方程，獲取節(jié)點(diǎn)電壓的不平衡量 DetaU=zeros(2*n2,1)。Jacbi(2*h,2*k)=0。Jacbi(2*h1,2*k)=real(Y(i,j))*real(B2(i,3))+imag(Y(i,j))*imag(B2(i,3))。Jacbi(2*h,2*k)=2*real(B2(i,3))。Jacbi(2*h1,2*k)=real(Y(i,j))*real(B2(i,3))+imag(Y(i,j))*imag(B2(i,3))+real(I(h,1))。for j=1:n if j~=isb k=k+1。end end end end end k=0。Jacbi(2*h,2*k)=Jacbi(2*h1,2*k1)。Jacbi(2*h1,2*k)=real(Y(i,j))*real(B2(i,3))+imag(Y(i,j))*imag(B2(i,3))。Jacbi(2*h,2*k)=Jacbi(2*h1,2*k1)2*imag(I(h,1))。Jacbi(2*h1,2*k)=real(Y(i,j))*real(B2(i,3))+imag(Y(i,j))*imag(B2(i,3))+real(I(h,1))。for j=1:n if j~=isb k=k+1。k=0。end end I %創(chuàng)建Jacbi(雅可比矩陣)Jacbi=zeros(2*n2)。for i=1:n if i~=isb h=h+1。end end DetaS %創(chuàng)建I，用于存儲(chǔ)節(jié)點(diǎn)電流參數(shù) i=zeros(n1,1)。DetaS(2*h1,1)=real(B2(i,2))OrgS(2*h1,1)。B2(i,6)==3 h=h+1。end end t=0。DetaS(2*h1,1)=real(B2(i,2))OrgS(2*h1,1)。for i=1:n %對(duì)PQ節(jié)點(diǎn)的處理 if i~=isbamp。PVU(t,1)=B2(i,3)。t=0。OrgS(2*h,1)=OrgS(2*h,1)+imag(B2(i,3))*(real(Y(i,j))*real(B2(j,3))imag(Y(i,j))*imag(B2(j,3)))real(B2(i,3))*(real(Y(i,j))*imag(B2(j,3))+imag(Y(i,j))*real(B2(j,3)))。B2(i,6)==3 h=h+1。OrgS(2*h,1)=OrgS(2*h,1)+imag(B2(i,3))*(real(Y(i,j))*real(B2(j,3))imag(Y(i,j))*imag(B2(j,3)))real(B2(i,3))*(real(Y(i,j))*imag(B2(j,3))+imag(Y(i,j))*real(B2(j,3)))。B2(i,6)==2 h=h+1。j=0。DetaS=zeros(2*n2,1)。Y(q,q)=Y(q,q)+1/(B1(i,5)^2*B1(i,3))。Y(q,p)=Y(p,q)。q=B1(i,2)。Y(q,q)=Y(q,q)+1/B1(i,3)+*B1(i,4)。Y(q,p)=Y(p,q)。q=B1(i,2)。Times=1。)。X=input(39。請(qǐng)輸入節(jié)點(diǎn)參數(shù):B2=39。)。B1=input(39。請(qǐng)輸入誤差精度:pr=39。)。isb=input(39。請(qǐng)輸入支路數(shù):n1=39。)。n=input(39。%X為節(jié)點(diǎn)號(hào)和對(duì)地參數(shù)矩陣。%第五烈為含變壓器支路的變壓器的變比%第六列為變壓器是否是否含有變壓器的參數(shù)，其中“1”為含有變壓器，%“0”為不含有變壓器。%第三列為支路的串列阻抗參數(shù)。第一篇：基于MATLAB的電力系統(tǒng)潮流計(jì)算基于MATLAB的電力系統(tǒng)潮流計(jì)算%簡單潮流計(jì)算的小程序，相關(guān)的原始數(shù)據(jù)數(shù)據(jù)數(shù)據(jù)輸入格式如下：%B1是支路參數(shù)矩陣，第一列和第二列是節(jié)點(diǎn)編號(hào)。節(jié)點(diǎn)編號(hào)由小到大編寫 %對(duì)于含有變壓器的支路，第一列為低壓側(cè)節(jié)點(diǎn)編號(hào)，第二列為高壓側(cè)節(jié)點(diǎn) %編號(hào)，將變壓器的串聯(lián)阻抗置于低壓側(cè)處理。%第四列為支路的對(duì)地導(dǎo)納參數(shù)。%B2為節(jié)點(diǎn)參數(shù)矩陣，其中第一列為節(jié)點(diǎn)注入發(fā)電功率參數(shù)；第二列為節(jié)點(diǎn) %負(fù)荷功率參數(shù)；第三列為節(jié)點(diǎn)電壓參數(shù)；第六列為節(jié)點(diǎn)類型參數(shù)，其中 %“1”為平衡節(jié)點(diǎn)，“2”為PQ節(jié)點(diǎn)，“3”為PV節(jié)點(diǎn)參數(shù)。其中第一列為節(jié)點(diǎn)編號(hào)，第二列為節(jié)點(diǎn)對(duì)地 %參數(shù)。請(qǐng)輸入節(jié)點(diǎn)數(shù):n=39。n1=input(39。)。請(qǐng)輸入平衡節(jié)點(diǎn)號(hào):isb=39。pr=input(39。)。請(qǐng)輸入支路參數(shù):B1=39。B2=input(39。)。節(jié)點(diǎn)號(hào)和對(duì)地參數(shù):X=39。Y=zeros(n)。%置迭代次數(shù)為初始值 %創(chuàng)建節(jié)點(diǎn)導(dǎo)納矩陣 for i=1:n1 if B1(i,6)==0 %不含變壓器的支路 p=B1(i,1)。Y(p,q)=Y(p,q)1/B1(i,3)。Y(p,p)=Y(p,p)+1/B1(i,3)+*B1(i,4)。else %含有變壓器的支路 p=B1(i,1)。Y(p,q)=Y(p,q)1/(B1(i,3)*B1(i,5))。Y(p,p)=Y(p,p)+1/B1(i,3)。end end Y OrgS=zeros(2*n2,1)。%將OrgS、DetaS初始化%創(chuàng)建OrgS，用于存儲(chǔ)初始功率參數(shù) h=0。for i=1:n %對(duì)PQ節(jié)點(diǎn)的處理 if i~=isbamp。for j=1:n OrgS(2*h1,1)=OrgS(2*h1,1)+real(B2(i,3))*(real(Y(i,j))*real(B2(j,3))imag(Y(i,j))*imag(B2(j,3)))+imag(B2(i,3))*(real(Y(i,j))*imag(B2(j,3))+imag(Y(i,j))*real(B2(j,3)))。end end end for i=1:n %對(duì)PV節(jié)點(diǎn)的處理，注意這時(shí)不可再將h初始化為0 if i~=isbamp。for j=1:n OrgS(2*h1,1)=OrgS(2*h1,1)+real(B2(i,3))*(real(Y(i,j))*real(B2(j,3))imag(Y(i,j))*imag(B2(j,3)))+imag(B2(i,3))*(real(Y(i,j))*imag(B2(j,3))+imag(Y(i,j))*real(B2(j,3)))。end end end OrgS %創(chuàng)建PVU 用于存儲(chǔ)PV節(jié)點(diǎn)的初始電壓 PVU=zeros(nh1,1)。for i=1:n if B2(i,6)==3 t=t+1。end end PVU %創(chuàng)建DetaS，用于存儲(chǔ)有功功率、無功功率和電壓幅值的不平衡量 h=0。B2(i,6)==2 h=h+1。DetaS(2*h,1)=imag(B2(i,2))OrgS(2*h,1)。for i=1:n %對(duì)PV節(jié)點(diǎn)的處理，注意這時(shí)不可再將h初始化為0 if i~=isbamp。t=t+1。DetaS(2*h,1)=real(PVU(t,1))^2+imag(PVU(t,1))^2real(B2(i,3))^2imag(B2(i,3))^2。h=0。I(h,1)=(OrgS(2*h1,1)OrgS(2*h,1)*sqrt(1))/conj(B2(i,3))。h=0。for i=1:n %對(duì)PQ節(jié)點(diǎn)的處理 if B2(i,6)==2 h=h+1。if i==j %對(duì)角元素的處理Jacbi(2*h1,2*k1)=imag(Y(i,j))*real(B2(i,3))+real(Y(i,j))*imag(B2(i,3))+imag(I(h,1))。Jacbi(2*h,2*k1)=Jacbi(2*h1,2*k)+2*real(I(h,1))。else %非對(duì)角元素的處理Jacbi(2*h1,2*k1)=imag(Y(i,j))*real(B2(i,3))+real(Y(i,j))*imag(B2(i,3))。Jacbi(2*h,2*k1)=Jacbi(2*h1,2*k)。end if k==(n1)%將用于內(nèi)循環(huán)的指針置于初始值，以確保雅可比矩陣換行k=0。for i=1:n %對(duì)PV節(jié)點(diǎn)的處理 if B2(i,6)==3 h=h+1。if i==j %對(duì)角元素的處理Jacbi(2*h1,2*k1)=imag(Y(i,j))*real(B2(i,3))+real(Y(i,j))*imag(B2(i,3))+imag(I(h,1))。Jacbi(2*h,2*k1)=2*imag(B2(i,3))。else %非對(duì)角元素的處理Jacbi(2*h1,2*k1)=imag(Y(i,j))*real(B2(i,3))+real(Y(i,j))*imag(B2(i,3))。Jacbi(2*h,2*k1)=0。end if k==(n1)%將用于內(nèi)循環(huán)的指針置于初始值，以確保雅可