【正文】 力集中在算法的構(gòu)思而不是編程上。不同領(lǐng)域、不同層次的用戶通過相應(yīng)工具的學(xué)習(xí)和應(yīng)用，可以方便地進(jìn)行計(jì)算、分析及設(shè)計(jì)工作。原始數(shù)據(jù)輸入格式的設(shè)計(jì)，主要應(yīng)從使用的角度出發(fā)，原則是簡單明了，便于修改。 基礎(chǔ)內(nèi)容電力系統(tǒng)潮流計(jì)算是研究電力系統(tǒng)穩(wěn)定運(yùn)行情況的一種計(jì)算，它根據(jù)給定的運(yùn)行條件及系統(tǒng)接線情況確定電力系統(tǒng)各部分的運(yùn)行狀態(tài)。用計(jì)算機(jī)進(jìn)行潮流計(jì)算的主要步驟是：建立數(shù)學(xué)模型；確定計(jì)算方法；制定框圖并編制程序；上機(jī)計(jì)算并對(duì)計(jì)算結(jié)果進(jìn)行分析等四步。運(yùn)用計(jì)算機(jī)研究潮流計(jì)算問題，對(duì)數(shù)學(xué)模型的建立和解題方法的確定又必須解決好收斂性、計(jì)算速度和內(nèi)存量等方面的要求。根據(jù)電力系統(tǒng)的特點(diǎn)，對(duì)極坐標(biāo)形式的牛頓拉夫遜的修正方程加以改進(jìn)和簡化發(fā)展形成了PQ分解法。本設(shè)計(jì)以牛頓——拉夫遜直角坐標(biāo)法利用計(jì)算機(jī)求解簡單網(wǎng)絡(luò)的潮流分布為主其主要步驟如下【2】：1． 形成節(jié)點(diǎn)導(dǎo)納矩陣Y。445c。4． 將各節(jié)點(diǎn)的電壓初值代入公式441a；441b；449a；449c求修正方程的系數(shù)矩陣——雅可比矩陣的各個(gè)元熟H（0）、N（0）、J（0）、L（0）以及R（0）、S（0）。6． 計(jì)算各節(jié)點(diǎn)的電壓新值U（i）.7． 運(yùn)用各節(jié)點(diǎn)電壓的新值自地三步開始進(jìn)入下一次跌代。第二章 計(jì)算內(nèi)容解：形成節(jié)點(diǎn)導(dǎo)納矩陣YB由上圖可得相應(yīng)的節(jié)點(diǎn)導(dǎo)納矩陣Y=[ 5+15i + 0 0。 + +5i 10+30i 0。0 + 0 + ]由于節(jié)點(diǎn)1為平衡節(jié)點(diǎn)，其他的節(jié)點(diǎn)均為PQ節(jié)點(diǎn)，系數(shù)矩陣，階數(shù)相同。他們就由導(dǎo)納矩陣的虛部部分中除第一行第一列外的各個(gè)元素所組成計(jì)算各節(jié)點(diǎn)有功功率不平衡量取=，=0；====。節(jié)點(diǎn)總數(shù)為：39。N=5disp(39。)。PQ節(jié)點(diǎn)為：39。JD=[2,3,4,5]Y=[ 5+15i + 0 0。 + +5i 10+30i 0。0 + 0 + ]G=real(Y)。B1=B(2:5, 2:5)。 b1=inv(B1)。dt(2)=0。 dt(4)=0。u(2)=。u(4)=。p(2)=。p(4)=。q(2)=。q(4)=。k=0。N1=4。 dt(1)=0。 for m=2:5 for n=1:5 pt(n)=u(m)*u(n)*(G(m,n)*cos(dt(m)dt(n))+B(m,n)*sin(dt(m)dt(n)))。各節(jié)點(diǎn)有功率不平衡量為：39。 u2=inv(u1(2:5,2:5))。 u4=u3(1:4,:)。 dp2=dp1(1:4,:)。 dtt=dtt(:,1:1)。各節(jié)點(diǎn)電壓的相位角dt為：39。 end disp(39。) qq(m)=q(m)sum(qt) end for m=2:5 qq1=flipud(rot90(qq(2:5)))。ut=b2*qq2/u4。各節(jié)點(diǎn)電壓的修正量：39。各節(jié)點(diǎn)電壓的大小為：39。endfor n=1:5 I(n)=Y(1,n)*U(n)。各節(jié)點(diǎn)功率s為：39。 endenddisp(39。)Swucha=max(abs(ut))k=k+1end、設(shè)計(jì)結(jié)果：節(jié)點(diǎn)總數(shù)為：N = 5平衡節(jié)點(diǎn)為：ans = 1PQ節(jié)點(diǎn)為： JD = 2 3 4 5Y = + + 0 0 + + + + + + + 0 0 + + + 0 + 0 + 各節(jié)點(diǎn)有功率不平衡量為：dp =0 各節(jié)點(diǎn)有功率不平衡量為：dp = 0 各節(jié)點(diǎn)有功率不平衡量為：dp = 0 各節(jié)點(diǎn)有功率不平衡量為：dp = 0 dtt = 各節(jié)點(diǎn)電壓的相位角dt為：dt =0 0 0 0dtt = 各節(jié)點(diǎn)電壓的相位角dt為：dt = 0 0 0dtt = 各節(jié)點(diǎn)電壓的相位角dt為：dt = 0 0dtt = 各節(jié)點(diǎn)電壓的相位角dt為：dt = 0 各節(jié)點(diǎn)無功率不平衡量為：qq = 0 各節(jié)點(diǎn)無功率不平衡量為：qq = 0 各節(jié)點(diǎn)無功率不平衡量為：qq = 0 各節(jié)點(diǎn)無功率不平衡量為：qq =0 各節(jié)點(diǎn)電壓的修正量：ut = 0 0 0 0 0 0 0 0 0 0 0 0各節(jié)點(diǎn)電壓的大小為：u = 各節(jié)點(diǎn)電壓的修正量：ut = 0 0 0 0 0 0 0 0 0 0 0 0各節(jié)點(diǎn)電壓的大小為：u = 各節(jié)點(diǎn)電壓的修正量：ut = 0 0 0 0 0 0 0 0 0 0 0 0各節(jié)點(diǎn)電壓的大小為：u = 各節(jié)點(diǎn)電壓的修正量：ut = 0 0 0 0 0 0 0 0 0 0 0 0各節(jié)點(diǎn)電壓的大小為：u = 各節(jié)點(diǎn)功率s為：S1 = + 各節(jié)點(diǎn)功率sij為：S = 0 + + 0 0 0 + + + 0 0 0 + 0 + 0 0 0 k = 1各節(jié)點(diǎn)有功率不平衡量為：dp = 0 各節(jié)點(diǎn)有功率不平衡量為：dp =0 各節(jié)點(diǎn)有功率不平衡量為：dp = 0 各節(jié)點(diǎn)有功率不平衡量為：dp = 0 dtt = 各節(jié)點(diǎn)電壓的相位角dt為：dt = 0 dtt = 各節(jié)點(diǎn)電壓的相位角dt為：dt = 0 dtt = 各節(jié)點(diǎn)電壓的相位角dt為：dt = 0 dtt = 各節(jié)點(diǎn)電壓的相位角dt為：dt = 0 各節(jié)點(diǎn)無功率不平衡量為：qq = 0 各節(jié)點(diǎn)無功率不平衡量為：qq = 0 各節(jié)點(diǎn)無功率不平衡量為：qq = 0 各節(jié)點(diǎn)無功率不平衡量為：qq =0 各節(jié)點(diǎn)電壓的修正量：ut = 0 0 0 0 0 0 0 0 0 0 0 0各節(jié)點(diǎn)電壓的大小為：u = 各節(jié)點(diǎn)電壓的修正量：ut = 0 0 0 0 0 0 0 0 0 0 0 0各節(jié)點(diǎn)電壓的大小為：u = 各節(jié)點(diǎn)電壓的修正量：ut = 0 0 0