测量的磁滞回线数据使用matlab拟合的方法.zip

clc clear all Hmin=0;Hmax=0;H1=0;H2=0; %Hmin为极小值，Hmax为极大值，H1为前一点，H2为后一点 n=1;b=0;z_zhe=zeros(1,2);B=zeros(1,200); for t=0:0.0001:0.1 H=30*sin(2*pi*50*t); %一次电流表达式 % H=20*sin(2*pi*50*t); % H=5*sin(2*pi*50*t); H2=H; if H2>=30 %判断集磁环有无进入饱和 B(n)=0.002*H2+1.44; %饱和段基本磁化曲线 b=1; n=n+1; continue; elseif H2<=-30 B(n)=0.002*H2-1.44; b=-1; n=n+1; continue; end if Hmin==0&Hmax==0&b==0 %判断是否在初始磁化曲线上 if H2>H1||H2==0 B(n)=BH_initial(H2); H1=H2; n=n+1; continue; else Hmax=H1;Bmax=B(n-1); end elseif Hmin==0&Hmax==0&b==1 %判断是否在饱和磁滞回线的下降支上 B(n)=BH_limit(H2); n=n+1; H1=H2; continue; elseif Hmin==0&Hmax==0&b==-1 %判断是否在饱和磁滞回线的上升支上 B(n)=-BH_limit(-H2); n=n+1; H1=H2; continue; end if H2>H1 %判断是否出现转折 if z_zhe(1)==0 z_zhe(1)=1; else z_zhe(2)=1; end elseif H2<H1 if z_zhe(1)==0 z_zhe(1)=2; else z_zhe(2)=2; end end if z_zhe(1)+z_zhe(2)==3 %出现转折进行赋值 if z_zhe(1)==1 Hmax=H1;Bmax=B1; else Hmin=H1;Bmin=B1; end z_zhe(1)=z_zhe(2); z_zhe(2)=0; end if H2>H1 %求解次磁滞回环上的磁化曲线 B(n)=-BH_limit(-H2)-BH_limit(Hmin)+Bmin+2*F(-Hmin)*F(H2); B1=B(n-1); H1=H2; elseif H2<H1 B(n)=BH_limit(H2)+BH_limit(-Hmax)+Bmax-2*F(Hmax)*F(-H2); B1=B(n-1); H1=H2; end n=n+1; end t=0:0.0001:0.1; H=60*sin(2*pi*50*t); %B=H/(0.03/(0.1*pi*4*pi*0.0000001)); %plot(t,B); plot(H,B); %hold on %plot(t,H,'red');