用M语言编写绘制基本磁化曲线和磁路未饱和时的磁化电流曲线的MATLAB 程序

在MATLAB中，M语言是一种用于数值计算和图形化的脚本语言，对于绘制磁化曲线和磁路未饱和时的磁化电流曲线，你可以创建一些函数来模拟这个过程。以下是一个简单的示例，假设我们有一个描述磁场与电流关系的函数magnetic_function：

function [B, I] = magnetic_curve(Hsat, I0, B0, n, I_range)
% Hsat: 饱和场强度
% I0: 起始电流
% B0: 初始磁通密度
% n: 磁化曲线指数
% I_range: 指定的电流范围

H = linspace(0, Hsat, 100); % 创建磁场范围
I = I0 * (H ./ Hsat).^n; % 根据磁化曲线方程计算电流

% 如果磁路未饱和，则取I在指定范围内的部分
if I_range(2) > Hsat
    saturated_index = find(H > Hsat);
    I(saturated_index:end) = [];
end

% 绘制磁化曲线
figure;
plot(H, I, 'LineWidth', 2);
xlabel('磁场 (安培/米)');
ylabel('电流 (安培)');
title('磁化曲线');
hold on;

% 如果有磁化电流曲线需求，同样处理
if nargin > 4
    plot(I_range, zeros(size(I_range)), 'r--', 'LineWidth', 2);
    hold off;
    legend('磁化曲线', '磁路未饱和范围');
end

运行此程序前，你需要先确定`Hsat`, `I0`, `B0`, `n`以及`I_range`的具体值。这个脚本首先计算了从0到`Hsat`的磁场对应的电流值，然后绘制了磁化曲线。如果指定了`I_range`，则会显示该范围内的线作为未饱和区的指示。