【正文】 。 另外，矩陣的數(shù)組運(yùn)算中還有冪運(yùn)算（運(yùn)算符為 .^ ）、指數(shù)運(yùn)算（ exp）、對(duì)數(shù)運(yùn)算（ log）、和開方運(yùn)算（ sqrt）等。有了“對(duì)應(yīng)元素”的規(guī)定，數(shù)組的運(yùn)算實(shí)質(zhì)上就是針對(duì)數(shù)組內(nèi)部的每個(gè)元素進(jìn)行的。矩陣的冪運(yùn)算與數(shù)組的冪運(yùn)算有很大的區(qū)別。 邏輯關(guān)系運(yùn)算 邏輯運(yùn)算是 MATLAB 中數(shù)組運(yùn)算所特有的一種運(yùn)算形式，也是幾乎所有的高級(jí)語(yǔ)言普遍適用的一種運(yùn)算。 26 — 拉夫遜法潮流計(jì)算程序 %本程序的 功能是用 牛頓 拉夫遜法進(jìn)行潮流計(jì)算 n=input(39。請(qǐng)輸入節(jié)點(diǎn)數(shù)： n=39。)。 n1=input(39。請(qǐng)輸入支路數(shù)： n1=39。)。 isb=input(39。請(qǐng)輸入平衡母線節(jié)點(diǎn)號(hào)： isb=39。)。 pr=input(39。請(qǐng)輸入誤差精度： pr=39。)。 B1=input(39。請(qǐng)輸入由支路參數(shù)形成的矩陣： B1=39。)。 B2=input(39。請(qǐng)輸入各節(jié)點(diǎn)參數(shù)形成的矩陣： B2=39。)。 X=input(39。請(qǐng)輸入由節(jié)點(diǎn)參數(shù)形成的矩陣： X=39。)。 Y=zeros(n)。e=zeros(1,n)。f=zeros(1,n)。V=seros(1,n)。O=zeros(1,n)。S1=zeros(n1)。 for i=1:n if X(i,2)~=0。 p=X(i,1)。 Y(p,p)=1./X(i,2)。 end end for i=1:n1 if B1(i,6)==0 p=B1(i,1)。q=B1(i,2)。 else p=B1(i,2)。q=B1(i,1)。 end Y(p,q)=Y(p,q)1./(B1(i,3)*B1(i,5)。 Y(p,q)=Y(p,q)。 Y(p,q)=Y(q,q)+1./(B1(i,3)*B1(i,5)^2)+B1(i,4)./2。 Y(p,p)=Y(p,p)+1./B1(i,3)+B1(i,4)./2。 end %求導(dǎo)納矩陣 G=real(Y)。B=imag(Y)。 for i=1:n e(i)=real(B2(i,3))。 f(i)=imag(B2(i,3))。 V(i)=B2(i,4)。 end for i=1:n S(i)=B2(i,1)B2(i,2)。 B(i,i)=B(i,i)+B2(i,5)。 end P=rea(S)。Q=imag(S)。 27 ICT1=0。IT2=1。NO=2*n。N=NO+1。a=0。 while IT2~=0 IT2=0。a=a+1。 for i=1:n。 if i~=isb C(i)=0。 D(i)=0。 for j1=1:n C(i)=C(i)+G(i,j1)*e(j1)B(i,j1)*f(j1)。 D(i)=D(i)+G(i,j1)*f(j1)+B(i,j1)*e(j1)。 end P1=C(i)*e(i)+f(i)*D(i)。 Q1=f(i)*C(i)D(i)*e(i)。 %求 39。P,Q39。 V2=e(i)^2+f(i)^2。 if B2(i,6)~=3 DP=P(i)P1。 DQ=Q(i)Q1。 for j1=1:n if j1~=isbamp。j1~=i X1=G(i,j1)*e(i)B(i,j1)*f(i)。 X2=B(i,j1)*e(i)G(i,j1)*f(i)。 X3=X2。 X4=X1。 p=2*i1。q=2*j11。J(p,q)=X3。J(p,N)=DQ。m=p+1。 J(m,q)=X1。J(m,N)=DP。q=q+1。J(p,q)=X4。J(m,q)=X2。 end end else DP=P(i)P1。 DV=V(i)~2V2。 for j1=1:n if j1~=isbamp。j1~=i X1=G(i,j1)*e(i)B(i,j1)*f(i)。 X2=B(i,j1)*e(i)G(i,j1)*f(i)。 X5=0。 X6=0。 p=2*i1。q=2*j11。J(p,q)=X5。J(p,N)=DV。m=p+1。 J(m,q)=X1。J(m,N)=DP。q=q+1。J(p,q)=X6。J(m,q)=X2。 elseif j1==iamp。j1~=isb X1=C(i)G(i,i)*e(i)B(i,i)*f(i)。 X2=D(i)+B(i,i)*e(i)G(i,i)*f(i)。 X5=2*e(i)。 X6=2*f(i)。 p=2*i1。q=2*j11。J(p,q)=X5。J(p,N)=DV。m=p+1。 28 J(m,q)=X1。J(m,N)=DP。q=q+1。J(p,q)=X6。J(m,q)=X2。 end end end end end %求雅可比矩陣 for k=3:N0 k1=k+1。N1=N。 for k2=k1:N1 J(k,k2)=J(k,k2)./J(k,k)。 end J(k,k)=1。 if k~=3 k4=k1。 for k3=3:k4 for k2=k1:N1 J(k3,k2)=J(k3,k2)J(k3,k)*J(k,k2)。 end J(k3,k)=0。 end end for k3=k1:N0 for k2=k1:N1 J(k3,k2)=J(k3,k2)J(k3,k)*J(k,k2)。 end J(k3,k)=0。 end end end for k=3:2:N01 L=(k+1)./2。 e(L)=e(L)J(k,N)。 k1=k+1。 f(L)=f(L)J(k1,N)。 end for k=3:N0 DET=abs(J(k,N))。 if DET=pr IT2=IT2+1 end end ICT2(a)=IT2 ICT1=ICT1+1。 for k=1:n 29 dy(k)=sqrt(e(k)^2+f(k)^2)。 end for i=1:n Dy(k)=sqrt(e(k)^2+f(k)^2)。 end for i=1:n Dy(ICT1,i)=dy(i)。 end end %用高斯消去法解“ w=J*V” disp(39。迭代次數(shù) 39。)。 disp(ICT1)。 disp(39。沒(méi)有達(dá)到精度要求的個(gè)數(shù) 39。)。 disp(ICT2)。 for k=1:n V(k)=sqrt(e(k)^2+f(k)^2)。 O(k)=atan(f(k)./e(k))*180./pi。 end E=e+f*j。 disp(39。各節(jié)點(diǎn)的實(shí)際電壓標(biāo)么值 E為（節(jié)點(diǎn)號(hào)從小到大的排列）： 39。)。 disp(E)。 disp(39。各節(jié)點(diǎn)的電壓大小 V為（節(jié)點(diǎn)號(hào)從小到大的排列）： 39。)。 disp(V)。 disp(39。各節(jié)點(diǎn)的電壓相角 O為（節(jié)點(diǎn)號(hào)從小到大的排列）： 39。)。 disp(O)。 for p=1:n C(p)=0。 for q=1:n C(p)=C(p)+conj(Y(p,q))*conj(E(q))。 end S(p)=E(p)*C(p)。 end disp(39。各節(jié)點(diǎn)的功率 S為（節(jié)點(diǎn)號(hào)從小到大排列）：‘）； disp(S)。 disp(39。各條支路的首端功率 Si為（順序同您輸入 B1時(shí)一樣）：‘）； for i=1:n1 if B1（ i,6） ==0 p=B1(i,1)。q=B1(i,2)。 else p=B1(i,2)。q=B1(i,1)。 end Si(p,q)=E(p)*(conj(E(p))*conj(B1(i,4)./2)+(conj(E(p)*B1(i,5))conj(E(q)))*conj(1./(B1(i,3)*B1(i,5))))。 disp(Si())。 end disp(39。各條支路的末端功率 Sj為（順序同您的輸入 B1時(shí)一樣）：‘）； 30 for i=1:n1 if B1(i,6)==0 p=B1(i,1)。q=B1(i,2)。 else p=B1(i,2)。q=B1(i,1)。 end Sj(q,p)=E(q)*(conj(E(q))*conj(B1(i,4)./2)+(xonj(E(q)./B1(i,5))conj(E(p)))*xonj(1./(B1(i,3)*B1(i,5))))。 disp(Sj(q,p))。 end disp(39。各條支路的功率損耗 DS為（順序同您輸入 B1時(shí)一樣）：’）； for i=1:n1 if B1(i,6)==0 p=B1(i,1)。q=B1(i,2)。 else p=B1(i,2)。q=B1(i,1)。 end DS(i)=Si(p,q)+Sj(q,p)。 disp(DS(i))。 end for i=1:ICT1 Cs(i)=i。 end disp(39。以下是每次迭代后各節(jié)點(diǎn)的電壓值（如圖所示）‘）； plot(Cs,Dy),xlabel(39。迭代次數(shù)’）， ylabel(39。電壓 39。),title(39。電壓迭代次數(shù)曲線 39。)。 31 總結(jié) 在電力系統(tǒng)調(diào)度運(yùn)行的多個(gè)領(lǐng)域都涉及到電網(wǎng)潮流計(jì)算。潮流是確定電力網(wǎng)絡(luò)運(yùn)行狀態(tài)的基本因素，潮流問(wèn)題是研究電力系統(tǒng)穩(wěn)態(tài)問(wèn)題的基礎(chǔ)和前提。 潮流計(jì)算靈活性和方便性的要求，對(duì)數(shù)字計(jì)算機(jī)的應(yīng)用也是一個(gè)很關(guān)鍵的問(wèn)題。過(guò)去在很長(zhǎng)時(shí)間內(nèi)，電力系統(tǒng)潮流計(jì)算是借助于交流臺(tái) 進(jìn)行的。交流臺(tái)模擬了電力系統(tǒng)，因此在交流計(jì)算臺(tái)上計(jì)算潮流時(shí)，計(jì)算人員可以隨時(shí)監(jiān)視系統(tǒng)各部分運(yùn)行狀態(tài)是否滿足要求，如發(fā)現(xiàn)某些部分運(yùn)行不合理，則可以立即進(jìn)行調(diào)整。這樣，計(jì)算的過(guò)程就相當(dāng)于運(yùn)算人員丟系統(tǒng)進(jìn)行操作 .調(diào)整的過(guò)程，非常直觀，物理概念也很清楚。當(dāng)利用數(shù)字計(jì)算機(jī)進(jìn)行潮流計(jì)算時(shí)，就失去了這種直觀性。 為了彌補(bǔ)這個(gè)缺點(diǎn)，潮流程序的編制必須盡可能使計(jì)算人員在計(jì)算機(jī)計(jì)算的過(guò)程中加強(qiáng)對(duì)計(jì)算機(jī)過(guò)程的監(jiān)視和控制，并便于作各種修改和調(diào)整。電力系統(tǒng)潮流計(jì)算問(wèn)題并不是單純的計(jì)算問(wèn)題，把它當(dāng)作一個(gè)運(yùn)行方式的調(diào)整問(wèn)題可能更為確切 。為了得到一個(gè)合理的運(yùn)行方式，往往需要不斷根據(jù)計(jì)算結(jié)果，修改原始數(shù)據(jù)。在這個(gè)意義上，我們?cè)诰幹瞥绷饔?jì)算程序時(shí)，對(duì)使用的方便性和靈活性必須予以足夠的重視。因此，除了要求計(jì)算方法盡可能適應(yīng)各種修改 .調(diào)整以外，還要注意輸入和輸出的方便性和靈活性，加強(qiáng)人機(jī)聯(lián)系，以便使計(jì)算人員能及時(shí)監(jiān)視計(jì)算過(guò)程并適當(dāng)?shù)乜刂朴?jì)算的進(jìn)行。