给出用MATLAB表示电偶极子电力线的代码

以下是用 MATLAB 表示电偶极子电力线的代码，其中包括绘制电力线的函数、电荷分布函数以及电场强度计算函数：

% 绘制电力线函数
function draw_electric_field_lines()
    % 定义网格
    [x, y] = meshgrid(-2:0.1:2);
    % 定义电场强度
    E = electric_field(x, y);
    % 绘制等势线和电力线
    contour(x, y, potential(x, y), 20, 'k');
    hold on;
    streamline(x, y, -E(:,:,1), -E(:,:,2), linspace(-2, 2, 30), linspace(-2, 2, 30));
    axis equal;
    xlabel('x');
    ylabel('y');
end

% 定义电荷分布函数
function rho = charge_density(x, y)
    % 定义电偶极子的位置和电荷量
    x1 = -1; y1 = 0; q1 = 1;
    x2 = 1; y2 = 0; q2 = -1;
    % 计算距离和电荷分布
    r1 = sqrt((x-x1).^2 + (y-y1).^2);
    r2 = sqrt((x-x2).^2 + (y-y2).^2);
    rho = q1./r1.^3 - q2./r2.^3;
end

% 定义电场强度计算函数
function E = electric_field(x, y)
    % 定义电荷密度
    rho = charge_density(x, y);
    % 计算电场强度
    [Ex, Ey] = gradient(potential(x, y));
    E(:,:,1) = -Ex;
    E(:,:,2) = -Ey;
    % 排除奇点
    E(isinf(E)) = NaN;
end

% 定义电势函数
function V = potential(x, y)
    % 定义电荷密度
    rho = charge_density(x, y);
    % 计算电势
    V = -1/(4*pi)*conv2(rho, 1./sqrt(x.^2+y.^2), 'same');
end

运行 `draw_electric_field_lines()` 函数即可绘制出电偶极子的电力线。