* Данный текст распознан в автоматическом режиме, поэтому может содержать ошибки
1148
Аналитическое и спектральное моделирование
Рис. 15.23. Баллистические траектории камня при массе 100 г
Сила Лоренца, действующая на движущуюся частицу: F=q*(E+[v,B]). Проек ции векторного произведения [v,B] на оси х, у, z:
[v,B]x=vy*Bz-vz*By [v,B]y=vz*Bx-vx*Bz [v,B]z=vx*By-vy*Bz
В соответствии с этим известные из курса физики дифференциальные уравне ния, описывающие траекторию полета частицы по осям х,у,z, имеют вид:
> sys:=diff(x(t),t$2)=q*(Ex+(diff(y(t),t)*Bz-diff(z(t),t)*By))/ massa,diff(y(t),t$2)=q*(Ey+(diff(z(t),t)*Bx – diff(x(t),t)*Bz))/ massa,diff(z(t),t$2)= q*(Ez+(diff(x(t),t)* By-diff(y(t),t)*Bx))/ massa;
Зададим исходные числовые данные (опустив размерности):
> q:=-1.6e-19:massa:=9.1e-31:V:=1e7:alpha:=80*Pi/180: > Vx:=V*cos(alpha):Vy:=V*sin(alpha):Ex:=0:Ey:=0:Ez:=0: Bx:=0.1:By:=0: Bz:=0:
