潮流計(jì)算的MATLAB源程序_第1頁(yè)
潮流計(jì)算的MATLAB源程序_第2頁(yè)
潮流計(jì)算的MATLAB源程序_第3頁(yè)
潮流計(jì)算的MATLAB源程序_第4頁(yè)
潮流計(jì)算的MATLAB源程序_第5頁(yè)
免費(fèi)預(yù)覽已結(jié)束，剩余2頁(yè)可下載查看

下載本文檔

版權(quán)說(shuō)明：本文檔由用戶提供并上傳，收益歸屬內(nèi)容提供方，若內(nèi)容存在侵權(quán)，請(qǐng)進(jìn)行舉報(bào)或認(rèn)領(lǐng)

文檔簡(jiǎn)介

1、%簡(jiǎn)單潮流計(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è)處理。%第三列為支路的串列阻抗參數(shù)。%第四列為支路的對(duì)地導(dǎo)納參數(shù)。%第五列為含變壓器支路的變壓器的變比%第六列為變壓器是否含有變壓器的參數(shù)，其中“1為含有變壓器,% “ 0);不含有變壓器。%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)

2、，“3為PV節(jié)點(diǎn)參數(shù)。1 .1' %X為節(jié)點(diǎn)號(hào)和對(duì)地參數(shù)矩陣。其中第一列為節(jié)點(diǎn)編號(hào)，第二列為節(jié)點(diǎn)對(duì)地%參數(shù)。n=input('請(qǐng)輸入節(jié)點(diǎn)數(shù) ：n=');n1=input('請(qǐng)輸入支路數(shù):n1=');isb=input('請(qǐng)輸入平衡節(jié)點(diǎn)號(hào)：isb=');1pr=input('請(qǐng)輸入誤差精度:pr=');B1=input('請(qǐng)輸入支路參數(shù):B1=');B2=input('請(qǐng)輸入節(jié)點(diǎn)參數(shù)舊2=');X=input('節(jié)點(diǎn)號(hào)和對(duì)地參數(shù) ：X=');Y=zeros(n);Times=1

3、; %置迭代次數(shù)為初始值%創(chuàng)建節(jié)點(diǎn)導(dǎo)納矩陣fori=1:n1ifB1(i,6)=0% 不含變壓器的支路p=B1(i,1);q=B1(i,2);Y(p,q尸Y(p,q)-1/B1(i,3);Y(q,p)=Y(p,q);Y(p,p尸丫(p,p)+1/B1(i,3)+0.5*B1(i,4);Y(q,q尸丫(q,q)+1/B1(i,3)+0.5*B1(i,4);else %含有變壓器的支路p=B1(i,1);q=B1(i,2);Y(p,q尸Y(p,q)-1/(B1(i,3)*B1(i,5);Y(q,p)=Y(p,q);Y(p,p)=Y(p,p)+1/B1(i,3);Y(q,q)=Y(q,q)+1/(B

4、1(i,52*B1(i,3);endendYOrgS=zeros(2*n-2,1);DetaS=zeros(2*n-2,1); % 將 OrgS、DetaS 初始化%創(chuàng)建OrgS,用于存儲(chǔ)初始功率參數(shù)h=0;j=0;fori=1:n %對(duì)PQ節(jié)點(diǎn)的處理ifi=isb&B2(i,6)=2h=h+1;forj=1:nOrgS(2*h-1,1)=OrgS(2*h-1,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)*(re al(Y(i,j)*imag(B2(j,3)+imag(Y(

5、i,j)*real(B2(j,3);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);endendip/ j 一產(chǎn) j :endfori=1:n %對(duì)PV節(jié)點(diǎn)的處理，注意這時(shí)不可再將h初始化為0ifi=isb&B2(i,6)=3h=h+1;forj=1:nOrgS(2*h-1,1)=OrgS(2*h-1,1)+real(B2(i,3)

6、*(real(Y(i,j)*real(B2(j,3)-imag(Y(i,j)*imag(B2(j,3)+imag(B2(i,3)*(re al(Y(i,j)*imag(B2(j,3)+imag(Y(i,j)*real(B2(j,3);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);I .,endendendOrgS%創(chuàng)建PVU用于存儲(chǔ)PV節(jié)

7、點(diǎn)的初始電壓PVU=zeros(n-h-1,1);t=0;fori=1:nifB2(i,6)=3t=t+1;PVU(t,1)=B2(i,3);endendPVU%創(chuàng)建DetaS,用于存儲(chǔ)有功功率、無(wú)功功率和電壓幅值的不平衡量歡迎共閱h=0;fori=1:n%對(duì)PQ節(jié)點(diǎn)的處理ifi=isb&B2(i,6)=2h=h+1;DetaS(2*h-1,1)=real(B2(i,2)-OrgS(2*h-1,1);DetaS(2*h,1)=imag(B2(i,2)-OrgS(2*h,1);endendt=0;fori=1:n %對(duì)PV節(jié)點(diǎn)的處理，注意這時(shí)不可再將h初始化為0ifi=isb&B

8、2(i,6)=3h=h+1;t=t+1;DetaS(2*h-1,1)=real(B2(i,2)-OrgS(2*h-1,1);DetaS(2*h,1)=real(PVU(t,1)A2+imag(PVU(t,1)A2-real(B2(i,3)A2-imag(B2(i,3)A2;endendDetaS%創(chuàng)建I,用于存儲(chǔ)節(jié)點(diǎn)電流參數(shù)i=zeros(n-1,1);h=0;. I I '1"fori=1:nifi=isbh=h+1;I(h,1)=(OrgS(2*h-1,1)-OrgS(2*h,1)*sqrt(-1)/conj(B2(i,3);endendI%創(chuàng)建Jacbi(雅可比矩陣)J

9、acbi=zeros(2*n-2);h=0；k=0;fori=1:n % 對(duì)PQ節(jié)點(diǎn)的處理,1ifB2(i,6)=2h=h+1;forj=1:nifj=isbk=k+1;ifi=j %對(duì)角元素的處理Jacbi(2*h-1,2*k-1)=-imag(Y(i,j)*real(B2(i,3)+real(Y(i,j)*imag(B2(i,3)+imag(I(h,1);Jacbi(2*h-1,2*k尸real(Y(i,j)*real(B2(i,3)+imag(Y(i,j)*imag(B2(i,3)+real(I(h,1);Jacbi(2*h,2*k-1)=-Jacbi(2*h-1,2*k)+2*real

10、(I(h,1);Jacbi(2*h,2*k尸Jacbi(2*h-1,2*k-1)-2*imag(I(h,1);else%非對(duì)角元素的處理Jacbi(2*h-1,2*k-1)=-imag(Y(i,j)*real(B2(i,3)+real(Y(i,j)*imag(B2(i,3);Jacbi(2*h-1,2*k尸real(Y(i,j)*real(B2(i,3)+imag(Y(i,j)*imag(B2(i,3);Jacbi(2*h,2*k-1)=-Jacbi(2*h-1,2*k);Jacbi(2*h,2*k尸Jacbi(2*h-1,2*k-1);endifk=(n1)%將用于內(nèi)循環(huán)的指針置于初始值，以

11、確保雅可比矩陣換行k=0;endendendendendk=0;fori=1:n % 對(duì)PV節(jié)點(diǎn)的處理ifB2(i,6)=3h=h+1;forj=1:nifj=isbk=k+1;ifi=j %對(duì)角元素的處理Jacbi(2*h-1,2*k-1)=-imag(Y(i,j)*real(B2(i,3)+real(Y(i,j)*imag(B2(i,3)+imag(I(h,1);I I 'L XJacbi(2*h-1,2*k尸real(Y(i,j)*real(B2(i,3)+imag(Y(i,j)*imag(B2(i,3)+real(I(h,1);Jacbi(2*h,2*k-1)=2*imag(B

12、2(i,3);Jacbi(2*h,2*k)=2*real(B2(i,3);|" f.Z' -f :else%非對(duì)角元素的處理Jacbi(2*h-1,2*k-1)=-imag(Y(i,j)*real(B2(i,3)+real(Y(i,j)*imag(B2(i,3);Jacbi(2*h-1,2*k尸real(Y(i,j)*real(B2(i,3)+imag(Y(i,j)*imag(B2(i,3);Jacbi(2*h,2*k-1)=0;Jacbi(2*h,2*k)=0;endifk=(n-i) %將用于內(nèi)循環(huán)的指針置于初始值，以確保雅可比矩陣換行k=0;endendendenden

13、d. . 'IJacbi%求解修正方程，獲取節(jié)點(diǎn)電壓的不平衡量DetaU=zeros(2*n-2,1);DetaU=inv(Jacbi)*DetaS;DetaU%修正節(jié)點(diǎn)電壓j=0；fori=1:n % 對(duì)PQ節(jié)點(diǎn)處理ifB2(i,6)=2j=j+1；B2(i,3)=B2(i,3)+DetaU(2*j,1)+DetaU(2*j-1,1)*sqrt(-1);endendfori=1:n % 對(duì)PV節(jié)點(diǎn)的處理ifB2(i,6)=3j=j+1;B2(i,3)=B2(i,3)+DetaU(2*j,1)+DetaU(2*j-1,1)*sqrt(-1);endendB2% 開始循環(huán) *whilea

14、bs(max(DetaU)>prOrgS=zeros(2*n-2,1); %!初始功率參數(shù)在迭代過程中是不累加的，所以在這里必須將其初始化為零矩陣h=0;j=0; fori=1:nifi=isb&B2(i,6)=2h=h+1;forj=1:nOrgS(2*h-1,1)=OrgS(2*h-1,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)*(re I I I, /' " . " fal(Y(i,j)*imag(B2(j,3)+imag(Y(i,

15、j)*real(B2(j,3);OrgS(2*h,1)=OrgS(2*h,1)+imag(B2(i,3)*(real(Y(i,j)*real(B2Q,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);endend 'l 'lendfori=1:nifi=isb&B2(i,6)=3h=h+1; I Xforj=1:nOrgS(2*h-1,1)=OrgS(2*h-1,1)+real(B2(i,3)*(real(Y(i,j)*real(B2(j

16、,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);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);endend * %、 -iendOrgS%創(chuàng)建DetaSh=0;fori=1:nifi=isb&B2

17、(i,6)=2h=h+1;DetaS(2*h-1,1)=real(B2(i,2)-OrgS(2*h-1,1);DetaS(2*h,1)=imag(B2(i,2)-OrgS(2*h,1);end endt=0；fori=1:nifi=isb&B2(i,6)=3h=h+1;t=t+1;DetaS(2*h-1,1)=real(B2(i,2)-OrgS(2*h-1,1);DetaS(2*h,1)=real(PVU(t,1)A2+imag(PVU(t,1)A2-real(B2(i,3)A2-imag(B2(i,3)A2;endendDetaS%創(chuàng)建Ii=zeros(n-1,1);h=0;fori

18、=1:nifi=isbh=h+1;I(h,1)=(OrgS(2*h-1,1)-OrgS(2*h,1)*sqrt(-1)/conj(B2(i,3);endend I I I%創(chuàng)建Jacbi. I I '1"Jacbi=zeros(2*n-2); |" f. 'A" "：:h=0;k=0;i ：fori=1:nifB2(i,6)=2h=h+1;forj=1:nifj=isbk=k+1;ifi=jJacbi(2*h-1,2*k-1)=-imag(Y(i,j)*real(B2(i,3)+real(Y(i,j)*imag(B2(i,3)+imag(

19、I(h,1);Jacbi(2*h-1,2*k尸real(Y(i,j)*real(B2(i,3)+imag(Y(i,j)*imag(B2(i,3)+real(I(h,1);Jacbi(2*h,2*k-1)=-Jacbi(2*h-1,2*k)+2*real(I(h,1);Jacbi(2*h,2*k尸Jacbi(2*h-1,2*k-1)-2*imag(I(h,1);elseJacbi(2*h-1,2*k-1)=-imag(Y(i,j)*real(B2(i,3)+real(Y(i,j)*imag(B2(i,3);Jacbi(2*h-1,2*k尸real(Y(i,j)*real(B2(i,3)+imag

20、(Y(i,j)*imag(B2(i,3);Jacbi(2*h,2*k-1)=-Jacbi(2*h-1,2*k);Jacbi(2*h,2*k尸Jacbi(2*h-1,2*k-1);endifk=(n-1)k=0;endendendendk=0;fori=1:nifB2(i,6)=3h=h+1;forj=1:nifj=isbk=k+1;ifi=jJacbi(2*h-1,2*k-1)=-imag(Y(i,j)*real(B2(i,3)+real(Y(i,j)*imag(B2(i,3)+imag(I(h,1);Jacbi(2*h-1,2*k尸real(Y(i,j)*real(B2(i,3)+imag(Y(i,j)*imag(B2(i,3)+real(I(h,1);Jacbi(2*h,2*k-1)=2*imag(B2(i,3);Jacbi(2*h,2*k)=2*real(B2(i,3);elseJacbi(2*h-1,2*k-1)=-imag(Y(i,j)*re

溫馨提示

  • 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ù)覽，若沒有圖紙預(yù)覽就沒有圖紙。
  • 4. 未經(jīng)權(quán)益所有人同意不得將文件中的內(nèi)容挪作商業(yè)或盈利用途。
  • 5. 人人文庫(kù)網(wǎng)僅提供信息存儲(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ì)自己和他人造成任何形式的傷害或損失。

最新文檔

評(píng)論

0/150

提交評(píng)論