Z>0半空间是相对介电常数为11的介质，Z<0的半空间是相对介电常数为2的介质 当电量为1C的点电荷放在介质分界面，用matlab绘制电场强度的空间三维分布。

好的，根据您的描述，电场强度的空间三维分布可以通过以下步骤和代码实现：

1. 定义空间范围和网格数量

对于 Z>0 和 Z<0 两个半空间，我们可以分别定义它们的范围和网格数量。假设 Z>0 的介质高度为 h1，Z<0 的介质高度为 h2，网格数量分别为 n1 和 n2。则可以定义如下：

h1 = 1; % Z>0 的介质高度
h2 = -1; % Z<0 的介质高度
n1 = 50; % Z>0 的网格数量
n2 = 50; % Z<0 的网格数量

2. 计算电场强度

根据点电荷在介质分界面处产生的电场公式，可以计算出空间中每个网格点的电场强度。假设点电荷位置为 [0, 0, 0]，电量为 1C，空气和介质的相对介电常数分别为 1 和 11（或 2），则可以定义如下的电场计算函数：

function E = calcElectricField(x, y, z, eps1, eps2)
% 计算点电荷在介质分界面处产生的电场强度
d1 = sqrt((x.^2 + y.^2 + (z-h1).^2).^3);
d2 = sqrt((x.^2 + y.^2 + (z-h2).^2).^3);
E = 1 ./ (4*pi*eps0) .* (1./d1 - 1./d2) .* (eps2-eps1) ./ (eps2+eps1);
end

其中，eps0 是真空介电常数。注意，这个函数的输入变量是 x、y、z 三个网格坐标，而不是网格编号。

3. 绘制空间分布图

最后，可以使用 Matlab 中的 surf 函数绘制电场强度的空间分布图。具体实现代码如下：

[x, y] = meshgrid(linspace(-1, 1, n1), linspace(-1, 1, n1));
z1 = h1 * ones(n1);
z2 = h2 * ones(n2);
E1 = calcElectricField(x, y, z1, 1, 11);
E2 = calcElectricField(x, y, z2, 1, 2);
figure;
surf(x, y, z1, E1); hold on;
surf(x, y, z2, E2);
colormap(jet);
colorbar;
xlabel('X'); ylabel('Y'); zlabel('Z');
view([45, 45]);
axis equal;
title('Electric Field Distribution');

这段代码首先用 meshgrid 函数生成 x、y 平面上的网格点坐标，然后计算出 z=0 和 z=h2 处的电场强度 E1 和 E2。最后以 surf 函数展示出来，其中 colormap 和 colorbar 用于将电场强度的大小映射到一种颜色渐变中。标题和坐标轴标签可根据需要进行修改。

希望这样的回答可以帮助到您。