【正文】 nj(B1(i,4)./2)+(conj(E(q)*B1(i,5))...conj(E(p)))*conj(1./(B1(i,3)*B1(i,5))))。,num2str(p),39。各條支路的功率損耗DS為(順序同您輸入B1時(shí)一致)：39。DS(39。39。%input(39。%input(39。+0。0+002。)。q=B1(i,2)。導(dǎo)納矩陣 Y=39。%i節(jié)點(diǎn)注入功率SGSLB(i,i)=B(i,i)+B2(i,5)。%迭代次數(shù)ICTa；不滿足收斂要求的節(jié)點(diǎn)數(shù)IT2 while IT2~=0% N0=2*n 雅可比矩陣的階數(shù)；N=N0+1擴(kuò)展列IT2=0。%節(jié)點(diǎn)功率Q計(jì)算fiΣ(Gij*ejBij*fj)eiΣ(Gij*fj+Bij*ej)%求i節(jié)點(diǎn)有功和無(wú)功功率P39。% dP/de=dQ/dfX2=B(i,j1)*e(i)G(i,j1)*f(i)。% X3=dQ/de J(p,N)=DQ節(jié)點(diǎn)無(wú)功功率差J(m,q)=X1。% dP/deX2=D(i)+B(i,i)*e(i)G(i,i)*f(i)。J(m,q)=X1。j1~=i%非平衡節(jié)點(diǎn)amp。J(p,N)=DV。j1~=isb %非平衡節(jié)點(diǎn)amp。J(p,N)=DV。N1=N。%消去運(yùn)算end%用當(dāng)前行K2列元素減去當(dāng)前行k列元素乘以第k行K2列元素J(k3,k)=0。if DET=pr%電壓偏差量是否滿足要求IT2=IT2+1。沒(méi)有達(dá)到精度要求的個(gè)數(shù)：39。)。disp(V)。%顯示各節(jié)點(diǎn)的電壓相角 for p=1:nC(p)=0。39。Siz(i)=Si(p,q)。,39。end disp(39。Sjy(i)=Sj(q,p)。)=39。)。,num2str(p),39。)。請(qǐng)輸入支路數(shù)：nl=39。請(qǐng)輸入誤差精度：pr=39。+01。00002。00002]。sida=zeros(1,n)。%非對(duì)角元Y(q,q)=Y(q,q)+1./(B1(i,3)*B1(i,5)^2)+B1(i,4)./2。f(i)=imag(B2(i,3))。N0=2*n。%Σ(Gij*ejBij*fj)D(i)=D(i)+G(i,j1)*f(j1)+B(i,j1)*e(j1)。%節(jié)點(diǎn)無(wú)功功率差%=============== 以上為除平衡節(jié)點(diǎn)外其它節(jié)點(diǎn)的功率計(jì)算 ================= %================= 求取Jacobi矩陣 ===================for j1=1:nif j1~=isbamp。J(p,q)=X3。% X4=dQ/df X2=dp/dfelseif j1==iamp。J(p,q)=X3。endendelse%=============== 下面是針對(duì)PV節(jié)點(diǎn)來(lái)求取Jacobi矩陣的元素 ===========DP=P(i)P1。p=2*i1。J(p,q)=X6。p=2*i1。J(p,q)=X6。%消去運(yùn)算（當(dāng)前行k列元素消為0）end%用當(dāng)前行K2列元素減去當(dāng)前行k列元素乘以第k行K2列元素J(k3,k)=0。%修改節(jié)點(diǎn)電壓實(shí)部k1=k+1。)。%計(jì)算各節(jié)點(diǎn)電壓的角度E(k)=e(k)+f(k)*j。disp(39。各節(jié)點(diǎn)的電壓相角sida為(節(jié)點(diǎn)號(hào)從小到大排列)：39。)。for i=1:nlp=B1(i,1)。ZF=[39。disp(39。if B1(i,6)==0Sj(q,p)=E(q)*(conj(E(q))*conj(B1(i,4)./2)+(conj(E(q)./B1(i,5))...conj(E(p)))*conj(1./(B1(i,3)*B1(i,5))))。,num2str(q),39。)。DDS(i)=DS(i)。disp(ZF)。節(jié)點(diǎn)號(hào)39。xlabel(39。subplot(2,2,1)。節(jié)點(diǎn)注入有功39。ylabel(39。xlabel(39。bar(Q1)?？偨Y(jié)通過(guò)本次課程設(shè)計(jì)讓我有復(fù)習(xí)了一次潮流計(jì)算的相關(guān)知識(shí)，跟家清晰了什么事潮流計(jì)算以及潮流計(jì)算的在電力系統(tǒng)的重要性。而且對(duì)細(xì)節(jié)的把握也得到了鍛煉，做題的精細(xì)程度得到了提高。%X為節(jié)點(diǎn)號(hào)和對(duì)地參數(shù)矩陣。isb=input(39。)。Times=1。)。B1=input(39。請(qǐng)輸入支路數(shù):n1=39。%第五烈為含變壓器支路的變壓器的變比%第六列為變壓器是否是否含有變壓器的參數(shù)，其中“1”為含有變壓器，%“0”為不含有變壓器。不過(guò)在計(jì)算潮流計(jì)算的過(guò)程中卻對(duì)以往學(xué)過(guò)的電力系統(tǒng)分析的相關(guān)知識(shí)進(jìn)行了一次較為深入的復(fù)習(xí)。)。grid on。Q1=imag(Siz)。節(jié)點(diǎn)號(hào)39。)。grid on。subplot(1,2,2)。plot(V)。)=39。DS(i)=Si(p,q)+Sj(q,p)。disp(39。ZF=[39。for i=1:nlp=B1(i,1)。,num2str(SSi(p,q))]。enddisp(Si(p,q))。各條支路的首端功率Si為(順序同您輸入B1時(shí)一致)：39。%計(jì)算各節(jié)點(diǎn)的功率 S = 電壓 X 注入電流的共軛值 end disp(39。)。39。for k=1:nV(k)=sqrt(e(k)^2+f(k)^2)。%迭代次數(shù) end %用高斯消去法解”w=J*V disp(39。%當(dāng)前行第3列元素已消為0endendend %====上面是用線性變換方式高斯消去法將Jacobi矩陣化成單位矩陣=====for k=3:2:N01L=(k+1)./2。% 對(duì)角元規(guī)格化K行K列對(duì)角元素賦1%==================== 回代運(yùn)算=======================================if k~=3% 不是第三行k 3k4=k1。J(m,N)=DP。% dP/dfX5=2*e(i)。J(m,N)=DP。% dP/dfX5=0。J(m,N)=DP。% dQ/dfp=2*i1。% X1=dP/de J(m,N)=DP節(jié)點(diǎn)有功功率差J(p,q)=X4。% X1=dP/de X4=dQ/dfp=2*i1。%電壓模平方%========= 以下針對(duì)非PV節(jié)點(diǎn)來(lái)求取功率差及Jacobi矩陣元素 =========if B2(i,6)~=3%非PV節(jié)點(diǎn)DP=P(i)P1。D(i)=0。%分解出各節(jié)點(diǎn)注入的有功和無(wú)功功率 ICT1=0。B=imag(Y)。endY(p,q)=Y(p,q)1./(B1(i,3)*B1(i,5))。f=zeros(1,n)。0+002。00002。+0。pr=。nl=12。disp(39。ZF=[39。end disp(39。,39。Sjy(i)=Sj(q,p)。39。S(39。q=B1(i,2)。disp(S)。)。各節(jié)點(diǎn)的電壓大小V為(節(jié)點(diǎn)號(hào)從小到大排列)：39。%將各節(jié)點(diǎn)電壓用復(fù)數(shù)表示 end %=============== 計(jì)算各輸出量 =========================== disp(39。disp(ICT1)。f(L)=f(L)J(k1,N)。%當(dāng)前行第k列元素已消為0endif k==N0%若已到最后一行break。% PV節(jié)點(diǎn)有功誤差J(m,q)=X2。q=2*j11。% PV節(jié)點(diǎn)有功誤差J(m,q)=X2。q=2*j11。% PV節(jié)點(diǎn)有功誤差DV=V(i)^2V2。J(p,N)=DQ。j1~=isb %非平衡節(jié)點(diǎn)amp。J(p,N)=DQ。j1~=i%非平衡節(jié)點(diǎn)amp。%Σ(Gij*fj+Bij*ej)endP1=C(i)*e(i)+f(i)*D(i)。N=N0+1。V(i)=B2(i,4)。%對(duì)角元K側(cè)Y(p,p)=Y(p,p)+1./B1(i,3)+B1(i,4)./2。S1=zeros(nl)。% B1矩陣：支路首端號(hào)；末端號(hào)；支路阻抗；支路對(duì)地電納 %支路的變比；支路首端處于K側(cè)為1，1側(cè)為0 % B2矩陣：該節(jié)點(diǎn)發(fā)電機(jī)功率；該節(jié)點(diǎn)負(fù)荷功率；節(jié)點(diǎn)電壓初始值 %PV節(jié)點(diǎn)電壓V的給定值；節(jié)點(diǎn)所接的無(wú)功補(bǔ)償設(shè)備的容量 %節(jié)點(diǎn)分類(lèi)標(biāo)號(hào)：1為平衡節(jié)點(diǎn)（應(yīng)為1號(hào)節(jié)點(diǎn)）；2為PQ節(jié)點(diǎn)； %3為PV節(jié)點(diǎn)；%input(39。0+002。+01。)。)。disp(ZF)。DDS(i)=DS(i)。)。,num2str(q),39。if B1(i,6)==0Sj(q,p)=E(q)*(conj(E(q))*conj(B1(i,4)./2)+(conj(E(q)./B1(i,5))...conj(E(p)))*conj(1./(B1(i,3)*B1(i,5))))。disp(39。ZF=[39。for i=1:nlp=B1(i,1)。)。各節(jié)點(diǎn)的電壓相角sida為(節(jié)點(diǎn)號(hào)從小到大排列)：39。disp(39。%計(jì)算各節(jié)點(diǎn)電壓的角度E(k)=e(k)+f(k)*j。)。%修改節(jié)點(diǎn)電壓實(shí)部k1=k+1。%消去運(yùn)算（當(dāng)前行k列元素消為0）end%用當(dāng)前行K2列元素減去當(dāng)前行k列元素乘以第k行K14 列元素J(k3,k)=0。J(p,q)=X6。p=2*i1。J(p,q)=X6。p=2*i1。endendelse%=============== 下面是針對(duì)PV節(jié)點(diǎn)來(lái)求取Jacobi矩陣的元素 ===========DP=P(i)P1。J(p,q)=X3。% X4=dQ/df X2=dp/dfelseif j1==iamp。J(p,q)=X3。%節(jié)點(diǎn)無(wú)功功率差%=============== 以上為除平衡節(jié)點(diǎn)外其它節(jié)點(diǎn)的功率計(jì)算 ================= %================= 求取Jacobi矩陣 ===================for j1=1:nif j1~=isbamp。%Σ(Gij*ejBij*fj)D(i)=D(i)+G(i,j1)*f(j1)+B(i,j1)*e(j1)。N0=2*n。f(i)=imag(B2(i,3))。%非對(duì)角元Y(q,q)=Y(q,q)+1./(B1(i,3)*B1(i,5)^2)+B1(i,4)./2。sida=zeros(1,n)。00002]。0+002。+0。pr=。nl=11。B2為節(jié)點(diǎn)參數(shù)矩陣，其中第一列為節(jié)點(diǎn)注入發(fā)電功率參數(shù)；第二列為節(jié)點(diǎn)負(fù)荷功率參數(shù)；第三列為節(jié)點(diǎn)電壓參數(shù)；第六列為節(jié)點(diǎn)類(lèi)型參數(shù)，其中“1”為平衡節(jié)點(diǎn)，“2”為PQ節(jié)點(diǎn)，“3”為PV節(jié)點(diǎn)參數(shù)。+180。(+)=(+)MVA 電壓分布：Se1=SL3+DSL3=(+)+(+)=(+)180。(+)=()MVAST4=ZZ+Z*L3*L3*Se+T4(VbVN)Z+ZL3180。16DVL1====VbDVL1=：功率分布：SL2=ZZ+Z*T3*L2*Sd+T3(VbVN)Z+ZL2****180。(+)=(+)MVA35Sh=15+DST7+DS07+jDQL5=(+)+=180。20=180。x0180。l4=180。180。180。l2=180。l1=180。10=180。180。180。[Vs(13)%+Vs(23)%Vs(12)%]=%=21Vs1%180。SN)=(+)MVA100DS03=4180。10)XT31Vs%180。(Dp0+j變壓器T3：（四并聯(lián)）I0%180。10)21Vs%180。103DS01=Dp0+j變壓器T2：（雙并聯(lián)）I0%180。11023RT1=180。12=36MWQ3=36180。I%=%===63kW239。xd========“d在Simulink中按照電力系統(tǒng)原型選擇元件進(jìn)行建模。Dp0=44kW239。I%=%=,d112180。在加人測(cè)量模塊，并對(duì)各元件的參數(shù)進(jìn)行設(shè)置后，用measurement和sink中的儀器可以觀察各元件的電壓、電流、功率的大小。同時(shí)，為了實(shí)時(shí)監(jiān)控電力系統(tǒng)的運(yùn)行狀態(tài)，也需要進(jìn)行大量而快速的潮流計(jì)算。本位就是運(yùn)用潮流計(jì)算具體分析，并有MATLAB仿真。MATLAB simulation目 錄 任務(wù)提出與方案論證....................................................................................................................................2 2 總體設(shè)計(jì)........................................................................................................................................................3 .............................................................................................................................3 .............................................................................................................................3 .............................................................................................................................3 3 詳細(xì)設(shè)計(jì)...............