箍缩磁滞回线matlab程序

在MATLAB中创建一个模拟磁滞回线（Hysteresis Loop）的程序，通常涉及到绘制B-H曲线，其中B代表磁场强度（Magnetic Flux Density），H代表磁场梯度（Magnetic Field Intensity）。以下是编写一个简单程序的基本步骤：

% 定义磁化常数、矫顽力和饱和磁化强度
Ms = 800; % 饱和磁化强度 (A/m)
Hc = 500; % 矫顽力 (A/m)
BHmax = Ms * Hc; % 磁滞回线的最大值

% 创建磁场梯度范围
H_range = -Hc:Hc; 

% 计算相应的磁场强度
B = zeros(size(H_range));
for i = 1:length(H_range)
    if abs(H_range(i)) < Hc
        B(i) = Ms * H_range(i); % 当在磁滞回线范围内时
    else
        B(i) = sign(H_range(i)) * BHmax; % 达到饱和点后的恒定值
    end
end

% 绘制磁滞回线
figure;
plot(B, H_range, 'b-o');
xlabel('磁场强度 (B)');
ylabel('磁场梯度 (H)');
title('磁滞回线示例');
hold on;
grid on;

% 添加回线标注
text(BHmax, -Hc, ['饱和点 (B = ', num2str(BHmax), ')'], 'HorizontalAlignment', 'right');
text(-BHmax, Hc, ['饱和点 (B = ', num2str(-BHmax), ')'], 'HorizontalAlignment', 'left');