MATLAB 实现双电荷电场分布计算与可视化：考虑传播衰减与时间因素

'); 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; ######[AI写代码神器 | 1099点数解答 | 2024-11-11 19:42:27]