matlab用球画出电流密度分布

好的，以下是matlab代码。需要注意的是，这里的电流密度分布式仅适用于球内电场沿球心到某点的方向分布的情况。

a = 1; % 球半径
P = 1; % 极化强度
mu = 1; % 磁导率
N = 50; % 离散化数目

% 构造球体网格
[x,y,z] = sphere(N-1);
x = a*x;
y = a*y;
z = a*z;

% 计算电流密度
r = sqrt(x.^2 + y.^2 + z.^2);
Jx = (P/mu)*x./r;
Jy = (P/mu)*y./r;
Jz = (P/mu)*z./r;
J = cat(3,Jx,Jy,Jz);

% 画图
quiver3(x,y,z,J(:,:,1),J(:,:,2),J(:,:,3));
axis equal;

运行以上代码后，即可得到球内电流密度分布的可视化结果。如果需要考虑电场沿其他方向分布的情况，需要对电流密度分布式进行相应修改。