【正文】 有功功率差J(p,q)=X4。% dQ/dfp=2*i1。J(m,N)=DP。% dP/dfX5=0。J(m,N)=DP。% dP/dfX5=2*e(i)。J(m,N)=DP。% 對(duì)角元規(guī)格化K行K列對(duì)角元素賦1%==================== 回代運(yùn)算=======================================if k~=3% 不是第三行k 3k4=k1。%當(dāng)前行第3列元素已消為0endendend %====上面是用線(xiàn)性變換方式高斯消去法將Jacobi矩陣化成單位矩陣=====for k=3:2:N01L=(k+1)./2。%迭代次數(shù) end %用高斯消去法解”w=J*V“ disp(39。for k=1:nV(k)=sqrt(e(k)^2+f(k)^2)。39。)。%計(jì)算各節(jié)點(diǎn)的功率 S = 電壓 X 注入電流的共軛值 end disp(39。各條支路的首端功率Si為(順序同您輸入B1時(shí)一致)：39。enddisp(Si(p,q))。,num2str(SSi(p,q))]。for i=1:nlp=B1(i,1)。ZF=[39。disp(39。DS(i)=Si(p,q)+Sj(q,p)。)=39。%input(39。%input(39。+01。+01。0+002。+0。)。)。end%本程序的功能是用牛頓——拉夫遜法進(jìn)行12節(jié)點(diǎn)潮流計(jì)算 %本程序的功能是用牛頓——拉夫遜法進(jìn)行潮流計(jì)算 clear。,39。for i=1:nlp=B1(i,1)。,num2str(SSj(q,p))]。enddisp(Sj(q,p))。各條支路的末端功率Sj為(順序同您輸入B1時(shí)一致)：39。,num2str(q),39。elseSi(p,q)=E(p)*(conj(E(p))*conj(B1(i,4)./2)+(conj(E(p)./B1(i,5))...conj(E(q)))*conj(1./(B1(i,3)*B1(i,5))))。)。for q=1:nC(p)=C(p)+conj(Y(p,q))*conj(E(q))。%顯示各節(jié)點(diǎn)的電壓大小V的模值 disp(39。disp(E)。)。%不滿(mǎn)足要求的節(jié)點(diǎn)數(shù)加1endendICT2(a)=IT2。%當(dāng)前行第k列元素已消為0endelse%是第三行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)。% N=N0+1 即 N=2*n+1擴(kuò)展列△P、△Q 或 △Ufor k2=k1:N1% 從k+1列的Jacobi元素到擴(kuò)展列的△P、△Q 或 △UJ(k,k2)=J(k,k2)./J(k,k)。% PV節(jié)點(diǎn)電壓誤差m=p+1。對(duì)角元X1=C(i)G(i,i)*e(i)B(i,i)*f(i)。% PV節(jié)點(diǎn)電壓誤差m=p+1。非對(duì)角元X1=G(i,j1)*e(i)B(i,j1)*f(i)。q=q+1。% dP/dfX3=D(i)+B(i,i)*e(i)G(i,i)*f(i)。J(m,N)=DP。% dP/df=dQ/deX3=X2。,Q39。a=a+1。%i節(jié)點(diǎn)無(wú)功補(bǔ)償量end %===================== P=real(S)。)。else%左節(jié)點(diǎn)處于K側(cè)p=B1(i,2)。Y=zeros(n)。00002。+01。請(qǐng)輸入誤差精度：pr=39。請(qǐng)輸入支路數(shù)：nl=39。)。,num2str(p),39。)。)=39。Sjy(i)=Sj(q,p)。end disp(39。,39。Siz(i)=Si(p,q)。39。%顯示各節(jié)點(diǎn)的電壓相角 for p=1:nC(p)=0。disp(V)。)。沒(méi)有達(dá)到精度要求的個(gè)數(shù)：39。if DET=pr%電壓偏差量是否滿(mǎn)足要求IT2=IT2+1。%消去運(yùn)算end%用當(dāng)前行K2列元素減去當(dāng)前行k列元素乘以第k行K2列元素J(k3,k)=0。N1=N。J(p,N)=DV。j1~=isb %非平衡節(jié)點(diǎn)amp。J(p,N)=DV。j1~=i%非平衡節(jié)點(diǎn)amp。J(m,q)=X1。% dP/deX2=D(i)+B(i,i)*e(i)G(i,i)*f(i)。% X3=dQ/de J(p,N)=DQ節(jié)點(diǎn)無(wú)功功率差J(m,q)=X1。% dP/de=dQ/dfX2=B(i,j1)*e(i)G(i,j1)*f(i)。%節(jié)點(diǎn)功率Q計(jì)算fiΣ(Gij*ejBij*fj)eiΣ(Gij*fj+Bij*ej)%求i節(jié)點(diǎn)有功和無(wú)功功率P39。%迭代次數(shù)ICTa；不滿(mǎn)足收斂要求的節(jié)點(diǎn)數(shù)IT2 while IT2~=0% N0=2*n 雅可比矩陣的階數(shù)；N=N0+1擴(kuò)展列IT2=0。%i節(jié)點(diǎn)注入功率SGSLB(i,i)=B(i,i)+B2(i,5)。導(dǎo)納矩陣 Y=39。q=B1(i,2)。)。00002。+0。)。)。%本程序的功能是用牛頓——拉夫遜法進(jìn)行11節(jié)點(diǎn)潮流計(jì)算 clear。()==VgDVT62==180。ST6=ST4()180。+180。VN=()180。()=()+=180。120DST3DS03jQL2DSL2=(+)MVA3：DST4P2+Q262+272=(RT4+jXT4)=(+)=(+)MVAVN21102P2++=(RL3+jXL3)=(+j14)=(+)MVA22VN11039。20=4W 22BL5=2180。104S DQ12L4=2BL4VN=11RL5=180。10SDQ1L3=2B2L3VN==r0180。l3=2180。10180。10180。SN)=(+)MVA100再次是傳輸線(xiàn)參數(shù)計(jì)算，5條傳輸線(xiàn)的具體計(jì)算如下。180。10)=(+)MVA 100 RT7XT721Dps180。10=%180。180。180。11023RT3=180。180。11023=180。10)2Vs%180。tan()=31MVar發(fā)電機(jī)G5：P5=2180。發(fā)電機(jī)G1：P1=4180。建立系統(tǒng)模型，并設(shè)置好參數(shù)以后，就可以在Simulink環(huán)境下進(jìn)行仿真運(yùn)行。10MVA10kVGGG3180。15MW”236。Dp0==73kW237。是電力系統(tǒng)研究人員長(zhǎng)期研究的一個(gè)課題。MATLAB simulation目 錄 任務(wù)提出與方案論證....................................................................................................................................2 2 總體設(shè)計(jì)........................................................................................................................................................3 .............................................................................................................................3 .............................................................................................................................3 .............................................................................................................................3 3 詳細(xì)設(shè)計(jì)........................................................................................................................................................4 計(jì)算前提............................................................................................................................................4 .............................................................................................................................................7 4設(shè)計(jì)圖及源程序...........................................................................................................................................11 .......................................................................................................................................11 ...............................................................................................................................11 5 總結(jié).............................................................................................................................................................31 參考文獻(xiàn)..........................................................................................................................................................32 任務(wù)提出與方案論證潮流計(jì)算是在給定電力系統(tǒng)網(wǎng)絡(luò)結(jié)構(gòu)、參數(shù)和決定系統(tǒng)運(yùn)行狀態(tài)的邊界條件的情況下確定系統(tǒng)穩(wěn)態(tài)運(yùn)行狀態(tài)的一種基本方法,是電力系統(tǒng)規(guī)劃和運(yùn)營(yíng)中不可缺少的一個(gè)重要組成部分。fp1=fopen(”“,”w+“)。float ds[N][N]。else goto top。iguass(1,N1,y1,B,x1)。} getch()。g[j][j]+=lg。jHeadnode1。2 0。1 3 + 0 1 0。DetaU=inv(Jacbi)*DetaS。Jacbi(2*h,2*k1)=2*imag(B2(i,3))。Jacbi(2*h,2*k1)=Jacbi(2*h1,2*k)。for i=1:n if B2(i,6)==2 h=h+1。DetaS(2*h,1)=real(PVU(t,1))^2+imag(PVU(t,1))^2real(B2(i,3))^2imag(B2(i,3))^2。B2(i,6)==2 h=h+1。for j=1:n OrgS(2*h1,1)=OrgS(2*h1,1)+real(B2(i,3))*(real(Y(i,j))*real(B2(j,3))imag(Y(i,j))*imag(B2(j,3)))+imag(B2(i,3))*(real(Y(i,j))*imag(B2(j,3))+imag(Y(i,j))*real(B2(j,3)))。B2(i,3)=B2(i,3)+DetaU(2*j,1)+DetaU(2*j1,1)*sqrt(1)。Jacbi(2*h1,2*k)=real(Y(i,j))*real(B2(i,3))+imag(Y(i,j))*imag(B2(i,3))。end end end end end k=0。Jacbi(2*h1,2*k)=real(Y(i,j))*real(B2(i,3))+imag(Y(i,j))*imag(B2(i,3))+real(I(h,1))。for i=1:n if i~=isb h=h+1。end end t=0。t=0。B2(i,6)==2 h=h+1。Y(q,p)=Y(p,q)。q=B1(i,2)。請(qǐng)輸入節(jié)點(diǎn)參數(shù):B2=39。)。n=input(39。第一篇：基于MATLAB的電力系統(tǒng)潮流計(jì)算基于MATLAB的電力系統(tǒng)潮流計(jì)算%簡(jiǎn)單潮流計(jì)算的小程序，相關(guān)的原始數(shù)據(jù)數(shù)據(jù)數(shù)據(jù)輸入格式如下：%B1是支路參數(shù)矩陣，第一列和第二列是節(jié)點(diǎn)編號(hào)。其中第一列為節(jié)點(diǎn)編號(hào)，第二列為節(jié)點(diǎn)對(duì)地 %參數(shù)。請(qǐng)輸入平衡節(jié)點(diǎn)號(hào):isb=39。B2=input(39。%置迭代次數(shù)為初始值 %創(chuàng)建節(jié)點(diǎn)導(dǎo)納矩陣 for i=1:n1 if B1(i,6)==0 %不含變壓器的支路 p=B1(i,1)。Y(p,q)=Y(p,q)1/(B1(i,3)*B1(i,5))。for i=1:n %對(duì)PQ節(jié)點(diǎn)的處理 if i~=isbamp。end end end OrgS %創(chuàng)建PVU 用于存儲(chǔ)PV節(jié)點(diǎn)的初始電壓 PVU=zeros(nh1,1)。DetaS(2*h,1