лр2
.doc
МИНОБРНАУКИ РОССИИ
Санкт-Петербургский государственный
электротехнический университет
«ЛЭТИ» им. В.И. Ульянова (Ленина)
Кафедра РТ
отчет
по лабораторной работе №2
по дисциплине «ПСМВЭКБ »
Преподаватель |
|
Синев А.Е. |
Санкт-Петербург
2021
Цель работы.
Составить программу для расчета траекторий в электрическом поле. Составить программу для расчета траектории движения одиночного заряда вблизи “точечного магнитного заряда”
m0=9.1*10^(-31);
q=1.6*10^(-19);
Bz=2;
Ex=2;
Ez=5;
y0=10;
x0=10;
z0=1;
Y= zeros(50);
X= zeros(50);
Z= zeros(50);
A=(m0/(q*Bz*Bz)*sqrt(x0*x0+(Ex/Bz+y0)^2));
fi=atan(x0*Bz/(Ex+y0*Bz));
%Xc=(mo/(q*Bz*Bz))*(Ex+y0*Bz);
for t=1:100
X(t)=-A*cos((q*Bz/m0)*t+fi)+m0/(q*Bz*Bz)*(Ex+y0*Bz);
Y(t)=A*sin(q*Bz*t/m0+fi)-(Ex*t)/Bz-m0*x0/(q*Bz);
%Xc(t)=m0/(q*Bz*Bz)+(q*Bz^2)*x0;
%Yc(t)=-(Ex*t/Bz+(q*Bz^2)*x0);
Z(t)=(q*Ez/2*m0)*t^2+z0*Bz*t;
end
plot3(X,Y,Z);
Вывод программы
Рис.Траеттория частицы в магнитном и электрическом поле.
Программа для расчета траектории движения одиночного заряда вблизи “точечного магнитного заряда”
m0=9.1*10^(-31);
e=1.6*10^(-19);
v=10^6;
u=4*10^-7*pi;
fi=0.01;
teta0=50;
teta=0:0.1:36000;
rmin=e*u*tan(deg2rad(fi))/(m0*v);
r=rmin./(cos((deg2rad(teta-teta0))*sin(deg2rad(fi))));
x=r*cos(deg2rad(fi)).*sin(deg2rad(teta));
y=r*sin(deg2rad(fi)).*cos(deg2rad(teta));
z=r.*cos(deg2rad(fi));
plot3(x,y,z);
Вывод программы
Рис .4 Траектория движения одиночного зарадя
Вывод: В данной работе были расчитаны траекткория в магнитном и электрическом поле и траектория возле “точечного магнитного заряда”. Во второй задаче не хватает входных данных. Нужно задать или rmin, или R, или угол ν. В программе угол ν практически подобран