用matlab绘制无限长导线的磁场分布

无限长导线的磁场分布可以使用比奥萨伐尔-拉普拉斯定理求解。以下是MATLAB代码示例：

clear all
close all
clc

% 定义常数
u_0 = 4*pi*1e-7; % 真空磁导率
I = 1; % 电流强度
a = 0.01; % 导线半径

% 定义网格
x = linspace(-0.1,0.1,50);
y = linspace(-0.1,0.1,50);
[X,Y] = meshgrid(x,y);

% 计算磁场分布
R = sqrt(X.^2+Y.^2);
theta = atan2(Y,X);

B_theta = (u_0*I)/(2*pi*R).*sin(theta).*(R>a);
B_z = (u_0*I)/(2*pi)*log(R./a);

Bx = -B_theta.*sin(theta);
By = B_theta.*cos(theta);
Bz = B_z;

% 绘制磁场图像
quiver3(X,Y,zeros(size(X)),Bx,By,Bz);
xlabel('x');
ylabel('y');
zlabel('z');
title('无限长导线磁场分布');

该代码会生成一个三维矢量场图像，其中箭头的长度和方向表示磁场强度和方向，箭头的位置表示空间中的坐标。在该图像中，磁场沿着导线轴向延伸，沿着导线周围的环向分布。