【正文】 pQ(2*k)*q。p=p+g*eb*f。e=Q(2*j1,1)。for j=1:1:A(1,1)g=real(Y(k,j))。m=0。p=0。W(2*k1,x)=C(k,3)Q(2*k)*m+Q(2*k1,1)*n。n=n+g*f+b*e。q=q+g*f+b*e。f=Q(2*j,1)。b=imag(Y(k,j))。n=0。q=0。for x=1:1:50%設置迭代次數(shù)為50次for i=1:1:A(7,1)k=C(i,1)。)。fprintf(fid,39。====================節(jié)點電壓V=============================n39。)。,39。endfid=fopen(39。endQ(2*A(3,1)1,1)=A(4,1)。elseQ(2*i1,1)=1。for i=1:1:(A(1,1)1)if C(i,3)0Q(2*i1,1)=C(i,3)。Y(n,m)=Y(m,n)。Y(n,n)=Y(n,n)+1/(B(a,3)+j*B(a,4))j*B(a,5)。Y(n,m)=Y(m,n)。Y(n,n)=Y(n,n)+1/(B(a,3)+j*B(a,4))/(B(a,5)^2)。n=B(a,2)。Y=zeros(A(1,1))。C=C39。fclose(fid)。%f39。,[5,A(2,1)])。%讀8個數(shù) B=fscanf(fid,39。%f39。)。,39。三、程序框圖四、實驗過程實驗程序及說明 clear fid=fopen(39。通過自己設計電力系統(tǒng)計算程序加深對電力系統(tǒng)分析的理解，同時加強計算機實際應用能力。39。endendend end %%%============fclose(fid)。+%39。,im)。if imfprintf(fid,39。%39。,i,j,i,j)。if(re~=0||im~=0)fprintf(fid,39。re=real(k)。)。end %%%==========%%%========== %輸出打印支路電流 fprintf(fid,39。n39。,im)。elseif im=0fprintf(fid,39。%39。im=imag(k)。%39。re=real(k)。,i,i)。%%%==========%%%========== %打印輸出節(jié)點電壓 for i=1:1:5fprintf(fid,39。nVn39。,im)。elseif im=0fprintf(fid,39。%39。im=imag(If)。%39。re=real(If)。If=39。)。endendfprintf(fid,39。+%39。,im)。if imfprintf(fid,39。,re)。fprintf(fid,39。for i=1:1:5for j=1:1:5k=Z(i,j)。Z_matrixn39。)。endendfprintf(fid,39。+%39。,im)。if imfprintf(fid,39。,re)。fprintf(fid,39。%%%========= %輸出導納陣 for i=1:1:5for j=1:1:5k=Y(i,j)。Y_matrixn39。)。,39。end end %%%=========fid=fopen(39。if B(a,5)0I(m,n)=(B(a,5)*V(1,m)V(1,n))/(B(a,3)+j*B(a,4))。end %%%========== %計算各支路電流 for a=1:1:A(2,1)m=B(a,1)。If=1/Z(3,3)。end end %%%%%========= Y(2,2)=Y(2,2)+j*。Y(m,n)=1/(B(a,3)+j*B(a,4))。elseif B(a,5)Y(m,m)=Y(m,m)+1/(B(a,3)+j*B(a,4))j*B(a,5)。Y(m,n)=1/(B(a,3)+j*B(a,4))/B(a,5)。if B(a,5)0Y(m,m)=Y(m,m)+1/(B(a,3)+j*B(a,4))。% 形成55的0陣%%%%%====== %導納陣元素計算 for a=1:1:A(2,1)m=B(a,1)。B=B39。,[5,A(2,1)])。B=fscanf(fid,39。%d39。)。,39。二、實驗內(nèi)容編寫數(shù)據(jù)輸入、形成導納陣程序電力系統(tǒng)短路計算實用公式的計算程序及編制和調(diào)試三、實驗過程程序代碼及說明%清屏 clear %讀數(shù)據(jù)fid=fopen(39。第一篇：哈工大電力系統(tǒng)分析短路及潮流計算實驗上機程序上機實驗實驗一節(jié)點導納陣生成和短路電流計算 實驗二簡單系統(tǒng)的牛頓法潮流計算姓名：班級：141班 學號：實驗說明：本程序的電路結構來自翁增銀、何仰贊主編的《電力系統(tǒng)分析》的例題實驗一節(jié)點導納陣生成和短路電流計算一、實驗目的根據(jù)所給的電力系統(tǒng)，編制短路電流計算程序，通過計算機進行調(diào)試，最后完成一個切實可行的電力系統(tǒng)計算應用程序。通過自己設計電力系統(tǒng)計算程序加深對電力系統(tǒng)分析的理解，同時加強計算機實際應用能力。39。r39。A=fscanf(fid,39。,2)。%f39。fclose(fid)。%求節(jié)點導納矩陣Y=zeros(A(1,1))。n=B(a,2)。Y(n,n)=Y(n,n)+1/(B(a,3)+j*B(a,4))/(B(a,5)^2)。Y(n,m)=Y(m,n)。Y(n,n)=Y(n,n)+1/(B(a,3)+j*B(a,4))j*B(a,5)。Y(n,m)=Y(m,n)。Z=inv(Y)。%接金屬短路，Zf=0 for i=1:1:5V(1,i)=1Z(i,3)/Z(3,3)。n=B(a,2)。elseif B(a,5)I(m,n)=(V(1,m)V(1,n))/(B(a,3)+j*B(a,4))。39。w39。fprintf(fid,39。)。re=real(k)。%39。im=imag(k)。%39。elseif im=0fprintf(fid,39。,im)。n39。end %%%============%%%============ %輸出阻抗陣，導納陣求逆 fprintf(fid,39。)。re=real(k)。%39。im=imag(k)。%39。elseif im=0fprintf(fid,39。,im)。n39。end %%%===========%%%=========== %打印輸出短路電流fprintf(fid,39。)。fprintf(fid,39。,re)。if imfprintf(fid,39。,im)。+%39。end fprintf(fid,39。)。Note%d:V%d=39。k=V(1,i)。fprintf(fid,39。,re)。if imfprintf(fid,39。,im)。+%39。endfprintf(fid,39。)。Ibrn39。for i=1:1:5for j=1:1:4k=I(i,j)。im=imag(k)。Branch%d%d:I%d%d=39。fprintf(fid,39。,re)。%39。elseif im=0fprintf(fid,39。,im)。%關閉文件附：39。程序輸出結果 節(jié)點導納陣：節(jié)點阻抗陣：短路電流：If= + 節(jié)點電壓：各支路電流：四、實驗總結這是我的第一次上機實驗，感覺稍微有點難，主要還是在工具軟件—C語言或者MATLAB的運用上，但是我相信，以后學習中，我會努力掌握的，這是我把理論應用于實際中的必要橋梁！實驗二簡單系統(tǒng)的牛頓法潮流計算一、實驗目的根據(jù)所給的電力系統(tǒng)，編制牛頓法潮流計算程序，通過計算機進行調(diào)試，最后完成一個切實可行的電力系統(tǒng)計算應用程序。二、實驗內(nèi)容電力系統(tǒng)潮流計算的計算程序設計及編制和調(diào)試。39。r39。%打開輸入數(shù)據(jù) A=fscanf(fid,39。,8)。%f39。C=fscanf(fid,39。,[3,(A(1,1)1)])。B=B39。B(2,5)=1/B(2,5)。%得44的0陣%%%================= %與前一實驗同法求導納陣 for a=1:1:A(2,1)m=B(a,1)。if B(a,5)0Y(m,m)=Y(m,m)+1/(B(a,3)+j*B(a,4))。Y(m,n)=1/(B(a,3)+j*B(a,4))/B(a,5)。else if B(a,5)Y(m,m)=Y(m,m)+1/(B(a,3)+j*B(a,4))j*B(a,5)。Y(m,n)=1/(B(a,3)+j*B(a,4))。endend end %%%===================%%%=================== %節(jié)點電壓賦初值Q=zeros(2*A(1,1),1)。Q(2*i,1)=0。Q(2*i,1)=0。Q(2*A(3,1),1)=0。39。w39。fprintf(fid,39。)。迭代計數(shù)tt V1=e1+jf1tttt V2=e2+jf2ttttV3=e3+jf3ttn39。%%%=======================求W陣=================================== W=zeros(2*(A(1,1)1),50)。p=0。m=0。for j=1:1:A(1,1)g=real(Y(k,j))。e=Q(2*j1,1)。p=p+g*eb*f。m=m+g*eb*f。endW(2*k,x)=C(k,2)Q(2*k1,1)*pQ(2*k)*q。End %%%=========================PQ節(jié)點=================================%%%=========================PV節(jié)點================================= for l=1:1:A(8,1)k=C(l+A(7,1),1)。q=0。n=0。b=imag(Y(k,j))。f=Q(2*j,1)。q=q+g*f+b*e。W(2*k1,x)=(C(k,3))^2(Q(2*k1,1)^2+Q(2*k)^2)。for i=1:1:2*(A(1,1)1)Max=max(abs(W(i,x)),Max)。end %%%===================%%%=================== %求雅克比行列式 for i=1:1:A(1,1)1%===========================================for j=1:1:A(1,1)1%===========================i~=jif i~=jk=C(i,1)。b=imag(Y(k,j))。f=Q(2*k,1)。J(2*k1,2*j)=b*eg*f。J(2*k,2*j1)=J(2*k1,2*j)。J(2*k,2*j)=0。p=0。m=0。g=real(Y(k,k))。J(2*k1,2*j1)=g*Q(2*i1,1)b*Q(2*i,1)。if k~=3J(2*k,2*j1)=b*Q(2*i1,1)g*Q(2*i,1)。endfor r=1:1:A(1,1)g=real(Y(k,r))。e=Q(2*r1,1)。p=p+g*eb*f。endJ(2*k1,2*j1)=p+J(2*k1,2*j1)。if i~=3J(2*k,2*j1)=q+J(2*k,2*j1)。elseJ(2*k,2*j1)=2*Q(2*i1,1)。endendend end %%%=========