Практическая информатика. Практическая работа № 9. Практическая работа 9 Матричный метод решения
Скачать 23.45 Kb.
|
Практическая работа № 9 Матричный метод решения Задание 1.Обозначим через А — матрицу коэффициентов при неизвестных; X — матрицу-столбец неизвестных; B - матрицу-столбец свободных членов: Вектор B: BT=(-4,17,3) С учетом этих обозначений данная система уравнений принимает следующую матричную форму: А*Х = B. Если матрица А — невырожденная (ее определитель отличен от нуля, то она имеет обратную матрицу А-1. Умножив обе части уравнения на А-1, получим: А-1*А*Х = А-1*B, А-1*А=Е. Это равенство называется матричной записью решения системы линейных уравнений. Для нахождения решения системы уравнений необходимо вычислить обратную матрицу А-1. Система будет иметь решение, если определитель матрицы A отличен от нуля. Найдем главный определитель. ∆=2•(1•3-(-1•(-1)))-3•(-7•3-(-1•1))+1•(-7•(-1)-1•1)=70 Итак, определитель 70 ≠ 0, поэтому продолжаем решение. Для этого найдем обратную матрицу через алгебраические дополнения. Пусть имеем невырожденную матрицу А: = Тогда: = где Aij — алгебраическое дополнение элемента aij в определителе матрицы А, которое является произведением (—1)i+j на минор (определитель) n-1 порядка, полученный вычеркиванием i-й строки и j-го столбца в определителе матрицы А. Транспонированная матрица к матрице A имеет вид: Вычисляем алгебраические дополнения. ∆1,1=(1•3-(-1•(-1)))=2 ∆1,2=-(-7•3-1•(-1))=20 ∆1,3=(-7•(-1)-1•1)=6 ∆2,1=-(3•3-(-1•1))=-10 ∆2,2=(2•3-1•1)=5 ∆2,3=-(2•(-1)-1•3)=5 ∆3,1=(3•(-1)-1•1)=-4 ∆3,2=-(2•(-1)-(-7•1))=-5 ∆3,3=(2•1-(-7•3))=23 Из полученных алгебраических дополнений составим присоединенную матрицу C: Вычислим обратную матрицу: Вектор результатов X X=A-1 • B XT=(5,2,0) x1=350 / 70=5 x2=140 / 70=2 x3=0 / 70=0 Проверка. 2•5-7•2+1•0=-4 3•5+1•2-1•0=17 1•5-1•2+3•0=3 Задание 2. Вектор B: BT=(0,-21,-4) Найдем главный определитель. ∆=1•(-5•1-1•(-6))-3•(-4•1-1•(-2))+3•(-4•(-6)-(-5•(-2)))=49 Итак, определитель 49 ≠ 0, поэтому продолжаем решение. Для этого найдем обратную матрицу через алгебраические дополнения. Транспонированная матрица к матрице A имеет вид: Вычисляем алгебраические дополнения. ∆1,1=(-5•1-(-6•1))=1 ∆1,2=-(-4•1-(-2•1))=2 ∆1,3=(-4•(-6)-(-2•(-5)))=14 ∆2,1=-(3•1-(-6•3))=-21 ∆2,2=(1•1-(-2•3))=7 ∆2,3=-(1•(-6)-(-2•3))=0 ∆3,1=(3•1-(-5•3))=18 ∆3,2=-(1•1-(-4•3))=-13 ∆3,3=(1•(-5)-(-4•3))=7 Из полученных алгебраических дополнений составим присоединенную матрицу C: Вычислим обратную матрицу: Вектор результатов X X=A-1 • B = XT=(-2,-3,5) x1=-98 / 49=-2 x2=-147 / 49=-3 x3=245 / 49=5 Проверка. 1•(-2)-4•(-3)-2•5=0 3•(-2)-5•(-3)-6•5=-21 3•(-2)+1•(-3)+1•5=-4 Задание 3. Вектор B: BT=(1,-7,0) Найдем главный определитель. ∆=1•(-3•(-2)-1•(-1))-2•(2•(-2)-1•(-3))+4•(2•(-1)-(-3•(-3)))=-35 Итак, определитель -35 ≠ 0, поэтому продолжаем решение. Для этого найдем обратную матрицу через алгебраические дополнения. Транспонированная матрица к матрице A имеет вид: Вычисляем алгебраические дополнения. ∆1,1=(-3•(-2)-(-1•1))=7 ∆1,2=-(2•(-2)-(-3•1))=1 ∆1,3=(2•(-1)-(-3•(-3)))=-11 ∆2,1=-(2•(-2)-(-1•4))=0 ∆2,2=(1•(-2)-(-3•4))=10 ∆2,3=-(1•(-1)-(-3•2))=-5 ∆3,1=(2•1-(-3•4))=14 ∆3,2=-(1•1-2•4)=7 ∆3,3=(1•(-3)-2•2)=-7 Из полученных алгебраических дополнений составим присоединенную матрицу C: Вычислим обратную матрицу: Вектор результатов X X=A-1 • B XT=(0,2,1) x1=0 / (-35)=0 x2=-70 / (-35)=2 x3=-35 / (-35)=1 Проверка. 1•0+2•2-3•1=1 2•0-3•2-1•1=-7 4•0+1•2-2•1=0 |