【正文】 r i=1:5for j=i:5y(j,i)=y(i,j)。y(1,5)=1/(+)。% 輸入原始數(shù)據(jù)，求節(jié)點(diǎn)導(dǎo)納矩陣y(1,2)=1/(+)。為了便于用迭代法解方程組，需要將上述功率方程改寫成功率平衡方程，并對功率平衡方程求偏導(dǎo)，得出對應(yīng)的雅可比矩陣，給未知節(jié)點(diǎn)賦電壓初值，將初值帶入功率平衡方程，得到功率不平衡量，這樣由功率不平衡量、雅可比矩陣、節(jié)點(diǎn)電壓不平衡量（未知的）構(gòu)成了誤差方程，解誤差方程，得到節(jié)點(diǎn)電壓不平衡量，節(jié)點(diǎn)電壓加上節(jié)點(diǎn)電壓不平衡量構(gòu)成節(jié)點(diǎn)電壓新的初值，將新的初值帶入原來的功率平衡方程，并重新形成雅可比矩陣，然后計(jì)算新的電壓不平衡量，這樣不斷迭代，不斷修正，一般迭代三到五次就能收斂。計(jì)算精度要求各節(jié)點(diǎn)電壓修正量不大于106。潮流計(jì)算的計(jì)算機(jī)算法包含高斯—賽德爾迭代法、牛頓拉夫遜法和P—Q分解法等，其中牛拉法計(jì)算原理較簡單、計(jì)算過程也不復(fù)雜，而且由于人們引入泰勒級(jí)數(shù)和非線性代數(shù)方程等在算法里從而進(jìn)一步提高了算法的收斂性和計(jì)算速度。i=0。} } if(L==0)break。for(i=0。int i,j,k。float BB[3]={4,6,5}。重點(diǎn)學(xué)習(xí)如何聲明和調(diào)用子程序。for(i=1。39。} for(i=0。sL[i].R=d1。for(i=0。sB[i].GenQ=d4。isB[i].Num=i1。d1=d2=d3=d4=d5=d6=ps=qs=。struct Bus sB[NBUS]。int i,j,k,l,h。float X[NBUS]。struct Shunt { int Num,NumI。struct Bus { int Num。以《電力系統(tǒng)穩(wěn)態(tài)分析》教材124頁例3－6題的數(shù)據(jù)存儲(chǔ)格式。2． 將我提供的程序拷貝到d:lx目錄下，并將Options――Directories中的輸出目錄改為d:lx，就可以運(yùn)行。7． 文件的讀寫目錄用“”表示。3． Tc數(shù)組或結(jié)構(gòu)體從0開始計(jì)數(shù)。4． 對Tc環(huán)境的掌握：打開文件，保存文件（F2），另存文件，斷點(diǎn)的設(shè)置（Ctrl+F8），調(diào)試變量的設(shè)置（Ctrl+F7），單步調(diào)試（F7），路徑的設(shè)置。各種條件判斷語句，循環(huán)語句的使用。掌握“*.txt”文件的打開，讀寫，關(guān)閉的函數(shù)。學(xué)會(huì)看幫助（F1），在調(diào)試過程中結(jié)束程序運(yùn)行（Ctrl+F2）5． 《電力系統(tǒng)穩(wěn)態(tài)分析》教材中第四章涉及到的基本知識(shí)。4． 循環(huán)語句只對它后面的第一語句起作用，如果想對多個(gè)語句起作用，則用{}括起來。8． 建議用《電力系統(tǒng)穩(wěn)態(tài)分析》書中的例43 來校驗(yàn)結(jié)果是否正確。（我提供的四個(gè)文件，；導(dǎo)納陣形成文件需要的頭文件,，潮流計(jì)算的原始數(shù)據(jù)；,用高斯消去法解線方程組。5,7,0 1,2 2,0 3,0 4,0 5,0 1 1 3 2 3 4 3 4 5 4 2 5 電力系統(tǒng)潮流上機(jī)課程設(shè)計(jì)教學(xué)指導(dǎo)書 1 2 6 2 3 7 2 4 4． 節(jié)點(diǎn)導(dǎo)納矩陣的形成 struct Line { int Num,NumI,NumJ。float Volt,Phase,GenP,GenQ,LoadP,LoadQ。float G,B。int L。int i1,i2,i3,kp,kq。struct Shunt sSH[NBUS]。e=。sB[i].Volt=d1。sB[i].LoadP=d5,sB[i].LoadQ=d6。isL[i].Num=i1。sL[i].X=d2。i電力系統(tǒng)潮流上機(jī)課程設(shè)計(jì)教學(xué)指導(dǎo)書/*Make Y Matrix*/ for(i=1。 n”)。)fprintf(fp,“Y(%3d,%3d)=(%,%)n”,i,j,YG[i][j],YB[i][j])。在潮流計(jì)算中要把高斯消去法子程序放到潮流計(jì)算主程序后面。int NN=3。float d,t,x[NBUS]。id){ /*在一行中找到一個(gè)最大值賦值d,并用JS[K]記住這個(gè)最大值所在的列號(hào)*/ d=fabs(a[k][j])。for(j=k+1。i){ t=。同時(shí)基于MATLAB的計(jì)算機(jī)算法以雙精度類型進(jìn)行數(shù)據(jù)的存儲(chǔ)和運(yùn)算, 數(shù)據(jù)精確度高，能進(jìn)行潮流計(jì)算中的各種矩陣運(yùn)算，使得傳統(tǒng)潮流計(jì)算方法更加優(yōu)化。二 牛頓拉夫遜法潮流計(jì)算 1 基本原理牛頓法是取近似解x(k)之后，在這個(gè)基礎(chǔ)上，找到比x(k)更接近的方程的根，一步步地迭代，找到盡可能接近方程根的近似根。2 基本步驟和設(shè)計(jì)流程圖形成了雅克比矩陣并建立了修正方程式，運(yùn)用牛頓拉夫遜法計(jì)算潮流的核心問題已經(jīng)解決，已有可能列出基本計(jì)算步驟并編制流程圖。y(4,5)=0。y(2,3)=1/(+)。endendY=0。B=imag(Y)。S(4)=。% 賦初值，U為節(jié)點(diǎn)電壓的幅值，a為節(jié)點(diǎn)電壓的相位角U=ones(1,5)。x2=ones(8,1)。M(i,j)=0。endend% 求有功、無功功率不平衡量for i=1:4for j=1:5oP(i)=oP(i)U(i)*U(j)*(G(i,j)*cos(a(i)a(j))+B(i,j)*sin(a(i)a(j)))。endx2=[oP,oQ]39。M(i,j)=N(i,j)。L(i,i)=L(i,i)U(i)*U(j)*(G(i,j)*sin(a(i)a(j))B(i,j)*cos(a(i)a(j)))endendN(i,i)=N(i,i)2*(U(i))^2*G(i,i)。% x1為所求△x的列向量% 求節(jié)點(diǎn)電壓新值，準(zhǔn)備下一次迭代for i=1:4oa(i)=x1(i)。endk=k+1。endS(5)=P(5)+Q(5)*sqrt(1)。利用x(*)=x(k+1)+DX(k+1)進(jìn)行多次迭代，通過迭代判據(jù)得到所需要的精度值即準(zhǔn)確值x(*)。電力系統(tǒng)潮流的計(jì)算和分析是電力系統(tǒng)運(yùn)行和規(guī)劃工作的基礎(chǔ)。在一些應(yīng)用場所，如輸電網(wǎng)絡(luò)中，只要考慮的是電力系統(tǒng)中有功功率的分布，而不需要計(jì)算各個(gè)節(jié)點(diǎn)電壓幅值，且要計(jì)算速度要快，這勢必要對潮流計(jì)算進(jìn)行簡化處理，本節(jié)實(shí)驗(yàn)就是研究直流潮流計(jì)算，編程與調(diào)試，獲得電力系統(tǒng)中各支路的有功分布，為進(jìn)一步進(jìn)行電力系統(tǒng)分析作準(zhǔn)備。雖然這些方法都具有一定的精度，但計(jì)算量較大，這顯然不適應(yīng)形成電網(wǎng)規(guī)劃方案時(shí)多次而反復(fù)的潮流計(jì)算要求。相當(dāng)于注入的有功功率。j206。i,j185。x239。237。式（4）為直流潮流方程，因?yàn)楹雎粤私拥刂?，同時(shí)忽略了支路電阻，所以沒有有功功率損耗。三、直流潮流算法計(jì)算步驟1）選擇平衡節(jié)點(diǎn)。15）通過式（2）的潮流計(jì)算公式求出各條支路的功率情況。實(shí)驗(yàn)數(shù)據(jù)如圖所示三母線電力系統(tǒng)中，支路電抗和節(jié)點(diǎn)注入的攻入如圖所示，編寫程序，求個(gè)各條支路的有功潮流分布。潮流計(jì)算可以選擇電力系統(tǒng)的運(yùn)行方式，便于定期對電力系統(tǒng)中的元件進(jìn)行檢修。④輸電線路資料：發(fā)電廠和變電所之間的輸電線路的電壓等級(jí)及長度標(biāo)于圖中，*106s。（；電壓35kv要求調(diào)整范圍在3536之間）。并且，將節(jié)點(diǎn)1設(shè)為平衡節(jié)點(diǎn)，將節(jié)點(diǎn)10設(shè)為pv節(jié)點(diǎn)，其余節(jié)點(diǎn)設(shè)為pq節(jié)點(diǎn)。由matlab編程調(diào)節(jié)后，可得到表4的發(fā)電廠電壓和變壓器分接頭電壓得標(biāo)幺值。因此我們可以看出，經(jīng)過調(diào)節(jié)后，節(jié)點(diǎn)9點(diǎn)電壓已經(jīng)滿足了系統(tǒng)的要求。 對比由上面的三種方法簡單地比較，我們可以看出，在同一個(gè)電力系統(tǒng)中，用不同的方法進(jìn)行潮流計(jì)算，所得