最新南昌大學電力系統(tǒng)分析實驗報告(共17頁)

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

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

3、網(wǎng)系統(tǒng)，S1=S2=2+j2(MV·A)，S3=S4=S5=4+j4(MV·A)，假定所有負荷均為恒功率負荷，節(jié)點1為參考節(jié)點相角為0o。計算電網(wǎng)個節(jié)點電壓及支路功率分布，收斂條件為。5、 實驗程序% 配電網(wǎng)前推回代潮流計算程序 % 使用IEEE 33節(jié)點配電系統(tǒng)作為算例，可實現(xiàn)弱環(huán)網(wǎng)情況下的潮流計算 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;

4、7 8 0.7114 0.2351; 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

5、 0.9373; 3 23 0.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.

6、00; 5 60.00 30.00; 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

7、 50.00; 24 420.00 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é)點，負荷 UB = 12.66; % 電壓基準 kVSB = 10; % 功率基準 MVAZB = UB2/SB; % 阻抗基準 ohm BranchData(:,3,4) = BranchData(:,3,4) / ZB; % 阻抗標幺

8、化NodeData(:,2,3) = NodeData(:,2,3) / SB / 1000;% 功率標幺化 NN = 33; % 節(jié)點數(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

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

10、/ V); % 節(jié)點注入電流 for n = 1: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

11、st=BranchData(i,1); en=BranchData(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實驗五 輸電網(wǎng)潮流計算實驗一、實驗?zāi)康?本實驗針對高壓輸電網(wǎng)絡(luò)具有多環(huán)、多電源的結(jié)構(gòu)特性，要求根據(jù)功率平衡方程的泰勒展開式近似的求取系統(tǒng)的狀態(tài)變量，使系統(tǒng)潮流計算能夠由計算機自行完成，即根據(jù)已知的電力網(wǎng)的數(shù)學模型（節(jié)點導納矩陣）及各節(jié)點參數(shù)，由計算程序運行完成該電力系統(tǒng)的潮流計算。通過實驗教學加深學生對復雜電力系統(tǒng)潮流計算計算方法的理解，學會運用電力

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

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

14、');nl=4;%input('請輸入支路數(shù)：nl=');isb=1;%input('請輸入平衡母線節(jié)點號：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é)點參數(shù)形成的矩陣：B2=');B2=0 0 1.05 1.05 0 1; 0 0.

15、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(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é)點處于1側(cè) p=B1(i,1);q=B1(i,2); else %左節(jié)點處于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

16、(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,3)+B1(i,4)./2; %對角元1側(cè) end%求導納矩陣disp('導納矩陣 Y=');disp(Y)%-G=real(Y);B=imag(Y);%分解出導納陣的實部和虛部 for i=1:n%給定各節(jié)點初始電壓的實部和虛部 e(i)=real(B2(i,3); f(i)=imag(B2(i,3); V(i)=B2(i,4);%PV節(jié)點電壓給定模值 endfor i=1:n%給定各節(jié)點注入功率 S(i)=B2(i,1)-B2

17、(i,2); %i節(jié)點注入功率SG-SL B(i,i)=B(i,i)+B2(i,5);%i節(jié)點無功補償量 end%-P=real(S);Q=imag(S); %分解出各節(jié)點注入的有功和無功功率ICT1=0;IT2=1;N0=2*n;N=N0+1;a=0; %迭代次數(shù)ICT1、a；不滿足收斂要求的節(jié)點數(shù)IT2while IT2=0 % N0=2*n 雅可比矩陣的階數(shù)；N=N0+1擴展列 IT2=0;a=a+1; for i=1:n if i=isb%非平衡節(jié)點 C(i)=0;D(i)=0; for j1=1:n C(i)=C(i)+G(i,j1)*e(j1)-B(i,j1)*f(j1);%(Gi

18、j*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)+f(i)*D(i);%節(jié)點功率P計算ei(Gij*ej-Bij*fj)+fi(Gij*fj+Bij*ej) Q1=C(i)*f(i)-e(i)*D(i);%節(jié)點功率Q計算fi(Gij*ej-Bij*fj)-ei(Gij*fj+Bij*ej)%求i節(jié)點有功和無功功率P',Q'的計算值 V2=e(i)2+f(i)2;%電壓模平方 %- 以下針對非PV節(jié)點來求取功率差及Jacobi矩陣元素 - if B2(i,6)=

19、3%非PV節(jié)點 DP=P(i)-P1;%節(jié)點有功功率差 DQ=Q(i)-Q1; %節(jié)點無功功率差 %- 以上為除平衡節(jié)點外其它節(jié)點的功率計算 -%-求取Jacobi矩陣 - for j1=1:n if j1=isb&j1=i%非平衡節(jié)點&非對角元 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)=

20、X3;J(p,N)=DQ;m=p+1; % X3=dQ/de J(p,N)=DQ節(jié)點無功功率差 J(m,q)=X1;J(m,N)=DP;q=q+1; % X1=dP/de J(m,N)=DP節(jié)點有功功率差 J(p,q)=X4;J(m,q)=X2; % X4=dQ/df X2=dp/df elseif j1=i&j1=isb%非平衡節(jié)點&對角元 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

21、/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;%擴展列Q m=p+1; J(m,q)=X1;q=q+1;J(p,q)=X4;J(m,N)=DP;%擴展列P J(m,q)=X2; end end else%-下面是針對PV節(jié)點來求取Jacobi矩陣的元素 - DP=P(i)-P1;% PV節(jié)點有功誤差 DV=V(i)2-V2;% PV節(jié)點電壓誤差 for j1=1:n if j1=isb&j1=i%非平衡節(jié)點&非對角元 X1=-G(i,j1)*e(i)-B(i,j

22、1)*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é)點電壓誤差 m=p+1; J(m,q)=X1;J(m,N)=DP;q=q+1;J(p,q)=X6; % PV節(jié)點有功誤差 J(m,q)=X2; elseif j1=i&j1=isb %非平衡節(jié)點&對角元 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

23、 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é)點電壓誤差 m=p+1; J(m,q)=X1;J(m,N)=DP;q=q+1;J(p,q)=X6; % PV節(jié)點有功誤差 J(m,q)=X2; end end end end end%- 以上為求雅可比矩陣的各個元素及擴展列的功率差或電壓差 - for k=3:N0 % N0=2*n （從第三行開始，第一、二行是平衡節(jié)點） k1=k+1;N1=N; % N=N0+1 即 N=2*n+1擴展列P、Q 或 U for k2=k1:N1% 從k+1列的Jacobi元

24、素到擴展列的P、Q 或 U J(k,k2)=J(k,k2)./J(k,k);% 用K行K列對角元素去除K行K列后的非對角元素進行規(guī)格化 end J(k,k)=1; % 對角元規(guī)格化K行K列對角元素賦1 %- 回代運算 - if k=3 % 不是第三行 k > 3 k4=k-1; for k3=3:k4% 用k3行從第三行開始到當前行的前一行k4行消去 for k2=k1:N1% k3行后各行上三角元素 J(k3,k2)=J(k3,k2)-J(k3,k)*J(k,k2);%消去運算（當前行k列元素消為0） end %用當前行K2列元素減去當前行k列元素乘以第k行K2列元素 J(k3,k)=

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

26、N1 %從第四列到2n+1列（即擴展列） J(k3,k2)=J(k3,k2)-J(k3,k)*J(k,k2);%消去運算（當前行3列元素消為0） end %用當前行K2列元素減去當前行3列元素乘以第三行K2列元素 J(k3,k)=0; %當前行第3列元素已消為0 end end end%-上面是用線性變換方式高斯消去法將Jacobi矩陣化成單位矩陣- for k=3:2:N0-1 L=(k+1)./2; e(L)=e(L)-J(k,N); %修改節(jié)點電壓實部 k1=k+1; f(L)=f(L)-J(k1,N); %修改節(jié)點電壓虛部 end %-修改節(jié)點電壓- for k=3:N0 DET=ab

27、s(J(k,N); if DET>=pr %電壓偏差量是否滿足要求 IT2=IT2+1; %不滿足要求的節(jié)點數(shù)加1 end end ICT2(a)=IT2; %不滿足要求的節(jié)點數(shù) ICT1=ICT1+1; %迭代次數(shù)end%用高斯消去法解"w=-J*V"disp('迭代次數(shù)：');disp(ICT1);disp('沒有達到精度要求的個數(shù)：');disp(ICT2);for k=1:n V(k)=sqrt(e(k)2+f(k)2); %計算各節(jié)點電壓的模值 sida(k)=atan(f(k)./e(k)*180./pi; %計算各節(jié)點電壓

28、的角度 E(k)=e(k)+f(k)*j; %將各節(jié)點電壓用復數(shù)表示end%- 計算各輸出量 -disp('各節(jié)點的實際電壓標幺值E為(節(jié)點號從小到大排列)：');disp(E); %顯示各節(jié)點的實際電壓標幺值E用復數(shù)表示disp('-');disp('各節(jié)點的電壓大小V為(節(jié)點號從小到大排列)：');disp(V); %顯示各節(jié)點的電壓大小V的模值disp('-');disp('各節(jié)點的電壓相角sida為(節(jié)點號從小到大排列)：');disp(sida); %顯示各節(jié)點的電壓相角for p=1:n C(p)=0;

29、 for q=1:n C(p)=C(p)+conj(Y(p,q)*conj(E(q); %計算各節(jié)點的注入電流的共軛值 end S(p)=E(p)*C(p); %計算各節(jié)點的功率 S = 電壓 X 注入電流的共軛值enddisp('各節(jié)點的功率S為(節(jié)點號從小到大排列)：');disp(S); %顯示各節(jié)點的注入功率 五、實驗結(jié)果 在上圖所示的簡單電力系統(tǒng)中，系統(tǒng)中節(jié)點1、2為PQ節(jié)點，節(jié)點3為PV節(jié)點，節(jié)點4為平衡節(jié)點，已給定P1s+jQ1s=-0.30-j0.18 P2s+jQ2s=-0.55-j0.13 P3s=0.5 V3s=1.10 V4s=1.050°容許

30、誤差=10-5節(jié)點導納矩陣：各節(jié)點電壓： 節(jié)點 e f v 1. 0.984637 -0.008596 0.984675 -0.500172 2. 0.958690 -0.108387 0.964798 -6.450306 3. 1.092415 0.128955 1.100000 6.732347 4. 1.050000 0.000000 1.050000 0.000000各節(jié)點功率： 節(jié)點 P Q 1 -0.300000 -0.180000 2 0.550000 -0.130000 3 0.500000 -0.551305 4 0.367883 0.264698導納矩陣 Y=1.0421

31、- 8.2429i-0.5882 + 2.3529i0.0000 + 3.6667i-0.4539 + 1.8911i-0.5882 + 2.3529i1.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é)點的實際電壓標幺值E為(節(jié)點號從小到大排列)： 0.9854 -

32、 0.0485i 0.9691 - 0.0759i 1.0968 + 0.0844i 1.0500 + 0.0000i各節(jié)點的電壓大小e為(節(jié)點號從小到大排列)： 0.9854 0.9691 1.0968 1.0500各節(jié)點的電壓大小f為(節(jié)點號從小到大排列)： -0.0485 -0.0759 0.0844 0各節(jié)點的電壓大小V為(節(jié)點號從小到大排列)： 0.9866 0.9721 1.1000 1.0500各節(jié)點的電壓相角sida為(節(jié)點號從小到大排列)： -2.8167 -4.4784 4.4018 0各節(jié)點的功率S為(節(jié)點號從小到大排列)： -0.5500 - 0.1300i -0.30

33、00 - 0.1800i 0.5000 + 0.0857i 0.3595 + 0.2343i 實驗六 水火電廠間有功功率負荷的經(jīng)濟分配一、實驗?zāi)康?本實驗通過對水火電廠間有功功率負荷的經(jīng)濟分配的計算機程序編制與調(diào)試，獲得能計算有功負荷經(jīng)濟分配的計算機程序，即根據(jù)已知的水火電廠的耗量特性、發(fā)電用水量及各時間段的有功負荷功率，由計算機程序自動完成各時段各電廠的有功出力的計算。通過實驗教學加深學生對電力系統(tǒng)經(jīng)濟運行概念的理解，學會運用數(shù)學模型進行水火電廠間的有功功率負荷的經(jīng)濟分配計算，熟悉各種常用應(yīng)用軟件，熟悉硬件設(shè)備的使用方法，加強編制調(diào)試計算機程序的能力，提高工程計算的能力，學習如何將理論知識和

34、實際工程問題結(jié)合起來。二、實驗器材 計算機、軟件（已安裝，包括各類編程軟件C語言、C+、VB、VC等、應(yīng)用軟件MATLAB等）、移動存儲設(shè)備（學生自備，軟盤、U盤等）三、實驗內(nèi)容 編制調(diào)試水、火電廠間有功功率負荷經(jīng)濟分配的計算機程序。程序要求根據(jù)已知參數(shù)，完成水火電廠間有功功率負荷的經(jīng)濟分配計算，要求計算出各電廠滿足經(jīng)濟運行要求所分配的有功負荷等參數(shù)。1.熟悉電力系統(tǒng)經(jīng)濟運行的計算方法，按照計算方法編制程序。2.將事先編制好的計算程序原代碼由自備移動存儲設(shè)備導入計算機。3.在相應(yīng)的編程環(huán)境下對程序進行組織調(diào)試。4.應(yīng)用計算例題驗證程序的計算效果。5.對調(diào)試正確的計算程序進行存儲、打印。6.完成

35、本次實驗的實驗報告。4、 實驗數(shù)據(jù) 一個火電廠和一個水電廠并聯(lián)運行?；痣姀S的燃料消耗特性為： 水電廠的耗水量特性為： 水電廠的給定日用水量為。系統(tǒng)的日負荷變化如下： 時 負荷為350MW，時，負荷為700MW，時 負荷為 500MW。火電廠容量為600MW，水電廠容量為450MW。試確定水、火電廠間的功率經(jīng)濟分配。5、 實驗程序%程序說明：ph變量是多行三列的變量，每一列分別對應(yīng)書中PH1 PH2 PH3；行列數(shù)等于迭代次數(shù)%pt變量是多行三列的變量，每一列分別對應(yīng)書中PT1 PT2 PT3；行列數(shù)等于迭代次數(shù)%W變量是一行多列向量，代表書中的全日發(fā)電耗水量，列數(shù)等于迭代次數(shù)clc;clear

36、;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');%火電消耗方程求導sym_dw=char(diff(sym_w,'sym_ph');%水電消耗方程求導sym_mw=double(sym_df) double('=sym_y*(') double(s

37、ym_dw) double(')'); %建立等式方程sym_ph,sym_pt=solve('sym_ph+sym_pt-sym_pld=0',char(sym_mw),'sym_ph','sym_pt');%求解方程組format long %用水量誤差值設(shè)定e=50;%實際全日發(fā)電耗水量wd=1.5e7;%y值的初值y(1)=0.5;%負荷功率pld=350 700 500;%實際負荷功率初值w(1)=0;%迭代次數(shù)k=1;%y每次迭代增量add_num=1e-4;ri=8 10 6;while (abs(w(k)-wd)>e) && k<=3000%公式ph(k,:)=(0.4-0.8.*y(k)+0.0007.*pld)./(0.003.*y(k)+0.0007);ph(k,:)=double(subs(sym_ph,sym_pld,sym_y,pld,y(k);% 公式pt(k,:)=(0.8.*y(k)-0.4+0.003.