用matlab写电磁场的代码

以下是用MATLAB编写电磁场的代码示例：

% 定义常数 mu0 = 4pi10^-7; % 真空磁导率 epsilon0 = 8.8541878176*10^-12; % 真空介电常数

% 定义电荷和电流 q = 1.6*10^-19; % 电子电荷量 i = 1; % 电流强度

% 定义空间范围 [x, y, z] = meshgrid(-1:0.01:1, -1:0.01:1, -1:0.01:1);

% 计算电磁场 r = sqrt(x.^2 + y.^2 + z.^2); % 距离 Bx = mu0q/(4pir.^3).((3x.z)./(r.^2)-x); % 磁场x方向分量 By = mu0q/(4pir.^3).((3y.z)./(r.^2)-y); % 磁场y方向分量 Bz = mu0q/(4pir.^3).((2z.^2-x.^2-y.^2)./(r.^2)); % 磁场z方向分量 Ex = (1/(4piepsilon0))(q./(r.^2)).x./r; % 电场x方向分量 Ey = (1/(4piepsilon0))(q./(r.^2)).y./r; % 电场y方向分量 Ez = (1/(4piepsilon0))(q./(r.^2)).*z./r; % 电场z方向分量

% 绘制电磁场 quiver3(x, y, z, Bx, By, Bz); % 绘制磁场 hold on; quiver3(x, y, z, Ex, Ey, Ez); % 绘制电场 axis equal; xlabel('x'); ylabel('y'); zlabel('z'); title('Electric and Magnetic Fields of a Point Charge')