在 MATLAB 中绘制三维电偶极子的电势分布图和电场线图，可以按照以下步骤进行。此代码会创建一个三维网格，计算电势和电场，然后绘制相应的图形。

在MATLAB中，要绘制三维电偶极子的电势分布图和电场线图，你需要使用PDE Toolbox或者`surf`、`quiver3`等基本绘图函数结合数值积分或物理公式。这里提供一个简化的示例，我们将假设你知道电偶极子的位置和方向，以及电常数（通常为真空中的8.854e-12 F/m）。首先，你需要定义电偶极子的参数：

% 假设电偶极子位置 (x, y, z) 和方向 (dx, dy, dz)
polarization = [0.1; 0; 0]; % 电偶极矩向量
position = [0; 0; 0]; % 电偶极子中心点

% 电场强度公式 E = (k * p * r^(-3)) * n
% 其中 k 是电常数，p 是电偶极矩，r 是距离，n 是单位方向向量
electric_field_strength = @(r, n) (8.854e-12 * polarization ./ norm(polarization).^3) .* r;

接下来，我们创建一个三维网格并计算电场值：

% 创建一个立方体网格，这里设置大小为[-1, 1] x [-1, 1] x [-1, 1]
[x, y, z] = meshgrid(linspace(-1, 1, 50), linspace(-1, 1, 50), linspace(-1, 1, 50));

% 计算每个网格点到电偶极子的距离和单位方向向量
distances = sqrt(sum((repmat(position, size(x), 1) - [x(:), y(:), z(:)]) .^ 2, 2));
directions = [(x(:) - position(1)) / distances, (y(:) - position(2)) / distances, (z(:) - position(3)) / distances];

% 计算电场值
electric_fields = electric_field_strength(distances, directions);

现在我们可以用`surf`画出电势分布图：

% 将电场转换为电势 (电场强度乘以距离)
potentials = -electric_fields .* distances;

% 绘制电势分布
figure;
surf(x, y, z, reshape(potentials, size(x)));
xlabel('X');
ylabel('Y');
zlabel('Z');
title('Electric Potential Distribution');

对于电场线，我们可以使用`quiver3`：

% 使用quiver3绘制电场线
figure;
hold on;
quiver3(x(:), y(:), z(:), electric_fields(:, 1), electric_fields(:, 2), electric_fields(:, 3), 'Color', 'red');
quiver3(x(position(1)), y(position(2)), z(position(3)), polarization(1), polarization(2), polarization(3), 'Color', 'blue', 'LineWidth', 2);
view(3); % 设置视图
xlabel('X');
ylabel('Y');
zlabel('Z');
title('Electric Field Lines');
hold off;

请注意，这只是一个基础示例，实际应用中可能需要根据具体物理模型和边界条件调整。此外，如果你的电偶极子是非对称的，你可能还需要考虑使用矢量场图（`vectorfield3d`）或其他更复杂的绘图方法。