matlab磁滞回线

以下是MATLAB绘制磁滞回线的步骤和代码：

1.磁滞模型推导

磁滞回线是描述铁磁性材料在外加磁场作用下磁化状态随磁场变化的曲线。铁磁性材料的磁滞回线可以用Jiles-Atherton模型来描述，该模型是一种基于微观物理机制的宏观磁滞模型。该模型的方程如下：

$$\frac{dM}{dH}=\frac{1}{\mu_0}\frac{d}{dM}\left(\frac{W}{V}\right)$$

其中，$M$是磁化强度，$H$是磁场强度，$\mu_0$是真空磁导率，$W$是磁能密度，$V$是体积。

2.MATLAB绘制磁滞模型

2.1无应力磁滞回线

以下是MATLAB绘制无应力磁滞回线的代码：

clc;clear;
tspan=[-20000,20000];
M0=-1.54.*10.^6;
[H,M]=ode45('odeM',tspan,M0);
tspan=[20000,-20000];
M1=1.54.*10.^6;
[Hf,Mf]=ode45('odeM_',tspan,M1);
figure(1);plot(H,M); hold on;figure(1);plot(Hf,Mf);

2.2有应力磁滞回线

以下是MATLAB绘制有应力磁滞回线的代码：

clc;clear;
tspan=[-20000,20000];
M0=-1.54.*10.^6;
[H,M]=ode45('odeM',tspan,M0);
tspan=[20000,-20000];
M1=1.54.*10.^6;
[Hf,Mf]=ode45('odeM_',tspan,M1);
figure(1);plot(H,M); hold on;figure(1);plot(Hf,Mf);

3.结果展示

运行以上代码后，可以得到无应力磁滞回线和有应力磁滞回线的图像。