【正文】 endtemp=line。i{int statemp,endtemp。coutt(bus[i].)。i{coutComplex::PrintfComplex(bus[i].V)。bus[i].V=Complex::Rec2Polar(bus[i].V)。i{bus[i]. = bus[i]. + bus[i].。}coutBus::JisuanNodeScal(X,bus,N)。dd=Complex::CaddC(aa,cc)。bb=Complex::subComplex(Complex::getconj(bus[endtemp1].V), Complex::getconj(bus[statemp1].V))。line[i].stoe=Complex::productComplex(bus[statemp1].V,dd)。cc=Complex::productComplex(bb , Complex::getconj(line[i].Y))。aa=Complex::productComplex(Complex::getconj(bus[statemp1].V), B)。statemp=line[i].start。=0。i{int statemp,endtemp。goto LOOP。//修正節(jié)點電壓// Bus::PrintfNodeV(bus,N)。gauss::gauss_slove(JacAug,x,NN)。//whether converbence看迭代是否結束if(icon==1){coutBus::JisuanJacAug(JacAug,X,bus,N)。//計算節(jié)點功率差值Bus::PrintfNodeScal(X,bus,N)。//計算節(jié)點功率Bus::JisuanNodeScal(X,bus,N)。LOOP:Bus::JisuanNodeI(X,bus,N)。x=new double[NN]。for(i=0。double **JacAug。//計算節(jié)點導納矩陣Bus::PrintfNodeDnz(X,N)。for(i=0。//輸出結點參數Complex **X。}Bus::ScanfBusData(bus)。ibus[i].=0。//輸出支路參數Bus *bus=new Bus[N]。//動態(tài)分配支路結構體Line::ScanfLineData(line)。N)==0){return 0。//輸入支路個數if((Lamp。//i作為整個程序的循環(huán)變量int N=Bus::ScanfBusNo()。謝謝include “” include include include include“” include“” include“” using namespace std。2007，35(1)：8184．[11] 方思立，劉增煌，孟慶和．大型汽輪發(fā)電機自并勵勵 磁系統的應用條件【J]．中國電力，1994，27(12)：61．63． FANG Si．Ii，LIU Zenghu鋤g，MENG Qin爭hc．m application conditions of large turbine generator selfexcitation system【J】．Electric Powef，1994，27(12)： 61．63．[12]梁小冰，黃方能．利用EMTDC進行長持續(xù)時間過程 的仿真研究【J】．電網技術，2002，26(9)：55．57． LIANG Xiaobing，HUANG Fan爭眥ng．How to cany out simulalion of long dul‘ation processes by use of EMTDC【J】．Power System 11echnology，2002，26(9)： 5557．[13]王卉，陳楷，彭哲，等．數字仿真技術在電力系統中 的應用及常用的幾種數字仿真工具【J】．繼電器，2004，32(21)：7l一75．wANG Hui，CHEN Kai，PENG zhe，et a1．Application of digital simulation眥hniques棚d severaJ simulation tools in power system[J】．Relay，2004，32(21)：71End2實例驗證 【例題】設有一系統網絡結線見圖1，各支路阻抗和各節(jié)點功率均已以標幺值標示于圖1中，其中節(jié)點2連接的是發(fā)電廠，設節(jié)點1電壓保持U1＝，試計算其中的潮流分布，請輸入節(jié)點數：n=5 請輸入支路數：m=7 請輸入平衡母線的節(jié)點號：ph=l 請輸入支路信息：BI=[ l 2 0．02+0．06i O l 1；1 3 0．08+0．24i 0 1 2；2 3 0．06+0．18i 0 l 3： 2 4 0．06+0．18i O l 4： 2 5 0．04+0．12i 0 l 5： 3 4 0．01+0．03i 0 l 6： 4 5 0．08+0．24i O 1 7] 請輸入節(jié)點信息：B2=[ 0 0 1．06 3；0．2+0．20i 0 1 1；一O．45一O．15i 0 l l；一0．40．05i 0 l 1；一0．6—0．1i 0 1 l] 請輸入節(jié)點號及對地阻抗： A=[l 0；2 0；3 0；4 0；5 O ] 請輸入修正值：ip=0．000 0l參考文獻[1]陳珩．電力系統穩(wěn)定分析[M]．北京：中國電力出版社，2002：139—187．[2]鄭阿奇．MATLAB實用教程[M]．北京：電子工業(yè)出版社，2005：1243．[3] 束洪春，孫士云，等．云電送粵交商流混聯系統全過 程動態(tài)電壓研究[J】．中國電力，2008，4l(10)：l4． SHU Hong—ch吼，SUN Shiyun，et a1．Research on fun prc39。J(:,no)=[]。J(no:)=[]。J(no:)=[]。end J(m,q)=X3。end J(p,q)=X1。P=2*i1。X4= C(i)+G(i,j)*f(i)B(i,j)*e(i)。X3= C(i)+G(i,j)*e(i)B(i,j)*f(i)。X2=C(i)+G(i,j)*e(i)+B(i,j)*f(i)。X1=D(i)+G(i,j)*f(i)B(i,j)*e(i)。Else if j=amp。J(m,q)=X4。J(p,q)=X2。m=p+1。J(p,q)=X1。q=2*j1。p=2*i1。X4=X1。X3=X2。X2=G(i,j)*e(i)+B(i,j)*f(i)。X1=G(i)*f(i)B(i,j)*e(i)。ifJ=amp。case 1 J(p,q)=X1。case 3 P=2*i1。%該子程序是用來求取jacci矩陣for i=1:n X3=0。End End case 2 Function for j=1:n J=jacci(Y,G,B,PQ,e,f,V,C,D,B2,n,ph,no)X1=G(i,j)*f(i)B(i,j)*e(i)。Df(k)=DX(i)。[x,y]=size(DX)。End ElseDX=[DX1,0,0]。DX(no+1)=0。x1=length(DX1)。endfunction [De,Df]=hxf(J,Df,ph,n,no)%該子函數是為求取De Df DX=JDf。end end end Df=Df’。p=p+1。else p=2*i1。Df(p)=P(i)P1。V1=e(i)^2+f(i)^2。endP1=C(i)*e(i)+D(i)*f(i)。For j=i:nC(i)=C(i)+G(i,j)*e(j)B(i,j)*f(j)。type resultm function [C,D,Df]=xxf(G,B,e,f,P,Q,n,B2,ph,V,no)%該子程序是用來求取Df for i=1:nIfi=phC(i)=0。end result2=[B1(:,6),B1(:,1),B1(:,2),real(b’),imag(b’),real(c’),imag(c’), real(b’+c’),imag(b’+c’)]。b(i)=v(B1(i1))*a(i)j*B1(i4)*v(B1(i))^2/2。resulte1=[A(:,1),real(v),imag(v),V,jd,real(S’)，imag(S’),real(B2(:1)),imag(B2(:1)),real(B2(:2)),imag(B2(:,2))]。V=abs(v)。end S(ph)=sum(hh)*v(ph)。end v=e+f*j。e=e+De；f=f+Df。whilemax(abs(De))ipamp。[De,Di]=hxf(J,D,F,ph,n,no)。[C,D,DF]=xxf(G,B,e,f,P,Q,n,B2,ph,V,no)。end Y(p,p)=Y(p,p)+1./(B1(i3)*B1(i5)^2+B1(i4)./2P=real(S)。V(i)=B2(i3)。For i=1:m S(i)=B2(i1)B2(i2)。for i=1：n End e(i)=real(B2(i3))。if A（i2）=0B=imag(Y)。e=zeros（1，n）；Y(p,q)=Y(p,q)。%它以矩陣形式存貯支路的情況，每行存貯一條支路 %第一列存貯支路的一個端點 %第二列存貯支路的另一個端點 %第三列存貯支路阻抗%第四列存貯支路的對地導納%第五列存貯變壓器的變比，注意支路為1 %第六列存貯支路的序號B2=input(‘n請輸入節(jié)點信息：B2=’)； %第一列為電源側的功率 %第二列為負荷側的功率 %第三列為該點的電壓值%第四列為該點的類型：1為PQ，2為PV節(jié)點，3為平衡節(jié)點 A=input（‘n請輸入節(jié)點號及對地阻抗：A=’）； ip=input（‘n請輸入修正值：ip=’）； %ip為修正值)?！尽俊尽縁unction tisco %這是一個電力系統潮流計算的程序 n=input（‘n請輸入節(jié)點數：n=’）； m=input(‘請輸入支路數：m=’)。在電網諧波量較高的地點，要逐步安裝諧波警報指示器，以便進一步分析諧波情況，并采取措施，保證電力設備安全運行。諧波電壓和電流應選取5次測量接近數值的算術平均值。即需測量或計算1111111119次諧波電壓和諧波電流。在該點測量諧波電壓，并在向用戶供電的線路的送電端測量諧波電流。為此，電力部門和用戶所安裝的電力電容器組，應根據實際存在的諧波情況，采取加裝串聯電抗器等措施，保證電力設備安全運行。在新的非線性用電設備接入電網前后，均要進行現場測量，檢查諧波電流、電壓正弦波形畸變率是否符合本規(guī)定。第三章 諧波管理第十七條 各級電力部門對電網的諧波情況，應定期進行測量分析。還應檢驗并采取措施防止發(fā)生諧振和諧波放大。用戶注入電網的諧波電流允許值 表2 見附表中(2)之3。第十四條 單相或三相不對稱非線性用電設備接入電網，按第十條和第十一條的規(guī)定執(zhí)行。注入電網的短時間諧波電流，不屬于本規(guī)定第九條和第十二條的限制范圍。但是，在公共連接點接入多個換流設備或交流調壓裝置，總容量超過表3的規(guī)定值時，應按第十條、第十一條和附錄中(二)的有關規(guī)定辦理。第十一條 用戶與電網連接點原有的總電壓正弦波形畸變率已超過表1規(guī)定的極限值的75%，或新增非線性用電設備向電網注入致諧波電流超過表2規(guī)定的允許值時，應由電力部門經過核算，確定該用電設備接入電網的技術措施和供電方案。任一用戶注入電網連接點的各次諧波電流均不得超過表2規(guī)定的允許值。表1 電網電壓正弦畸形畸變率極限值(相電壓)%，但經過測量和計算，下一級電網的電壓正弦形畸變率未超過表1規(guī)定的極限值時，110千伏電網電壓正弦波形畸變率可限制在2%以內。第七條 用戶用電設備對諧波電壓的要求較本規(guī)定的電壓正弦波疇變率極限更嚴格時，由用戶自行采取限制諧波電壓的措施。但如對方的國家標準或企業(yè)標準的全部或部分規(guī)定比本規(guī)定嚴格，則應按對方較嚴格的規(guī)定執(zhí)行。如用戶的非線性用電設備接入電網，增加或改變了電網的諧波值及其分布，特別是使與電網連接點的諧波電壓、電流升高，用戶必須采取措施，把諧波電流限制在允許的范圍內，方能接入電網運行。所需投資和設備由非線性用電設備的所屬單位負責。第三條 電網原有的諧波超過本規(guī)定的電壓正弦波形畸變率極限值時，應查明諧波源并采取措施，把電壓正弦波形畸變率限制在規(guī)定的極限值以內。為保證向國民經濟各部門提供質量合格的50赫茲電能，必須對各種非線性用電設備注入電網的諧波電流加以限制，以保證電網和用戶用電設備的安全經濟運行，特制訂本規(guī)定。第三篇：SD 1261984 電力系統諧波管理暫行規(guī)定電力系統諧波管理暫行規(guī)定SD126～84 第一章 總則第一條 電力系統中的諧波主要是治金、化工、電氣化鐵路等換流設備及其他非線性用電設備產生的。利用x(*)=x(k+1)+DX(k+1)進行多次迭代，通過迭代判據得到所需要的精度值即準確值x(*)。四運行結果節(jié)點導納矩陣經過五次迭代后的雅克比矩陣迭代次數以及節(jié)點電壓的幅值和相角（弧度數）節(jié)點注入功率和電流五 結果分析在這次學習和實際操作過程里：首先，對電力系統分析中潮流計算的部分特別是潮流計算的計算機算法中的牛頓拉夫遜法進行深入的研讀，弄明白了其原理、計算過程、公式推導以及設計流程。endS(5)=P(5)+Q(5)*sqrt(1)。for j=1:5P(i)=U(i)*U(j)*(G(i,j)*cos(a(i)a(j))+B(i,j)*sin(a(i)a(j)))+P(i)。endk=k+1。endfor i=1:4a(i)=a(i)+oa(i)。% x1為所求△x的列向量% 求節(jié)點電壓新值，準備下一次迭代for i=1:4oa(i)=x1(i