【正文】 ction dIm = transtatea(t,Im,K)R1=。L1=。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)。% 文件名：% 目的：仿真電流互感器的暫態(tài)過程%%%%%%%%%% 初始化%%%%%%%%%%%%clear。 %%清除變量%磁滯回線起始點B0a=0。I0a=0。B0b=0。I0b=0。B0c=0。I0c=0。B0d=0。I0d=0。%變量：磁通Ba=0。Bb=0。Bc=0。Bd=0。K=[。0]。 %變量：勵磁電感options=[]。%RK法參數(shù)設定，此處為默認值deltaT=。 %采樣時間間隔Tmax=。 %仿真時間tspan = [0,deltaT]。 %RK法時間寬度yzero = [0。0。0。0]。 %RK法初始值%勵磁電流變化率，象征磁滯回線的改換方向的點DI0a=0。DI0b=0。DI0c=0。DI0d=0。%%%%%%%%%%%%%%%%仿真開始，是個循環(huán)過程%%%%%%%%%%5 for i=1:Tmax/deltaT% 四階RK解微分方程,注意此為定步長RK[t,num_y] =rk4fixed(@transtatea,tspan,yzero,2,K)。% 更新參數(shù)，為下次RK做準備tspan=[t(length(t),1),t(length(t),1)+deltaT]。tt(i)=[t(length(t),1)]。yzero=num_y(2,:)。yzero=yzero39。Im(i,:)=num_y(2,:)。%取得勵磁電流變化率DI1a=diff(num_y(:,1))./diff(t)。DI1b=diff(num_y(:,2))./diff(t)。DI1c=diff(num_y(:,3))./diff(t)。DI1d=diff(num_y(:,4))./diff(t)。%取得當前狀態(tài)的勵磁電感值[Ka,B0a,I0a,Ba]=getLm(Ba,Im(i,1),B0a,I0a,DI1a,DI0a,i)。[Kb,B0b,I0b,Bb]=getLm(Bb,Im(i,2),B0b,I0b,DI1b,DI0b,i)。[Kc,B0c,I0c,Bc]=getLm(Bc,Im(i,3),B0c,I0c,DI1c,DI0c,i)。[Kd,B0d,I0d,Bd]=getLm(Bd,Im(i,4),B0d,I0d,DI1d,DI0d,i)。%為計算勵磁電感更新參數(shù)DI0a=DI1a。DI0b=DI1b。DI0c=DI1c。DI0d=DI1d。%輸出中間變量，電感和磁通K=[Ka。Kb。Kc。0]。end。%文件名： 功能:按磁滯回線計算動態(tài)電感可以使用167。%文件名： 功能：四階定步長龍格－庫塔法?？梢允褂?67。%文件：% 目的：得到系統(tǒng)微分方程 dIm=inv(L)*(Ur*IM),此文件適合于仿真涌流過程function dIm = transtatea(t,Im,K)R1=。L1=。Um=100000。dIm=(inv([(L1+K(1)),0,0,L1。0,(L1+K(2)),0,L1。0,0,(L1+K(3)),L1。L1,L1,L1,6*L1]))*(Um*[sin(314*t)。sin(314*t+pi*2/3)。sin(314*t+pi*4/3)。0][R1,0,0,R1。0,R1,0,R1。0,0,R1,R1。R1,R1,R1,6*R1]*Im)。%文件名： 目的：仿真電流互感器的暫態(tài)過程 可以使用167。%文件名： 功能:按磁滯回線計算動態(tài)電感可以使用167。%文件名： 功能：四階定步長龍格－庫塔法?？梢允褂?67。%文件：% 目的：得到系統(tǒng)微分方程 dIm=inv(L)*(Ur*IM),此文件適合于仿真涌流過程function dIm = transtatea(t,Im,K)R1=。L1=。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)。% 文件名： 目的：仿真電流互感器的暫態(tài)過程 可以使用167。%文件名： 功能:按磁滯回線計算動態(tài)電感可以使用167。%文件名： 功能：四階定步長龍格－庫塔法?？梢允褂?67。%文件：% 目的：得到系統(tǒng)微分方程 dIm=inv(L)*(Ur*IM),此文件適合于仿真涌流過程function dIm = transtatea(t,Im,K)R1=。L1=。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)。% 文件名： 目的：仿真電流互感器的暫態(tài)過程 可以使用167。%文件名：%功能:按磁滯回線計算動態(tài)電感%(I0,B0):回環(huán)開始點function [L,B0,I0,B]=getLm(B,I,B0,I0,DI1,DI0,begin)alfa=。beta=20。gama=2*。 C=。slope=9。Hz=。Bz=380。% 處理初始剩磁if (begin==1)Hzero=I0。Bzero=B0。 if(DI10) B=(alfa*B0(gama*I0pi/2))/(alfa*(atan(beta*(I0+C))+pi/2))*(atan(beta*(I+C))+pi/2)+1 /alfa*(gama*Ipi/2)。 L=(alfa*B0(gama*I0pi/2))/(alfa*(atan(beta*(I0+C))+pi/2))*beta/(1+(beta*(I+C))*(beta*(I+C)))+gama/alfa。 elseB=(alfa*B0(gama*I0+pi/2))/(alfa*(atan(beta*(I0C))pi/2))*(atan(beta*(IC))pi/2)+1/alfa*(gama*I+pi/2)。L=(alfa*B0(gama*I0+pi/2))/(alfa*(atan(beta*(I0C))pi/2))*beta/(1+(beta*(IC))*(beta*(IC)))+gama/alfa。 end。% 結束處理起始磁通else% 處理回環(huán)時的狀態(tài) if (IHz) if(DI0*DI10) I0=I。 B0=B。 end B=Bz+slope*(IHz)。 L=slope。 elseif(IHz) if(DI0*DI10) I0=I。 B0=B。 end B=Bz+slope*(I+Hz)。 L=slope。 else if(DI0*DI10) I0=I。 B0=B。 if(DI10)B=(alfa*B0(gama*I0+pi/2))/(alfa*(atan(beta*(I0C))pi/2))*(atan(beta*(IC))pi/2)+1/alfa*(gama*I+pi/2)。L=(alfa*B0(gama*I0+pi/2))/(alfa*(atan(beta*(I0C))pi/2))*beta/(1+(beta*(IC))*(beta*(IC)))+gama/alfa。 elseB=(alfa*B0(gama*I0pi/2))/(alfa*(atan(beta*(I0+C))+pi/2))*(atan(beta*(I+C))+pi/2)+1/alfa*(gama*Ipi/2)。L=(alfa*B0(gama*I0pi/2))/(alfa*(atan(beta*(I0+C))+pi/2))*beta/(1+(beta*(I+C))*(beta*(I+C)))+gama/alfa。 end else if(DI10)B=(alfa*B0(gama*I0pi/2))/(alfa*(atan(beta*(I0+C))+pi/2))*(atan(beta*(I+C))+pi/2)+1/alfa*(gama*Ipi/2)。L=(alfa*B0(gama*I0pi/2))/(alfa*(atan(beta*(I0+C))+pi/2))*beta/(1+(beta*(I+C))*(beta*(I+C)))+gama/alfa。 elseB=(alfa*B0(gama*I0+pi/2))/(alfa*(atan(beta*(I0C))pi/2))*(atan(beta*(IC))pi/2)+1/alfa*(gama*I+pi/2)。L=(alfa*B0(gama*I0+pi/2))/(alfa*(atan(beta*(I0C))pi/2))*beta/(1+(beta*(IC))*(beta*(IC)))+gama/alfa。 end。 end。 end。%結束處理回環(huán)時狀態(tài)end。%文件名： 功能：四階定步長龍格－庫塔法?？梢允褂?67。%文件：目的：得到系統(tǒng)微分方程 dIm=inv(L)*(Ur*IM),此文件適合于仿真涌流過程可以使用167。 % 文件名： 目的：仿真電流互感器的暫態(tài)過程 可以使用167。%文件名： 功能:按磁滯回線計算動態(tài)電感可以使用167。%文件名： 功能：四階定步長龍格－庫塔法?？梢允褂?67。%文件：目的：得到系統(tǒng)微分方程 dIm=inv(L)*(Ur*IM),此文件適合于仿真涌流過程可以使用167。 % 文件名： 目的：仿真電流互感器的暫態(tài)過程 可以使用167。%文件名： 功能:按磁滯回線計算動態(tài)電感可以使用167。%文件名： 功能：四階定步長龍格－庫塔法?？梢允褂?67。%文件名：目的：得到系統(tǒng)微分方程 dIm=inv(L)*(Ur*IM),此文件適合于仿真涌流過程可以使用167。 % 文件名： 目的：仿真電流互感器的暫態(tài)過程 可以使用167。%文件名： 功能:按磁滯回線計算動態(tài)電感可以使用167。%文件名： 功能：四階定步長龍格－庫塔法?？梢允褂?67。%文件：目的：得到系統(tǒng)微分方程 dIm=inv(L)*(Ur*IM),此文件適合于仿真涌流過程可以使用167。 56