三相变压器建模及仿真研究 %文件名：rk4fixed.m
功能：四阶定步长龙格－库塔法。 可以使用§1中的rk4fixed.m程序。
%文件：trastate.m
% 目的：得到系统微分方程 dIm=inv(L)*(U-r*IM),此文件适合于仿真涌流过程
function dIm = transtatea(t,Im,K)
R1=0.7925; L1=0.07336; Um=100000;
dIm=(inv([(L1+K(1)),0,0;0,(L1+K(2)),0;0,0,(L1+K(3))]))*(Um*[sin(314*t);sin(314*t+pi
*2/3);sin(314*t+pi*4/3)]-[R1,0,0;0,R1,0;0,0,R1]*Im);
附录D.在Yy0接线方式下两段反正切函数拟和极限磁滞回环的程序
% 文件名：transa.m
目的：仿真电流互感器的暂态过程 可以使用§1中的transa.m程序。
%文件名：getLm.m
功能:按磁滞回线计算动态电感 可以使用§1中的getLm.m程序。
%文件名：rk4fixed.m
功能：四阶定步长龙格－库塔法。 可以使用§1中的rk4fixed.m程序。
35
三相变压器建模及仿真研究 %文件：trastate.m
% 目的：得到系统微分方程 dIm=inv(L)*(U-r*IM),此文件适合于仿真涌流过程
function dIm = transtatea(t,Im,K)
R1=0.7925; L1=0.07336; Um=100000;
dIm=(inv([(3*L1+2*K(1)),-K(2),-K(3);-K(1),(3*L1+2*K(2)),-K(3);-K(1),-K(2),(3*L1
+2*K(3))]))*(3*Um*[sin(314*t);sin(314*t+pi*2/3);sin(314*t+pi*4/3)]-[3*R1,0,0;0,3*R1,0;0,0,3*R1]*Im);
附录E.在Yd11接线方式下两段反正切函数加两段直线拟和极限磁滞回环的程序
% 文件名：transa.m
目的：仿真电流互感器的暂态过程 可以使用§1中的transa.m程序。
%文件名：getLm.m
%功能:按磁滞回线计算动态电感 %(I0,B0):回环开始点
function [L,B0,I0,B]=getLm(B,I,B0,I0,DI1,DI0,begin)
alfa=0.00377; beta=20;
gama=2*0.016867; C=0.2227; slope=9; Hz=1.2; Bz=380;
% 处理初始剩磁 if (begin==1) Hzero=I0; Bzero=B0;
36
三相变压器建模及仿真研究 if(DI1<0)
B=(alfa*B0-(gama*I0-pi/2))/(alfa*(atan(beta*(I0+C))+pi/2))*(atan(beta*(I+C))+pi/2)+1
/alfa*(gama*I-pi/2);
L=(alfa*B0-(gama*I0-pi/2))/(alfa*(atan(beta*(I0+C))+pi/2))*beta/(1+(beta*(I+C))*(bet
a*(I+C)))+gama/alfa;
else
B=(alfa*B0-(gama*I0+pi/2))/(alfa*(atan(beta*(I0-C))-pi/2))*(atan(beta*(I-C))-pi/2)
+1/alfa*(gama*I+pi/2);
L=(alfa*B0-(gama*I0+pi/2))/(alfa*(atan(beta*(I0-C))-pi/2))*beta/(1+(beta*(I-C))*(
beta*(I-C)))+gama/alfa;
end;
% 结束处理起始磁通 else
% 处理回环时的状态 if (I>Hz)
if(DI0*DI1<0) I0=I; B0=B; end
B=Bz+slope*(I-Hz); L=slope; elseif(I<-Hz)
if(DI0*DI1<0) I0=I; B0=B; end
B=-Bz+slope*(I+Hz); L=slope; else
if(DI0*DI1<0) I0=I; B0=B; if(DI1>0)
B=(alfa*B0-(gama*I0+pi/2))/(alfa*(atan(beta*(I0-C))-pi/2))*(atan(beta*(
I-C))-pi/2)+1/alfa*(gama*I+pi/2);
L=(alfa*B0-(gama*I0+pi/2))/(alfa*(atan(beta*(I0-C))-pi/2))*beta/(1+(bet
a*(I-C))*(beta*(I-C)))+gama/alfa;
else
B=(alfa*B0-(gama*I0-pi/2))/(alfa*(atan(beta*(I0+C))+pi/2))*(atan(beta*(
I+C))+pi/2)+1/alfa*(gama*I-pi/2);
L=(alfa*B0-(gama*I0-pi/2))/(alfa*(atan(beta*(I0+C))+pi/2))*beta/(1+(bet
37
三相变压器建模及仿真研究 a*(I+C))*(beta*(I+C)))+gama/alfa;
end else
if(DI1<0)
B=(alfa*B0-(gama*I0-pi/2))/(alfa*(atan(beta*(I0+C))+pi/2))*(atan(beta*(
I+C))+pi/2)+1/alfa*(gama*I-pi/2);
L=(alfa*B0-(gama*I0-pi/2))/(alfa*(atan(beta*(I0+C))+pi/2))*beta/(1+(bet
a*(I+C))*(beta*(I+C)))+gama/alfa;
else
B=(alfa*B0-(gama*I0+pi/2))/(alfa*(atan(beta*(I0-C))-pi/2))*(atan(beta*(
I-C))-pi/2)+1/alfa*(gama*I+pi/2);
L=(alfa*B0-(gama*I0+pi/2))/(alfa*(atan(beta*(I0-C))-pi/2))*beta/(1+(bet
a*(I-C))*(beta*(I-C)))+gama/alfa;
end; end; end;
%结束处理回环时状态 end;
%文件名：rk4fixed.m
功能：四阶定步长龙格－库塔法。 可以使用§1中的rk4fixed.m程序。
%文件：trastate.m
目的：得到系统微分方程 dIm=inv(L)*(U-r*IM),此文件适合于仿真涌流过程 可以使用§1中的trastate.m程序。
附录F.在Ynd11接线方式下两段反正切函数加两段直线拟和极限磁滞回环的程序
% 文件名：transa.m
目的：仿真电流互感器的暂态过程 可以使用§2中的transa.m程序。
38
三相变压器建模及仿真研究
%文件名：getLm.m
功能:按磁滞回线计算动态电感 可以使用§5中的getLm.m程序。
%文件名：rk4fixed.m
功能：四阶定步长龙格－库塔法。 可以使用§2中的rk4fixed.m程序。
%文件：trastate.m
目的：得到系统微分方程 dIm=inv(L)*(U-r*IM),此文件适合于仿真涌流过程 可以使用§2中的trastate.m程序。
附录G.在Yny0接线方式下两段反正切函数加两段直线拟和极限磁滞回环的程序
% 文件名：transa.m
目的：仿真电流互感器的暂态过程 可以使用§3中的transa.m程序。
%文件名：getLm.m
功能:按磁滞回线计算动态电感 可以使用§5中的getLm.m程序。
%文件名：rk4fixed.m
功能：四阶定步长龙格－库塔法。
39
三相变压器建模及仿真研究 可以使用§3中的rk4fixed.m程序。
%文件名：trastate.m
目的：得到系统微分方程 dIm=inv(L)*(U-r*IM),此文件适合于仿真涌流过程 可以使用§3中的trastate.m程序。
附录H.在Yy0接线方式下两段反正切函数加两段直线拟和极限磁滞回环的程序
% 文件名：transa.m
目的：仿真电流互感器的暂态过程 可以使用§4中的transa.m程序。
%文件名：getLm.m
功能:按磁滞回线计算动态电感 可以使用§5中的getLm.m程序。
%文件名：rk4fixed.m
功能：四阶定步长龙格－库塔法。 可以使用§4中的rk4fixed.m程序。
%文件：trastate.m
目的：得到系统微分方程 dIm=inv(L)*(U-r*IM),此文件适合于仿真涌流过程 可以使用§4中的trastate.m程序。
40