基于Matlab的电偶极子辐射仿真

✅作者简介：热爱科研的Matlab仿真开发者，修心和技术同步精进，matlab项目合作可私信。

🍎个人主页：Matlab科研工作室

🍊个人信条：格物致知。

更多Matlab仿真内容点击👇

智能优化算法       神经网络预测       雷达通信       无线传感器        电力系统

信号处理              图像处理               路径规划       元胞自动机        无人机 

⛄ 内容介绍

采用了镜像法等方法对电偶极子及其产生的静电场,电磁场及辐射等做了较系统和深入的分析,研究,使分析方便,简化,推出的结论有一定实际指导意义.

⛄ 代码

%电偶极子的仿真

%程序仿真里面，绘制了对称振子的二维和三维方向图，同时以动画的形式展示了电磁波的辐射过程

clear;clc;close;

​theta0=linspace(0,2*pi,100);

phi0=linspace(0,2*pi,200);

[theta,phi]=meshgrid(theta0,phi0);

lambda=5;%波长

k=2*pi/lambda;%波数k

w=pi/24;%角频率w

Pattern_2D=sin(theta0);%用来绘制二维极坐标

Pattern_3D=sin(theta);%用来绘制三维方向图

figure(1);% 1.绘制3维立体方向图

[x,y,z]=sph2cart(phi,0.5*pi-theta,Pattern_3D);%球坐标形式3维方向图

mesh(x,y,z);

title('电偶极子3维方向图仿真');

xlabel('x');ylabel('y');zlabel('z');

figure(2)%2.绘制二维极坐标方向图

polarplot(theta0,abs(Pattern_2D));%极坐标形式二维方向图

title('电偶极子2维方向图仿真');

pax=gca;

pax.ThetaDir='clockwise';

pax.ThetaZeroLocation='top';

hold on;

rmax=10;

X=rmax*sin(theta0)'*cos(phi0);

Y=rmax*sin(theta0)'*sin(phi0);

t=0;%时间

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

fmat=moviein(48);

figure(3)%3.动态演示电磁波的辐射动态过程

title('电偶极子辐射动态仿真图');

for t=1:1:48

T=(X./R).^2.*(cos(w*t-k*R)-k*R.*sin(w*t-k.*R))./(k*R);

a=-1:0.4:1;

h= polar([0 2*pi], [0 10]); %产生在极坐标下的一条直线  

delete(h);   %删除上述直线，但留下坐标轴   

hold on

contour(X,Y,T,a,'b'); 

hold on

yy=cos(w*t);

plot(0,yy,'or','markerfacecolor','r','markersize',10);

text(-0.2,yy,'+','color','w');

plot(0,-yy,'ob','markerfacecolor','b','markersize',10);%电荷展示的设置

text(-0.1,-yy,'-','color','w');

hold off

fmat(t)=getframe(gcf);

end

movie(gcf,fmat,10,24);

⛄ 运行结果

⛄ 参考文献

[1]吕宽州, 姜俊. 电偶极子的场及辐射[J]. 河南科学, 2003, 21(5):3.

⛳️ 代码获取关注我

❤️部分理论引用网络文献，若有侵权联系博主删除

❤️ 关注我领取海量matlab电子书和数学建模资料