Результат выполнения программы
При заданных начальных значениях:
Iw=100; { Круговой ток, А*виток }
R1=0.9E-2; R2=1.6E-2; { Внутренний и внешний радиусы катушки, м }
h=0.5E-2; { Толщина катушки, м }
V0=1E7; { Начальная скорость электрона, м/с }
Nz=5; Nr=7; { Число витков катушки вдоль осей z и r }
m0=1.2566371E-6; { Магнитная постоянная }
b0=-1.6E-19/9.1E-31; { Удельный заряд электрона }
Ib=Iw/(Nr*Nz); { Ток в одном витке катушки, A }
m=1; { Магнитная проницаемость среды }
r0=0.2E-2; { Начальное расстояние электрона от оси r }
z0=-3*h; zmax=3*h; { Начальное и конечное положения электрона по оси z }
dt=(zmax-z0)/(750*V0); { Шаг по времени }
path='C:\Results.dat'; { Имя файла с результатами }
k=20000; { Масштаб }
r[1]:=r0; dr[1]:=0; { Задание начального положения и скорости по r }
z[1]:=z0; dz[1]:=V0; { Задание начального положения и скорости по z }
q[1]:=0; dq[1]:=0; { Задание начального угла и угловой скорости }
получен следующий результат:
Результаты расчета, сохраненные в файле, показывают, что на всем протяжении полета электрона модуль его скорости не изменяется, что соответствует теории.
При уменьшении начальной скорости электрона получен следующий результат:
V0=2E5; { Начальная скорость электрона, м/с }
Модуль скорости электрона также не изменялся.
Литература
Сысун В.И. «Математическое моделирование объектов физической электроники»
Курсков С.Ю. Методические указания к лабораторным работам по курсу «Математическое моделирование физических объектов»