用 MATLAB绘制静电场线图

静电场线图是在三维空间中绘制等势线，其代码实现可以用 MATLAB 的 `contourslice` 函数。以下是一个静电场线图的 MATLAB 代码示例：

% 设置电荷和观测点的参数
k = 9e9; % 库仑常数
q = 1e-9; % 电荷电量，单位是库仑
xq = 0; % 电荷所在的 x 坐标，单位是米
yq = 0; % 电荷所在的 y 坐标，单位是米
zq = 0; % 电荷所在的 z 坐标，单位是米
xmin = -1; % 绘图区域的 x 最小值，单位是米
xmax = 1; % 绘图区域的 x 最大值，单位是米
ymin = -1; % 绘图区域的 y 最小值，单位是米
ymax = 1; % 绘图区域的 y 最大值，单位是米
zmin = -1; % 绘图区域的 z 最小值，单位是米
zmax = 1; % 绘图区域的 z 最大值，单位是米

% 设置绘图参数
numContours = 10; % 等势面数量
numLines = 30; % 字段线数量
sliceSpacing = 0.1; % 切片间距

% 计算电势
[X,Y,Z] = meshgrid(xmin:0.01:xmax, ymin:0.01:ymax, zmin:0.01:zmax);
V = k*q./sqrt((X-xq).^2 + (Y-yq).^2 + (Z-zq).^2);

% 绘制静电场线图
figure;
contourslice(X, Y, Z, V, [], [], linspace(zmin, zmax, numLines), numContours);
xlabel('x (m)');
ylabel('y (m)');
zlabel('z (m)');
title('静电场线图');
axis tight;

该代码将绘制一个以原点为中心的点电荷产生的静电场线图。你可以根据需要调整电荷和观测点的参数、绘图区域、等势面数量、字段线数量和切片间距。