matlab绘图

半径为a的均匀细圆环，环上带4πεq的电荷

(1)求空间的电势分布

(2)求空间的电场分布

要求：以等势线、等势面表示电势；以箭头

、场线表示电场。

 

% 设置参数

a = 1; % 圆环半径

q = 1; % 电荷量

epsilon = 8.854e-12; % 真空介电常数

 

% 设置网格点

N = 100; % 网格数量

x = linspace(-3*a, 3*a, N); % x坐标范围

y = linspace(-3*a, 3*a, N); % y坐标范围

z = linspace(-3*a, 3*a, N); % z坐标范围

[X, Y, Z] = meshgrid(x, y, z); % 构建网格

 

% 计算电势

V = zeros(N, N, N); % 初始化电势矩阵

for i = 1:N

    for j = 1:N

        for k = 1:N

            r = sqrt(X(i,j,k)^2 + Y(i,j,k)^2 + Z(i,j,k)^2); % 点(r,theta,z)到原点的距离

            theta = atan2(Y(i,j,k), X(i,j,k)); % 点(r,theta,z)的极角

            phi = acos(Z(i,j,k)/r); % 点(r,theta,z)与z轴的夹角

            V(i,j,k) = (q/(4*pi*epsilon)) * cos(phi)/r; % 电势公式

        end

    end

end

 

% 绘制等势线和等势面

figure;

contourslice(X, Y, Z, V, [], [], [-2*a, -a, 0, a, 2*a], 'cubic');

xlabel('x'); ylabel('y'); zlabel('z');

title('等势线和等势面');

 

% 计算电场

[E_x, E_y, E_z] = gradient(-V, x, y, z);

 

% 绘制电场箭头

figure;

quiver3(X, Y, Z, E_x, E_y, E_z);

xlabel('x'); ylabel('y'); zlabel('z');

title('电场分布');

 

% 绘制电场场线

figure;

streamslice(X, Y, Z, E_x, E_y, E_z);

xlabel('x'); ylabel('y'); zlabel('z');

title('电场分布');