优化代码% 定义电偶极子位置和方向 r = [0 0 0]; p = [1 0 0]; % 确定绘制坐标系 [x,y,z] = meshgrid(-2:0.2:2,-2:0.2:2,-2:0.2:2); % 计算电场分布 e0 = 8.854e-12; k = 1/(4*pi*e0); r1 = [x(:)-r(1), y(:)-r(2), z(:)-r(3)]; r2 = sqrt(r1(:,1).^2+r1(:,2).^2+r1(:,3).^2); E = k*(3*(r1*p')*p-r1)./(r2.^5); % 绘制电力线 figure; quiver3(x,y,z,E(:,1),E(:,2),E(:,3),0.5); axis equal; xlabel('X'); ylabel('Y'); zlabel('Z');

可以对代码进行如下优化：

1. 减少重复计算：可以将 `r1(:,1).^2+r1(:,2).^2+r1(:,3).^2` 这个计算结果提前计算出来，避免重复计算。

2. 合并常数计算：可以将 `k*3` 这个常数计算提前计算出来，避免重复计算。

3. 使用向量化运算：可以使用向量化运算，避免循环和条件判断等造成的性能损耗。

优化后的代码如下：

% 定义电偶极子位置和方向
r = [0 0 0];
p = [1 0 0];

% 确定绘制坐标系
[x,y,z] = meshgrid(-2:0.2:2,-2:0.2:2,-2:0.2:2);

% 计算电场分布
e0 = 8.854e-12;
k = 1/(4*pi*e0);
r1 = [x(:)-r(1), y(:)-r(2), z(:)-r(3)];
r2 = sqrt(sum(r1.^2, 2));
c = k*3;
E = c*(r1*p')*p-r1./(r2.^5);

% 绘制电力线
figure;
quiver3(x,y,z,E(:,1),E(:,2),E(:,3),0.5);
axis equal;
xlabel('X');
ylabel('Y');
zlabel('Z');

通过减少重复计算、合并常数计算和使用向量化运算等方法，可以提高代码的运行效率。