【正文】 sign)。 fscanf(rf,%d%d,amp。n1,amp。m0)。 } fscanf(rf,%d%d,amp。 if(rf==NULL) { fprintf(wf,不能打開文件 !)。 /*正零序結(jié)點數(shù)支路數(shù) */ rf=fopen(d:\,r)。 /*控制循環(huán)變量 */ int kp,sign。 /*正 ,零序支路阻抗 */ int n11[M1],n12[M1],n01[M0],n02[M0],zl1[M1],zl0[M0]。 /*結(jié)點序相電壓 ,支路序相電流*/ float UP1,ZP1,ZP2,ZP0,IF[3],IFX[3]。 /*結(jié)點導(dǎo)納 ,阻抗矩陣 */ float I0[N1]。 void prog5()。 void prog3()。 /*輸入 ,輸出文件指針 */ main() { void prog1()。 附錄 1 開 始 數(shù)組、變量的定義 從文件輸入原始數(shù)據(jù) 形成正、零序網(wǎng)的導(dǎo)納矩陣，并輸出 對導(dǎo)納矩陣求逆形成阻抗矩陣，并輸出 計算各節(jié)點電壓的正常分量，并輸出 故障類型 FType FType=1 FType=2 FType=3 FType=4 K(3)故障 K(1)故障 K(2)故障 K(1， 1)故障 各節(jié)點序電壓 各節(jié)點序電流 各電氣量的相量 輸出計算結(jié)果 結(jié) 束 附錄 2 短 路 計 算 程 序 的 源 程 序 清 單 include include define N1 5 /*正序結(jié)點數(shù) */ define M1 9 /*正序支路數(shù) */ define N0 5 /*零序結(jié)點數(shù) */ define M0 9 /*零序支路數(shù) */ float SR[3][3],SI[3][3]。解算過程雖不一樣但結(jié)果是基本吻合的。 運行結(jié)果分析 經(jīng)核對，程序的運行結(jié)果與指導(dǎo)書和手算結(jié)果基本一致，故本程序能正確地實現(xiàn)故障計算的功能。 程序使用說明 本程序使用簡單，只需按規(guī)定格式輸入電力系統(tǒng)的網(wǎng)絡(luò)信息到一個名為 “ ”的文本文件，再運行可執(zhí)行文件后按屏幕提示輸入即可。 程序的主要功能 （ 1） 形成網(wǎng)絡(luò)的節(jié)點導(dǎo)納矩陣和節(jié)點阻抗矩陣； （ 2） 可以計算各種短路故障： A） 故障口的序電流和相電流； B） 網(wǎng)絡(luò)中各支路的序電流 C） 網(wǎng)絡(luò)中各節(jié)點的 序電壓和相電壓。 程序設(shè)計的總體思想 系統(tǒng)的原始數(shù)據(jù)以文件的形式輸入，程序根據(jù)這些信息形成導(dǎo)納矩陣，再對之求逆形成阻抗矩陣，由阻抗矩陣求得故障口三序電流，進(jìn)而求得各節(jié)點的三序電壓，并由之求得各支路三序電流，最后將序分量合成相分量，完成一個完整的故障電氣量計算。然而，在電力系統(tǒng)故障計算中，節(jié)點的注入電流通常為已知量，用節(jié)點阻抗方程進(jìn)行故障計算比較直觀，解算簡單且 物理意義明確。 Uc ∠ 176。 Uc ∠ 176。 Uc ∠ 137176。 Uc ∠ 150176。 Uc= aUa1+ a178。Ub= a178。 ∠ 106176。 ∠110176。 ∠ 254176。 ∠250176。 ∠ 120176。 ∠120176。 ∠ 240176。 ∠240176。Ic= aIa1+178。 Ib= a178。 化簡序網(wǎng)絡(luò) 對系統(tǒng)網(wǎng)絡(luò)進(jìn)行“星 三角”變換，得下圖： CD正序網(wǎng)絡(luò)化簡圖A零序網(wǎng)絡(luò)化簡圖AD C 計算電流分支系數(shù) 在特定故障點加入單位電流，求出各支路電流，該電流即為支路的電流分支系數(shù)。 （三）、 手工計算系統(tǒng)故障時各電氣量 為了便于驗證程序運行的正確性，在編制程序之前先對系統(tǒng)進(jìn)行手工的故障計算，確定故障時系統(tǒng)的各個電氣量。由假設(shè) 2可知，正、負(fù)序網(wǎng)絡(luò)相同，故在