Этот пример находит среднюю степень двойки излучения, привлекающую заряды, перемещающиеся в эллиптическую орбиту (электрический диполь).
Два противоположных заряда, e1
и e2
, сформируйте электрический диполь. Массами заряженных частиц является m1
и m2
, соответственно. Для общего центра массового m1*r1 + m2*r2 = 0
, где r1
и r2
векторы расстояния к заряженным частицам. Расстоянием между заряженными частицами является r = r1 - r2
.
syms m1 m2 e1 e2 r1 r2 r [r1,r2] = solve(m1*r1 + m2*r2 == 0, r == r1 - r2, r1, r2)
r1 =
r2 =
Найдите дипольный момент этой системы:
d = e1*r1 + e2*r2; simplify(d)
ans =
Согласно формуле Larmor, общая степень, излученная в модуле времени, , или, в терминах расстояния между заряженными частицами, . Здесь точка означает производную времени. Закон Кулона позволяет вам найти значения ускорения в терминах уменьшенной массы системы, , и продукт зарядов частиц, .
alpha = sym('alpha'); syms m c m = m1*m2/(m1 + m2); r2 = -alpha/(m*r^2); J = simplify(subs(2/(3*c^3)*d^2, r, r2))
J =
Главная полуось a и эксцентриситет из эллиптической орбиты даны следующими выражениями, где E
общая орбитальная энергия, и угловой момент.
syms E L phi a = alpha/(2*E)
a =
eccentricity = sqrt(1-2*E*L^2/(m*alpha^2))
eccentricity =
Уравнение эллиптической орбиты, , позволяет вам выразить расстояние r
в терминах угла phi
.
r = a*(1 - eccentricity^2)/(1 + eccentricity*cos(phi));
Средние заряженные частицы степени двойки излучения, перемещающиеся в эллиптическую орбиту, являются интегралом степени излучения по одному полному циклу движения, нормированного периодом движения, . Период движений T
T = 2*pi*sqrt(m*a^3/alpha);
Изменение переменной интегрирования t
к phi
, вы получаете следующий результат. Используйте simplify
функция, чтобы получить более короткий результат интегрирования. Здесь, используйте subs
оценивать J
.
J = subs(J); Javg = simplify(1/T*int(J*m*r^2/L, phi, 0, 2*pi))
Javg =
Оцените среднюю степень излучения электрического диполя с одной частицей, намного более тяжелой, чем, m1>>m2
. Для этого вычислите предел выражения для степени излучения, приняв тот m1
стремится к бесконечности.
limJ = limit(Javg, m1, Inf); simplify(limJ)
ans =