【正文】 [13]黃紹平，[J].繼電器，2004，卷期頁. [14][M]，北京：機(jī)械工業(yè)出版社，2006. [15]萬凱，劉會(huì)金，陳允平．考慮鐵磁特性的有載調(diào)壓變壓器無暫態(tài)切換過程的理論研究[J].中國電機(jī)工程學(xué)報(bào)，2003，卷期頁.[16] 連小洲，[J].江西電力職業(yè)技術(shù)學(xué)院學(xué)報(bào)，2005，卷期頁.[17]張誠，廖勇，劉刃，[MJ].重慶：大學(xué)學(xué)報(bào)（自然科學(xué)版），2003，26（10）：6366.[18] Greene J D, Gross C A. Noniinear modeiing of transformer[J]. IEEE Trans. Ind. Appi.,1998，卷期頁.[19] Cavaiio A, Setoia R, Vasca F. Using MATLAB,SIMULINK and control system toolbox[M]. New Jersey: Prentice Haii,1994.[20] Niemela V A，Owen H A Jr， coupled secondarismodel for multiwinding transformers with parameter values calculated from shortcircllit impedances．PESC’90，SanAntonio，TX，USA，1990．附 錄：function I=KLBH(xh,h,T)fai=[,,0,,]。i=[,,,]。py=csapi(fai,i)。faiz=,fais=2,iz=,is=4。k=(faisfaiz)/(isiz),b1=faizk*iz,b2=faiz+k*iz。iR=inline(39。(faib1)/k39。,39。k39。,39。b139。,39。fai39。)。iL=inline(39。(faib2)/k39。,39。k39。,39。b239。,39。fai39。)。n=1for t=0:h:T if abs(xh(n))=faiz I(n)=fnval(py,xh(n))。 end if xh(n)faiz I(n)=iR(k,b1,xh(n))。 end if xh(n)faiz I(n)=iL(k,b2,xh(n))。 end n=n+1。end：U=1。fai1=[,,,]。i1=[,,,]。fai2=[,,,]。i2=[,,,]。fais=[],faiz=[]。is=[],iz=[]。pp1=csapi(i1,fai1)。fnplt(pp1,39。xk39。)。hold onpp2=csapi(i2,fai2)。fnplt(pp2,39。xk39。)。hold offgrid ontitle(39。磁滯回線39。)xlabel(39。勵(lì)磁電流i/.39。)ylabel(39。磁通:Q/.39。)axis([ ])h=pi/50,wt=0:h:4*pi。f=U*sin(wt)。figure(2)plot(wt,f)grid ontitle(39。磁通波形39。)ylabel(39。磁通:Q/.39。)xlabel(39。wt/rad39。)axis([0 14 ])py1=csapi(fai1,i1)。py2=csapi(fai2,i2)。f=diff(f)。n=1。for wt=h:h:4*pi if f(n)=0 i0(n)=fnval(py2,U*sin(wt)) end if f(n)0 i0(n)=fnval(py1,U*sin(wt)) end n=n+1。endwt=h:h:4*pi。figure(3)plot(wt,i0,39。*39。)grid ontitle(39?？紤]磁滯時(shí)的勵(lì)磁電流波形39。)ylabel(39。電流:i/.39。)xlabel(39。wt/rad39。)axis([0 14 ])：c2=,s2=,xk=,rk=。c1=,s1=。c0=1,s0=0。x=0::。du2=x*(rk*c2+xk*s2)。du1=x*(rk*c1+xk*s1)。du0=x*(rk*c0+xk*s0)。u0=1du0。u1=1du1。u2=1du2。plot(x,u2,39。x39。,x,u1,39。o39。,x,u0,x,1,39。39。)xlabel(39。負(fù)載系數(shù)B39。)。ylabel(39。u2/.39。)。title(39。變壓器的外特性39。)。legend(39。容性負(fù)載39。,39。感性負(fù)載39。,39。純電阻負(fù)載39。,3)axis([0 ])：f=inline(39。x*sn*c/(x*sn*c+p0+x.^2*pkn)39。,39。sn39。,39。c39。,39。p039。,39。pkn39。,39。x39。)。c=,sn=1,p0=,pkn=。n=1。for x=0:: h(n)=f(sn,c,p0,pkn,x)。 n=n+1。endx=0::plot(x,h)axis([0 0 1])xlabel(39。負(fù)載系數(shù)B39。)ylabel(39。變壓器效率n39。)title(39。效率特性曲線39。)：h=,T=,t=0:h:T。i=*sin(100*pi*t)+*sin(500*pi*t)+*sin(700*pi*t)。i2=*sin(100*pi*t)+*sin(500*pi*t)+*sin(700*pi*t)*sin(300*pi*t)。i3=*sin(300*pi*t)。figure(1)plot(t,i,t,i2,39。r39。,t,i3,39。r39。)xlabel(39。時(shí)間:t/s39。)。ylabel(39。電流:i/.39。)。title(39。近似正弦波加入三次諧波39。)。legend(39。近似正弦波39。,39。加入三次諧波后的波形39。,39。三次諧波39。)grid onfai=[,,,]。I=[,,,]。pp=csapi(I,fai)。figure(2)ct=fnval(pp,i)。ct2=fnval(pp,i2)。plot(t,ct,39。39。,t,ct2,39。r39。)xlabel(39。時(shí)間:t/s39。)。ylabel(39。磁通Q/.39。)。legend(39。原磁通波形39。,39。加入三次諧波時(shí)的磁通波形39。)grid ont=h:h:T。e=diff(ct)。e2=diff(ct2)。t=h:h:T。figure(3)plot(t,e/,39。39。,t,e2/,39。k39。)xlabel(39。時(shí)間:t/s39。)。ylabel(39。電壓:U/.39。)。legend(39。原電壓波形39。,39。加入三次諧波時(shí)的電壓波形39。)。grid on：h=,T=,rad=0,t=0:h:T。U1m=1,N1=1,f=50,w=2*pi*f。L1=,R1=。u=U1m/N1.*sin(w*t+rad)。X=。faiM=1。fai1=faiM.*cos(w*t+rad),fai2=(X+faiM).*cos(rad)*exp(R1/L1*t)。xh=fai1+fai2。plot(t,u)title(39。電壓波形39。)。xlabel(39。時(shí)間:t/s39。)。ylabel(39。電壓:u/.39。)。axis([0 2 2])。grid onfigure(2)plot(t,fai1,39。39。,t,fai2,39。:39。,t,xh,39。39。)title(39。磁通波形39。)。xlabel(39。時(shí)間:t/s39。)。ylabel(39。磁通:Q/.39。)。legend(39。強(qiáng)制分量Q~39。,39。自由分量Q=39。,39。Q~+Q=39。)。grid oni=KLBH(xh,h,T)。figure(3)plot(t,i)title(39。勵(lì)磁電流波形39。)。xlabel(39。時(shí)間:t/s39。)。ylabel(39。電流i/.39。)。grid on：L1=,R1=,w=100*pi,a=pi/2。B=R1/L1。h=,T=。ft=inline(39。B*X+310*sin(100*pi*t+a)39。,39。B39。,39。t39。,39。X39。,39。a39。)。X(1)=。n=1。for t=h:h:T k1=h*ft(B,(n1)*h,X(n),a)。 k2=h*ft(B,(n1)*h+h/2,X(n)+k1/2,a)。 k3=h*ft(B,(n1)*h+h/2,X(n)+k2/2,a)。 k4=h*ft(B,(n1)*h,X(n)+k3,a)。 X(n+1)=X(n)+(k1+2*k2+2*k3+k4)/6。 n=n+1。endt=0:h:Tplot(t,X)grid onxlabel(39。時(shí)間:t/s39。)ylabel(39。磁通Q/.39。)i=KLBH(X,h,T)。figure(2)plot(t,i)xlabel(39。時(shí)間:t/s39。)。ylabel(39。電流:i/.39。)grid on：T=,h=,t=0:h:T。f=50,w=2*pi*f,a=pi。Rk=,Lk=。i1=sin(w*t+aatan(w*Lk/Rk)),i2=sin(aatan(w*Lk/Rk))*exp(Rk/Lk*t)。i=i1+i2。plot(t,i,39。39。,t,i1,39。39。,t,i2,39。.39。)。legend(39。短路電流39。,39。強(qiáng)制分量39。,39。自由分量39。,3)。xlabel(39。時(shí)間:t/s39。)。ylabel(39。電流:i/.39。)。title(39。突發(fā)短路電流波形39。)grid on56