基于matlab的电偶极子辐射代码

以下是一个基于Matlab的电偶极子辐射代码示例：

% 定义常数
c = 3e8; % 真空中光速
mu0 = pi*4e-7; % 真空磁导率
eps0 = 1/(mu0*c^2); % 真空介电常数

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

% 定义观察点的位置
theta = linspace(0, pi, 181); % 观察点的极角
phi = linspace(0, 2*pi, 361); % 观察点的方位角
[theta, phi] = meshgrid(theta, phi);
r = [sin(theta).*cos(phi), sin(theta).*sin(phi), cos(theta)]; % 观察点的位置

% 计算电偶极子辐射场的幅度和相位
k = 2*pi/(c/10); % 波矢量
kr = k*sqrt(sum((r-r0).^2, 2)); % 距离因子
E = (1/(4*pi*eps0))*exp(-1i*kr)./kr.*dot(p, (r-r0), 2); % 电场
E = reshape(E, size(theta)); % 调整电场的形状

% 绘制电偶极子辐射场的三维图像
figure;
surf(r(:, :, 1), r(:, :, 2), r(:, :, 3), abs(E));
xlabel('x');
ylabel('y');
zlabel('z');
title('Electric Dipole Radiation Field');

这个代码计算了一个沿x轴方向的电偶极子在原点处的辐射场，并在球形网格上绘制了这个辐射场的三维图像。通过调整电偶极子的方向和位置，可以计算任意方向和位置的电偶极子辐射场。