clear all; clc; % 常量定义 epsilon0 = 8.854187817e-12; % 真空介电常数 (f/m) mu0 = 4*pi*1e-7; % 真空磁导率 (h/m) q1 = single(1e-7); % 每个电荷的电量 (库伦) r1 = single(0.1); r2 = single(0.02); d = single(0.16); c = single(299792458); % 光速 (m/s) % 时间参数 t = linspace(0, 1, 2000); % 时间向量 (秒)，减少时间步长 % 创建空间网格，降低分辨率 [x, y, z_grid] = meshgrid(linspace(-10, 10, 44), linspace(-10, 10, 44), linspace(-10, 10, 44)); % 初始化电场数组 e_total_x = zeros(size(x), 'single'); e_total_y = zeros(size(y), 'single

'); e_total_z = zeros(size(z_grid), 'single'); % 计算各电荷贡献电场 e1_x = (1/(4*pi*epsilon0)) * (q1 * (x - d/2) ./ ((x - d/2).^2 + y.^2 + z_grid.^2).^1.5); e1_y = (1/(4*pi*epsilon0)) * (q1 * y ./ ((x - d/2).^2 + y.^2 + z_grid.^2).^1.5); e1_z = (1/(4*pi*epsilon0)) * (q1 * z_grid ./ ((x - d/2).^2 + y.^2 + z_grid.^2).^1.5); e2_x = (1/(4*pi*epsilon0)) * (q1 * (x + d/2) ./ ((x + d/2).^2 + y.^2 + z_grid.^2).^1.5); e2_y = (1/(4*pi*epsilon0)) * (q1 * y ./ ((x + d/2).^2 + y.^2 + z_grid.^2).^1.5); e2_z = (1/(4*pi*epsilon0)) * (q1 * z_grid ./ ((x + d/2).^2 + y.^2 + z_grid.^2).^1.5); % 计算电场传播时间 t_emit = d / c; % 计算电场传播距离 d_propagate = c * (t - t_emit); % 计算各电场在传播距离上的衰减 e1_propagate_x = e1_x .* exp(-1j * 2*pi * d_propagate / (c / sqrt(epsilon0 * mu0))); e1_propagate_y = e1_y .* exp(-1j * 2*pi * d_propagate / (c / sqrt(epsilon0 * mu0))); e1_propagate_z = e1_z .* exp(-1j * 2*pi * d_propagate / (c / sqrt(epsilon0 * mu0))); e2_propagate_x = e2_x .* exp(-1j * 2*pi * d_propagate / (c / sqrt(epsilon0 * mu0))); e2_propagate_y = e2_y .* exp(-1j * 2*pi * d_propagate / (c / sqrt(epsilon0 * mu0))); e2_propagate_z = e2_z .* exp(-1j * 2*pi * d_propagate / (c / sqrt(epsilon0 * mu0))); % 计算总电场 e_total_x = e1_propagate_x + e2_propagate_x; e_total_y = e1_propagate_y + e2_propagate_y; e_total_z = e1_propagate_z + e2_propagate_z; % 绘制电场分布 figure; quiver3(x(:), y(:), z_grid(:), real(e_total_x(:)), real(e_total_y(:)), real(e_total_z(:))); xlabel('X'); ylabel('Y'); zlabel('Z'); title('Electric Field Distribution'); xlim([-10, 10]); ylim([-10, 10]); zlim([-10, 10]); grid on; axis equal; [2024-11-11 19:42:27 | AI写代码神器 | 1099点数解答]