matlab实现的磁滞回线的源码,并附有中文解释.zip

%磁滞回线 %采用三次样条插值画磁滞回线 %通过磁滞回线画磁通为正弦时磁化电流波形 clear; %清除workspace中变量 H_data=[ 0.005 0.026 0.050 0.065 0.078 0.096 0.126 0.164 0.270 0.460... 0.529 0.471 0.341 0.203 0.074 0.000 -0.028 -0.049 -0.063 -0.077... -0.093 -0.124 -0.162 -0.270 -0.460 -0.530 -0.469 -0.338 -0.200 -0.070]; B_data=[ -0.174 -0.152 -0.123 -0.089 -0.051 -0.003 0.092 0.177 0.253 0.297... 0.308 0.298 0.279 0.252 0.217 0.182 0.157 0.126 0.092 0.055... 0.000 -0.087 -0.174 -0.249 -0.293 -0.303 -0.294 -0.274 -0.247 -0.212]; H=H_data;B=B_data; %磁滞回线的数据 H=H.*1e4;B=B.*5; subplot(2,2,1); Bx=-pi/2:0.01:3*pi/2; Bsin=1.5*sin(Bx); %计算正弦值，扩大1.5倍适应磁滞曲线 plot(Bx,Bsin); %画磁场密度正弦曲线 grid on; xlim([-pi/2,3*pi/2]); %限定X轴显示范围 xlabel('wt'); %X轴标号 ylabel('磁场密度B(T)'); %Y轴标号 subplot(2,2,2); hold on %保持图形 plot(H,B,'ro'); %画磁滞回线 B1=B(7:11);H1=H(7:11); %第1象限数据B>0数据 B2=B(11:20);H2=H(11:20); %第2象限数据B>0数据 B3=B(21:26);H3=H(21:26); %第3象限数据B>0数据 B4=[B(26:30),B(1:6)];H4=[H(26:30),H(1:6)]; %第4象限数据B>0数据 BB1=[B1,B4];HH1=[H1,H4]; %磁滞回线的下分支 XI1=-6290:10:6290; YI1=interp1(HH1,BB1,XI1,'spline'); %用3次样条插值计算 plot(XI1,YI1); %画插值后回线下分支 BB2=[B2,B3];HH2=[H2,H3]; %磁滞回线的上分支 XI2=-5290:10:5290; YI2=interp1(HH2,BB2,XI2,'spline'); %用3次样条插值计算 plot(XI2,YI2); %画插值后回线上分支 grid on; xlabel('磁化强度H(A/M)'); %X轴标号 ylabel('磁场密度B(T)'); %Y轴标号 M_X=1.39; % XI1=sin((-M_X:0.01:M_X)./M_X.*pi./2).*M_X; %磁通正弦变化 YI1=interp1(BB1,HH1,XI1,'spline'); %用3次样条插值计算 subplot(2,2,4); hold on; len_X=length(XI1); %计算先前波形的横坐标长度 X1=(1:len_X); X1=X1/len_X*pi; plot(X1,YI1); %画电流波形 XI2=sin((-M_X:1:M_X)./M_X.*pi./2).*M_X; YI2=interp1(BB1,HH1,XI1,'spline'); %用3次样条插值计算 X2=(1:len_X)+len_X-1; X2=X2/len_X*pi; %平移到先前波形的右侧 plot(X2,-YI2); %画电流波形 grid on; xlim([0,2*pi]); %限定X轴显示范围 xlabel('wt'); %X轴标号 ylabel('磁化电流I(安匝)'); %Y轴标号