電磁場仿真實驗【matlab】帶電粒子在磁場中的受力

實驗： 電磁波的反射干涉

實驗：靜電場的邊值問題

實驗：帶電粒子在磁場中的受力

實驗：電磁波的極化

• (1)電場強度和磁場強度都不爲零；

• (2)電場強度爲零，磁場強度不爲零；

• (3)電場強度不爲零，磁場強度爲零。

• 磁聚焦現象

• 磁鏡現象

Matlab主函數部分：

%不同情況，需要更改調用的函數
[t,w]=ode23(@ifun1,[0:0.1:20],[0,0.01,0,6,0,0.01],[ ]);
plot3(w(:,1),w(:,3),w(:,5));%繪製軌跡
grid on
title('b\neq0,e\neq0')%標題每次改一下

• 電場強度和磁場強度都不爲零；

function dwdt=ifun(t,w)%b e不等於0
q=1.6e-2;
e=1;
b=1;
m=0.02;%設置條件
dwdt=[w(2);
    q.*b./m.*w(4);
    w(4);
    (q.*e./m)-(q.*b./m).*w(2);
    w(6);
    0;];%建立微分方程
end

• 電場強度爲零，磁場強度不爲零；

function dwdt=ifun1(t,w)%b e不等於0
q=1.6e-2;
e=0;
b=1;
m=0.02;%設置條件
dwdt=[w(2);
    q.*b./m.*w(4);
    w(4);
    (q.*e./m)-(q.*b./m).*w(2);
    w(6);
    0;];%建立微分方程
end

• 電場強度不爲零，磁場強度爲零。

function dwdt=ifun2(t,w)%b e不等於0
q=1.6e-2;
e=1;
b=0;
m=0.02;%設置條件
dwdt=[w(2);
    q.*b./m.*w(4);
    w(4);
    (q.*e./m)-(q.*b./m).*w(2);
    w(6);
    0;];%建立微分方程
end

• 磁聚焦現象

t=0:0.01:2*pi;
a1=0.5.*(t-pi);
for m=[-16:2:10]*pi/180
    axis([0 7 -1 1 -0.4 0.4]);
    grid on;
    view(12,18);
    hold on;
plot3(cos(m).*t,0*sin(m).*cos(a1).^2,2*sin(m).*cos(a1).*sin(a1));	plot3(cos(m).*t,2*sin(m).*cos(a1).^2,2*sin(m).*cos(a1).*sin(a1));
end
xlabel('x');ylabel('y'); zlabel('z');
title('磁聚焦')

• 磁鏡現象

q=1.6e-19;
m=1.6725e-27;
v0=1e6;sita=pi/4;
vy=0,vz=v0*cos(sita),vx=v0*sin(sita);
c=[0,1.1,0,vx,vy,vz];
[t,y]=ode45('s4',[0:1e-7:1.5e-6],c,[],q,m)
plot3(y(:,1),y(:,2),y(:,3))