用MATLAB画出电偶极子的E面和H面方向图的代码
时间: 2023-11-21 13:06:19 浏览: 37
以下是用MATLAB画出电偶极子的E面和H面方向图的代码:
```matlab
% 定义电偶极子的位置和电荷大小
x1 = -1; y1 = 0; q1 = 1;
x2 = 1; y2 = 0; q2 = -1;
% 定义观察点位置
x0 = 0; y0 = 1;
% 计算电场和磁场大小
k = 9e9; % 库仑常数
r1 = sqrt((x0-x1).^2 + (y0-y1).^2);
r2 = sqrt((x0-x2).^2 + (y0-y2).^2);
E1 = k*q1*(x0-x1)./r1.^3 + k*q2*(x0-x2)./r2.^3;
E2 = k*q1*(y0-y1)./r1.^3 + k*q2*(y0-y2)./r2.^3;
B = k*(q1-q2)/(r1^2+r2^2);
% 计算E面和H面方向
theta = atan2(y0-y1, x0-x1);
phi = atan2(y2-y0, x2-x0);
Etheta = E1*cos(theta) + E2*sin(theta);
Ephi = -E1*sin(theta) + E2*cos(theta);
Htheta = -B*cos(phi);
Hphi = B*sin(phi);
% 绘制E面和H面方向图
figure;
quiver(x0, y0, Etheta, Ephi, 'r', 'LineWidth', 2);
hold on;
quiver(x0, y0, Htheta, Hphi, 'b', 'LineWidth', 2);
axis equal;
xlim([-2 2]);
ylim([-2 2]);
title('E and H plane of a dipole');
xlabel('x');
ylabel('y');
legend('E plane', 'H plane');
hold off;
```
运行上述代码,就可以得到一个电偶极子的E面和H面方向图像。其中,红色箭头表示电场方向,蓝色箭头表示磁场方向。可以看到,E面和H面方向图像在观察点处都垂直于电偶极子的方向,这是符合电磁理论的预期结果的。