南昌大學(xué)電力系統(tǒng)分析實(shí)驗(yàn)報(bào)告

1、本科生實(shí)驗(yàn)報(bào)告 課程名稱： 電力系統(tǒng)分析實(shí)驗(yàn) 專業(yè)班級(jí)： 電力系統(tǒng)124 班 姓 名： 學(xué) 號(hào)： 所在學(xué)期： 2014-2015-2 2015年 6 月 20 日實(shí)驗(yàn)四 配電網(wǎng)潮流計(jì)算實(shí)驗(yàn)一、實(shí)驗(yàn)?zāi)康?本實(shí)驗(yàn)針對配電網(wǎng)具有輻射式開環(huán)結(jié)構(gòu)的特性，基于前推回代的方法編制程序使系統(tǒng)潮流計(jì)算能夠由計(jì)算機(jī)自行完成，即已知配電網(wǎng)首端節(jié)點(diǎn)電壓和各末端節(jié)點(diǎn)的功率，由末端向首端計(jì)算網(wǎng)絡(luò)功率分布，再由首端向末端計(jì)算節(jié)點(diǎn)電壓分布，迭代直至收斂條件。通過實(shí)驗(yàn)教學(xué)加深學(xué)生對配電網(wǎng)潮流計(jì)算方法的理解，掌握樹狀輻射式配電網(wǎng)潮流計(jì)算的基本算法。熟悉各種常用應(yīng)用軟件，熟悉硬件設(shè)備的使用方法，加強(qiáng)編制調(diào)試計(jì)算機(jī)程序的能力，提高工

2、程計(jì)算的能力，學(xué)習(xí)如何將理論知識(shí)和實(shí)際工程問題結(jié)合起來。二、實(shí)驗(yàn)器材 計(jì)算機(jī)、軟件（已安裝，包括各類編程軟件C語言、C+、VB、VC等、應(yīng)用軟件MATLAB等）、移動(dòng)存儲(chǔ)設(shè)備（學(xué)生自備，軟盤、U盤等）三、實(shí)驗(yàn)內(nèi)容 編制調(diào)試配電網(wǎng)潮流計(jì)算的計(jì)算機(jī)程序。程序要求根據(jù)已知的配電網(wǎng)參數(shù)，完成該配電系統(tǒng)的潮流計(jì)算，要求計(jì)算出節(jié)點(diǎn)電壓、功率等參數(shù)。1.先編制好的電力系統(tǒng)潮流計(jì)算的計(jì)算程序原代碼由自備移動(dòng)存儲(chǔ)設(shè)備導(dǎo)入計(jì)算機(jī)。2.應(yīng)的編程環(huán)境下對程序進(jìn)行組織調(diào)試。3.應(yīng)用計(jì)算例題驗(yàn)證程序的計(jì)算效果。4.對調(diào)試正確的計(jì)算程序進(jìn)行存儲(chǔ)、打印。5.完成本次實(shí)驗(yàn)的實(shí)驗(yàn)報(bào)告。四、實(shí)驗(yàn)數(shù)據(jù)如圖所示一個(gè)5節(jié)點(diǎn)的配電網(wǎng)系統(tǒng)，

3、S1=S2=2+j2(MVA)，S3=S4=S5=4+j4(MVA)，假定所有負(fù)荷均為恒功率負(fù)荷，節(jié)點(diǎn)1為參考節(jié)點(diǎn)相角為0o。計(jì)算電網(wǎng)個(gè)節(jié)點(diǎn)電壓及支路功率分布，收斂條件為。5、 實(shí)驗(yàn)程序% 配電網(wǎng)前推回代潮流計(jì)算程序 % 使用IEEE 33節(jié)點(diǎn)配電系統(tǒng)作為算例，可實(shí)現(xiàn)弱環(huán)網(wǎng)情況下的潮流計(jì)算 countnum=0;BranchData = 1 2 0.0922 0.0470; 2 3 0.4930 0.2511; 3 4 0.3660 0.1864; 4 5 0.3811 0.1941; 5 6 0.8190 0.7070; 6 7 0.1872 0.6188; 7 8 0.7114 0.235

4、1; 8 9 1.0300 0.7400; 9 10 1.0440 0.7400; 10 11 0.1966 0.0650; 11 12 0.3744 0.1238; 12 13 1.4680 1.1550; 13 14 0.5416 0.7129; 14 15 0.5910 0.5260; 15 16 0.7463 0.5450; 16 17 1.2890 1.7210; 17 18 0.7320 0.5740; 2 19 0.1640 0.1565; 19 20 1.5042 1.3554; 20 21 0.4095 0.4784; 21 22 0.7089 0.9373; 3 23 0.

5、4512 0.3083; 23 24 0.8980 0.7091; 24 25 0.8960 0.7011; 6 26 0.2030 0.1034; 26 27 0.2842 0.1447; 27 28 1.0590 0.9337; 28 29 0.8042 0.7006; 29 30 0.5075 0.2585; 30 31 0.9744 0.9630; 31 32 0.3105 0.3619; 32 33 0.3410 0.5302; % 支路，阻抗 NodeData = 2 100.00 60.00; 3 90.00 40.00; 4 120.00 80.00; 5 60.00 30.0

6、0; 6 60.00 20.00; 7 200.00 100.00; 8 200.00 100.00; 9 60.00 20.00; 10 60.00 20.00; 11 45.00 30.00; 12 60.00 35.00; 13 60.00 35.00; 14 120.00 80.00; 15 60.00 10.00; 16 60.00 20.00; 17 60.00 20.00; 18 90.00 40.00; 19 90.00 40.00; 20 90.00 40.00; 21 90.00 40.00; 22 90.00 40.00; 23 90.00 50.00; 24 420.0

7、0 200.00; 25 420.00 200.00; 26 60.00 25.00; 27 60.00 25.00; 28 60.00 20.00; 29 120.00 70.00; 30 200.00 600.00; 31 150.00 70.00; 32 210.00 100.00; 33 60.00 40.00; % 節(jié)點(diǎn)，負(fù)荷 UB = 12.66; % 電壓基準(zhǔn) kVSB = 10; % 功率基準(zhǔn) MVAZB = UB2/SB; % 阻抗基準(zhǔn) ohm BranchData(:,3,4) = BranchData(:,3,4) / ZB; % 阻抗標(biāo)幺化NodeData(:,2,3)

8、 = NodeData(:,2,3) / SB / 1000;% 功率標(biāo)幺化 NN = 33; % 節(jié)點(diǎn)數(shù) A0 = zeros(NN);for n = 1:NN-1 A0(BranchData(n,1),BranchData(n,2) = 1;end % 形成 A0 AssociatedMatrix=0; for n=2:NN-1 AssociatedMatrix(n,n)=1; temp=BranchData(n-1,1); AssociatedMatrix(n,1:n-1)=AssociatedMatrix(temp,1:n-1);end A0T = A0; % 形成 A0 的轉(zhuǎn)置S =

9、 0;-NodeData(:,2) - i*NodeData(:,3); % 形成 SZL = 0;BranchData(:,3) + i*BranchData(:,4); % 形成 ZL V = ones(NN,1);V(1) = 1; % 各個(gè)節(jié)點(diǎn)電壓賦初值 IL(NN,1) = -conj(S(NN) / V(NN); % 最末支路電流賦初值 Delta = 1; % 收斂判據(jù)賦初值TempV = V; % 賦初值，用于記憶上次迭代結(jié)果 while Delta 1e-8 countnum=countnum+1; IN = conj(S ./ V); % 節(jié)點(diǎn)注入電流 for n = 1:

10、NN-1 IL(NN-n) = A0(NN-n,NN-n+1:end) * IL(NN-n+1:end) - IN(NN-n); end % 電流回代過程 for n = 2:NN V(n) = A0T(n,1:n-1) * V(1:n-1) - ZL(n) * IL(n); end % 電壓前推過程 Delta = max(abs(V-TempV); % 更新收斂判據(jù) TempV = V; % 記憶迭代結(jié)果end Vangle(:,1)=abs(V);Vangle(:,2)=angle(V)/3.1415*180; for i=1:NN-1 st=BranchData(i,1); en=Br

11、anchData(i,2); Sij(i,1)=V(st)*conj(V(st)-V(en)/ZL(i+1); Sji(i,1)=V(en)*conj(V(en)-V(st)/ZL(i+1);End實(shí)驗(yàn)五 輸電網(wǎng)潮流計(jì)算實(shí)驗(yàn)一、實(shí)驗(yàn)?zāi)康?本實(shí)驗(yàn)針對高壓輸電網(wǎng)絡(luò)具有多環(huán)、多電源的結(jié)構(gòu)特性，要求根據(jù)功率平衡方程的泰勒展開式近似的求取系統(tǒng)的狀態(tài)變量，使系統(tǒng)潮流計(jì)算能夠由計(jì)算機(jī)自行完成，即根據(jù)已知的電力網(wǎng)的數(shù)學(xué)模型（節(jié)點(diǎn)導(dǎo)納矩陣）及各節(jié)點(diǎn)參數(shù)，由計(jì)算程序運(yùn)行完成該電力系統(tǒng)的潮流計(jì)算。通過實(shí)驗(yàn)教學(xué)加深學(xué)生對復(fù)雜電力系統(tǒng)潮流計(jì)算計(jì)算方法的理解，學(xué)會(huì)運(yùn)用電力系統(tǒng)的數(shù)學(xué)模型，掌握潮流計(jì)算的過程及其特點(diǎn)，熟悉各

12、種常用應(yīng)用軟件，熟悉硬件設(shè)備的使用方法，加強(qiáng)編制調(diào)試計(jì)算機(jī)程序的能力，提高工程計(jì)算的能力，學(xué)習(xí)如何將理論知識(shí)和實(shí)際工程問題結(jié)合起來。二、實(shí)驗(yàn)器材 計(jì)算機(jī)、軟件（已安裝，包括各類編程軟件C語言、C+、VB、VC等、應(yīng)用軟件MATLAB等）、移動(dòng)存儲(chǔ)設(shè)備（學(xué)生自備，軟盤、U盤等）三、實(shí)驗(yàn)內(nèi)容 編制調(diào)試電力系統(tǒng)潮流計(jì)算的計(jì)算機(jī)程序。程序要求根據(jù)已知的電網(wǎng)的數(shù)學(xué)模型（節(jié)點(diǎn)導(dǎo)納矩陣）及各節(jié)點(diǎn)參數(shù)，完成該電力系統(tǒng)的潮流計(jì)算，要求計(jì)算出節(jié)點(diǎn)電壓、功率等參數(shù)。1.采用牛頓拉夫遜法編制程序編制潮流程序。2.先編制好的電力系統(tǒng)潮流計(jì)算的計(jì)算程序原代碼由自備移動(dòng)存儲(chǔ)設(shè)備導(dǎo)入計(jì)算機(jī)。3.應(yīng)的編程環(huán)境下對程序進(jìn)行組織調(diào)

13、試。4.應(yīng)用計(jì)算例題驗(yàn)證程序的計(jì)算效果。5.對調(diào)試正確的計(jì)算程序進(jìn)行存儲(chǔ)、打印。6.完成本次實(shí)驗(yàn)的實(shí)驗(yàn)報(bào)告。4、 實(shí)驗(yàn)程序%牛頓拉夫遜法進(jìn)行潮流計(jì)算% B1矩陣：1、支路首端號(hào)；2、末端號(hào)；3、支路阻抗；4、支路對地電納% 5、支路的變比；6、支路首端處于K側(cè)為1，1側(cè)為0% B2矩陣：1、該節(jié)點(diǎn)發(fā)電機(jī)功率；2、該節(jié)點(diǎn)負(fù)荷功率；3、節(jié)點(diǎn)電壓初始值% 4、PV節(jié)點(diǎn)電壓V的給定值；5、節(jié)點(diǎn)所接的無功補(bǔ)償設(shè)備的容量% 6、節(jié)點(diǎn)分類標(biāo)號(hào)：1為平衡節(jié)點(diǎn)；2為PQ節(jié)點(diǎn)；% 3為PV節(jié)點(diǎn)；clear;clc;n=4;%input(請輸入節(jié)點(diǎn)數(shù)：n=);nl=4;%input(請輸入支路數(shù)：nl=);isb=

14、1;%input(請輸入平衡母線節(jié)點(diǎn)號(hào)：isb=);pr=0.00001;%input(請輸入誤差精度：pr=);B1=2 4 0.1+0.4i 0.01528*2i 1 0 ; 3 4 0.3i 0 1./1.1 0 ; 1 4 0.12+0.5i 0.01920*2i 1 0; 1 2 0.08+0.4i 0.01413*2i 1 0 ; %input(請輸入各節(jié)點(diǎn)參數(shù)形成的矩陣：B2=);B2=0 0 1.05 1.05 0 1; 0 0.3+0.18i 1 0 0 2; 0.5 0 1.1 1.1 0 3; 0 0.55+0.13i 1 0 0 2;Y=zeros(n);e=zeros

15、(1,n);f=zeros(1,n);V=zeros(1,n);sida=zeros(1,n);S1=zeros(nl);% % %-for i=1:nl%支路數(shù) if B1(i,6)=0%左節(jié)點(diǎn)處于1側(cè) p=B1(i,1);q=B1(i,2); else %左節(jié)點(diǎn)處于K側(cè) p=B1(i,2);q=B1(i,1); end Y(p,q)=Y(p,q)-1./(B1(i,3)*B1(i,5);%非對角元 Y(q,p)=Y(p,q); %非對角元 Y(q,q)=Y(q,q)+1./(B1(i,3)*B1(i,5)2)+B1(i,4)./2;%對角元K側(cè) Y(p,p)=Y(p,p)+1./B1(i,

16、3)+B1(i,4)./2; %對角元1側(cè) end%求導(dǎo)納矩陣disp(導(dǎo)納矩陣 Y=);disp(Y)%-G=real(Y);B=imag(Y);%分解出導(dǎo)納陣的實(shí)部和虛部 for i=1:n%給定各節(jié)點(diǎn)初始電壓的實(shí)部和虛部 e(i)=real(B2(i,3); f(i)=imag(B2(i,3); V(i)=B2(i,4);%PV節(jié)點(diǎn)電壓給定模值 endfor i=1:n%給定各節(jié)點(diǎn)注入功率 S(i)=B2(i,1)-B2(i,2); %i節(jié)點(diǎn)注入功率SG-SL B(i,i)=B(i,i)+B2(i,5);%i節(jié)點(diǎn)無功補(bǔ)償量 end%-P=real(S);Q=imag(S); %分解出各節(jié)

17、點(diǎn)注入的有功和無功功率ICT1=0;IT2=1;N0=2*n;N=N0+1;a=0; %迭代次數(shù)ICT1、a；不滿足收斂要求的節(jié)點(diǎn)數(shù)IT2while IT2=0 % N0=2*n 雅可比矩陣的階數(shù)；N=N0+1擴(kuò)展列 IT2=0;a=a+1; for i=1:n if i=isb%非平衡節(jié)點(diǎn) C(i)=0;D(i)=0; for j1=1:n C(i)=C(i)+G(i,j1)*e(j1)-B(i,j1)*f(j1);%(Gij*ej-Bij*fj) D(i)=D(i)+G(i,j1)*f(j1)+B(i,j1)*e(j1);%(Gij*fj+Bij*ej) end P1=C(i)*e(i)+

18、f(i)*D(i);%節(jié)點(diǎn)功率P計(jì)算ei(Gij*ej-Bij*fj)+fi(Gij*fj+Bij*ej) Q1=C(i)*f(i)-e(i)*D(i);%節(jié)點(diǎn)功率Q計(jì)算fi(Gij*ej-Bij*fj)-ei(Gij*fj+Bij*ej)%求i節(jié)點(diǎn)有功和無功功率P,Q的計(jì)算值 V2=e(i)2+f(i)2;%電壓模平方 %- 以下針對非PV節(jié)點(diǎn)來求取功率差及Jacobi矩陣元素 - if B2(i,6)=3%非PV節(jié)點(diǎn) DP=P(i)-P1;%節(jié)點(diǎn)有功功率差 DQ=Q(i)-Q1; %節(jié)點(diǎn)無功功率差 %- 以上為除平衡節(jié)點(diǎn)外其它節(jié)點(diǎn)的功率計(jì)算 -%-求取Jacobi矩陣 - for j1=1

19、:n if j1=isb&j1=i%非平衡節(jié)點(diǎn)&非對角元 X1=-G(i,j1)*e(i)-B(i,j1)*f(i);% dP/de=-dQ/df X2=B(i,j1)*e(i)-G(i,j1)*f(i);% dP/df=dQ/de X3=X2; % X2=dp/df X3=dQ/de X4=-X1; % X1=dP/de X4=dQ/df p=2*i-1;q=2*j1-1; J(p,q)=X3;J(p,N)=DQ;m=p+1; % X3=dQ/de J(p,N)=DQ節(jié)點(diǎn)無功功率差 J(m,q)=X1;J(m,N)=DP;q=q+1; % X1=dP/de J(m,N)=DP節(jié)點(diǎn)有功功率差

20、J(p,q)=X4;J(m,q)=X2; % X4=dQ/df X2=dp/df elseif j1=i&j1=isb%非平衡節(jié)點(diǎn)&對角元 X1=-C(i)-G(i,i)*e(i)-B(i,i)*f(i);% dP/de X2=-D(i)+B(i,i)*e(i)-G(i,i)*f(i);% dP/df X3=D(i)+B(i,i)*e(i)-G(i,i)*f(i); % dQ/de X4=-C(i)+G(i,i)*e(i)+B(i,i)*f(i);% dQ/df p=2*i-1;q=2*j1-1;J(p,q)=X3;J(p,N)=DQ;%擴(kuò)展列Q m=p+1; J(m,q)=X1;q=q+1;

21、J(p,q)=X4;J(m,N)=DP;%擴(kuò)展列P J(m,q)=X2; end end else%-下面是針對PV節(jié)點(diǎn)來求取Jacobi矩陣的元素 - DP=P(i)-P1;% PV節(jié)點(diǎn)有功誤差 DV=V(i)2-V2;% PV節(jié)點(diǎn)電壓誤差 for j1=1:n if j1=isb&j1=i%非平衡節(jié)點(diǎn)&非對角元 X1=-G(i,j1)*e(i)-B(i,j1)*f(i); % dP/de X2=B(i,j1)*e(i)-G(i,j1)*f(i); % dP/df X5=0;X6=0; p=2*i-1;q=2*j1-1;J(p,q)=X5;J(p,N)=DV; % PV節(jié)點(diǎn)電壓誤差 m=p+

22、1; J(m,q)=X1;J(m,N)=DP;q=q+1;J(p,q)=X6; % PV節(jié)點(diǎn)有功誤差 J(m,q)=X2; elseif j1=i&j1=isb %非平衡節(jié)點(diǎn)&對角元 X1=-C(i)-G(i,i)*e(i)-B(i,i)*f(i);% dP/de X2=-D(i)+B(i,i)*e(i)-G(i,i)*f(i);% dP/df X5=-2*e(i); X6=-2*f(i); p=2*i-1;q=2*j1-1;J(p,q)=X5;J(p,N)=DV; % PV節(jié)點(diǎn)電壓誤差 m=p+1; J(m,q)=X1;J(m,N)=DP;q=q+1;J(p,q)=X6; % PV節(jié)點(diǎn)有功誤

23、差 J(m,q)=X2; end end end end end%- 以上為求雅可比矩陣的各個(gè)元素及擴(kuò)展列的功率差或電壓差 - for k=3:N0 % N0=2*n （從第三行開始，第一、二行是平衡節(jié)點(diǎn)） k1=k+1;N1=N; % N=N0+1 即 N=2*n+1擴(kuò)展列P、Q 或 U for k2=k1:N1% 從k+1列的Jacobi元素到擴(kuò)展列的P、Q 或 U J(k,k2)=J(k,k2)./J(k,k);% 用K行K列對角元素去除K行K列后的非對角元素進(jìn)行規(guī)格化 end J(k,k)=1; % 對角元規(guī)格化K行K列對角元素賦1 %- 回代運(yùn)算 - if k=3 % 不是第三行 k

24、 3 k4=k-1; for k3=3:k4% 用k3行從第三行開始到當(dāng)前行的前一行k4行消去 for k2=k1:N1% k3行后各行上三角元素 J(k3,k2)=J(k3,k2)-J(k3,k)*J(k,k2);%消去運(yùn)算（當(dāng)前行k列元素消為0） end %用當(dāng)前行K2列元素減去當(dāng)前行k列元素乘以第k行K2列元素 J(k3,k)=0; %當(dāng)前行第k列元素已消為0 end if k=N0 %若已到最后一行 break; end %- 前代運(yùn)算 - for k3=k1:N0 % 從k+1行到2*n最后一行 for k2=k1:N1 % 從k+1列到擴(kuò)展列消去k+1行后各行下三角元素 J(k3,

25、k2)=J(k3,k2)-J(k3,k)*J(k,k2);%消去運(yùn)算 end %用當(dāng)前行K2列元素減去當(dāng)前行k列元素乘以第k行K2列元素 J(k3,k)=0; %當(dāng)前行第k列元素已消為0 end else %是第三行k=3 %- 第三行k=3的前代運(yùn)算 - for k3=k1:N0 %從第四行到2n行（最后一行） for k2=k1:N1 %從第四列到2n+1列（即擴(kuò)展列） J(k3,k2)=J(k3,k2)-J(k3,k)*J(k,k2);%消去運(yùn)算（當(dāng)前行3列元素消為0） end %用當(dāng)前行K2列元素減去當(dāng)前行3列元素乘以第三行K2列元素 J(k3,k)=0; %當(dāng)前行第3列元素已消為0

26、end end end%-上面是用線性變換方式高斯消去法將Jacobi矩陣化成單位矩陣- for k=3:2:N0-1 L=(k+1)./2; e(L)=e(L)-J(k,N); %修改節(jié)點(diǎn)電壓實(shí)部 k1=k+1; f(L)=f(L)-J(k1,N); %修改節(jié)點(diǎn)電壓虛部 end %-修改節(jié)點(diǎn)電壓- for k=3:N0 DET=abs(J(k,N); if DET=pr %電壓偏差量是否滿足要求 IT2=IT2+1; %不滿足要求的節(jié)點(diǎn)數(shù)加1 end end ICT2(a)=IT2; %不滿足要求的節(jié)點(diǎn)數(shù) ICT1=ICT1+1; %迭代次數(shù)end%用高斯消去法解w=-J*Vdisp(迭代次

27、數(shù)：);disp(ICT1);disp(沒有達(dá)到精度要求的個(gè)數(shù)：);disp(ICT2);for k=1:n V(k)=sqrt(e(k)2+f(k)2); %計(jì)算各節(jié)點(diǎn)電壓的模值 sida(k)=atan(f(k)./e(k)*180./pi; %計(jì)算各節(jié)點(diǎn)電壓的角度 E(k)=e(k)+f(k)*j; %將各節(jié)點(diǎn)電壓用復(fù)數(shù)表示end%- 計(jì)算各輸出量 -disp(各節(jié)點(diǎn)的實(shí)際電壓標(biāo)幺值E為(節(jié)點(diǎn)號(hào)從小到大排列)：);disp(E); %顯示各節(jié)點(diǎn)的實(shí)際電壓標(biāo)幺值E用復(fù)數(shù)表示disp(-);disp(各節(jié)點(diǎn)的電壓大小V為(節(jié)點(diǎn)號(hào)從小到大排列)：);disp(V); %顯示各節(jié)點(diǎn)的電壓大小V的

28、模值disp(-);disp(各節(jié)點(diǎn)的電壓相角sida為(節(jié)點(diǎn)號(hào)從小到大排列)：);disp(sida); %顯示各節(jié)點(diǎn)的電壓相角for p=1:n C(p)=0; for q=1:n C(p)=C(p)+conj(Y(p,q)*conj(E(q); %計(jì)算各節(jié)點(diǎn)的注入電流的共軛值 end S(p)=E(p)*C(p); %計(jì)算各節(jié)點(diǎn)的功率 S = 電壓 X 注入電流的共軛值enddisp(各節(jié)點(diǎn)的功率S為(節(jié)點(diǎn)號(hào)從小到大排列)：);disp(S); %顯示各節(jié)點(diǎn)的注入功率 五、實(shí)驗(yàn)結(jié)果 在上圖所示的簡單電力系統(tǒng)中，系統(tǒng)中節(jié)點(diǎn)1、2為PQ節(jié)點(diǎn)，節(jié)點(diǎn)3為PV節(jié)點(diǎn)，節(jié)點(diǎn)4為平衡節(jié)點(diǎn)，已給定P1s+

29、jQ1s=-0.30-j0.18 P2s+jQ2s=-0.55-j0.13 P3s=0.5 V3s=1.10 V4s=1.050容許誤差=10-5節(jié)點(diǎn)導(dǎo)納矩陣：各節(jié)點(diǎn)電壓： 節(jié)點(diǎn) e f v 1. 0. -0. 0. -0. 2. 0. -0. 0. -6. 3. 1. 0. 1. 6. 4. 1. 0. 1. 0.各節(jié)點(diǎn)功率： 節(jié)點(diǎn) P Q 1 -0. -0. 2 0. -0. 3 0. -0. 4 0. 0.導(dǎo)納矩陣 Y=1.0421 - 8.2429i-0.5882 + 2.3529i0.0000 + 3.6667i-0.4539 + 1.8911i-0.5882 + 2.3529i1.

30、0690 - 4.7274i 0.0000 + 0.0000i-0.4808 + 2.4038i0.0000 + 3.6667i0.0000 + 0.0000i 0.0000 - 3.3333i0.0000 + 0.0000i-0.4539 + 1.8911i-0.4808 + 2.4038i 0.0000 + 0.0000i0.9346 - 4.2616i 迭代次數(shù)： 4各節(jié)點(diǎn)的實(shí)際電壓標(biāo)幺值E為(節(jié)點(diǎn)號(hào)從小到大排列)： 0.9854 - 0.0485i 0.9691 - 0.0759i 1.0968 + 0.0844i 1.0500 + 0.0000i各節(jié)點(diǎn)的電壓大小e為(節(jié)點(diǎn)號(hào)從小到大排

31、列)： 0.9854 0.9691 1.0968 1.0500各節(jié)點(diǎn)的電壓大小f為(節(jié)點(diǎn)號(hào)從小到大排列)： -0.0485 -0.0759 0.0844 0各節(jié)點(diǎn)的電壓大小V為(節(jié)點(diǎn)號(hào)從小到大排列)： 0.9866 0.9721 1.1000 1.0500各節(jié)點(diǎn)的電壓相角sida為(節(jié)點(diǎn)號(hào)從小到大排列)： -2.8167 -4.4784 4.4018 0各節(jié)點(diǎn)的功率S為(節(jié)點(diǎn)號(hào)從小到大排列)： -0.5500 - 0.1300i -0.3000 - 0.1800i 0.5000 + 0.0857i 0.3595 + 0.2343i 實(shí)驗(yàn)六 水火電廠間有功功率負(fù)荷的經(jīng)濟(jì)分配一、實(shí)驗(yàn)?zāi)康?本實(shí)驗(yàn)通

32、過對水火電廠間有功功率負(fù)荷的經(jīng)濟(jì)分配的計(jì)算機(jī)程序編制與調(diào)試，獲得能計(jì)算有功負(fù)荷經(jīng)濟(jì)分配的計(jì)算機(jī)程序，即根據(jù)已知的水火電廠的耗量特性、發(fā)電用水量及各時(shí)間段的有功負(fù)荷功率，由計(jì)算機(jī)程序自動(dòng)完成各時(shí)段各電廠的有功出力的計(jì)算。通過實(shí)驗(yàn)教學(xué)加深學(xué)生對電力系統(tǒng)經(jīng)濟(jì)運(yùn)行概念的理解，學(xué)會(huì)運(yùn)用數(shù)學(xué)模型進(jìn)行水火電廠間的有功功率負(fù)荷的經(jīng)濟(jì)分配計(jì)算，熟悉各種常用應(yīng)用軟件，熟悉硬件設(shè)備的使用方法，加強(qiáng)編制調(diào)試計(jì)算機(jī)程序的能力，提高工程計(jì)算的能力，學(xué)習(xí)如何將理論知識(shí)和實(shí)際工程問題結(jié)合起來。二、實(shí)驗(yàn)器材 計(jì)算機(jī)、軟件（已安裝，包括各類編程軟件C語言、C+、VB、VC等、應(yīng)用軟件MATLAB等）、移動(dòng)存儲(chǔ)設(shè)備（學(xué)生自備，軟盤

33、、U盤等）三、實(shí)驗(yàn)內(nèi)容 編制調(diào)試水、火電廠間有功功率負(fù)荷經(jīng)濟(jì)分配的計(jì)算機(jī)程序。程序要求根據(jù)已知參數(shù)，完成水火電廠間有功功率負(fù)荷的經(jīng)濟(jì)分配計(jì)算，要求計(jì)算出各電廠滿足經(jīng)濟(jì)運(yùn)行要求所分配的有功負(fù)荷等參數(shù)。1.熟悉電力系統(tǒng)經(jīng)濟(jì)運(yùn)行的計(jì)算方法，按照計(jì)算方法編制程序。2.將事先編制好的計(jì)算程序原代碼由自備移動(dòng)存儲(chǔ)設(shè)備導(dǎo)入計(jì)算機(jī)。3.在相應(yīng)的編程環(huán)境下對程序進(jìn)行組織調(diào)試。4.應(yīng)用計(jì)算例題驗(yàn)證程序的計(jì)算效果。5.對調(diào)試正確的計(jì)算程序進(jìn)行存儲(chǔ)、打印。6.完成本次實(shí)驗(yàn)的實(shí)驗(yàn)報(bào)告。4、 實(shí)驗(yàn)數(shù)據(jù) 一個(gè)火電廠和一個(gè)水電廠并聯(lián)運(yùn)行?；痣姀S的燃料消耗特性為： 水電廠的耗水量特性為： 水電廠的給定日用水量為。系統(tǒng)的日負(fù)荷變

34、化如下： 時(shí) 負(fù)荷為350MW，時(shí)，負(fù)荷為700MW，時(shí) 負(fù)荷為 500MW。火電廠容量為600MW，水電廠容量為450MW。試確定水、火電廠間的功率經(jīng)濟(jì)分配。5、 實(shí)驗(yàn)程序%程序說明：ph變量是多行三列的變量，每一列分別對應(yīng)書中PH1 PH2 PH3；行列數(shù)等于迭代次數(shù)%pt變量是多行三列的變量，每一列分別對應(yīng)書中PT1 PT2 PT3；行列數(shù)等于迭代次數(shù)%W變量是一行多列向量，代表書中的全日發(fā)電耗水量，列數(shù)等于迭代次數(shù)clc;clear;syms sym_pt sym_ph sym_f sym_w sym_pld sym_y %定義變量sym_f=3+0.4.*sym_pt+0.00035.*sym_pt2; %火電消耗方程sym_w=2+0.8.*sym_ph+1.5.*0.001.*sym_ph2;%水電消耗方程sym_df=char(diff(sym_f,sym_pt);%火電消耗方程求導(dǎo)sym_dw=char(diff(sym_w,sym_ph);%水電消耗方程求導(dǎo)sym_mw=double(sym_d