Применение магнитного поля типа «магнитный кокон»
В биологии и медицине.
УДК 577.23, УДК 577.24
Медицина. {14.04.2017 03.00 МСК} Теория, требующая совершенно безопасной и относительно дешёвой экспериментальной проверки, предсказывает, что при помощи особой конфигурации такого магнитного поля, становится возможным лечить болезни суставов, кровеносную, лимфатическую и нервную системы организма за счёт глубинного эндофореза. Я, также, предполагаю, что в определённых условиях такое поле способно содействовать омоложению организма и значительному увеличению продолжительности жизни за счёт значительного продления молодости и зрелости. Известно, что живые организмы это коллоидные растворы с упорядоченным и целенаправленным движением материи (атомов и молекул, взаимодействующих при помощи электрического и магнитного полей) и информации с помощью энергетики. Об этом знают биофизики и специалисты в области информации. Так вот, при помощи такого поля, можно направленно управлять энергетикой и процессами, происходящими в организме на молекулярном и атомарном уровнях, а именно, сшивать повреждённые участки белков, РНК и ДНК и механизма, компенсирующего укорочение теломер (теломеразы).
Приложение 1
© Евгений Александрович Григорьев с помощью программы MathCad
Январь 1997 года
Впервые в науке - расчет магнитного поля внутри системы проводников,
эквивалентной по электромагнитным свойствам полому замкнутому проводнику.
Рассчитывается напряженность магнитного поля вдоль координатной оси X.
Поверхностный ток Is, текущий по полому замкнутому проводнику между его полюсами А и В, можно представить в виде "нитей тока" i [И.Е. Тамм, Основы теории электричества, стр. 140], токов, текущих по математическим меридианам сферы от B к А.
Это значительно упрощает программу расчета магнитного поля внутри системы проводников, эквивалентной по электромагнитным свойствам (в необходимом приближении) сферическому полому замкнутому проводнику (ПЗП ).
Закон Био-Савара-Лапласа
в векторной форме
Видно, что в пределах точности машинного счёта, МП внутри ПЗП равно нулю.
Приложение 2
Программа расчёта интенсивности магнитного поля внутри сферического полого замкнутого проводника вдоль координаты Х на языке программирования Quick BASIC
1995 год
'SFERA
SCREEN 12
TIMER ON
LOCATE 12, 1
COLOR 15
PRINT "Calculation of a magnetostatic field inside a Spherical Conductors System, which";
LOCATE 14, 2
PRINT "is equivalent by electromagnetic properties to a Closed and Hollow Conductor."
ON TIMER(7) GOSUB 1
DO WHILE INKEY$ = "": LOOP
1
CLS
Nl = 150 'number of lines with current in Spherical Conductor System
Nds = 50 'number of elements of current in line of SCS
DIM SG(Nds + 1), CG(Nds + 1)
C$ = " +###.##"
Pi = 3.14159265359#
I = 10 'SCS line current
IL = I * Nl 'current of Linear Conductor
Ndl = Nds 'number of elements of current in Linear Conductor
Rs = 20
L = Rs 'length of Linear Conductor strength
dl = L / Ndl
dg = 2 * Pi / Nds
ds = Rs * dg
FOR j = 0 TO Nds - 1
SG(j) = SIN(j * dg + dg / 2)
CG(j) = COS(j * dg + dg / 2)
NEXT j
LINE (500, 10)-(500, 60), 9
LOCATE 3, 64
COLOR 9
PRINT "LC"
CIRCLE (500, 110), 50, 15
LINE (500, 110)-(550, 110), 8
LOCATE 7, 66
COLOR 8
PRINT "Rs"
COLOR 5
LINE (20, 40)-(320, 40) 'X-axis
LINE (20, 105)-(320, 105) 'Y-axis
LINE (20, 165)-(320, 165) 'Z-axis
LINE (20, 230)-(320, 230) 'SCS module-axis
LINE (20, 300)-(320, 300) 'Linear current module-axis
LINE (320, 10)-(320, 350), 8, , &HFF00
LOCATE 23, 40
COLOR 8
PRINT "Rs (radius of SCS,"
LOCATE 24, 30
PRINT "along which calculates a magnetic field)"
FOR X = 0 TO Rs - 1 STEP 1
HVx = 0
HVy = 0
HVz = 0
Hs = 0
HCx = 0
HCy = 0
HCz = 0
HL = 0
FOR v = 0 TO Nds - 1
HUx = 0
HUy = 0
HUz = 0
FOR u = 0 TO Nl - 1
Rx = X - Rs * SG(v) * COS(u)
Ry = -Rs * SG(v) * SIN(u)
Rz = -Rs * CG(v)
R3 = (SQR(Rx * Rx + Ry * Ry + Rz * Rz)) ^ 3
dsx = -ds * CG(v) * COS(u)
dsy = -ds * CG(v) * SIN(u)
dsz = ds * SG(v)
Hsx = (dsy * Rz - Ry * dsz) / R3
Hsy = (dsz * Rx - Rz * dsx) / R3
Hsz = (dsx * Ry - Rx * dsy) / R3
HUx = HUx + Hsx
HUy = HUy + Hsy
HUz = HUz + Hsz
NEXT u
HVx = HVx + HUx
HVy = HVy + HUy
HVz = HVz + HUz
Hs = SQR(HVx ^ 2 + HVy ^ 2 + HVz ^ 2)
RLx = X
'RLy = 0
RLz = Rs + dl * v + dl / 2
RL3 = (SQR(RLx * RLx + RLz * RLz)) ^ 3
'dlx = 0
'dly = 0
dlz = dl
'Hlx = (dly * RLz - RLy * dlz) / RL3
Hly = (dlz * RLx - RLz * dlx) / RL3
'Hlz = (dlx * RLy - RLx * dly) / RL3
'HCx = HCx + Hlx
HCy = HCy + Hly
'HCz = HCz + Hlz
'HL = SQR(HCx ^ 2 + HCy ^ 2 + HCz ^ 2)
HL = HCy
NEXT v
LOCATE 2, 1
COLOR 14
PRINT "Hsx ="; USING C$; I * HVx
PSET (X * 300 / Rs + 20, HVx + 40)
COLOR 10
LOCATE 6, 1
PRINT "Hsy ="; USING C$; I * HVy
PSET (X * 300 / Rs + 20, HVy + 105)
COLOR 11
LOCATE 10, 1
PRINT "Hsz ="; USING C$; I * HVz
PSET (X * 300 / Rs + 20, HVz + 165)
COLOR 15
LOCATE 14, 1
PRINT "Hs ="; USING C$; I * Hs
PSET (X * 300 / Rs + 20, Hs + 230)
LOCATE 14, 25
PRINT "Magnetic intensity in Spherical Conductors System (SCS)"
'LOCATE 14, 25
'PRINT "in Spherical Conductors System (SCS)"
COLOR 9
LOCATE 18, 1
PRINT "HL ="; USING C$; HL * IL * 10
PSET (X * 300 / Rs + 20, HL * IL * 10 + 300)
LOCATE 18, 25
PRINT "Magnetic intensity from Linear Conductor current (LC)"
NEXT X
FOR X = Rs TO Rs * 2 STEP 1
HCx = 0
HCy = 0
HCz = 0
HL = 0
FOR v = 0 TO Nds - 1
RLx = X
'RLy = 0
RLz = Rs + dl * v + dl / 2
RL3 = (SQR(RLx * RLx + RLz * RLz)) ^ 3
'dlx = 0
'dly = 0
dlz = dl
'Hlx = (dly * RLz - RLy * dlz) / RL3
Hly = (dlz * RLx - RLz * dlx) / RL3
'Hlz = (dlx * RLy - RLx * dly) / RL3
'HCx = HCx + Hlx
HCy = HCy + Hly
'HCz = HCz + Hlz
'HL = SQR(HCx ^ 2 + HCy ^ 2 + HCz ^ 2)
HL = HCy
NEXT v
PSET (X * 300 / Rs + 20, HL * IL * 10 + 300)
NEXT X
STOP
COLOR 15
END
Приложение 3
Движение одной заряженной частицы в магнитном поле, создаваемом двумя другими движущимися частицами и между ними.
SCREEN 12
TIMER ON
COLOR 3
LOCATE 2, 12
PRINT " Motion of one charged corpuscle (OMC)"
LOCATE 4, 12
PRINT " in a magnetic field (MF) such as a Magnetic; Cocoon,"
LOCATE 6, 12
PRINT "which is created by two others mobile charged corpuscles (TMC)"
LOCATE 8, 25
PRINT "Etmc = 5 eV; Eomc = 50 eV"
ON TIMER(5) GOSUB 1
DO WHILE INKEY$ = "": LOOP
1
CLS
C$ = " +#.#####^^^^"
Pi = 3.141593
Mo = 4 * Pi * .0000001
K = Mo / (4 * Pi) 'SI
M = 1.67E-27 * 2 '[kg]
Q = 1.6E-19 '[C]
LOCATE 6, 2
COLOR 9
INPUT "distance between TMC [sm](10 - 50)="; A0
a = A0 * .01 '[m] distance between TMC
Xk1 = 0
Yk1 = 0
Xk2 = Xk1 + a
Yk2 = 0
CIRCLE (Xk1 * 1000 + 50, Yk1 + 240), 10, 9
CIRCLE (Xk2 * 1000 + 50, Yk2 + 240), 10, 9
Ek = 5 '[eV] energy of TMC
TVk = 0 'angle of moving of TMC
Vk = SQR(2 * Ek * ABS(Q) / M) 'velocity of TMC
Vkx = Vk * COS(TVk)
Vky = 0
Vkz = 0
LOCATE 1, 60
PRINT "V tmc ="; USING C$; Vk
COLOR 13
Eq = 50 '[eV] energy of entry of OMC
LOCATE 8, 2
INPUT "angle of entry of OMC (0 - 360)="; TVq
Vq = SQR(2 * Eq * ABS(Q) / M) 'velocity of OMC
Vqx = Vq * COS(TVq)
Vqy = Vq * SIN(TVq)
Vqz = 0
LOCATE 3, 60
PRINT "V omc ="; USING C$; Vq
dlq = .0001# 'step of moving of OMC
Xq = Xk1 + .5 * a 'coordinate began movements of OMC
LOCATE 10, 2
INPUT "Coordinate Y began movements of OMC [sm](-25/+25)="; Yq0
COLOR 8
Yq = Yq0 * .01
Zq = 0
Fx = 0
Fy = 0
Fz = 0
FOR t = 0 TO 100000
dt = dlq / Vq
Rx1 = Xq - Xk1
Rx2 = Xq - Xk2
Ry1 = Yq - Yk1
Ry2 = Yq - Yk2 'Rz1 = 0 Rz2 = 0
R1 = SQR(Rx1 * Rx1 + Ry1 * Ry1)
R2 = SQR(Rx2 * Rx2 + Ry2 * Ry2)
'B = K * Q * [V x R] / R^3 MF, created by TMC [Tl]
'Љ—1
'Bx1 = K * Q * (Vky * Rz1 - Ry1 * Vkz) / R1 ^ 3 == 0
'By1 = K * Q * (Vkz * Rx1 - Rz1 * Vkx) / R1 ^ 3 == 0
Bz1 = K * Q * (Vkx * Ry1 - Rx1 * Vky) / R1 ^ 3
'Љ—2
'Bx2 = K * Q * (Vky * Rz2 - Ry2 * Vkz) / R2 ^ 3 == 0
'By2 = K * Q * (Vkz * Rx2 - Rz2 * Vkx) / R2 ^ 3 == 0
Bz2 = K * Q * (Vkx * Ry2 - Rx2 * Vky) / R2 ^ 3
'Lorentz force which influence on OMC FL = Q*[V x B]
'Px = (Vy * (Hz1 + Hz2) - (Hy1 + Hy2) * Vz)
'Py = (Vz * (Hx1 + Hx2) - (Hz1 + Hz2) * Vx)
'Pz = (Vx * (Hy1 + Hy2) - (Hx1 + Hx2) * Vy)
Fx = Vqy * (Bz1 + Bz2)
Fy = -(Bz1 + Bz2) * Vqx
'Fz = (Vqx * (By1 + By2) - (Bx1 + Bx2) * Vqy) ==0
'F = M*dV/dt = FL = Q*[V x Bk] ==> dV=PX*dt/M
dVqx = Fx * dt / M
dVqy = Fy * dt / M
dVqz = 0 'dVz = Fz * dt / M
Vqxt = Vqx + dVqx
Vqyt = Vqy + dVqy
Vqzt = 0 'Vqzt = Vqz + dVqz
Vqt = SQR(Vqxt ^ 2 + Vqyt ^ 2) '+ Vqzt ^ 2)
kV = Vq / Vqt
'kV - because velocity of OMC must be constant
Vqx = Vqxt * kV
Vqy = Vqyt * kV
Vqz = 0 'Vqz = Vqzt * kV
V = SQR(Vqx ^ 2 + Vqy ^ 2) '+ Vqz ^ 2)
dXq = Vqx * dt
dYq = Vqy * dt
dZq = 0 'dZq = Vz * dt
Xq = Xq + dXq
Yq = Yq + dYq
Zq = 0 'Zq = Zq + dZq
PSET (Xq * 1000 + 50, Yq * 1000 + 240), 13
NEXT t
END
Приложение 4
Движение заряженной частицы в магнитном поле типа «Магнитный кокон».
SCREEN 12
TIMER ON
COLOR 3
LOCATE 3, 17
PRINT " Motion of one deuteron in a magnetic field "
LOCATE 5, 17
PRINT " such as a Magnetic Cocoon "
LOCATE 7, 17
PRINT " Ed = 500 eV"
COLOR 15
ON TIMER(3) Gosub 1
DO While INKEY$ = "": LOOP
1
CLS
Dim X1(80), X2(80)
C$ = " +#.#####^^^^"
Pi = 3.141593
Mo = 4 * Pi * .0000001
M = 1.67E-27 * 2
Q = 1.6E-19
E = 500 '[eV]
2 CLS
LOCATE 2, 10
INPUT "ANGLE (50 - 175)="; AN
If 175 < AN Goto 2
If AN < 50 Goto 2
TV = AN * Pi / 180
V = Sqr(2 * E * Abs(Q) / M)
Vx = V * Cos(TV)
Vy = V * Sin(TV)
Vz = 0
LOCATE 2, 60
PRINT "V ="; USING C$; V
dlq = .0002
'dt = dlq / V
A = .04
YL1 = 0
YL2 = 0
Xv = 0
Yv = .001
Zv = 0
L = .08
DL = .001
N = L / DL
IL = 3000000
K = Mo / (4 * Pi)
For I = 0 To N - 1
X1(I) = -I * DL - DL / 2 - A
X2(I) = I * DL + DL / 2 + A
Next I
For t = 0 To 3000
dt = dlq / V
PXx = 0
PXy = 0
PXz = 0
For I = 0 To N - 1
Rx1 = Xv - X1(I)
Rx2 = Xv - X2(I)
Ux = DL
Uy = 0
Uz = 0
Ry1 = Yv + YL1
Ry2 = Yv + YL2
Rz = Zv
R1 = Sqr(Rx1 * Rx1 + Ry1 * Ry1 + Rz * Rz)
R2 = Sqr(Rx2 * Rx2 + Ry2 * Ry2 + Rz * Rz)
Hy1 = (-Rz * Ux) / R1 ^ 3
Hz1 = (Ux * Ry1) / R1 ^ 3
Hy2 = (-Rz * Ux) / R2 ^ 3
Hz2 = (Ux * Ry2) / R2 ^ 3
Px = (Vy * (Hz1 + Hz2) - (Hy1 + Hy2) * Vz)
Py = -(Hz1 + Hz2) * Vx
Pz = Vx * (Hy1 + Hy2)
PXx = PXx + Px
PXy = PXy + Py
PXz = PXz + Pz
Next I
dVx = Q * K * IL * PXx * dt / M
dVy = Q * K * IL * PXy * dt / M
dVz = Q * K * IL * PXz * dt / M
Vxt = Vx + dVx
Vyt = Vy + dVy
Vzt = Vz + dVz
Vt = Sqr(Vxt ^ 2 + Vyt ^ 2 + Vzt ^ 2)
kV = V / Vt
'kV - the particle must have a constant value modulo of velocity
Vx = Vxt * kV
Vy = Vyt * kV
Vz = Vzt * kV
V = Sqr(Vx ^ 2 + Vy ^ 2 + Vz ^ 2)
dXv = Vx * dt
dYv = Vy * dt
dZv = Vz * dt
Xv = Xv + dXv
Yv = Yv + dYv
Zv = Zv + dZv
PSET (Xv * 10000 + 320, Yv * 10000 + 240), 12
Next t
End
Приложение 5
Движение двух заряженных частиц в магнитном поле типа «Магнитный кокон» с учётом Кулоновского взаимодействия.
SCREEN 12
TIMER ON
COLOR 3
LOCATE 3, 17
PRINT " Motion of two deuterones in a magnetic field "
LOCATE 5, 17
PRINT " such as a Magnetic Cocoon "
COLOR 15
ON TIMER(3) GOSUB 1
DO WHILE INKEY$ = "": LOOP
1
CLS
DIM XP1(150), XP2(150)
C$ = " #.###^^^^"
Pi = 3.141593
Eo = 8.85E-12
Mo = 4 * Pi * .0000001
M = 1.67E-27 * 2
Q = 1.6E-19
E = 700! '[eV]
J = 8E+08
LOCATE 3, 10
INPUT "ANGLE (30 - 170)="; AN
IF 170 < AN GOTO 1
IF AN < 30 GOTO 1
TV = AN * Pi / 180
V = SQR(2 * E * ABS(Q) / M)
Vx1 = V * COS(TV)
Vy1 = V * SIN(TV)
Vz1 = 0
Vx2 = Vx1
Vy2 = -Vy1
dlq = .0001
A = .007
Xv1 = 0
Yv1 = .0003
Xv2 = 0
Yv2 = -.0003
L = .25
DL = .01
N = L / DL
Du = .06
S = Pi * Du ^ 2 / 4
IL = S * J
K = Mo / (4 * Pi)
FOR I = 0 TO N - 1
XP1(I) = -I * DL - DL / 2 - A
XP2(I) = I * DL + DL / 2 + A
NEXT I
FOR t = 0 TO 8000
IL = S * J
K1 = 5000000!
dt = dlq / V
LOCATE 2, 65
PRINT "t ="; USING C$; t * dt
PXx1 = 0
PXy1 = 0
PXx2 = 0
PXy2 = 0
PCx1 = 0
PCy1 = 0
PCx2 = 0
PCy2 = 0
FOR I = 0 TO N - 1
Rx11 = Xv1 - XP1(I)
Rx12 = Xv1 - XP2(I)
Rx21 = Xv2 - XP1(I)
Rx22 = Xv2 - XP2(I)
Ux = DL
REM Uy = 0 Uz = 0
Ry1 = Yv1
Ry2 = Yv2
R11 = SQR(Rx11 * Rx11 + Ry1 * Ry1)
R12 = SQR(Rx12 * Rx12 + Ry1 * Ry1)
R21 = SQR(Rx21 * Rx21 + Ry2 * Ry2)
R22 = SQR(Rx22 * Rx22 + Ry2 * Ry2)
RGx1 = Xv1 - Xv2
RGy1 = Yv1 - Yv2
RG1 = SQR(RGx1 ^ 2 + RGy1 ^ 2)
RGx2 = -RGx1
RGy2 = -RGy1
RG2 = RG1
REM Hx1 = (Uy * Rz - Ry * Uz) / R1 ^ 3
REM Hy1 = (Uz * Rx1 - Rz * Ux) / R1 ^ 3
REM Hz1 = (Ux * Ry - Rx1 * Uy) / R1 ^ 3
Hz11 = (Ux * Ry1) / R11 ^ 3
Hz12 = (Ux * Ry1) / R12 ^ 3
Hz21 = (Ux * Ry2) / R21 ^ 3
Hz22 = (Ux * Ry2) / R22 ^ 3
REM Px = (Vy * (Hz1 + Hz2) - (Hy1 + Hy2) * Vz)
REM Py = (Vz * (Hx1 + Hx2) - (Hz1 + Hz2) * Vx)
REM Pz = (Vx * (Hy1 + Hy2) - (Hx1 + Hx2) * Vy)
PGx1 = RGx1 / (ABS(RG1 ^ 3)) * K1
PGy1 = RGy1 / (ABS(RG1 ^ 3)) * K1
PGx2 = RGx2 / (ABS(RG2 ^ 3)) * K1
PGy2 = RGy2 / (ABS(RG2 ^ 3)) * K1
'Coefficient "K1" imitates group of particles
Px1 = Vy1 * (Hz11 + Hz12)
Py1 = -(Hz11 + Hz12) * Vx1
Px2 = Vy2 * (Hz21 + Hz22)
Py2 = -(Hz21 + Hz22) * Vx2
PCx1 = PCx1 + PGx1
PCy1 = PCy1 + PGy1
PCx2 = PCx2 + PGx2
PCy2 = PCy2 + PGy2
PXx1 = PXx1 + Px1
PXy1 = PXy1 + Py1
PXx2 = PXx2 + Px2
PXy2 = PXy2 + Py2
NEXT I
dVx1 = Q * K * IL * PXx1 * dt / M
dVy1 = Q * K * IL * PXy1 * dt / M
dVx2 = Q * K * IL * PXx2 * dt / M
dVy2 = Q * K * IL * PXy2 * dt / M
Vxt1 = Vx1 + dVx1
Vyt1 = Vy1 + dVy1
Vt1 = SQR(Vxt1 ^ 2 + Vyt1 ^ 2)
kV1 = V / Vt1
Vxt2 = Vx2 + dVx2
Vyt2 = Vy2 + dVy2
Vt2 = SQR(Vxt2 ^ 2 + Vyt2 ^ 2)
kV2 = V / Vt2
'kV - the particle must have a constant value modulo of velocity
dVCx1 = Q ^ 2 * PCx1 * dt / M / (4 * Pi * Eo)
dVCy1 = Q ^ 2 * PCy1 * dt / M / (4 * Pi * Eo)
dVCx2 = Q ^ 2 * PCx2 * dt / M / (4 * Pi * Eo)
dVCy2 = Q ^ 2 * PCy2 * dt / M / (4 * Pi * Eo)
Vx1 = Vxt1 * kV1 + dVCx1
Vy1 = Vyt1 * kV1 + dVCy1
Vx2 = Vxt2 * kV2 + dVCx2
Vy2 = Vyt2 * kV2 + dVCy2
Vt = SQR(Vx1 ^ 2 + Vy1 ^ 2 + Vz1 ^ 2)
REM LOCATE 23, 20
REM PRINT "Vt ="; USING C$; Vt
dXv1 = Vx1 * dt
dYv1 = Vy1 * dt
dXv2 = Vx2 * dt
dYv2 = Vy2 * dt
Xv1 = Xv1 + dXv1
Yv1 = Yv1 + dYv1
Xv2 = Xv2 + dXv2
Yv2 = Yv2 + dYv2
PSET (Xv1 * 70000 + 320, Yv1 * 60000 + 240), 9
PSET (Xv2 * 70000 + 320, Yv2 * 60000 + 240), 12
IF t = 2000 THEN CLS 1
IF t = 5000 THEN CLS 1
IF t = 6500 THEN CLS 1
NEXT t
LOCATE 29, 5
INPUT "repeat under other angle - 1 "; U1
IF U1 = 1 THEN GOTO 1
IF U1 <> 1 THEN END
Приложение 6
Движение двух заряженных частиц в нарастающем магнитном поле типа «Магнитный кокон» - «Сферический пинч».
SCREEN 12
TIMER ON
COLOR 3
LOCATE 3, 12
PRINT " Motion of two deuterones"
LOCATE 5, 12
PRINT " in a increasing magnetic field such as a Magnetic Cocoon "
LOCATE 7, 12
PRINT " (Spherical pinch)"
COLOR 15
ON TIMER(4) GOSUB 1
DO WHILE INKEY$ = "": LOOP
1
CLS
DIM XP1(80), XP2(80)
C$ = " ##.###^^^^"
Pi = 3.141593
Eo = 8.85E-12
Mo = 4 * Pi * .0000001
M = 1.67E-27 * 2
Q = 1.6E-19
E = 10000 '[eV]
J = 7E+10
LOCATE 2, 10
INPUT "ANGLE (60 - 300)="; AN
IF 300 < AN GOTO 1
IF AN < 60 GOTO 1
TV = AN * Pi / 180
V = SQR(2 * E * ABS(Q) / M)
Vx1 = V * COS(TV)
Vy1 = V * SIN(TV)
Vz1 = 0
Vx2 = Vx1
Vy2 = -Vy1
dlq = .0003
A = .015
Xv1 = 0
Yv1 = .001
Xv2 = 0
Yv2 = -.001
L = .08
DL = .001
N = L / DL
Du = .01
s = Pi * Du ^ 2 / 4
IL = s * J
K = Mo / (4 * Pi)
FOR i = 0 TO N - 1
XP1(i) = -i * DL - DL / 2 - A
XP2(i) = i * DL + DL / 2 + A
'PSET ((XP1(I) * 2000 + 320 - A * 5000), 240)
'PSET ((XP2(I) * 2000 + 320 + A * 5000), 240)
NEXT i
tt = 0
FOR t = 0 TO 8000
dt = dlq / V
tt = tt + dt
LOCATE 2, 65
PRINT "t ="; USING C$; tt
LOCATE 4, 65
PRINT "I ="; USING C$; IL
PXx1 = 0
PXy1 = 0
PXx2 = 0
PXy2 = 0
PCx1 = 0
PCy1 = 0
PCx2 = 0
PCy2 = 0
FOR i = 0 TO N - 1
Rx11 = Xv1 - XP1(i)
Rx12 = Xv1 - XP2(i)
Rx21 = Xv2 - XP1(i)
Rx22 = Xv2 - XP2(i)
Ux = DL
REM Uy = 0 Uz = 0
Ry1 = Yv1
Ry2 = Yv2
R11 = SQR(Rx11 * Rx11 + Ry1 * Ry1)
R12 = SQR(Rx12 * Rx12 + Ry1 * Ry1)
R21 = SQR(Rx21 * Rx21 + Ry2 * Ry2)
R22 = SQR(Rx22 * Rx22 + Ry2 * Ry2)
RGx1 = Xv1 - Xv2
RGy1 = Yv1 - Yv2
RG1 = SQR(RGx1 ^ 2 + RGy1 ^ 2)
RGx2 = -RGx1
RGy2 = -RGy1
RG2 = RG1
REM Hx1 = (Uy * Rz - Ry * Uz) / R1 ^ 3
REM Hy1 = (Uz * Rx1 - Rz * Ux) / R1 ^ 3
REM Hz1 = (Ux * Ry - Rx1 * Uy) / R1 ^ 3
Hz11 = (Ux * Ry1) / R11 ^ 3
Hz12 = (Ux * Ry1) / R12 ^ 3
Hz21 = (Ux * Ry2) / R21 ^ 3
Hz22 = (Ux * Ry2) / R22 ^ 3
REM Px = (Vy * (Hz1 + Hz2) - (Hy1 + Hy2) * Vz)
REM Py = (Vz * (Hx1 + Hx2) - (Hz1 + Hz2) * Vx)
REM Pz = (Vx * (Hy1 + Hy2) - (Hx1 + Hx2) * Vy)
PGx1 = RGx1 / (ABS(RG1 ^ 3))
PGy1 = RGy1 / (ABS(RG1 ^ 3))
PGx2 = RGx2 / (ABS(RG2 ^ 3))
PGy2 = RGy2 / (ABS(RG2 ^ 3))
Px1 = Vy1 * (Hz11 + Hz12)
Py1 = -(Hz11 + Hz12) * Vx1
Px2 = Vy2 * (Hz21 + Hz22)
Py2 = -(Hz21 + Hz22) * Vx2
PCx1 = PCx1 + PGx1
PCy1 = PCy1 + PGy1
PCx2 = PCx2 + PGx2
PCy2 = PCy2 + PGy2
PXx1 = PXx1 + Px1
PXy1 = PXy1 + Py1
PXx2 = PXx2 + Px2
PXy2 = PXy2 + Py2
NEXT i
dVx1 = Q * K * IL * PXx1 * dt / M
dVy1 = Q * K * IL * PXy1 * dt / M
dVx2 = Q * K * IL * PXx2 * dt / M
dVy2 = Q * K * IL * PXy2 * dt / M
Vxt1 = Vx1 + dVx1
Vyt1 = Vy1 + dVy1
Vt1 = SQR(Vxt1 ^ 2 + Vyt1 ^ 2)
kV1 = V / Vt1
Vxt2 = Vx2 + dVx2
Vyt2 = Vy2 + dVy2
Vt2 = SQR(Vxt2 ^ 2 + Vyt2 ^ 2)
kV2 = V / Vt2
'kV - the particle must have a constant value modulo of velocity
kq = 3000 'Magnification of distance of a Coulomb interaction
'imitates group of particles
dVCx1 = (kq * Q) ^ 2 * PCx1 * dt / M / (4 * Pi * Eo)
dVCy1 = (kq * Q) ^ 2 * PCy1 * dt / M / (4 * Pi * Eo)
dVCx2 = (kq * Q) ^ 2 * PCx2 * dt / M / (4 * Pi * Eo)
dVCy2 = (kq * Q) ^ 2 * PCy2 * dt / M / (4 * Pi * Eo)
Vx1 = Vxt1 * kV1 + dVCx1
Vy1 = Vyt1 * kV1 + dVCy1
Vx2 = Vxt2 * kV2 + dVCx2
Vy2 = Vyt2 * kV2 + dVCy2
Vt = SQR(Vx1 ^ 2 + Vy1 ^ 2 + Vz1 ^ 2)
REM LOCATE 23, 20
REM PRINT "Vt ="; USING C$; Vt
dXv1 = Vx1 * dt
dYv1 = Vy1 * dt
dXv2 = Vx2 * dt
dYv2 = Vy2 * dt
Xv1 = Xv1 + dXv1
Yv1 = Yv1 + dYv1
Xv2 = Xv2 + dXv2
Yv2 = Yv2 + dYv2
PSET (Xv1 * 25000 + 320, Yv1 * 25000 + 240), 9
PSET (Xv2 * 25000 + 320, Yv2 * 25000 + 240), 4
Y12 = Yv1 - Yv2
IF Y12 < .0000000000002# THEN dlq = .00000000000001#
IF Y12 > .0000000000002# THEN dlq = .00000000000002#
IF Y12 > .000000000002# THEN dlq = .0000000000002#
IF Y12 > .00000000002# THEN dlq = .000000000002#
IF Y12 > .0000000002# THEN dlq = .00000000002#
IF Y12 > .000000002# THEN dlq = .0000000002#
IF Y12 > .00000002# THEN dlq = .000000002#
IF Y12 > .0000002# THEN dlq = .00000002#
IF Y12 > .000002 THEN dlq = .0000002#
IF Y12 > .00002 THEN dlq = .000004#
IF Y12 > .0002 THEN dlq = .00004#
IF Y12 > .002 THEN dlq = .0001#
IF t = 7000 THEN CLS 1
J = J + 8E+07
IL = s * J
NEXT t
LOCATE 29, 5
INPUT "Repeat with other angle - 1 "; U1
IF U1 = 1 THEN GOTO 1
IF U1 <> 1 THEN END
Приложение 7
Движение электрона в магнитном поле тора с полоидальным током – «Ондулятор».
SCREEN 12
COLOR 3
LOCATE 5, 5
PRINT " Motion of an electron in a magnetic field "
LOCATE 7, 5
PRINT " such as a Magnetic Cocoon (undulator) "
COLOR 15
C$ = " +#.#####^^^^"
Pi = 3.14159265359#
M = 9.1E-31
Q = 1.6E-19
jc = 6E+08
dt = 3E-10
Uq = 400!
Du = .001
TV = 30 * Pi / 180
K = 4 * Pi * .0000001 * (-Q) * (jc * Du ^ 2 / 4) / (4 * Pi)
V = SQR(2 * Q * Uq / M)
Vx = V * COS(TV)
Vy = V * SIN(TV)
Vz = 0
LOCATE 3, 60
PRINT "V= "; USING C$; V
A = .03
X1 = -.05
X2 = .05
dX = (X2 - X1) / 100
dG = Pi / 50
Xv = -.32
Yv = 0
Zv = 0
PSET (Xv * 1200 + 420, Yv * 4000 + 300)
DIM X(100)
FOR i = 0 TO 99
X(i) = X1 + i * dX + dX / 2
PSET (X(i) * 1200 + 420, A * 4000 + 300)
PSET (X(i) * 1200 + 420, -A * 4000 + 300)
NEXT i
PGx = 0
PGy = 0
PGz = 0
FOR t = 1 TO 100
FOR j = 0 TO 99
CG = COS(j * dG)
SG = SIN(j * dG)
PXx = 0
PXy = 0
PXz = 0
Ry = -A * CG + Yv
Rz = -A * SG + Zv
Ry2 = Ry * Ry
Rz2 = Rz * Rz
Ux = dX 'Uy = 0 Uz = 0
FOR i = 0 TO 99
Rx = -X(i) + Xv 'R = SQR(Rx * Rx + Ry2 + Rz2)
R3 = (SQR(Rx * Rx + Ry2 + Rz2)) ^ 3
Hx = (Vy * Rz - Vz * Ry) / R3
Hy = (Vz * Rx - Vx * Rz) / R3
Hz = (Vx * Ry - Vy * Rx) / R3
'Px = Uy * Hz - Uz * Hy
Py = -Ux * Hz 'Py = Uz * Hx - Ux * Hz
Pz = Ux * Hy 'Pz = Ux * Hy - Uy * Hx
'PXx = PXx + Px
PXy = PXy - Py
PXz = PXz - Pz
NEXT i
'PGx = PGx + PXx
PGy = PGy + PXy
PGz = PGz + PXz
NEXT j
'dVx = K * PGx * dt / M
dVy = K * PGy * dt / M
dVz = K * PGz * dt / M
Vxt = Vx + dVx
Vyt = Vy + dVy
Vzt = Vz + dVz
Vt = SQR(Vxt ^ 2 + Vyt ^ 2 + Vzt ^ 2)
kV = V / Vt
'kV - the particle must have a constant value modulo of velocity
Vx = Vxt * kV
Vy = Vyt * kV
Vz = Vzt * kV
V = SQR(Vx ^ 2 + Vy ^ 2 + Vz ^ 2)
dXv = Vx * dt
dYv = Vy * dt
dZv = Vz * dt
Xv = Xv + dXv
Yv = Yv + dYv
Zv = Zv + dZv
PSET (Xv * 1200 + 420, Yv * 4000 + 300), 12
NEXT t
Приложение 8
© Григорьев Евгений Александрович. 2004 год.
Дата добавления: 2018-09-20; просмотров: 365; Мы поможем в написании вашей работы! |
Мы поможем в написании ваших работ!