用MTALAB绘制铁磁性材料磁滞回线

1.磁滞模型推导

磁滞模型采用了J-A模型进行磁滞回线的绘制。J-A模型最终可推导出M-H的微分关系，根据此关系可利用MATLAB的数值方法解微分方程，最终绘制磁滞回线。

2.MATLAB绘制磁滞模型

2.1无应力磁滞回线

正向磁化磁滞模型，写成函数形式

function dM=odeM(H,M)
Ms=1.585.*10.^6;c=0.7;alpha=0.001;k=1000;a=600;
He=H+alpha.*M;
Man=Ms.*(coth(He./a)-a./He);
dM=((c.*Ms./a).*(1-(coth(He./a)).^2+(a./He).^2)+(1-c).*(Man-M)./(k*(1-c)-alpha.*(Man-M))).*(1./(1-alpha.*c));

反向磁化磁滞模型，写成函数形式

function dM=odeM_(H,M)
Ms=1.585.*10.^6;c=0.7;alpha=0.001;k=1000;a=600;
He=H+alpha.*M;
Man=Ms.*(coth(He./a)-a./He);
dM=((c.*Ms./a).*(1-