Небесная механика - один из разделов астрономии.
Разделы астрономии:
- общая астрономия
- космогония
- космология
- астрометрия (геодезическая, сферическая, оптическая)
- небесная механика (теоретическая астрономия)
По определению Дубошина предметом небесной механики является изучение поступательного, вращательного и деформационного движений небесных тел как искусственного, так и естественного происхождения, под действием различных сил: сил гравитации, сил отталкивания, сил сопротивления среды, сил электромагнитных и др.
Всякое сложное движение тела можно разложить на 3 составляющих:
- поступательное движение
- деформационное движение.
Пример 1:
Движение Земли в солнечной системе можно разложить на поступательное движение центра масс Земли по орбите вокруг Солнца, вращательное движение Земли вокруг центра масс и деформационное движение Земли под действием приливов в оболочках Земли, движение литосферных плит и тепловое расширение и сжатие Земли.
Пример 2:
Полет футбольного мяча.
В курсе небесной механики будем изучать только поступательное движение центра масс.
Основные задачи небесной механики:
Все множество задач делится на две группы:
1. Определение орбит небесных тел.
2. Вычисление эфемерид небесных тел.
1 группа задач: по наблюдаемому (измеряемому) видимому движению небесного тела определяются параметры траектории движения, такие, которые во время движения остаются постоянными величинами. Эти параметры называют элементами орбиты и их определение и составляет предмет первой группы задач.
Наблюдать (измерять) – означает регистрировать с помощью аппаратуры различные кинематические движения: расстояние, направление (углы), скорость изменения кинематических величин (радиальная, угловая).
2 группа задач: осуществляется прогнозирование движения различных небесных объектов с помощью построения теоретической модели движения и параметров орбиты, полученных из решения 1-й группы задач.
Прогнозирование движения тождественно понятию вычислений эфемериды небесного объекта. Эфемерида может быть задана в двух видах: либо в аналитическом (совокупность формул, явно выражающих зависимость координат векторов положения и скорости небесного тела от элементов орбиты и времени), другой вид эфемерид табличный, когда координаты векторов в положении скорости спутников задаются в виде таблицы.
координаты вектора положения | координаты вектора скорости | |||||
Примеры таких таблиц есть в астрономических ежегодниках.
Методы небесной механики
Все множество методов решения двух данных задач делится на две группы:
1. аналитические методы
2. численные методы
В аналитических методах небесной механики математическая модель движения небесного объекта представляется в явном виде, когда любая координата (x, y,z, ẋ,ẏ,ż) представляется как аналитическая функция:
x = x(а, е, i, , ω, tπ, t)
ż= ż(а, е, i, Ω, ω, tπ, t) ,
где а, е, i, Ω, ω, tπ – элементы орбиты небесного тела, а t – время.
В численных методах математическая модель движения представляется в дифференциальном виде, то есть в виде дифференциальных уравнений движений, когда искомые функции x(t), y(t), z(t), ẋ(t), ẏ(t), ż(t), находятся под знаком оператора дифференцирования, и отыскание (нахождение) их решений осуществляется путем численного интегрирования дифференциальных уравнений одним из методов (Рунге, Кутты, Адамса и др.).
Так, например, движение ИСЗ очень просто и точно описывается дифференциальными уравнениями:
d2x/dt2=-μx/r3+Φx(x, y, z, t)
)
d2 z/dt2=-μz/r3+Φz(x, y, z, t) .
В формулах (1.2) искомые функции x (t), y (t), z (t) находятся под знаком оператора двойного дифференцирования, под знаком квадратного корня в
третьей степени r3=( √x2+y2+z2 )3 и под знаком возмущающих функций Φx, Φy, Φz.
Если результатом аналитического метода являются формулы (1.1), то результатом численных методов небесной механики является таблица в виде эфемериды, описанной выше.
Аналитические методы позволяют качественно охарактеризовать геометрическую кинематику и динамику движения небесного тела, но точность аналитических методов всякий раз отстает от запросов практики. Прогресс измерительной техники в точности измерений идет более быстрыми темпами, чем повышается точность аналитических методов. Численные же методы наоборот, не позволяют получить качественную картину движения небесного тела, но дают возможность прогнозировать движение методически с неограниченной точностью. Предел точности ставят сами измерения и адекватность дифференциального уравнения реальному движению.
РАЗДЕЛ 1. НЕВОЗМУЩЕННОЕ ДВИЖЕНИЕ СПУТНИКА
Тема 2. Постановка задачи невозмущенного движения спутника
Тема 2.1. Понятие динамической системы, используемой для решения геодезических задач.
В космической геодезии исследуется состояние реальной динамической системы (рис. 2.1): 
Рис. 2.1. Реальная динамическая система
А – Земля,
S - искусственный спутник Земли,
- планеты,
Fa – сила притяжения спутника S Землей А,
F - Солнцем,
F - Луной,
F - планетами,
Fатм – сопротивление движения S плотными слоями атмосферы,
Fсв д – световое давление (солнечная радиация или солнечный ветер),
Fтепл д – тепловое давление Земли, тепловая радиация,
Fотраж – отраженное световое давление от Земли и Луны,
Fэл магн – электромагнитные силы взаимодействия магнитного поля Земли с электромагнитным полем спутника,
и другие силы.
Для выделения главной составляющей движения спутника реальное ЭДС заменяется моделью, в которой делается 4 следующих основных допущения:
1 допущение: в динамической системе рассматриваются только силы гравитации (притяжения), а все остальные силы отбрасываются, так как они по величине примерно в 1000 раз меньше главной силы гравитации
(Fʘ=0; Fатм=0; Fтепл=0).
2 допущение: из всех тел, создающих гравитацию, рассматриваются только два тела: центральное тело, масса которого (М) во много раз превышает массу спутника (в космической геодезии это Земля). Второе тело, движение которого измеряется и используется в дальнейшем, называется спутником (m)
m<<M.
3 допущение: несферическая структура гравитационного поля Земли (ГПЗ) заменяется сферической моделью ГПЗ.
Термин несферическое ГПЗ эквивалентен термину нецентральное ГПЗ (рис. 2.2).
Термин сферическое ГПЗ эквивалентен термину центральное ГПЗ
(рис. 2.3).
Рис. 2.2 – нецентральное ГПЗ
В нецентральном ГПЗ силовые линии, пространственные кривые не пересекаются в одной точке. Касательные к силовым линиям в точке нахождения спутника также не пересекаются в одной точке.
На рисунке 2.3 покажем три эквивалентных модели центрального ГПЗ.
а) б) в)
Сферически симметричная однородный материальная точка
Модель ГПЗ шар
Рис. 2.3. Три модели центрального ГПЗ
В случае сферической симметрии ГПЗ плотность (ρ) вещества внутри Земли есть функция только одной координаты r – расстояния от внутренней элементарной массы до точки О
ρ = ρ(r). (2.1.)
Тогда как в случае несферической структуры ГПЗ плотность ρ в земле – есть функция трех координат.
ρ = ρ(r,φ,λ) , (2.2.)
где φ и λ – сферические координаты.
В случае ГПЗ в виде шара ρ – есть величина постоянная
ρ =const.
В случае ГПЗ, создаваемого материальной точкой, понятие плотности теряет смысл, в этом случае пользуются понятием материальной точки – это геометрическая точка, снабженная массой М всего тела А.
Для Земли обычно в качестве геометрической точки выбирают центр масс О и снабжают его массой
4 допущение: поскольку размеры спутника во много раз меньше, чем расстояние до Земли, то несферическое гравитационное поле, создаваемое телом спутника, заменяется сферическим, это значит, что спутник заменяется на материальную точку S с массой m.
Таким образом, мы приходим к следующей идеальной модели динамической системы (рис.2.4):
Рис. 2.4. Идеальная модель динамической системы (ДС)
Изучение движения двух материальных точек называется называется задачей двух тел или невозмущенным движением.
В задаче двух тел действуют только две силы взаимного притяжения
(Fа и Fs), с которыми Земля и спутник притягивают друг друга, больше никакие другие силы не рассматриваются. Если к этим двум силам прибавить любую другую из рассматриваемых выше сил, такое движение называется возмущенным движением.
Тема 3. Понятие универсальной инерциальной системы координат
Мы будем изучать решение задачи двух тел в инерциальной системе координат, так как в ней наиболее просто составляется уравнение движения спутника. Инерциальной системой называется такая система, которая движется в пространстве поступательно, равномерно, прямолинейно, с постоянной ориентировкой осей относительно внешнего пространства.
Рис. 3.1. Поступательная прямолинейная равномерная ИСК
Оси координат должны быть неподвижны относительно других имеющихся тел. В природе все тела движутся криволинейно, неравномерно и вращаются. Создать строгую абстрактную ИСК невозможно, но так как на практике все задачи имеют ограниченную точность решения на ограниченном отрезке времени, вместо ИСК пользуются понятием квазиинерциальных систем координат. Квазиинерциальная система координат на ограниченном отрезке времени движется практически прямолинейно, равномерно и почти не меняет своей ориентировки.
Пример 1: Квази ИСК, используемая в космической геодезии – это система координат, у которой начало координат помещается в центр масс Земли. Ось аппликат направляется по оси вращения Земли, а ось абсцисс по линии узлов орбиты Земли, то есть по линии пересечения плоскости орбиты Земли с плоскостью экватора Земли (рис. 3.2).
Рис. 3.2. Квазиинерциальная система координат Оxyz c началом О в центре масс Земли, осью аппликат Оz по оси вращения, и осью абсцисс по линии узлов орбиты Земли
На рис.3.2 точка γ – восходящий узел орбиты Земли или точка весеннего равноденствия.
На ограниченном отрезке времени, на котором ищется решение конкретной практической задачи, кусок траектории орбиты Земли от t1 до t2 можно рассматривать прямой линией и движение Земли происходит равномерно и прямолинейно.
Тема 4. Составление дифференциального уравнения невозмущенного движения спутника
Тема 4.1. Постановка задачи.
Даны две материальные точки: М – центральное тело и m – спутник, находящийся в движении под действием только силы взаимного притяжения (рис. 4.1.).
Движение двух тел (М и m) изучается относительно ИСК Оxyz.
Рис. 4.1. Неограниченная задача двух тел
Требуется описать движение двух материальных точек М и m относительно ИСК Оxyz в дифференциальной форме (то есть в виде дифференциальных уравнений).
Тема 4.2. Три закона Ньютона, используемые для составления дифференциальных уравнений (ДУ).
Первый закон – закон всемирного тяготения:
Две материальные точки притягиваются друг к другу с силой, пропорциональной их массам и обратно пропорциональной квадрату расстояния r между ними:
F=fMm/r
Второй закон – закон инерции: сила, действующая на тело, пропорциональна массе тела m и его ускорению ω:
F= mω .
Третий закон – закон равенства действия противодействию (4.1)
|Fm| = |FM| или Fm = - FM. (4.3)
Тема 4.3. Векторная форма ДУ невозмущенного движения спутника в неограниченной задаче двух тел.
Запишем силу притяжения спутника FM центральным телом M в векторной форме. Для этого используем понятие орта и модуля любого вектора
а=аоа, (4.4)
где ао - орт вектора а, а=|а| - модуль вектора а.
Применим (4.4) к векторам Fm и FM
Fm = Fmо Fm, FM = FMо FM. (4.5)
Используя (4.1), получаем, что
Fm= FmоfMm/r2, FM= FMо fMm/r
Введем радиусы - векторы:
ρ, характеризующий движение спутника m относительно инерциальной системы Oxyz ; и а, характеризующий движение центрального тела М относительно ИСК.
Введем вектор r, характеризующий взаимное положение двух тел m и M.
Из (рис. 1.1) следует, что орт вектора rо= FMо и rо= - Fm,
rо = r/ r. (4.7)
С учетом (4.7) перепишем (4.5) и (4.6) так:
|
FM = (fMm/r2)(r/ r) .
Используем второй закон Ньютона, имеем:
Fm= mωm ; FM = MωM. (4.8*)
Найдем ускорение спутника ωm и центрального тела ωM. Из треугольника ОmM имеем: сумма векторов в замкнутой фигуре равна нулю, то есть
r =ρ – а . (4.9)
Дифференцируя дважды (4.8) по времени t, получаем:
r(t) = ρ(t) – а(t) , (4.9)
ŕ(t) = ρ˙(t) – а˙(t) , (4.10)
ŕ˙(t) = ˙ρ˙(t) – ˙а˙(t) , (4.11)
где точками сверху обозначены производные по времени t.
ŕ(t) =(d/dt)(r(t)) ;
ŕ˙(t)= d/dt(d/dt(r(t)))= (d2F(t))/dt2.
Из курса физики известно, что первая производная по времени t от радиуса-вектора r(t) есть вектор скорости:
ŕ(t) = V(t) , (4.12)
а вторая производная
ŕ˙(t) = ω(t) = d/dt(V(t)
из (рис. 4.1.) имеем
ωm =˙ρ˙ ; ωM = ˙а˙ . (4.14)
Подставляя в (4.8.1) получим
Fm=m˙ρ˙ , FM =М˙а˙ . (4.15)
Ньютон высказал гипотезу, что приравнивая (4.8) и (4.15), получаем:
|
М˙а˙= (fMm/r2)(r/|r|) ,
поскольку М≠0, m≠0 можно сократить в (4.16)
˙ρ˙= - (fMm/r2)(r/|r|) , ˙а˙= (fMm/r2)(r/|r|
Выражение (4.17) и есть дифференциальное уравнение движения спутника m и центрального тела М относительно ИСК Oxyz в неограниченной задаче двух тел.
Тема 4.4. Дифференциальное уравнение невозмущенного движения в ограниченных задачах двух тел.
Покажем, что модуль ускорения движения центрального тела |˙а˙| намного меньше модуля вектора ускорения спутника |˙ρ˙|. Вычислим модуль |˙ρ˙|:
|˙ρ˙|=|- (fM/r2)(r/|r|)| = (fM/r2)(r/|r|) = fM/r2 .
Вычислим |˙а˙|:
|˙а˙|=|- (fm/r2)(r/|r|)| = fm/r2.
Беря отношение, получаем:
|˙ρ˙|/|˙а˙|=fM/r2 ; fm/r2=М/m
или
|˙а˙|/|˙ρ˙|= m/М. (4.18)
Отсюда видно, если М<< m, то |˙ρ˙| >> |˙а˙|.
Если применить задачу двух тел к ИСЗ, масса которого 1 тонна, и поделить на массу Земли, то отношение будет миллиардным.
Вывод: начало ИСК из точки О можно перенести в точку М(центр масс центрального тела). И тогда мы получаем ограниченную задачу двух тел. (рис.4.2.).
Рис. 4.2. Ограниченная задача двух тел.
Получим ДУ для ограниченной задачи двух тел.
Вычтем левые и правые части в выражениях:
˙ρ˙- ˙а˙ = (fM/r2)(-r/|r|) - (fm/r2)(r/|r|) ,
˙ρ˙ -˙а˙ = (f(M+m)/r2)(-r/|r|) .
Задача двух тел или ˙ρ˙- ˙а˙= ŕ˙, то
ŕ˙=((fM(1+m/М))/r2) / (-r/|r|
Обычно в ограниченной задаче членом m/М≈0 пренебрегают, тогда
ŕ˙=(fM/r2)(-r/|r|) , (4.20)
где вводят обозначение
μ=fM, (4.21)
величина μ – гравитационный параметр тела М. Коэффициент f – универсальная постоянная тяготения, которая справедлива для Земли.
Тема 4.5. Дифференциальные уравнения невозмущенного движения спутника в координатной форме.
Мы получили ДУ невозмущенного движения в векторной форме относительно ИСК Oxyz в векторном виде.
ŕ˙+ (μ/ r2)(r/|r|)=
Рис. 4.2. Невозмущенное движение спутника относительно ИСК
Чтобы преобразовать векторную форму уравнения (4.19) в координатную, надо спроектировать это равенство на оси ИСК Oxyz. Проектирование любого вектора на заданное направление осуществляется с помощью скалярного произведения двух векторов.
Из векторной алгебры скалярное произведение – есть скаляр (число), равный произведению модуля этих векторов (сомножителей) на cos угла между ними
аb , (4.22)
(а,b
В одном тексте используется только какое-либо одно обозначение. Мы будем использовать первое (4.22)
аb=|а||b|cos(а,b) , (4.24)
или так: аb=|а||b|cos(α) , (4.25)
где α=(а˄,b) – угол между а и b.
Рассмотрим геометрический смысл скалярного произведения (рис. 4.3):
Рис. 4.3. Геометрия скалярного произведения двух векторов а и b
Перепишем (4.25),как
аb=а(bоb)= b[аbо]= b[аcosα]= аb=bапр на b , (4.26)
аbо= апр на b . (4.27)
Из векторной алгебры: любой вектор можно представить его разложением по координатным осям (а в более общем случае в виде разложения по двум, трем и более направлениям).
Во всякой исходной СК (в нашем случае ИСК Oxyz) вводятся по умолчанию орты осей координат i,j,k с координатами:
i={1,0,0},j={0,1,0},k={0,0,1} . (4.28)
Теперь, чтобы спроектировать равенство (4.19) на оси координат, надо представить каждый вектор, входящий в него r, ˙r˙ в виде разложения по координатам:
r=x
x=ОS=ri
y=ОP=rj (4.29)
z=ОU=rk
r=x+y+z
или
r = xi+yj+zk . (4.30)
Рис. 4.4 Разложение r по координатным осям
Если даны координаты какого-либо вектора
r= {x, y,z} задание вектора прямоугольными координатами
r=xi+yj+zk – разложение вектора.
Дифференцируя дважды равенство (4.30) получаем:
r˙=x˙i + y˙j + z˙k + xdi/dt + ydj/dt + zdk/dt,
но так как Oxyz есть ИСК (i=const, j=const, k=const)
di/dt = dj/dt = dk/dt =0,
r˙ = x˙i + y˙j + z˙k . (4.31)
Вектор r˙ есть вектор скорости с координатами x˙ y˙ z˙
r˙= {x˙,y˙,z˙}=V . (4.32)
где V – другое обозначение вектора скорости
r˙˙= (d/dt)i = x˙˙i y˙˙j+ z˙˙k , (4.33)
или
r˙˙= {x˙˙, y˙˙, z˙˙} , (4.34)
x˙˙, y˙˙, z˙˙ - координаты вектора ускорения спутника.
Теперь, подставляя в дифференциальное уравнение (4.19) выражение (4.34) и (4.31) получаем:
x˙˙i + y˙˙j + z˙˙k= - (μ/ r3)( xi + yj + zk
Проектируя (4.35) на оси Оx, затем на ось Oy и, наконец, на ось Oz путем умножения скалярного на орт i, затем на j и наконец на k, получаем
x˙˙(i i)+y˙˙(j i)+z˙˙(k i)= - (μ/ r3)( x(ii)+y(ji)+z(ki)),
x˙˙=- (μ/ r3) x.
Поступая аналогично, получаем два других равенства:
y˙˙=- (μ/ r3) y
z˙˙=- (μ/ r3) z .
Выражение (4.36) – дифференциальное уравнение невозмущенного движения в ИСК в координатной форме.
Тема 5. Постановка задачи интегрирования дифференциальных уравнений невозмущенного движения спутника
Мы получили математическую модель движения спутника в наиболее простом виде (в виде ДУ)
r˙˙+-(μ/ r3)r=
Это уравнение не дает возможности получить качественную характеристику траектории, необходимо его разрешить относительно вектора положения r
r=r(t) или x=x(t), y(t), z(t
Искомые функции входят в дифференциальное уравнение под знаком дифференцирования дважды:
(d2/dt2)(r(t))=˙r˙ ,
так же находится степень 3/2
r3=(√x2+y2+z2 )3.
Это нелинейное дифференциальное уравнение.
Из двух способов решения ДУ (4.19) мы рассмотрим аналитический метод решения уравнения. В этом методе ищутся явные функции вектора положения r={x, y, z} и вектора скорости r˙={x, y,z} через начальные условия движения движения спутника
|
rо˙=r˙(tо) или xо˙=x˙(tо), yо˙=y˙(tо), zо˙=z˙(tо) .
Итак, требуется получить формулы
xо = x(μ, rо,iо, tо, t)
yо = y(μ, rо,iо, tо, t)
zо = z(μ, rо,iо, tо, t) (5.3)
zо˙ = z˙ (μ, rо,iо, tо, t) .
Тема 6. Интеграл площадей.
В этом разделе мы получим один из, так называемых, первых интегралов.
r˙˙+(μ/r3) r = 0 (6.1)≡(4.19)
Из раздела дифференциального и интегрального исчисления в математике термин «проинтегрировать выражение» f(t), то есть найти
∫f(t)dt,
означает найти такую функцию F(t), называемую первообразной, производная которой по t дает подынтегральное выражение f(t), то есть
F(t)= ∫f(t)dt, (6.2)
(dF/ dt)= f(t).
Преобразуем (6.1) к виду, удобному для интегрирования, то есть удобному для получения первообразной функции; для этого умножим векторно слева равенство (6.1) на текущий радиус - вектор спутника.
r = r(t),
r× r˙˙ + (μ/ r3)( r×r)= r×0 . (6.3)
По свойству модуля векторного произведения имеем:
|а × b |=а bsin(а˄b),
| r × r |= r rsin0о =0.
Тогда получаем вектор
r × r˙˙=0. (6.4)
Ищем интеграл по t:
∫(r × r˙˙)dt= ∫0dt. (6.5)
Подберем первообразную к подынтегральному выражению r × r˙˙.
Пусть первообразная
r × r˙. (6.6)
Проверим, так ли это, дифференцируя по t:
(d/dt)(r × r˙) = r˙× r˙+ r × r˙˙= r × r˙˙.
Подставляя (6.7) в (6.5) получаем
∫(d/dt)(r × r˙)dt = с
Или
r × r ˙= с . (6.8)
Выражение (6.8) в небесной механике носит название интеграла площадей.
Векторная константа с – производная постоянная интегрирования и носит название векторной константы площадей.
с определяется по начальным условиям движения спутника в положении rо= r(tо) и в скорости r˙о= r˙(tо) , которая должна быть известна.
Интеграл площадей справедлив для любого момента времени t
r(t) × r˙( t) = с , (6.8)
значит он справедлив для момента tо
с = rо × rо˙ . (6.9)
Теперь с получает конкретные численные координаты. Мы получили интеграл площадей в векторной форме (6.8).
Интеграл площадей в координатной форме.
Здесь надо поступать так же, как это делалось в переходе от векторной формы Оу (6.1) к координатной форме (см раздел 5).
Для этого надо каждый вектор, входящий в векторное равенство (6.8), представить его разложение по координатам ( в проекции на оси координат Оx, Оy, Оz):
с = сxi + сyj + сzk
r = xi + yj + zk (6.10)
r˙= x˙i + y˙j + z˙k ,
с = {сх, сy, сz}
r = { x, y, z } (6.10)
r˙= { x˙,y˙,z˙} ,
и подставим (6.10) в (6.8)
(x i+ y j+ z k) ×( x˙i+ y˙j+ z˙k) = сxi + сyj + сzk (6.11)
xx˙i × i + y x˙j × i + z x˙k × i +… =с x i+…
…+ xy˙i × j + y˙j × j + zy˙k × j +… +с y j+… (6.12)
…+ xz˙i × k + z˙yj × k + zz˙k × k +с z k
j × k = + i; j × i= - k; k × i= j .
- yx˙k + zx˙j + xy˙k - zy˙i – xzj + yz˙i = сxi + сyj + сzk . (6.13)
Приравниваем множители при одинаковых ортах слева и справа от знака равенства, получаем:
с x= y z ˙ - z y ˙
с y= z x ˙ - x z ˙ (6.14)
с z= x y ˙ - y x ˙ .
Правило циклической перестановки:
 |
Тема 7. Нормальное уравнение плоскости орбиты в векторной форме. Физический смысл векторной константы площадей
Умножим интеграл площадей (6.8)
r × r˙ = с (6.8) ≡(7.1)
скалярно на текущий радиус – вектор спутника r = r(t)
(r × r˙)r = сr . (7.2)
Левая часть равенства (7.2) есть смешанное произведение (векторно – скалярное) трех векторов:
(а × b) с = f – скаляр, (7.3)
Численно равный V параллелепипеда, построенного на этих векторах.
Отсюда следует, что если два вектора из трех равны друг другу и совпадают по направлению, то V = 0.
Смешанные произведения подчиняются циклической перестановке векторов:
(а × b) с = (b × с) а = - (с × b) а, (7.4)
(r × r˙) r = (r × r) r˙ = 0 r˙ = 0 = с r,
сr =
Покажем, что выражение (7.5) представляет собой уравнение плоскости.
Из аналитической геометрии известно, что всякое уравнение первой степени от трех переменных есть уравнение плоскости
Аx+Вy+Сz+D = 0 , (7.6)
А, В,С, D – производные константы.
x, y, z – координаты текущей точки в трехмерном пространстве Оxyz.
Перепишем скалярное уравнение плоскости (7.6) в векторном виде, для этого введем вектор N с тремя компонентами
N= {А, В, С} или N=Аi + Вj + Сk . (7.7)
|
Из за большого объема этот материал размещен на нескольких страницах:
1 2 3 4 5 |
