磁偶极子磁力线分布的matlab代码

以下是磁偶极子磁力线分布的 MATLAB 代码：

% 定义磁偶极子参数
m = 1; % 磁矩
a = 1; % 磁偶极子长度

% 定义网格范围和网格步长
x = linspace(-3, 3, 50);
y = linspace(-3, 3, 50);
z = linspace(-3, 3, 50);
[xx, yy, zz] = meshgrid(x, y, z);
dx = x(2) - x(1);
dy = y(2) - y(1);
dz = z(2) - z(1);

% 计算磁力线分布
Bx = zeros(size(xx));
By = zeros(size(yy));
Bz = zeros(size(zz));
for i = 1:numel(xx)
    r = sqrt(xx(i)^2 + yy(i)^2 + zz(i)^2);
    cos_theta = zz(i) / r;
    sin_theta = sqrt(xx(i)^2 + yy(i)^2) / r;
    cos_phi = xx(i) / sqrt(xx(i)^2 + yy(i)^2);
    sin_phi = yy(i) / sqrt(xx(i)^2 + yy(i)^2);
    Bx(i) = m * (2*cos_theta*sin_phi)/(r^3) * a^3;
    By(i) = m * (2*cos_theta*cos_phi)/(r^3) * a^3;
    Bz(i) = m * ((3*cos_theta^2-1)/(r^3)) * a^3;
end

% 绘制磁力线
figure;
streamline(xx, yy, zz, Bx, By, Bz);
axis equal;
xlabel('x');
ylabel('y');
zlabel('z');
title('磁偶极子磁力线分布');

这段代码可以绘制出磁偶极子的磁力线分布。需要注意的是，这里的磁偶极子是一个理想化的模型，实际情况下可能存在一些误差。