基于Matlab的电磁场与波“镜像电荷法”仿真——导体球外点电荷的几种情况

已于 2022-04-03 19:09:20 修改
阅读量6.4k 收藏 110
点赞数 19
分类专栏： 仿真 文章标签： matlab 开发语言 电磁学
于 2022-01-29 23:31:28 首次发布
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/qq_55710894/article/details/122745514
版权

目录

零、使用到的绘制电场线方向的原创函数

             drawEline函数 

1、目的

2、源代码

3、注意事项

一、导体球外点电荷的镜像法（导体球接地情况与不接地情况）

1.目的：

2.理论分析：

3.源代码与结果（分为接地与不接地两种子情况，每种情况又有二维三维两种仿真）：

3.1接地导体球外单个点电荷：

3.2不接地导体球外单个点电荷：

3.3接地导体球外，直径所在直线上左右对称放置两个正点电荷：

3.4接地导体球外，与球心连线夹角为90度放置两个正点电荷

3.5同心球壳，外壳接地，球壳内部有一点电荷(p49)

零、使用到的绘制电场线方向的原创函数

             drawEline函数 

作者：顾世杰，戴桢琦

1、目的

本次MATLAB仿真几乎每一个项目的二维图都需要画电场线，按照老师的要求，电场线头尾必须与电荷连接，因此我们选择使用streamline函数画电场线，由于streamline函数本身画出的线不带箭头，我们通过研究上一届学长学姐们的资料，研究并改良了其中画箭头的代码语句，将其做了必要的扩充，并封装成一个函数M文件。

2、源代码

function pipi=drawEline(h1,d1,d2)

posAxes = get(gca, 'Position');  %获取当前坐标区位置信息
posX = posAxes(1);
posY = posAxes(2);
width = posAxes(3);
height = posAxes(4);%分解位置信息
hold on;
limX = get(gca, 'Xlim');%获取x轴范围
limY = get(gca, 'Ylim');%获取y轴范围
minX = limX(1);
maxX = limX(2);
minY = limY(1);
maxY = limY(2);
x=[h1(1).XData(d1) h1(1).XData(d2)];
y=[h1(1).YData(d1) h1(1).YData(d2)];
xNew = posX + (x - minX) / (maxX - minX) * width;
yNew = posY + (y - minY) / (maxY - minY) * height;%通过坐标变换得到箭头的绝对坐标
pipi=annotation('arrow', xNew, yNew);

3、注意事项

该函数M文件必须与调用它的M文件放在同一个文件夹中，并且要与语句h1=streamline(…);配套使用

一、导体球外点电荷的镜像法（导体球接地情况与不接地情况）

                                                                  作者：顾世杰

1.目的：

通过Matlab中的相关函数公式，仿真课本上的导体球外镜像法，进一步理解镜像法的相关原理。

2.理论分析：

由真空中的库伦定律可知之间的力F满足：

其中

由电场强度E的定义可以得点电荷的电场强度为：

则点电荷形成的场内一点P处的电为为：

    

进而有公式：

     

通过使用①与②两个公式，即可使用Matlab计算出点电荷的电势以及电场强度，进而可以通过软件中的相关函数绘出电场线和等位线。

3.源代码与结果（分为接地与不接地两种子情况，每种情况又有二维三维两种仿真）：

3.1接地导体球外单个点电荷：

3.1.1计算公式：

其中电荷量：q=1e-6  

常系数：k= 8.9875e9

3.1.2二维图代码：

k = 8.9875e9;
R = 0.05;
d = 0.15;
b = R^2 / d;
q1 = 1e-6;
q2 = -R * q1 / d;

x = linspace(-0.3,0.3,60);
y = linspace(-0.16,0.16,60);
[X,Y] = meshgrid(x);

r1 = sqrt((X - d).^2 + Y.^2);
r2 = sqrt((X - b).^2 + Y.^2);
r3 = sqrt(X.^2 + Y.^2);
U = k * q1./r1 + k * q2./r2;

a = 0:pi/100:2*pi;
x = R * sin( a );
y = R * cos( a );
plot ( x,y,'k-');
hold on;
plot( b,0,'r.','MarkerSize',30);
plot( 0,0,'r.','MarkerSize',30);
plot( d,0,'k.','MarkerSize',30);
xlabel('x','fontsize',16);
ylabel('y','fontsize',16);

u = linspace(0,800000,201);
[C,H] = contour(X,Y,U,u,'k-');
contour(X,Y,U,u,'c-');
% clabel(C,H);
% grid on;

[Ex,Ey] = gradient(-U);
E = sqrt(Ex.^2 + Ey.^2);
Ex = Ex./E;
Ey = Ey./E;
th=linspace(0,2*pi,50);
x1=d+0.005*cos(th);
y1=0.005*sin(th);
for i=1:20
    h1 = streamline(X,Y,Ex,Ey,x1(i),y1(i));
    set(h1,'LineWidth',1);
%     pipi=drawEline(h1,120,121);
%     set(pipi,'HeadWidth',4,'HeadLength',4,'Color','k')；
end
for i=[21,30]
    h1 = streamline(X,Y,Ex,Ey,x1(i),y1(i));
    set(h1,'LineWidth',1);
%     pipi=drawEline(h1,90,91);
%     set(pipi,'HeadWidth',4,'HeadLength',4,'Color','k')
end
for i=22:29
    h1 = streamline(X,Y,Ex,Ey,x1(i),y1(i));
    set(h1,'LineWidth',1);
%     pipi=drawEline(h1,70,71);
%     set(pipi,'HeadWidth',4,'HeadLength',4,'Color','k')
end    
for i=31:50
    h1 = streamline(X,Y,Ex,Ey,x1(i),y1(i));
    set(h1,'LineWidth',1);
%     pipi=drawEline(h1,120,121);
%     set(pipi,'HeadWidth',7,'HeadLength',7,'Color','b')
end
fill(x,y,'w');
for i=1:20
    h1 = streamline(X,Y,Ex,Ey,x1(i),y1(i));
    set(h1,'LineWidth',1,'linestyle','--');
    pipi=drawEline(h1,120,121);
    set(pipi,'HeadWidth',4,'HeadLength',4,'Color','k')
end
for i=[21,30]
    h1 = streamline(X,Y,Ex,Ey,x1(i),y1(i));
    set(h1,'LineWidth',1,'linestyle','--');
    pipi=drawEline(h1,90,91);
    set(pipi,'HeadWidth',4,'HeadLength',4,'Color','k')
end
for i=22:29
    h1 = streamline(X,Y,Ex,Ey,x1(i),y1(i));
    set(h1,'LineWidth',1,'linestyle','--');
    pipi=drawEline(h1,70,71);
    set(pipi,'HeadWidth',4,'HeadLength',4,'Color','k')
end    
for i=31:50
    h1 = streamline(X,Y,Ex,Ey,x1(i),y1(i));
    set(h1,'LineWidth',1,'linestyle','--');
    pipi=drawEline(h1,120,121);
    set(pipi,'HeadWidth',4,'HeadLength',4,'Color','k')
end

plot([-0.05*cos(120),-0.05*cos(120)-0.01],[-0.05*sin(120),-0.05*sin(120)],'k','linewidth',2);
plot([-0.05*cos(120)-0.01,-0.05*cos(120)-0.01],[-0.05*sin(120),-0.05*sin(120)-0.01],'k','linewidth',2);
plot([-0.05*cos(120)-0.002,-0.05*cos(120)-0.018],[-0.05*sin(120)-0.01,-0.05*sin(120)-0.01],'k','linewidth',2);

% fill(x,y,'g');
plot( b,0,'r.','MarkerSize',30);
plot( 0,0,'y.','MarkerSize',5);
plot( d,0,'k.','MarkerSize',30);

title('不带电不接地导体球外点电荷的场以及电势分布','fontsize',15);

3.1.3二维图截图：

完善前

 完善后：

3.1.4三维图代码：

k = 8.9875e9;
R = 0.05;
d = 0.15;
b = R^2 / d;
q1 = 1e-6;
q2 = -R * q1 / d;
%题目中的电荷，半径等参数



x = linspace(-0.3,0.3,60);
y = linspace(-0.16,0.16,60);
[X,Y] = meshgrid(x);
%生成网格

r1 = sqrt((X - d).^2 + Y.^2);
r2 = sqrt((X - b).^2 + Y.^2);
r3 = sqrt(X.^2 + Y.^2);
U = k * q1./r1 + k * q2./r2;
u = linspace(0,400000,101);
[Ex,Ey] = gradient(-U);
E = sqrt(Ex.^2 + Ey.^2);
Ex = Ex./E;
Ey = Ey./E;
surf(U);
hold on
shading interp;
[c,h]=contour3(U,u);
colormap(jet);
h.EdgeColor='k';
h.LineStyle='--';
h1=streamslice(Ex,Ey);
set(h1,'Color',[0.1 0.9 0.8]);
%做出电场线的三维分布图

for j=1:length(h1)
    ui = interp2(U,h1(j).XData,h1(j).YData);
    h1(j).ZData = ui;
end
%标识出图像的流线

3.1.5三维图截图：

完善前

 完善后

3.2不接地导体球外单个点电荷：

3.2.1计算公式

其中电荷量：q=1e-6  

常系数：k= 8.9875e9

3.2.2二维图代码：

k = 8.9845e9;
R = 0.05;
d = 0.15;
b = R^2 / d;
q1 = 1e-6;
q2 = -R * q1 / d;
q3 = -q2;

x = linspace(-0.3,0.3,60);
y = linspace(-0.16,0.16,60);
[X,Y] = meshgrid(x);

r1 = sqrt((X - d).^2 + Y.^2);
r2 = sqrt((X - b).^2 + Y.^2);
r3 = sqrt(X.^2 + Y.^2);
U = k * q1./r1 + k * q2./r2 + k * q3./r3;

a = 0:pi/100:2*pi;
x = R * sin( a );
y = R * cos( a );
plot ( x,y,'c-');
hold on;
plot( b,0,'r.','MarkerSize',30);
plot( 0,0,'r.','MarkerSize',30);
plot( d,0,'r.','MarkerSize',30);
xlabel('x','fontsize',16);
ylabel('y','fontsize',16);

u = linspace(0,600000,251);
[C,H] = contour(X,Y,U,u,'k-');
contour(X,Y,U,u,'c-');
% clabel(C,H);
% grid on;

[Ex,Ey] = gradient(-U);
E = sqrt(Ex.^2 + Ey.^2);
Ex = Ex./E;
Ey = Ey./E;
th=linspace(0,2*pi,50);
x1=d+0.005*cos(th);
y1=0.005*sin(th);
x2=0+0.005*cos(th);
y2=0.005*sin(th);
x3=b+0.005*cos(th);
y3=0.005*sin(th);
for i=1:20
    h1 = streamline(X,Y,Ex,Ey,x1(i),y1(i));
    set(h1,'LineWidth',1);
%     pipi=drawEline(h1,120,121);
%     set(pipi,'HeadWidth',4,'HeadLength',4,'Color','b')
end
for i=[21,30]
    h1 = streamline(X,Y,Ex,Ey,x1(i),y1(i));
    set(h1,'LineWidth',1);
%     pipi=drawEline(h1,220,221);
%     set(pipi,'HeadWidth',4,'HeadLength',4,'Color','b')
end
for i=22:29
    h1 = streamline(X,Y,Ex,Ey,x1(i),y1(i));
    set(h1,'LineWidth',1);
%     pipi=drawEline(h1,40,41);
%     set(pipi,'HeadWidth',4,'HeadLength',4,'Color','b')
end    
for i=31:50
    h1 = streamline(X,Y,Ex,Ey,x1(i),y1(i));
    set(h1,'LineWidth',1);
%     pipi=drawEline(h1,120,121);
%     set(pipi,'HeadWidth',4,'HeadLength',4,'Color','b')
end

for i=19:32
    h2 = streamline(X,Y,Ex,Ey,x2(i),y2(i));
    set(h2,'LineWidth',1);
%     pipi=drawEline(h2,165,166);
%     set(pipi,'HeadWidth',4,'HeadLength',4,'Color','b')
end

fill(x,y,'w');
for i=1:20
    h1 = streamline(X,Y,Ex,Ey,x1(i),y1(i));
    set(h1,'LineWidth',1,'linestyle','--');
    pipi=drawEline(h1,120,121);
    set(pipi,'HeadWidth',4,'HeadLength',4,'Color','k')
end
for i=[21,30]
    h1 = streamline(X,Y,Ex,Ey,x1(i),y1(i));
    set(h1,'LineWidth',1,'linestyle','--');
    pipi=drawEline(h1,220,221);
    set(pipi,'HeadWidth',4,'HeadLength',4,'Color','k')
end
for i=22:29
    h1 = streamline(X,Y,Ex,Ey,x1(i),y1(i));
    set(h1,'LineWidth',1,'linestyle','--');
    pipi=drawEline(h1,40,41);
    set(pipi,'HeadWidth',4,'HeadLength',4,'Color','k')
end    
for i=31:50
    h1 = streamline(X,Y,Ex,Ey,x1(i),y1(i));
    set(h1,'LineWidth',1,'linestyle','--');
    pipi=drawEline(h1,120,121);
    set(pipi,'HeadWidth',4,'HeadLength',4,'Color','k')
end

for i=19:32
    h2 = streamline(X,Y,Ex,Ey,x2(i),y2(i));
    set(h2,'LineWidth',1,'linestyle','--');
    pipi=drawEline(h2,165,166);
    set(pipi,'HeadWidth',4,'HeadLength',4,'Color','k')
end
a = 0:pi/100:2*pi;
x = R * sin( a );
y = R * cos( a );
plot ( x,y,'c-','LineWidth',1);
plot( b,0,'g.','MarkerSize',30);
plot( 0,0,'g.','MarkerSize',30);
plot( d,0,'r.','MarkerSize',30);
title('镜像法——不接地导体球外点电荷','fontsize',15);

   3.2.3二维图截图：

完善前

 完善后

3.2.4三维图代码：

k = 8.9875e9;
k = 8.9875e9;
R = 0.05;
d = 0.15;
b = R^2 / d;
q1 = 1e-6;
q2 = -R * q1 / d;
q3 = -q2;
%题目中的电荷，半径等参数

x = linspace(-0.3,0.3,60);
y = linspace(-0.16,0.16,60);
[X,Y] = meshgrid(x);
%生成网格

r1 = sqrt((X - d).^2 + Y.^2);
r2 = sqrt((X - b).^2 + Y.^2);
r3 = sqrt(X.^2 + Y.^2);
U = k * q1./r1 + k * q2./r2 + k * q3./r3;
%由已知值计算出电势表达式

u = linspace(-2000000,2000000,1001);
[Ex,Ey] = gradient(-U);
E = sqrt(Ex.^2 + Ey.^2);
Ex = Ex./E;
Ey = Ey./E;
surf(U);
hold on
shading interp;

subplot(2,1,1)
[c,h]=contour3(U,u);
colormap(white);
h.EdgeColor='c';
% h.LineStyle='--';

subplot(2,1,2);
[c,h]=contour3(U,u);
colormap(white);
h.EdgeColor='w';
% % h.LineStyle='--';
h1=streamslice(Ex,Ey);
set(h1,'Color',[0.1 0.9 0.8]);
%做出电场线的三维分布图

for j=1:length(h1)
    ui = interp2(U,h1(j).XData,h1(j).YData);
    h1(j).ZData = ui;
end
%标识出图像的流线

3.2.5三维图截图：

完善前

 完善后

3.3接地导体球外，直径所在直线上左右对称放置两个正点电荷：

3.3.1计算公式：

其中电荷量：q=1e-6  

常系数：k= 8.9875e9

3.3.2二维图代码：

k = 8.9875e9;

R = 0.05;

d = 0.15;

b = R^2 / d;

q1 = 1e-6;

q2 = -R * q1 / d;

q3 = q1;

q4 = -R * q1 / d;



x = linspace(-0.3,0.3,60);

y = linspace(-0.16,0.16,60);

[X,Y] = meshgrid(x);



r1 = sqrt((X - d).^2 + Y.^2);

r2 = sqrt((X - b).^2 + Y.^2);

r3 = sqrt((X + d).^2 + Y.^2);

r4 = sqrt((X + b).^2 + Y.^2);

U = k * q1./r1 + k * q2./r2 + k * q3./r3 + k * q4./r4;



a = 0:pi/100:2*pi;

x = R * sin( a );

y = R * cos( a );

plot ( x,y,'k-');

hold on;

plot( b,0,'r.','MarkerSize',30);

plot( 0,0,'r.','MarkerSize',30);

plot( d,0,'k.','MarkerSize',30);

xlabel('x','fontsize',16);

ylabel('y','fontsize',16);



u = linspace(0,400000,101);

[C,H] = contour(X,Y,U,u,'k-');

contour(X,Y,U,u,'c-');

% clabel(C,H);

% grid on;



[Ex,Ey] = gradient(-U);

E = sqrt(Ex.^2 + Ey.^2);

Ex = Ex./E;

Ey = Ey./E;



th=linspace(0,2*pi,50);

x1=d+0.005*cos(th);

y1=0.005*sin(th);

for i=1:20

    h1 = streamline(X,Y,Ex,Ey,x1(i),y1(i));

    set(h1,'LineWidth',1);

    pipi=drawEline(h1,120,121);

    set(pipi,'HeadWidth',4,'HeadLength',4,'Color','k')

end

for i=[21,30]

    h1 = streamline(X,Y,Ex,Ey,x1(i),y1(i));

    set(h1,'LineWidth',1);

    pipi=drawEline(h1,90,91);

    set(pipi,'HeadWidth',4,'HeadLength',4,'Color','k')

end

for i=22:29

    h1 = streamline(X,Y,Ex,Ey,x1(i),y1(i));

    set(h1,'LineWidth',1);

    pipi=drawEline(h1,70,71);

    set(pipi,'HeadWidth',4,'HeadLength',4,'Color','k')

end    

for i=31:50

    h1 = streamline(X,Y,Ex,Ey,x1(i),y1(i));

    set(h1,'LineWidth',1);

    pipi=drawEline(h1,120,121);

    set(pipi,'HeadWidth',4,'HeadLength',4,'Color','k')

end



x1=-d-0.005*cos(th);

y1=0.005*sin(th);

for i=1:20

    h1 = streamline(X,Y,Ex,Ey,x1(i),y1(i));

    set(h1,'LineWidth',1);

    pipi=drawEline(h1,120,121);

    set(pipi,'HeadWidth',4,'HeadLength',4,'Color','k')

end

for i=[21,30]

    h1 = streamline(X,Y,Ex,Ey,x1(i),y1(i));

    set(h1,'LineWidth',1);

    pipi=drawEline(h1,90,91);

    set(pipi,'HeadWidth',4,'HeadLength',4,'Color','k')

end

for i=22:29

    h1 = streamline(X,Y,Ex,Ey,x1(i),y1(i));

    set(h1,'LineWidth',1);

    pipi=drawEline(h1,70,71);

    set(pipi,'HeadWidth',4,'HeadLength',4,'Color','k')

end    

for i=31:50

    h1 = streamline(X,Y,Ex,Ey,x1(i),y1(i));

    set(h1,'LineWidth',1);

    pipi=drawEline(h1,120,121);

    set(pipi,'HeadWidth',4,'HeadLength',4,'Color','k')

end



plot([-0.05*cos(120),-0.05*cos(120)-0.01],[-0.05*sin(120),-0.05*sin(120)],'k','linewidth',2);

plot([-0.05*cos(120)-0.01,-0.05*cos(120)-0.01],[-0.05*sin(120),-0.05*sin(120)-0.01],'k','linewidth',2);

plot([-0.05*cos(120)-0.002,-0.05*cos(120)-0.018],[-0.05*sin(120)-0.01,-0.05*sin(120)-0.01],'k','linewidth',2);



fill(x,y,'w');

th=linspace(0,2*pi,50);

x1=d+0.005*cos(th);

y1=0.005*sin(th);

for i=[1:20]

    h1 = streamline(X,Y,Ex,Ey,x1(i),y1(i));

    set(h1,'LineWidth',1,'linestyle','--');

    pipi=drawEline(h1,120,121);

    set(pipi,'HeadWidth',4,'HeadLength',4,'Color','k')

end

for i=[21,30]

    h1 = streamline(X,Y,Ex,Ey,x1(i),y1(i));

    set(h1,'LineWidth',1,'linestyle','--');

    pipi=drawEline(h1,90,91);

    set(pipi,'HeadWidth',4,'HeadLength',4,'Color','k')

end

for i=[22:29]

    h1 = streamline(X,Y,Ex,Ey,x1(i),y1(i));

    set(h1,'LineWidth',1,'linestyle','--');

    pipi=drawEline(h1,70,71);

    set(pipi,'HeadWidth',4,'HeadLength',4,'Color','k')

end    

for i=[31:50]

    h1 = streamline(X,Y,Ex,Ey,x1(i),y1(i));

    set(h1,'LineWidth',1,'linestyle','--');

    pipi=drawEline(h1,120,121);

    set(pipi,'HeadWidth',4,'HeadLength',4,'Color','k')

end

x1=-d-0.005*cos(th);

y1=0.005*sin(th);

for i=[1:20]

    h1 = streamline(X,Y,Ex,Ey,x1(i),y1(i));

    set(h1,'LineWidth',1,'linestyle','--');

    pipi=drawEline(h1,120,121);

    set(pipi,'HeadWidth',4,'HeadLength',4,'Color','k')

end

for i=[21,30]

    h1 = streamline(X,Y,Ex,Ey,x1(i),y1(i));

    set(h1,'LineWidth',1,'linestyle','--');

    pipi=drawEline(h1,90,91);

    set(pipi,'HeadWidth',4,'HeadLength',4,'Color','k')

end

for i=[22:29]

    h1 = streamline(X,Y,Ex,Ey,x1(i),y1(i));

    set(h1,'LineWidth',1,'linestyle','--');

    pipi=drawEline(h1,70,71);

    set(pipi,'HeadWidth',4,'HeadLength',4,'Color','k')

end    

for i=[31:50]

    h1 = streamline(X,Y,Ex,Ey,x1(i),y1(i));

    set(h1,'LineWidth',1,'linestyle','--');

    pipi=drawEline(h1,120,121);

    set(pipi,'HeadWidth',4,'HeadLength',4,'Color','k')

end

plot( b,0,'r.','MarkerSize',30);

plot( -b,0,'r.','MarkerSize',30);

plot( 0,0,'y.','MarkerSize',5);

plot( d,0,'k.','MarkerSize',30);

plot( -d,0,'k.','MarkerSize',30);

title('接地导体球外，直径所在直线上左右对称放置的两个正点电荷','fontsize',15);

3.3.3二维图截图：

完善前

 完善后

3.3.4 三维图代码：

k = 8.9875e9;

R = 0.05;

d = 0.15;

b = R^2 / d;

q1 = 1e-6;

q2 = -R * q1 / d;

q3 = q1;

q4 = -R * q1 / d;

%题目中的电荷，半径等参数



x = linspace(-0.3,0.3,60);

y = linspace(-0.16,0.16,60);

[X,Y] = meshgrid(x);

%生成网格



r1 = sqrt((X - d).^2 + Y.^2);

r2 = sqrt((X - b).^2 + Y.^2);

r3 = sqrt((X + d).^2 + Y.^2);

r4 = sqrt((X + b).^2 + Y.^2);

U = k * q1./r1 + k * q2./r2 + k * q3./r3 + k * q4./r4;

%由已知值计算出电势表达式



u = linspace(-2000000,2000000,1001);

[Ex,Ey] = gradient(-U);

E = sqrt(Ex.^2 + Ey.^2);

Ex = Ex./E;

Ey = Ey./E;

surf(U);

hold on

shading interp;



subplot(2,1,1);

[c,h]=contour3(U,u);

colormap(jet);

h.EdgeColor='c';



subplot(2,1,2);

[c,h]=contour3(U,u);

colormap(white);

h.EdgeColor='w';

% % h.LineStyle='--';

h1=streamslice(Ex,Ey);

set(h1,'Color',[0.1 0.9 0.8]);

% 做出电场线的三维分布图

for j=1:length(h1)

    ui = interp2(U,h1(j).XData,h1(j).YData);

    h1(j).ZData = ui;

end

%标识出图像的流线

3.3.5 三维图截图

完善前

 完善后

3.4接地导体球外，与球心连线夹角为90度放置两个正点电荷

3.4.1计算公式：

其中电荷量：q=1e-6  

常系数：k= 8.9875e9

3.4.2二维图代码：

k = 8.9875e9;

R = 0.05;

d = 0.15;

b = R^2 / d;

q1 = 1e-6;

q2 = -R * q1 / d;

q3 = q1;

q4 = -R * q1 / d;



x = linspace(-0.3,0.3,60);

y = linspace(-0.16,0.16,60);

[X,Y] = meshgrid(x);



r1 = sqrt((X - d).^2 + Y.^2);

r2 = sqrt((X - b).^2 + Y.^2);

r3 = sqrt(X.^2 + (Y - d).^2);

r4 = sqrt(X.^2 + (Y - b).^2);

U = k * q1./r1 + k * q2./r2 + k * q3./r3 + k * q4./r4;



a = 0:pi/100:2*pi;

x = R * sin( a );

y = R * cos( a );

plot ( x,y,'k-');

hold on;

plot( b,0,'r.','MarkerSize',30);

plot( 0,0,'r.','MarkerSize',30);

plot( d,0,'k.','MarkerSize',30);

xlabel('x','fontsize',16);

ylabel('y','fontsize',16);



u = linspace(0,400000,101);

[C,H] = contour(X,Y,U,u,'k-');

contour(X,Y,U,u,'c-');

% clabel(C,H);

% grid on;



[Ex,Ey] = gradient(-U);

E = sqrt(Ex.^2 + Ey.^2);

Ex = Ex./E;

Ey = Ey./E;



th=linspace(0,2*pi,50);

x1=d+0.005*cos(th);

y1=0.005*sin(th);

for i=1:50

    h1 = streamline(X,Y,Ex,Ey,x1(i),y1(i));

    set(h1,'LineWidth',1);

    pipi=drawEline(h1,70,71);

    set(pipi,'HeadWidth',4,'HeadLength',4,'Color','k')

end  



x1=0.005*cos(th);

y1=d+0.005*sin(th);

for i=1:50

    h1 = streamline(X,Y,Ex,Ey,x1(i),y1(i));

    set(h1,'LineWidth',1);

    pipi=drawEline(h1,70,71);

    set(pipi,'HeadWidth',4,'HeadLength',4,'Color','k')

end    



fill(x,y,'w');



x1=d+0.005*cos(th);

y1=0.005*sin(th);

for i=1:50

    h1 = streamline(X,Y,Ex,Ey,x1(i),y1(i));

    set(h1,'LineWidth',1,'linestyle','--');

    pipi=drawEline(h1,70,71);

    set(pipi,'HeadWidth',4,'HeadLength',4,'Color','k')

end  



x1=0.005*cos(th);

y1=d+0.005*sin(th);

for i=1:50

    h1 = streamline(X,Y,Ex,Ey,x1(i),y1(i));

    set(h1,'LineWidth',1,'linestyle','--');

    pipi=drawEline(h1,70,71);

    set(pipi,'HeadWidth',4,'HeadLength',4,'Color','k')

end    



plot([-0.05*cos(120),-0.05*cos(120)-0.01],[-0.05*sin(120),-0.05*sin(120)],'k','linewidth',2);

plot([-0.05*cos(120)-0.01,-0.05*cos(120)-0.01],[-0.05*sin(120),-0.05*sin(120)-0.01],'k','linewidth',2);

plot([-0.05*cos(120)-0.002,-0.05*cos(120)-0.018],[-0.05*sin(120)-0.01,-0.05*sin(120)-0.01],'k','linewidth',2);





plot( b,0,'r.','MarkerSize',30);

plot( 0,b,'r.','MarkerSize',30);

plot( 0,0,'y.','MarkerSize',5);

plot( d,0,'k.','MarkerSize',30);

plot( 0,d,'k.','MarkerSize',30);

title('接地导体球外，与球心连线夹角为90度放置两个电荷','fontsize',15);

3.4.3二维图截图：

完善前

 完善后

3.4.4 三维图代码：

k = 8.9875e9;

R = 0.05;

d = 0.15;

b = R^2 / d;

q1 = 1e-6;

q2 = -R * q1 / d;

q3 = q1;

q4 = -R * q1 / d;

%题目中的电荷，半径等参数



x = linspace(-0.3,0.3,60);

y = linspace(-0.16,0.16,60);

[X,Y] = meshgrid(x);

%生成网格



r1 = sqrt((X - d).^2 + Y.^2);

r2 = sqrt((X - b).^2 + Y.^2);

r3 = sqrt(X.^2 + (Y - d).^2);

r4 = sqrt(X.^2 + (Y - b).^2);

U = k * q1./r1 + k * q2./r2 + k * q3./r3 + k * q4./r4;

%由已知值计算出电势表达式



u = linspace(-2000000,2000000,1001);

[Ex,Ey] = gradient(-U);

E = sqrt(Ex.^2 + Ey.^2);

Ex = Ex./E;

Ey = Ey./E;

surf(U);

hold on

shading interp;



subplot(2,1,1);

[c,h]=contour3(U,u);

colormap(jet);

h.EdgeColor='c';



subplot(2,1,2);

[c,h]=contour3(U,u);

colormap(white);

h.EdgeColor='w';

% % h.LineStyle='--';

h1=streamslice(Ex,Ey);

set(h1,'Color',[0.1 0.9 0.8]);

% 做出电场线的三维分布图

for j=1:length(h1)

    ui = interp2(U,h1(j).XData,h1(j).YData);

    h1(j).ZData = ui;

end

%标识出图像的流线

3.4.5 三维图截图：

完善前

 完善后

3.5同心球壳，外壳接地，球壳内部有一点电荷(p49)

3.5.1计算公式：

3.5.2二维图代码：

k = 8.9875e9;

b = 0.05;

a = 0.03;

d = 0.01;

q = 1e-6;



x = linspace(-0.06,0.06,60);

y = linspace(-0.06,0.06,60);

[X,Y] = meshgrid(x,y);



r = sqrt(X.^2 + Y.^2);

r1 = sqrt((X - d).^2 + Y.^2);

r2 = sqrt((X - a^2/d).^2 + Y.^2);

U1 = k * q * (1./r - 1/b);

U2 = k * q * (1/a - 1/b) + k * q./r1 - k * a /d * q./r2;



c = linspace(0,2*pi,50);

x1 = a * sin( c );

y1 = a * cos( c );

plot ( x1,y1,'k-','LineWidth',3);

fill(x1,y1,'w','EdgeColor','w');

hold on;

x2 = b * sin( c );

y2 = b * cos( c );

plot ( x2,y2,'k-','LineWidth',3);

xlabel('x','fontsize',16);

ylabel('y','fontsize',16);



u = linspace(120000,2000000,51);

[C2,H2] = contour(X,Y,U2,u,'k-');

contour(X,Y,U2,u,'c-');

% clabel(C2,H2);

grid on;



 [Ex,Ey] = gradient(-U2);

 E = sqrt(Ex.^2 + Ey.^2);

 Ex = Ex./E;

 Ey = Ey./E;

 th=linspace(0,2*pi,50);

 x1=d + 0.001*cos(th);

 y1=0.001*sin(th);

 for i=43:50

    h1 = streamline(X,Y,Ex,Ey,x1(i),y1(i));

    set(h1,'LineWidth',1);

    pipi=drawEline(h1,60,61);

    set(pipi,'HeadWidth',4,'HeadLength',4,'Color','k')

end

for i=1:8

    h1 = streamline(X,Y,Ex,Ey,x1(i),y1(i));

    set(h1,'LineWidth',1);

    pipi=drawEline(h1,60,61);

    set(pipi,'HeadWidth',4,'HeadLength',4,'Color','k')

end

for i=9:42

    h1 = streamline(X,Y,Ex,Ey,x1(i),y1(i));

    set(h1,'LineWidth',1);

    pipi=drawEline(h1,80,81);

    set(pipi,'HeadWidth',4,'HeadLength',4,'Color','k')

end    



R = linspace(0.03,0.05,50);

c = linspace(0,2*pi,50);

for i=1:50

    if i~=[1,50]

X1 = R(i) * cos(c);

Y1 = R(i) * sin(c);

plot(X1,Y1,'w-','LineWidth',3);

hold on;

    else

X1 = R(i) * cos(c);

Y1 = R(i) * sin(c);

plot(X1,Y1,'k-','LineWidth',3);

hold on;

    end

end



for i=1:50

X2 = 0.03 * cos(c(i));

Y2 = 0.03 * sin(c(i));

X3 = 0.05 * cos(c(i));

Y3 = 0.05 * sin(c(i));

h1=plot([X2,X3],[Y2,Y3],'b','LineWidth',1);

set(h1,'LineWidth',1);

pipi=drawEline(h1,1,2);

set(pipi,'HeadWidth',4,'HeadLength',4,'Color','b')

hold on;

end



fill([0,-0.06,-0.06,-0.04],[0.06,0.06,0,0.035],'w','EdgeColor','w');

fill([0,0.06,0.06,0.04],[0.06,0.06,0,0.035],'w','EdgeColor','w');

fill([0,-0.06,-0.06,-0.04],[-0.06,-0.06,0,-0.035],'w','EdgeColor','w');

fill([0,0.06,0.06,0.04],[-0.06,-0.06,0,-0.035],'w','EdgeColor','w');

R = linspace(0.05,0.06,50);

c = linspace(0,2*pi,50);

for i=2:50

X1 = R(i) * cos(c);

Y1 = R(i) * sin(c);

plot(X1,Y1,'w-','LineWidth',3);



hold on;

end



u = linspace(0,120000,11);

[C1,H1] = contour(X,Y,U1,u,'k-');

contour(X,Y,U1,u,'c-');

% clabel(C1,H1);

grid on;



plot( d,0,'r.','MarkerSize',30);

plot([0.05*cos(120),0.05*cos(120)+0.01],[-0.05*sin(120),-0.05*sin(120)],'k','linewidth',2);

plot([0.05*cos(120)+0.01,0.05*cos(120)+0.01],[-0.05*sin(120),-0.05*sin(120)-0.01],'k','linewidth',2);

plot([0.05*cos(120)+0.002,0.05*cos(120)+0.018],[-0.05*sin(120)-0.01,-0.05*sin(120)-0.01],'k','linewidth',2);



title('镜像法——同心球壳，外壳接地，球壳内部有一点电荷 ','fontsize',15);

3.5.3二维图截图：

完善前

 完善后

、

3.5.4 三维图代码：

k = 8.9875e9;

b = 0.05;

a = 0.03;

d = 0.01;

q = 1e-6;

%题目中的电荷，半径等参数



x = linspace(-0.06,0.06,60);

y = linspace(-0.06,0.06,60);

[X,Y] = meshgrid(x,y);

%生成网格



r = sqrt(X.^2 + Y.^2);

r1 = sqrt((X - d).^2 + Y.^2);

r2 = sqrt((X - a^2/d).^2 + Y.^2);

U1 = k * q * (1./r - 1/b);

U = k * q * (1/a - 1/b) + k * q./r1 - k * a /d * q./r2;

%由已知值计算出电势表达式



u = linspace(0,20000000,1001);

[Ex,Ey] = gradient(-U);

E = sqrt(Ex.^2 + Ey.^2);

Ex = Ex./E;

Ey = Ey./E;

surf(U);

hold on

shading interp;



subplot(2,1,1);

[c,h]=contour3(U,u);

colormap(jet);

h.EdgeColor='c';



subplot(2,1,2);

[c,h]=contour3(U,u);

colormap(white);

h.EdgeColor='w';

% % h.LineStyle='--';

h1=streamslice(Ex,Ey);

set(h1,'Color',[0.1 0.9 0.8]);

% 做出电场线的三维分布图

for j=1:length(h1)

    ui = interp2(U,h1(j).XData,h1(j).YData);

    h1(j).ZData = ui;

end

%标识出图像的流线

3.5.5三维图截图：

完善前

完善后 

鸽呜顾
关注
  • 19
    点赞
  • 踩
  • 110
    收藏
    觉得还不错? 一键收藏
  • 打赏
    打赏
  • 15
    评论
专栏目录
MATLAB仿真镜像法.pdf
10-07
MATLAB仿真镜像法.pdf
一文! 解决镜像法，电轴法在电磁场中的应用
m0_52555753的博客
04-15 1614
镜像法成立的条件在于参见我写的这一篇博客所以我们只需要保证求解区域中拉普拉斯方程/泊松方程保持不变，同时边界条件不变，那么结果就是一致的。
基于Matlab的电磁场与波“镜像电荷法”仿真——点电荷位于半球形凸起导体附近
qq_55710894的博客
04-04 2332
点电荷位于半球形凸起导体附近，下半球实现了虚线绘制
用MATLAB画点电荷的电场线
08-21
西电电磁场大作业，画电荷构成的电场线和等势面。应该是可以应对了。
点电荷分布matlab仿真,利用Matlab-模拟点电荷电场的分布的仿真 实验报告.pdf
weixin_42436686的博客
03-16 1160
利用Matlab-模拟点电荷电场的分布的仿真 实验报告利用 Matlab模拟点电荷电场的分布一、实验目的：1．熟悉单个点电荷及一对点电荷的电场分布情况；2．学会使用 Matlab 进行数值计算，并绘出相应的图形；二、实验原理：根据库伦定律：在真空中，两个静止点电荷之间的作用力与这两个电荷的电量乘积成正比，与它们之间距离的平方成反比，作用力的方向在两个电荷的连线上，两电荷同号为斥力，异号为吸力，它们...
利用Matlab模拟点电荷系的电场线和等势面.pdf
07-03
利用Matlab模拟点电荷系的电场线和等势面.pdf
优化方法-最速下降法+SUMT外点法（MATLAB实现）
10-28
优化方法中最速下降法+SUMT外点法MATLAB实现，代码均为本人实现，如有需要请自行下载，下载后请根据自身需求调整方程，欢迎讨论。
outerpointmethod_外点法matlab程序实例_
10-03
人工智能算法中外点法matlab实例，示范作用
MATLAB的梯度法，内点法，外点法，罚函数，惩罚函数，线性梯度法，源程序，按照提示输入，可直接运行
07-12
MATLAB的梯度法，内点法，外点法，罚函数，惩罚函数，线性梯度法，源程序，按照提示输入，可直接运行
MATLAB的梯度法，内点法，外点法，罚函数，惩罚函数
01-30
MATLAB的梯度法，内点法，外点法，罚函数，惩罚函数，线性梯度法，源程序，按照提示输入，可直接运行-MATLAB' s gradient method, interior point method, outside the point of law, penalty function, penalty ...
球体的 Mie 电场模拟：模拟 3D 球体散射的电场-matlab开发
05-29
该例程不仅计算米氏散射系数，还计算所有三个场，入射电场、散射电场和内部电场，这些场要么是一系列对应于等距 phi 值的笛卡尔平面，要么是球体。 此代码在一系列波数和介电常数上执行这些模拟。 莎拉·帕奇（Sarah Patch）的原创节目。 这是 Nick Walter 的并行版本，它将所需的输出写入 .h5 或 .mat 文件。 如有问题，请通过 [email protected] 联系 Nick Walter。
利用MATLAB绘制点两个电荷电场线
05-15
利用MATLAB绘制点电荷电场线，软件 MATLAB，源于课程设计题目（2021.06），此为用软件画点电荷电场线，同种以及异种电荷，单个及2个电荷。 适用人群，学生课程设计参考。 原理，以及代码内部，可见于待会儿的文章。 突破点：场点坐标与点电荷坐标求得电场力向量的角度，又可通过电场力大小公式求得其大小，通过极坐标与直角坐标的转换，得到两个点电荷分别对场点的电场力的向量坐标，将其坐标存入数组。E1与E2的向量坐标求和，就是场点总电场力的向量坐标。这样就解决了判断分量是相加还是相减的问题！ 再根据后续坐标迭代公式，利用MATLAB中ploy函数，以及调整for循环，while循环语句的位置，实现了完整的电场线的绘制。
matlab绘制等量同/异号电荷电场线与电势图
11-07
matlab绘制等量同/异号电荷电场线与电势图 （matlab能力有限，绘制结果不够美观，还请见谅）
MATLAB 外点法优化设计代码
09-01
MATLAB上实现了外点惩罚函数法对问题进行优化求解，有兴趣的有需要的可以下载学习学习
惩罚函数法的matlab实现方法
02-08
工程优化设计中常常用到的惩罚函数法的Matlab实现的源程序代码 工程优化设计中常常用到的惩罚函数法的Matlab实现的源程序代码
基于Matlab的电磁场与波“镜像电荷法”仿真——平面导体劈的点电荷镜像
qq_55710894的博客
04-04 2519
平面导体劈的镜像电荷法仿真
matlab画一个点电荷电场线,matlab画点电荷电场线
weixin_31421613的博客
04-04 5411
实验一 利用 Matlab 模拟点电荷的电场分布一、实验目的: 1.熟悉点电荷的电场分布情况; 2.学会使用 Matlab 绘图二、实验原理 MATLAB 输入命令的方式有两种,一种就是在命令窗口中直接输入简单的语 句,这种方式适应于命令比较简单、且......电场强度、电势是描述静电场属性的重要物理量, 利用等势面和电场线可 以很好的描述静电场。但是电势分布是复杂抽象的,本文利用 Matlab ...
电磁场与电磁波（7）——唯一性定理的应用：镜像法
Lmhark的博客
11-05 6530
镜像法是用一组适当配置的电荷系统来替代要求解的给定系统。根据唯一性定理，只要新添加的电荷（镜像电荷）不进入原求解区域，以保证不改变原求解区域内的电荷分布，就能保证不改变原求解区域内的电位方程，并且原区域边界处的条件保持不变就可满足要求。镜像法实际上是用在求解区域外添加的电荷来替代原有边界对系统的影响，因此，在使用镜像法时，首先把原有的边界去掉，再在求解区域外适当添加若干电荷，调整电荷数量、电荷量和位置，使得在原有边界处的条件不变。这样根据唯一性定理，新系统的解在原求解区域中就是原系统的解。
约束优化_外点罚函数法_外罚函数_约束优化_
10-04
通过外点罚函数法求解约束问题中目标函数的最优解
外点法matlab实现
最新发布
04-25
外点法（Exterior Point Method）是一种用于求解凸优化问题的数值方法。它通过将优化问题转化为等式约束的问题，并在等式约束下寻找最优解。以下是外点法在MATLAB中的实现步骤： 1. 定义目标函数和约束条件：首先，需要定义目标函数和约束条件。目标函数是需要最小化或最大化的函数，约束条件是问题的限制条件。 2. 转化为等式约束问题：将原始的不等式约束问题转化为等式约束问题。这可以通过引入松弛变量或者拉格朗日乘子来实现。 3. 构建拉格朗日函数：根据等式约束问题，构建拉格朗日函数。拉格朗日函数是目标函数和约束条件的线性组合。 4. 求解拉格朗日对偶问题：将拉格朗日函数转化为对偶问题，并使用数值方法求解对偶问题。这可以通过使用MATLAB中的优化工具箱中的函数来实现，如fmincon()。 5. 恢复原始问题的最优解：根据对偶问题的最优解，恢复原始问题的最优解。这可以通过计算原始问题的最优解与约束条件的关系来实现。

“相关推荐”对你有帮助么？

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论 15
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

当前余额3.43元 前往充值 >
需支付：10.00元
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包

打赏作者

鸽呜顾

你的鼓励将是我创作的最大动力

¥1 ¥2 ¥4 ¥6 ¥10 ¥20
扫码支付：¥1
获取中
扫码支付

您的余额不足，请更换扫码支付或充值

打赏作者

实付元
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明：

1.余额是钱包充值的虚拟货币，按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载，可以购买VIP、付费专栏及课程。

余额充值