MATLAB 导体电流磁场分布模型 - 代码示例

以下是一个简单的示例代码，用于计算一个导体电流为 2 安培的周围磁场分布：

% 定义常数
mu0 = 4 * pi * 1e-7; % 真空中的磁导率
I = 2; % 导体电流

% 定义坐标范围和步长
x = linspace(-1, 1, 100);
y = linspace(-1, 1, 100);
z = linspace(-1, 1, 100);
dx = x(2) - x(1);
dy = y(2) - y(1);
dz = z(2) - z(1);

% 创建网格
[X, Y, Z] = meshgrid(x, y, z);

% 计算磁场
Bx = zeros(size(X));
By = zeros(size(Y));
Bz = zeros(size(Z));
for i = 1:numel(x)
    for j = 1:numel(y)
        for k = 1:numel(z)
            r = [x(i), y(j), z(k)];
            dl = [dx, 0, 0];
            dB = cross(dl, r) * I / (norm(r)^3);
            Bx(i,j,k) = Bx(i,j,k) + dB(1);
            By(i,j,k) = By(i,j,k) + dB(2);
            Bz(i,j,k) = Bz(i,j,k) + dB(3);
        end
    end
end
B = sqrt(Bx.^2 + By.^2 + Bz.^2);

% 绘制磁场图像
figure;
slice(X, Y, Z, B, [0], [0], [0]);
xlabel('x');
ylabel('y');
zlabel('z');
title('Magnetic Field Distribution');
colorbar;

此代码将创建一个 100x100x100 的网格，然后计算沿 x 轴上的导体电流 2A 所产生的磁场分布。最后，使用 slice 函数在原点处绘制磁场图像。可以通过更改 x，y 和 z 向量的范围和步长来调整网格的大小和分辨率，以获得更精细的结果。