用matlab模拟两个赫姆霍兹线圈靠近时的磁感应强度三维曲线

以下是用MATLAB模拟两个赫姆霍兹线圈靠近时的磁感应强度三维曲线的代码：

% 定义常量
mu = 4*pi*1e-7; % 真空磁导率
I = 1; % 电流强度
a = 0.1; % 线圈半径
d = 0.2; % 线圈间距

% 定义网格
[x,y,z] = meshgrid(-0.5:0.01:0.5,-0.5:0.01:0.5,-0.5:0.01:0.5);

% 计算磁感应强度
Bx = zeros(size(x));
By = zeros(size(y));
Bz = zeros(size(z));
for i = 1:numel(x)
    r1 = sqrt((x(i)+d/2)^2 + y(i)^2 + z(i)^2);
    r2 = sqrt((x(i)-d/2)^2 + y(i)^2 + z(i)^2);
    Bx(i) = mu*I*a^2/(2*(r1^3-r2^3))*(x(i)+d/2-x(i)+d/2);
    By(i) = mu*I*a^2/(2*(r1^3-r2^3))*(y(i)-y(i));
    Bz(i) = mu*I*a^2/(2*(r1^3-r2^3))*(z(i)-z(i));
end

% 画图
quiver3(x,y,z,Bx,By,Bz);
xlabel('X');
ylabel('Y');
zlabel('Z');
title('磁感应强度三维曲线');

运行上述代码，即可得到磁感应强度三维曲线的图像。其中，蓝色箭头表示磁感应强度的方向和大小。