用MATLAB画出旋转形成的面电流分布半径为a的均匀极化介质球以它的一条直径为轴匀速旋转，角速度为w，该轴平行于球的极化强度P，

好的，这是一个比较具体的物理问题，需要一些数学和物理知识的支持。我会尽力回答，但是可能会有一些细节问题需要您自己解决。

首先，我们可以将问题简化为一个二维情况，即一个半径为a的圆盘以其直径为轴匀速旋转。这个圆盘代表了极化介质球的一个截面，因为整个球对称，所以可以在这个截面上考虑问题。

接下来，我们需要确定电流分布的表达式。根据题意，这是一个均匀极化介质，所以电流密度应该与电场强度成正比，而电场强度又与极化强度成正比。因此，电流密度应该是一个常数乘以极化强度。具体来说，我们可以假设极化强度沿着圆盘的半径方向分布，即

$P(r) = Pr/a$

其中r为距离圆盘中心点的距离，P为圆盘中心点的极化强度。根据安培环路定理，电流密度应该是一个环绕圆盘中心点的环流密度，即

$J(r) = \frac{2P(r)}{a}$

其中系数2是因为我们只考虑了一个半圆，而整个圆盘应该是由两个半圆组成的。

现在我们需要确定电流元素的表达式。考虑到圆盘的旋转，我们可以假设电流元素在圆盘上匀布，然后沿着圆盘旋转。因此，电流元素可以表示为

$dI = J(r)dS$

其中dS为圆盘上的面积元素。在极坐标系下，这个面积元素可以表示为

$dS = r dr d\theta$

其中r为距离圆心的距离，$\theta$为与圆盘直径的夹角。因此，电流元素可以表示为

$dI = \frac{2P(r)}{a} r dr d\theta$

接下来，我们需要将电流元素投影到一个垂直于旋转轴的平面上，这个平面可以表示为

$x = r\cos\theta$

$y = r\sin\theta$

因此，电流元素可以表示为

$dI = \frac{2P(x)}{a} x dx dy$

现在我们可以用MATLAB来绘制电流分布了。我们可以使用一个二重循环来计算每个点的电流密度，然后将其投影到平面上，最后使用quiver函数来绘制电流密度矢量。具体代码如下：

% 圆盘半径
a = 1;

% 极化强度
P = 1;

% 角速度
w = 1;

% 生成网格点
[X,Y] = meshgrid(-a:0.1:a, -a:0.1:a);

% 计算每个点的电流密度
Jx = zeros(size(X));
Jy = zeros(size(X));
for i = 1:size(X,1)
    for j = 1:size(X,2)
        r = sqrt(X(i,j)^2 + Y(i,j)^2);
        if r < a
            P_r = P * r / a;
            J = 2 * P_r / a;
            Jx(i,j) = -J * sin(atan2(Y(i,j), X(i,j)) + w*t);
            Jy(i,j) = J * cos(atan2(Y(i,j), X(i,j)) + w*t);
        end
    end
end

% 绘制电流密度矢量
quiver(X, Y, Jx, Jy);

这个代码会生成一个电流分布的图像，其中电流密度矢量的长度和方向代表了电流密度的大小和方向。注意到我们在计算时使用了atan2函数来计算每个点的角度，这是因为MATLAB的atan函数只能返回$-\pi/2$到$\pi/2$的值，无法区分不同象限的角度。

希望这个回答能够帮助您解决问题，如果还有其他问题，请随时提出！