基于Matlab的电磁场与波“电轴法”仿真

已于 2022-04-05 18:55:44 修改
阅读量2.6k 收藏 58
点赞数 12
分类专栏： 仿真 文章标签： 电磁学 matlab
于 2022-04-04 14:45:46 首次发布
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/qq_55710894/article/details/123953451
版权

0、代码中用到的drawEline函数

在使用代码前记得把这个函数和仿真文件放到一个文件夹下哦

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);

1.目的：

  了解电轴法的等效原理，学会用MATLAB绘制电轴法的电场线与等势面。

2.理论分析：

 

3.源代码与结果

3.1两根相同半径轴线（电压已知，tao未知）：

3.1.1二维图源代码：

a = 3;b = 4;d = 5;u0 = 10;
e0 = 10^(-9)/(36*pi);
x = linspace(-10,10,500);
y = linspace(-10,10,500);
theta=linspace(0,2*pi,50);
xr1=b+a*cos(theta);
yr1=a*sin(theta);
xr2=-b+a*cos(theta);
yr2=a*sin(theta);
[xx, yy]=meshgrid(x, y);
tao = pi*e0*u0/(log((b+d-a)/(b-d+a)));
fai = tao*log(sqrt((xx-b).^2+yy.^2)./sqrt((xx+b).^2+yy.^2))/(2*pi*e0);
[Ex,Ey] = gradient(fai);
E = sqrt(Ex.^2+Ey.^2);
Ex = Ex./E;Ey = Ey./E;
cv=linspace(min(min(fai)),max(max(fai)),51);
contour(xx,yy,fai,cv,'r-');
title('\fontname{ ËÎÌå }\fontsize{11} 两半径相同电轴'),hold on
plot(xr1,yr1,'k','linewidth',2)
plot(xr2,yr2,'k','linewidth',2)
xlabel('x');ylabel('y'),hold off
th=linspace(0,2*pi,17);
x1=-b+0.01*cos(th);
y1=0.01*sin(th);
x2=b+0.01*cos(th);
y2=0.01*sin(th);
for i=1:17
    h1=streamline(xx,yy,-Ex,-Ey,x1(i),y1(i));  
        switch i
            case {1,2,9,16}
                drawEline(h1,1000,1001);
            case {3,15}
                drawEline(h1,1030,1031);
            case {4,14}
                drawEline(h1,1150,1151);
            case {5,13}
                drawEline(h1,1350,1351);                
            case {6,12}
                drawEline(h1,1650,1651);
            case {7,11}
                drawEline(h1,2700,2701);
            case {8,10}
                drawEline(h1,930,931);
            otherwise
                continue
        end
end
for i=[1,2,16]
    h1=streamline(xx,yy,Ex,Ey,x2(i),y2(i)); 
    switch i
        case {2,16}
            drawEline(h1,930,929);
        case 1
            drawEline(h1,1000,999);
        otherwise
            continue
    end
end


fill(xr2,yr2,'w')
fill(xr1,yr1,'w')
for i=1:17
    h1=streamline(xx,yy,-Ex,-Ey,x1(i),y1(i));  
        set(h1,'linestyle','--');
        switch i
            case {1,2,9,16}
                drawEline(h1,1000,1001);
            case {3,15}
                drawEline(h1,1030,1031);
            case {4,14}
                drawEline(h1,1150,1151);
            case {5,13}
                drawEline(h1,1350,1351);                
            case {6,12}
                drawEline(h1,1650,1651);
            case {7,11}
                drawEline(h1,2700,2701);
            case {8,10}
                drawEline(h1,930,931);
            otherwise
                continue
        end
end
for i=[1,2,16]
    h1=streamline(xx,yy,Ex,Ey,x2(i),y2(i)); 
    set(h1,'linestyle','--');
    switch i
        case {2,16}
            drawEline(h1,930,929);
        case 1
            drawEline(h1,1000,999);
        otherwise
            continue
    end
end

3.1.2二维图结果：

修改前

 修改后

 3.1.3三维图源代码：

x = linspace(-10,10,500);
y = linspace(-10,10,500);
[xx, yy]=meshgrid(x, y);
a = 3;b = 4;d = 5;u0 = 10;
e0 = 10^(-9)/(36*pi);
tao = pi*e0*u0/(log((b+d-a)/(b-d+a)));
fai = tao*log(sqrt((xx-b).^2+yy.^2)./sqrt((xx+b).^2+yy.^2))/(2*pi*e0);
[Ex,Ey] = gradient(fai);
E = sqrt(Ex.^2+Ey.^2);
Ex = Ex./E;Ey = Ey./E;
fai0 = tao*log(0.1/7.9)/(2*pi*e0);
F = linspace(-fai0,fai0,40);
surf(fai);
hold on
shading interp;

subplot(2,1,1);
[c,h]=contour3(fai,F);
h1=streamslice(Ex,Ey);
set(h1,'Color',[0.1 0.9 0.8]);

subplot(2,1,2);
for j=1:length(h1)
    ui = interp2(fai,h1(j).XData,h1(j).YData);
    h1(j).ZData = ui;
end
[c,h]=contour3(fai,F);
h.EdgeColor='k';
h.LineStyle='--';
colormap(white);

3.1.4三维图结果

修改前

 修改后

3.2两根半径不同轴线（非嵌套，tao已知）：

3.2.1计算公式：

 由以下关系式

 可得

 3.2.2二维图源代码：

a1 = 2.5;a2 = 4;d = 7;tao=2.5484e-10;
h1=(d^2+a1^2-a2^2)/(2*d);
h2=(d^2-a1^2+a2^2)/(2*d);
b=sqrt(h1^2-a1^2);
e0 = 10^(-9)/(36*pi);
x = linspace(-10,10,500);
y = linspace(-10,10,500);
theta=linspace(0,2*pi,50);
xr1=h2+a2*cos(theta);
yr1=a2*sin(theta);
xr2=-h1+a1*cos(theta);
yr2=a1*sin(theta);
[xx, yy]=meshgrid(x, y);
fai = tao*log(sqrt((xx-b).^2+yy.^2)./sqrt((xx+b).^2+yy.^2))/(2*pi*e0);
[Ex,Ey] = gradient(fai);
E = sqrt(Ex.^2+Ey.^2);
Ex = Ex./E;Ey = Ey./E;
cv=linspace(min(min(fai)),max(max(fai)),51);
contour(xx,yy,fai,cv,'r-');
title('\fontname{ ËÎÌå }\fontsize{11} 两根半径不同轴线（非嵌套） '),hold on
plot(xr1,yr1,'k','linewidth',2)
plot(xr2,yr2,'k','linewidth',2)
xlabel('x');ylabel('y'),hold off
th=linspace(0,2*pi,51);
x1=-b+0.01*cos(th);
y1=0.01*sin(th);
x2=b+0.01*cos(th);
y2=0.01*sin(th);
for i=[1,5,9,13,17,21,22,23,29,30,31,35,39,43,47]
    h1=streamline(xx,yy,-Ex,-Ey,x1(i),y1(i)); 
        switch i
            case 1
                drawEline(h1,350,351);
            case {5,47}
                drawEline(h1,360,361);
            case {9,43}
                drawEline(h1,388,389);                
            case {13,39}
                drawEline(h1,445,446);
            case {17,35}
                drawEline(h1,530,531);
            case {21,31}
                drawEline(h1,940,941);    
            case {22,30}
                drawEline(h1,1270,1271);
            case {23,29}
                drawEline(h1,1930,1931);
            otherwise
                continue
        end
end
 
for i=[24,25,26,27,28]
    h1=streamline(xx,yy,-Ex,-Ey,x1(i),y1(i));
    drawEline(h1,1250,1251); 
end
 
for i=[1,2,3,50,49]
    h1=streamline(xx,yy,Ex,Ey,x2(i),y2(i));
    drawEline(h1,1250,1249); 
end

3.2.3二维图结果：

（实在肝不动了，就不抠细节了）

 3.2.4三维图源代码：

a1 = 2.5;a2 = 4;d = 7;tao=2.5484e-10;
h1=(d^2+a1^2-a2^2)/(2*d);
h2=(d^2-a1^2+a2^2)/(2*d);
b=sqrt(h1^2-a1^2);
e0 = 10^(-9)/(36*pi);
x = linspace(-10,10,500);
y = linspace(-10,10,500);
[xx, yy]=meshgrid(x, y);
fai = tao*log(sqrt((xx-b).^2+yy.^2)./sqrt((xx+b).^2+yy.^2))/(2*pi*e0);
fai0 = tao*log(0.05/2.4878)/(2*pi*e0);
F = linspace(-fai0,fai0,60);
[Ex,Ey] = gradient(fai);
E = sqrt(Ex.^2+Ey.^2);
Ex = Ex./E;Ey = Ey./E;
surf(fai);
hold on
shading interp;
h1=streamslice(Ex,Ey);
set(h1,'Color',[0.1 0.9 0.8]);
for j=1:length(h1)
    ui = interp2(fai,h1(j).XData,h1(j).YData);
    h1(j).ZData = ui;
end
[c,h]=contour3(fai,F);
h.EdgeColor='k';
h.LineStyle='--';
colormap(jet);

3.2.5三维图结果：

3.3两根半径相同轴线（嵌套，tao已知）：

3.3.1计算公式：

 由以下关系式

可得

 3.3.2二维图源代码

a1 = 1.5;a2 = 4;d = 2;tao=2.5484e-10;
H1=(a2^2-a1^2-d^2)/(2*d);
H2=(a2^2-a1^2+d^2)/(2*d);
b=sqrt(H1^2-a1^2);
e0 = 10^(-9)/(36*pi);
x = linspace(-10,10,500);
y = linspace(-10,10,500);
theta=linspace(0,2*pi,50);
xr1=-H2+a2*cos(theta);
yr1=a2*sin(theta);
xr2=-H1+a1*cos(theta);
yr2=a1*sin(theta);
[xx, yy]=meshgrid(x, y);
fai = tao*log(sqrt((xx-b).^2+yy.^2)./sqrt((xx+b).^2+yy.^2))/(2*pi*e0);
[Ex,Ey] = gradient(fai);
E = sqrt(Ex.^2+Ey.^2);
Ex = Ex./E;Ey = Ey./E;
cv=linspace(min(min(fai)),max(max(fai)),51);
contour(xx,yy,fai,cv,'r-');
title('\fontname{ ËÎÌå }\fontsize{11} 两根半径不同轴线（嵌套）'),hold on
plot(xr1,yr1,'k','linewidth',2)
plot(xr2,yr2,'k','linewidth',2)
xlabel('x');ylabel('y'),hold off
th=linspace(0,2*pi,51);
x1=-b+0.01*cos(th);
y1=0.01*sin(th);
x2=b+0.01*cos(th);
y2=0.01*sin(th);
for i=[1,5,9,13,17,21,22,23,29,30,31,35,39,43,47]
    h1=streamline(xx,yy,-Ex,-Ey,x1(i),y1(i)); 
        switch i
            case {1,5,47}
                drawEline(h1,550,551);
            case {9,43}
                drawEline(h1,570,571);                
            case {13,39}
                drawEline(h1,680,681);
            case {17,35}
                drawEline(h1,950,951);
            case {21,31}
                drawEline(h1,1950,1951);    
            case {22,30}
                drawEline(h1,2570,2571);
            case {23,29}
                drawEline(h1,1230,1231);
            otherwise
                continue
        end
end
 
for i=[24,25,26,27,28]
    h1=streamline(xx,yy,-Ex,-Ey,x1(i),y1(i));
    drawEline(h1,1250,1251); 
end
 
for i=[1,2,3,4,48,49,50]
    h1=streamline(xx,yy,Ex,Ey,x2(i),y2(i));
    drawEline(h1,1250,1249); 
end

3.3.3二维图结果

 3.3.4三维图源代码

a1 = 1.5;a2 = 4;d = 2;tao=2.5484e-10;
H1=(a2^2-a1^2-d^2)/(2*d);
H2=(a2^2-a1^2+d^2)/(2*d);
b=sqrt(H1^2-a1^2);
e0 = 10^(-9)/(36*pi);
x = linspace(-10,10,500);
y = linspace(-10,10,500);
[xx, yy]=meshgrid(x, y);
fai = tao*log(sqrt((xx-b).^2+yy.^2)./sqrt((xx+b).^2+yy.^2))/(2*pi*e0);
fai0 = tao*log(0.1/3.7426)/(2*pi*e0);
F = linspace(-fai0,fai0,60);
[Ex,Ey] = gradient(fai);
E = sqrt(Ex.^2+Ey.^2);
Ex = Ex./E;Ey = Ey./E;
surf(fai);
hold on
shading interp;
h1=streamslice(Ex,Ey);
set(h1,'Color',[0.1 0.9 0.8]);
for j=1:length(h1)
    ui = interp2(fai,h1(j).XData,h1(j).YData);
    h1(j).ZData = ui;
end
[c,h]=contour3(fai,F);
h.EdgeColor='k';
h.LineStyle='--';
colormap(jet);

3.3.5三维图结果

鸽呜顾
关注
  • 12
    点赞
  • 踩
  • 58
    收藏
    觉得还不错? 一键收藏
  • 打赏
    打赏
  • 0
    评论
专栏目录
Matlab电磁场仿真
Green_hand00的博客
05-06 1603
matlab画图需要用到的常用函数：quiver( ), surf( ), gradient( ), 当然还有很重要的对meshgrid( )的理解
心电图导联及心电轴.doc
10-12
心电图导联及心电轴.doc
车床电主轴的发展历史现状以与趋势与主电轴工作原理毕业论文.doc
09-16
车床电主轴的发展历史现状以与趋势与主电轴工作原理毕业论文.doc
电磁场matlab仿真
06-24
对终端不同的三种驻波状态仿真，主要为短路、开路、负载为电感三种情况的电流仿真。
电磁场实验matlab仿真
11-11
．根据库伦定律，利用Matlab强大的绘图功能画出单个点电荷的电场分布情况，包括电力线和等势面。 2．根据库伦定律,利用Matlab强大的绘图功能画出一对点电荷的电场分布情况,包括电力线的分布和等势面。 3．实验内容1中，可以在正电荷和负电荷中任选一组画出其电场分布，实验内容2中，可以在一对正电荷，一对负电荷和一正一负一对电荷中选择一组画出其电场分布情况。
电磁场 电位 导体和电介质
04-04
电场强度 电介质 电轴法 电容 边值问题和唯一性
基于Matlab的电磁场与波“带电圆棒电场电势”仿真
qq_55710894的博客
04-04 5794
带电圆棒电场电势仿真
【全套完结】电磁场与电磁波实验-----全套Matlab仿真实验
热门推荐
胡毛毛的博客
05-13 2万+
目录实验一实验二更多相关文章点这里哦 实验一 用Matlab仿真点电荷电场强度和电位（势）求点电荷Q的电势强度和电势，点电荷的电场线和等势线是什么形状？ 实验二 更多相关文章点这里哦 通信工程(信息类，电子类，电气工程，自动化，计算机，软件工程，机电，等相关专业)全套学习指导 ============================================================================ 注： 创作不易，您的点赞、收藏加关注就是博主最大的更新动力 后续会对部分文章V
基于Matlab的电磁场与波“电磁波的入射、反射、透射”仿真—（可实现波形入射反射过程的仿真）
qq_55710894的博客
04-02 1万+
写在前面，本仿真与其他大多数仿真的改进点在于，本仿真实现了对于波形运动情况的仿真，并将大部分情况融合成为一个函数文件，可以通过改变参数实现不同的仿真 一、正入射的基础驻波实现 1.原理： 通过Matlab中的相关函数公式，仿真课本上电磁波的正入射过程中驻波的形成。代码中主要通过for循环对函数进行绘制 此代码与类似代码实现的效果并无不同 2.代码 w = 10; %角频率,可以改变周期和震荡速录 x = -3:0
电磁场实验 MATLAB仿真直导线的磁场分布
12-14
电磁场实验 直导线的磁场分布 MATLAB仿真 电磁场实验 直导线的磁场分布 MATLAB仿真 电磁场实验 直导线的磁场分布 MATLAB仿真
镜像法与电轴法静电场.ppt
01-26
镜像法与电轴法静电场.ppt
心电轴测定的新方法.pdf
07-08
论文参考资料
直导线的磁场分布matlab实验报告_磁场专题知识点汇总
weixin_32821533的博客
01-09 3803
磁现象磁性磁铁能吸引铁、钴、镍等物质的性质(吸铁性)磁体具有磁性的物质天然磁体、人造磁体磁极磁体上磁性最强的部分叫磁极。(磁体两端最强中间最弱)种类：水平面自由转动的磁体，指南的磁极叫南极(S)，指北的磁极叫北极(N)磁化使原来没有磁性的物体获得磁性的过程。磁铁之所以吸引铁钉是因为铁钉被磁化后，铁钉与磁铁的接触部分间形成 异名磁极，异名磁极相互吸引的结果。钢和软铁的磁化软铁被磁...
matlab模拟线圈电磁场,利用MATLAB的PDE工具箱对电场和磁场进行模拟
weixin_30694389的博客
03-18 898
中学物理Vol. 32 No.巧 疑的态度，认为弊大于利，笔者认为这些想法不无道理，只是大家的关注点应该在于如何更好发挥它的优势，而避免它的不足．在实践的过程中笔者有以下心得和体会，与大家分享： 3，1明确运用传感器实验是一种教学辅助手段 因为现在运用传感器实验还存在着种种不足之处：如器材新颖，学生好奇心强，难以集中学生的注意力；学生对传感器的原理不熟悉，学生对数据的由来存在疑虑．所以不能以为采用...
西安电子科技大学电磁场实验matlab实现
m0_56162460的博客
12-09 851
在非均匀磁场中，粒子发生旋转的半 径与磁感应强度成反比，图1示带电粒子在递增磁场的运动轨迹[3]。当带电粒子从弱磁场区 向强磁场区运动时，如图2恒有分力F1阻止带电粒子向强磁场区的运动，于是当带电粒子运动一段时间后，粒子速度减小至零，粒子像被反弹一样开始反向运动，这种现象被称为磁镜现象。（事实上上述的左下公式针对西电的要求进行了略微改变，也就是变量名和顺寻不一样，因为西电模板给的公式顺序就不一样，原理还是一样的，当然，这篇文章仿真用的公式没有写在文中，而是写在了仿真的代码里面）
【磁场】基于Matlab实现磁场正演
qq_59747472的博客
10-07 805
本文主要是利用MATLAB软件对磁性体的磁场进行正演模拟,了解磁性体的一般分布规律以及影响磁性体磁场的主要因素(如物性参数,走向等),掌握磁性体磁场正演的计算方法,为后续更为复杂的研究奠定一定的基础.
基于MATLAB的电力线路参数计算仿真
熟悉MATLAB/fpga/RK3588/LABVIEW/protus/单片机/dsp/Multisim
10-28 970
电力系统图由发电、变电、输电、配电和用电等环节组成的电能生产与消费系统。它的功能是将自然界的一次能源通过发电动力装置（主要包括锅炉、汽轮机、发电机及电厂辅助生产系统等）转化成电能，再经输、变电系统及配电系统将电能供应到各负荷中心，通过各种设备再转换成动力、热、光等不同形式的能量，为地区经济和人民生活服务。其中电力线路，主要分为输电线路和配电线路。输电线路一般电压等级较高，磁场强度大，击穿空气（电弧）距离长。如何计算电力线路的一些重要参数显得格外重要。
Matlab进阶绘图第55期—带填充纹理的横向单组柱状图
最新发布
qq_26447137的博客
05-20 64
Matlab进阶绘图第55期—带填充纹理的横向单组柱状图
MAST30013
LUCK_520的博客
05-14 926
MAST30013 – Techniques in Operations Research Semester 1, 2024 Group Project Submissions:Project reports are to be submitted before the (strict) due date: 11:59pm Friday of Week 11 (17 May, 2024). A complete submission will consist of aLaTex typeset r
matlab电磁场仿真作业
01-17
MATLAB是一种功能强大的数学软件和编程语言，可以用于电磁场的仿真。对于电磁场仿真作业，我会使用MATLAB来完成以下任务。 首先，我会使用MATLAB的数值计算能力来解决电磁场的方程。例如，麦克斯韦方程可以使用有限差分法或有限元法进行离散化，然后使用MATLAB进行数值求解。我会编写一个自定义的MATLAB函数来实现这些数值求解算法。 其次，我会使用MATLAB的绘图功能来可视化电磁场的分布。通过将电磁场方程的解表示为二维或三维图像，我可以帮助学生更好地理解电磁场的行为。我会使用MATLAB的绘图函数来绘制电磁场的等值线图、矢量图和流线图等。 此外，我还会使用MATLAB的数据分析和统计工具来分析电磁场的性能。例如，我可以计算电磁场的能量密度、功率流密度和辐射模式等。我还可以使用MATLAB的符号计算工具来分析电磁场的解析解。这些分析结果可以帮助学生进一步了解电磁场的相关性质。 最后，我会编写一个MATLAB脚本来自动化电磁场仿真的过程。通过使用MATLAB的脚本语言，我可以将多个任务组织在一起，并使用循环和条件语句来处理大量的仿真数据。这样可以节省学生的时间并提高仿真的效率。 综上所述，我会利用MATLAB的数值计算、绘图、数据分析和脚本编程功能来完成电磁场仿真作业。

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

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

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

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

打赏作者

鸽呜顾

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

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

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

打赏作者

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

抵扣说明：

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

余额充值