用matlab绘制无限长直导线的二维电场和磁场分布

假设无限长直导线的电流为I，直径为d。则该导线产生的磁场B为：

B = u0*I/(2*pi*r)

其中，u0是真空中的磁导率，r是距导线的距离。

电场E的计算比较复杂，需要使用椭圆积分。这里给出电场的表达式：

E = u0*I/(2*pi*d)*K(k)/r

其中，K(k)是第一类完全椭圆积分，k是椭圆的模，r是距导线的距离。

在MATLAB中，可以使用ellipk函数计算第一类完全椭圆积分。下面给出绘制电场和磁场分布的MATLAB代码：

% 定义常数
u0 = 4*pi*1e-7; % 真空中的磁导率
I = 1; % 电流
d = 0.01; % 直径

% 定义计算范围
x = linspace(-1, 1, 100);
y = linspace(-1, 1, 100);
[X,Y] = meshgrid(x,y);

% 计算电场和磁场
r = sqrt(X.^2 + Y.^2);
B = u0*I./(2*pi*r);
k = 2*sqrt(X.^2 + Y.^2)./d;
K = ellipk(k);
E = u0*I./(2*pi*d).*K./r;

% 绘制电场和磁场分布
figure;
subplot(1,2,1);
contourf(X,Y,E);
colorbar;
title('Electric field distribution');
xlabel('x (m)');
ylabel('y (m)');

subplot(1,2,2);
contourf(X,Y,B);
colorbar;
title('Magnetic field distribution');
xlabel('x (m)');
ylabel('y (m)');

执行以上代码，即可得到无限长直导线的电场和磁场分布图像。