磁偶极子matlab程序

以下是一个简单的磁偶极子模拟程序的 MATLAB 代码：

% 定义常量
mu_0 = 4 * pi * 1e-7;  % 真空中的磁导率
m = [0; 0; 1];  % 磁矩的方向（向上）
M = 1e6;  % 磁矩的大小

% 定义计算区域
x = linspace(-2, 2, 50);  % x方向坐标
y = linspace(-2, 2, 50);  % y方向坐标
z = linspace(-2,2, 50);  % z方向坐标
[X, Y, Z] = meshgrid(x, y, z);  % 组成网格

% 计算磁场
Bx = zeros(size(X));  % 初始化x方向的磁场为0
By = zeros(size(Y));  % 初始化y方向的磁场为0
Bz = zeros(size(Z));  % 初始化z方向的磁场为0
for i = 1:numel(x)
    for j = 1:numel(y)
        for k = 1:numel(z)
            r = [X(i,j,k); Y(i,j,k); Z(i,j,k)];  % 当前点到磁矩的位置向量
            R = norm(r);  % 当前点到磁矩的距离
            B = mu_0 / (4 * pi) * (3 * dot(m, r) * r / R^5 - m / R^3);  % 计算磁场
            Bx(i,j,k) = Bx(i,j,k) + B(1);  % 累加磁场分量
            By(i,j,k) = By(i,j,k) + B(2);
            Bz(i,j,k) = Bz(i,j,k) + B(3);
        end
    end
end

% 绘制磁场
figure;
quiver3(X, Y, Z, Bx, By, Bz);
title('磁偶极子模拟');
xlabel('x(m)');
ylabel('y(m)');
zlabel('z(m)');
axis equal;

这个程序计算了一个沿 $z$ 轴方向的磁矩的磁场分布，计算区域为 $[-2,2]\times[-2,2]\times[-2,2]$。你可以根据需要修改常量和计算区域。