Processing math: 100%

Теорема об умножении определителей

Теорема об умножении определителей. Определитель произведения двух квадратных матриц порядка n равен произведению определителей этих матриц: det(AB)=det(A)det(B) или полная формула: det(ki=1Ai)=ki=1detAi,Ai(P),i=1,,k.

Для доказательства рассмотрим случай k=2. Допустим заданы две матрицы A=aijMn(P) и B=bijMn(P). Воспользуемся вспомогательной блочной матрицей C=A0EB размера 2n×2n, определитель которой имеет вид: Δ=|a11a12a1n000a21a22a2n000an1an2ann000100b11b12b1n010b21b22b2n001bn1bn2bnn|
Вычислим Δ используя теорему Лапласа. Замечаем, что отличным от нуля будет только det(A). Следовательно, Δ=det(A)det(B). Теперь с помощью элементарных преобразований изменим Δ так, что в итоге получим определитель вида |ACEO|. Где C является произведением матриц A и B. Первый столбец умножим на b11 и прибавим к (n+1)-му столбцу, второй на элемент b21 и вновь прибавим к (n+1)-му столбцу. Так же обнулим остальные элементы матрицы B. Записав подробнее полученный определитель имеем: Δ=|a11a12a1nc11c12c1na21a22a2nc21c22c2nan1an2anncn1cn2cnn100000010000001000| Снова вычислим определитель Δ, разложением по последним n столбцам. В этом случае отличным от нуля минором nго порядка будет определитель матрицы C. Поэтому Δ=detCdet(E)=detC(1)n(1)S1+S2, где S1=2nk=n+1k, a S2=nk=1k. В результате получаем Δ=detC(1)2n(n2+n)=detC. Теперь, подставляя имеем доказательство теоремы: Δ=detC=det(AB)=det(A)det(B).

Замечание Известно, что произведение матриц в общем случае не коммутативно, т.е. ABBA. Но определитель это действительное число, а произведение действительных чисел коммутативно. Следовательно, det(AB)=detAdetB=detBdetA=det(BA)

Теорема об умножении определителей является следствием формулы Бине-Коши. Это теорема об определителе произведения прямоугольных матриц, в случае если это произведение дает квадратную матрицу. Справедлива для матриц с элементами любого коммутативного кольца.

Теорема (формула Бине-Коши). Пусть даны две матрицы A и B размеров (m×n) и (n×m) соответственно. Определитель матрицы равен нулю, если m>n, и равен сумме произведений всех соответствующих миноров m-го порядка мaтрицы A на соответствующие миноры m-го порядка матрицы B, если mn. Миноры матриц A и B одинакового порядка, равного наименьшему из чисел n и m, называются соответствующими друг другу, если они стоят в столбцах матрицы A и строках матрицы B с одинаковыми номерами: detAB=γ1<γ2<<γmAγ1<γ2<<γmBγ1<γ2<<γm,
где Aγ1<γ2<<γm — минор матрицы A, составленный из столбцов с номерами γ1<γ2<<γm, и Bγ1<γ2<<γm — минор матрицы B, составленный из строк с номерами γ1<γ2<<γm.

Допустим C=AB, cij=mγ=1aiγbγi. Значит detC=σ(1)σγ1a1γ1bγ1σ(1)γnanγnbγnσ(n)= =mγ1,,γn=1a1γ1annσ(1)σbγ1σ(1)bγnσ(n)=γ1,,γn=1a1γ1anγnBγ1γn. Минор Bγ1γn не равен нулю только в том случае, когда γ1,,γn попарно различны, значит и суммировать можно по парно различные номера γ1,,γn. Для любой перестановки τ этих номеров справедливо Bτ(γ1)τ(γn)=(1)τBγ1γn, из чего следует γ1,,γn=1a1γ1anγnBγ1γn=γ1<γ2<<γn(1)τa1τ(1)anτ(n)Bγ1γn= =γ1<γ2<<γmAγ1<γ2<<γmBγ1<γ2<<γm.

Примеры решения задач

Рассмотрим примеры решения задач связанных с рассмотренной теоремой. Читателю рекомендовано попытаться решить задачи самостоятельно, а затем сверить свое решение с приведенным ниже.

    1. Найти определитель произведения матриц: A=3418,B=2915

      Решение

      Находим определители данных матриц второго порядка: |3416|=18+4=14 и |2715|=107=3. По теореме об определителе произведения матриц получаем: det(AB)=det(A)det(B)=(14)(3)=42. Вычислим этот же определитель, находя произведение матриц: AB=|3416||2715|=|21423| Следовательно, det(AB)=46+4=42. Результаты совпадают.

    2. Найти определитель матрицы пятого порядка: M=12uvw34xyz003210025300342

      Решение

      Разобьём данную матрицу на 4 блока, M=ABOC где A=1234,
      B=uvwxyz, O=000000, C=321253342.
      Представим блочную матрицу как произведение (в справедливости этого представления можно убедиться, найдя произведение по правилам умножения блочных матриц). D=ABCD=E2OTOCE2BOE3AOTOE3, где E2,E3 — единичные матрицы соответствующих порядков.
      |AOTOE3|=detA=|A|, |E2OTOC|=detC=|C|.
      Матрица E2BOE3 — треугольная с единицами на главной диагонали, следовательно ее определитель равен 1 По теореме об определителе произведения получаем:
      |ABOC|=|E2OTOC| |E2BOE3| |AOTOE3|=|C|1|A|=|A||C| Найдем detA и detC. |1234|=2 |321253342|=15836+30+18=3. Подставляя, получаем, detM=23=6

    3. Представьте в виде определителя произведение определителей: |2111121111211112||4114||3113|
      Решение

      По теореме об определителе ступенчатой матрицы имеем:
      |4114||3113|=|4100140000310013| Предположим A=2111121111211112,B=4100140000310013,
      тогда AB=9644274455755515, по теореме об определителе произведения получаем искомый определитель det(AB)=|9644274455755515|.

Литература

  1. Белозеров Г.С. Конспект лекций по линейной алгебре.
  2. В.А. Ильин, Э.Г. Позняк. Линейная алгебра; 5-е изд., стереотипное. ФИЗМАТЛИТ. — 2002. С. 38-39
  3. А.И. Кострикин. Введение в алгебру. Основы алгебры С.138-139
  4. Курош А.Г. Курс высшей алгебры М.: Наука, 1968, С.93-95
  5. Фаддеев Д. К. Лекции по алгебре: Учебное пособие для вузов.— M.: Наука. Главная редакция физико-математической литературы, 1984.— 416 с. C. 130-134

Теорема об умножении определителей

Тест на знание темы «Теорема об умножении определителей».

Соответствие между действиями над операторами и действиями над их матрицами

Как известно, для любого линейного оператора можно определить матрицу этого оператора, при чем такая матрица будет единственной для заданной пары базисов (или одного базиса, в случае оператора из Ω(X), где (X,P) — линейное пространство). Тогда, действия над линейным операторами можно свести к операциям над их матрицами, заданными в фиксированных базисах.

Лемма. В фиксированных базисах, матрицей суммы операторов будет сумма матриц этих операторов.

Зададим два линейных пространства над одним и тем же полем (X,P) и (Y,P) и укажем их размерности, dimX=m, dimY=n. В пространстве X зададим базис e=e1,e2,,em, а в пространстве Yg=g1,g2,,gn.

Зададим линейный оператор AΩ(X,Y). Для оператора A можем записать систему:{Ae1=a11g1+a21g2++an1gn,Ae2=a12g2+a22g2++an2gn,Aem=a1mg1+a2mg2++anmgn.Или можем записать кратко, через сумму:Aej=ni=1aijgi, где j=¯1,m. Тогда, в базисах e и g матрица оператора A будет иметь вид: Age=(a11a12a1ma21a22a2man1an2anm).

Аналогично, зададим линейный оператор BΩ(X,Y). Для него можем записать систему:{Be1=b11g1+b21g2++bn1gn,Be2=b12g2+b22g2++bn2gn,Bem=b1mg1+b2mg2++bnmgn.Или можем записать кратко, через сумму:Bej=ni=1bijgi, где j=¯1,m. Тогда, в базисах e и g матрица оператора B будет иметь вид: Bge=(b11b12b1mb21b22b2mbn1bn2bnm).

Определим линейный оператор C=A+B, где CΩ(X,Y). Для оператора C можем записать систему:{Ce1=c11g1+c21g2++cn1gn,Ce2=c12g2+c22g2++cn2gn,Cem=c1mg1+c2mg2++cnmgn.Или можем записать кратко, через сумму:Cej=ni=1cijgi, где j=¯1,m. Тогда, в базисах e и g матрица оператора C будет иметь вид: Cge=(c11c12c1mc21c22c2mcn1cn2cnm).

Рассмотрим подробнее равенство. ni=1cijgi=Cej= (по определению оператора суммы) =(A+B)ej=Aej+Bej= (используя равенства для Aej и для Bej)=ni=1aijgi+ni=1bijgi=ni=1(aij+bij)gi.Следовательно, ni=1cijgi=ni=1(aij+bij)gi.

Таким образом, каждый элемент матрицы Cge представляет собой сумму соответствующих элементов матриц Age и Bge, что и означает, что Cge=Age+Bge.

Лемма. В фиксированных базисах, матрицей произведения оператора на число будет матрица этого оператора, умноженная на это число.

Зададим два линейных пространства над одним и тем же полем (X,P) и (Y,P) и укажем их размерности, dimX=m, dimY=n. В пространстве X зададим базис e=e1,e2,,em, а в пространстве Yg=g1,g2,,gn.

Зададим линейный оператор AΩ(X,Y). Для оператора A можем записать систему:{Ae1=a11g1+a21g2++an1gn,Ae2=a12g2+a22g2++an2gn,Aem=a1mg1+a2mg2++anmgn.Или можем записать кратко, через сумму:Aej=ni=1aijgi, где j=¯1,m. Тогда, в базисах e и g матрица оператора A будет иметь вид: Age=(a11a12a1ma21a22a2man1an2anm).

Определим линейный оператор C=λA, где CΩ(X,Y), λP. Для оператора C можем записать систему:{Ce1=c11g1+c21g2++cn1gn,Ce2=c12g2+c22g2++cn2gn,Cem=c1mg1+c2mg2++cnmgn.Или можем записать кратко, через сумму:Cej=ni=1cijgi, где j=¯1,m. Тогда, в базисах e и g матрица оператора C будет иметь вид: Cge=(c11c12c1mc21c22c2mcn1cn2cnm).

Рассмотрим подробнее равенство. ni=1cijgi=Cej= (по определению произведения оператора на число) =(λA)ej=λ(Aej)= (используя равенство для Aej)=λni=1aijgi=ni=1λaijgi.Следовательно, ni=1cijgi=ni=1λaijgi.

Таким образом, каждый элемент матрицы Cge представляет собой произведение числа λ на соответствующий элемент матрицы Age, что и означает, что Cge=λAge.

Лемма. В фиксированных базисах, матрицей произведения операторов будет произведение матриц этих операторов.

Зададим три линейных пространства над одним и тем же полем (X,P), (Y,P) и (Z,P) и укажем их размерности, dimX=m, dimY=n, dimZ=k. В пространстве X зададим базис e=e1,e2,,em, в пространстве Yg=g1,g2,,gn, а в пространстве Zt=t1,t2,,tk.

Зададим линейный оператор AΩ(X,Y). Для оператора A можем записать систему:{Ae1=a11g1+a21g2++an1gn,Ae2=a12g2+a22g2++an2gn,Aem=a1mg1+a2mg2++anmgn.Или можем записать кратко, через сумму:Aej=ni=1aijgi, где j=¯1,m. Тогда, в базисах e и g матрица оператора A будет иметь вид: Age=(a11a12a1ma21a22a2man1an2anm).

Аналогично, зададим линейный оператор BΩ(Y,Z). Для него можем записать систему:{Bg1=b11t1+b21t2++bk1tk,Bg2=b12t2+b22t2++bk2tk,Bgn=b1nt1+b2nt2++bkntk.Или можем записать кратко, через сумму:Bgi=kf=1bfitf, где i=¯1,n. Тогда, в базисах g и t матрица оператора B будет иметь вид: Btg=(b11b12b1nb21b22b2nbk1bk2bkn).

Определим линейный оператор C=BA, где CΩ(X,Z). Для оператора C можем записать систему:{Ce1=c11t1+c21t2++ck1tk,Ce2=c12t2+c22t2++ck2tk,Cem=c1mt1+c2mt2++ckmtk.Или можем записать кратко, через сумму:Cej=kd=1cdjtd, где j=¯1,m. Тогда, в базисах e и t матрица оператора C будет иметь вид: Cte=(c11c12c1kc21c22c2kck1ck2ckm).

Рассмотрим подробнее равенство. kd=1cdjtd=Cej= (по определению произведения операторов) =(BA)ej=B(Aej)= (используя равенство для Aej)=Bni=1aijgi=ni=1aijBgi=ni=1aij(Bgi)= (используя равенство для Bgi)=ni=1aijkf=1bfitf=ni=1kf=1aijbfitf==kf=1ni=1bfiaijtf=kf=1(ni=1bfiaij)tf.Следовательно, получили равенство: kd=1cdjtd=kf=1(ni=1bfiaij)tf, а так как d=¯1,k и f=¯1,k, то получаем следующее:cdj=ni=1bdiaij.

Таким образом, каждый элемент матрицы Cte, с индексами d и j равен сумме попарных произведений каждого элемента d-ой строки матрицы Btg на соответствующий элемент j-ого столбца матрицы Age. Это и означает, по определению произведения матриц, что Cte=BtgAge.

Примеры решения задач

  1. Пусть заданы два линейных оператора A(x1,x2,x3)=(x2+x3,2x1+x3,3x1x2+x3),B(x1,x2,x3)=(2x1x2x3,x12x2+x3,x1+x22x3)и базисe=(1,0,0),(0,1,0),(0,0,1).Найти матрицу суммы операторов C=A+B в базисе e.
    Решение
  2. Пусть задан оператор дифференцирования DΩ(R4[x]). Найти матрицу оператора F=2D (FΩ(R4[x])) в базисе e=1,x,x2,x3,x4.
    Решение
  3. Пусть заданы два линейных оператора A(x1,x2,x3)=(x1x2+x3,x3,x2),B(x1,x2,x3)=(2x1+3x2,x1,x2x3)и базисe=(1,0,1),(2,0,1),(1,1,0).Найти матрицу произведения операторов C=BA в базисе e.
    Решение
  4. Пусть заданы два линейных оператора A(x1,x2,x3)=(2x1x2,3x1+x3,2x22x3),B(x1,x2,x3)=(x1+x3,x2x1,3x2+x3)и базисe=(1,0,0),(0,1,0),(0,0,1).Найти матрицу оператора C=2BA+3A в базисе e.
    Решение
  5. Пусть заданы три линейных оператора A(x1,x2,x3)=(x1+x2+x3,2x1x2,3x2+x3),B(x1,x2,x3)=(2x23x3,x1+x3,2x13x2),C(x1,x2,x3)=(x1,x24x3,2x1+6x3)и базисe=(1,0,1),(1,1,0),(0,1,1).Найти матрицу оператора D=A25B+6C в базисе e.
    Решение

Соответствие между действиями над операторами и действиями над их матрицами

Тест на знание темы «Соответствие между действиями над операторами и действиями над их матрицами».

Смотрите также

  1. Воеводин В.В. Линейная алгебра 400 стр. М.: Наука, 1980, cтр. 194-196
  2. Личный конспект, составленный на основе лекций Белозерова Г.С.
  3. Проскуряков И.В. Сборник задач по линейной алгебре. 384 стр. М.: Наука, 1984, стр. 189-190

13.3 Матрица Якоби

Рассмотрим отображение f:ERm, где ERn. Оно состоит из m функций: f=(f1(x1,,xn),f2(x1,,xn),,fm(x1,,xn)), которые осуществляют отображение множества E из Rn в пространство Rm.

Предположим, что функции fk(x1,,xn), где k=¯1,m, дифференцируемы, то есть имеют частные производные по аргументам (x1,,xn):

f1x1,,fnxn,x=¯1,m.

Составим матрицу из этих частных производных по переменным x1,,xn

(f1x1f1x2f1xnf2x1f2x2f2xnfmx1fmx2fmxn)

Такая матрица называется матрицей Якоби.

Если m=n, то получаем квадратную матрицу, определитель которой называется определителем Якоби или якобианом Jf(x) и обозначается

Jf(x)=(f1,,fn)(x1,,xn)=|f1x1(x)f1x2(x)f1xn(x)f2x1(x)f2x2(x)f2xn(x)fnx1(x)fnx2(x)fnxn(x)|.

Замечание. Если все частные производные непрерывны, то и сам оределитель Якоби является непрерывной функцией.

Теорема. Якобиан тождественно равен нулю в некоторой области S:

D(f1,f2,,fn)D(x1,x2,,xn)0 при x=(x1,,xn)S

тогда и только тогда, когда между функциями f1,f2,,fn имеется функциональная зависимость в S, то есть существует функция G(y1,y2,,yn)0 такая, что

G(f1(x),f2(x),,fn(x))0 при всех x=(x1,,xn)S.

Пример 1.Являются ли функции функционально зависимыми?

{f1=x1+x2+x31;f2=x1x2+x1x3+x2x32;f3=x21+x22+x23+3.

Решение.

D(f1,f2,f3)D(x1,x2,x3)=|111x2+x3x1+x3x1+x22x12x22x3|=

=|111x1+x2+x3x1+x2+x3x1+x2+x32x12x22x3|0

Так как якобиан равен нулю, то эти функции функционально зависимы. Несложно найти эту зависимость:

(f1+1)22(f2+2)(f33)=0.

Пример 2. Для линейных функций f1=a11x1++a1nxnb1,,fm=am1x1+amnxnbm матрица Якоби будет матрицей коэффициентов при переменных:

Решение.

(a11a12a1nam1am2amn)

Если мы хотим разрешить систему f1=0,f2=0,,fn=0 относительно x1,,xn, то для случая m=n определитель Якоби

|a11a1nan1ann|

есть определитель системы и для её разрешимости он должен быть отличен от нуля.

Пример 3. Переход элементарной площади dS=dxdy от декартовых координат (x,y) к полярным координатам (r,ϕ):

Решение.

{x=rcos(ϕ);y=rsin(ϕ).

Матрица Якоби имеет вид:

J(r,ϕ)=(xrxϕyryϕ)=(cos(ϕ)rsin(ϕ)sin(ϕ)rcos(ϕ)).
Якобиан перехода от декартовых координат к полярным есть определитель матрицы Якоби:

J(r,ϕ)=detI(r,ϕ)=det(cos(ϕ)rsin(ϕ)sin(ϕ)rcos(ϕ)).

Таким образом, элемент площади при переходе от декартовых к полярным координатам будет выглядеть следующим образом:

dS=dxdy=J(r,ϕ)drdϕ=rdrdϕ.

Пример 4.Переход элементарного объёма dV=dx dy dz от декартовых координат (x,y,z) к сферическим координатам (r,θ,ϕ) :

Решение.

{x=rsin(θ)cos(ϕ);y=rsin(θ)sin(ϕ);z=rcos(θ).

Матрица Якоби имеет следующий вид: I(r,θ,ϕ)=(xrxθxϕyryθyϕzrzθzϕ)=

=(sin(θ)cos(ϕ)rcos(θ)cos(ϕ)rsin(θ)sin(ϕ)sin(θ)sin(ϕ)rcos(θ)sin(ϕ)rsin(θ)cos(ϕ)cos(θ)rsin(θ)0).

А якобиан перехода от декартовых координат к сферическим – есть определитель матрицы Якоби:

J(r,θ,ϕ)=detI(r,θ,ϕ) =

= |sin(θ)cos(ϕ)rcos(θ)cos(ϕ)rsin(θ)sin(ϕ)sin(θ)sin(ϕ)rcos(θ)sin(ϕ)rsin(θ)cos(ϕ)cos(θ)rsin(θ)0|=r2sin(θ).

Таким образом, элемент объёма при переходе от декартовых к сферическим координатам будет выглядеть следующим образом:

dV=dxdydz=J(r,θ,ϕ)drdθdϕ=r2sin(θ)drdθdϕ.

Матрица Якоби

Для закрепления пройденного материала предлагается пройти тест.

Список использованной литературы

  1. Коляда В.И., Кореновский А.А. Курс лекций по математическому анализу.-Одесса : Астропринт, 2009. стр.309-311
  2. Демидович Б.П. «Сборник задач и упражнений по математическому анализу» 13-е издание, 1997 М.: Изд-во Моск. ун-та, ЧеРо. №3990.
  3. Фихтенгольц Г.М. Курс дифференциального и интегрального исчисления: Том 1 / Г.М. Фихтенгольц – М.: Книга по Требованию, 2013. стр.455-456.