8024

Матричные операции и решения СЛАУ в MatLab.

Лабораторная работа

Информатика, кибернетика и программирование

Матричные операции и решения СЛАУ в MatLab. Индивидуальные задания ЗАДАНИЕ № 1. Найти определитель тремя способами: 1) методом треугольника, 2) по теореме Лапласа, раскладывая по элементам 1-ой строки, 3) по встроенной команде det(A). Сравнить все р...

Русский

2013-02-01

23.81 KB

21 чел.

Матричные операции и решения СЛАУ в MatLab.

Индивидуальные задания

ЗАДАНИЕ 1.

Найти определитель тремя способами: 1) методом треугольника, 2) по теореме

Лапласа, раскладывая по элементам 1-ой строки, 3) по встроенной команде det(A).

Сравнить все результаты между собой и в случае их несовпадения найти и исправить

свои ошибки.

4)

%Metod triugolnica

>>A=[[2 9 3];[1 -2 1];[3 7 2]];

>>d1=A(1,1)*A(2,2)*A(3,3)+A(1,2)*A(2,3)*A(3,1)+A(2,1)*A(3,2)*A(1,3);

>>d2=A(1,3)*A(2,2)*A(3,1)+A(2,1)*A(1,2)*A(3,3)+A(2,3)*A(3,2)*A(1,1);

>>Det=d1-d2

Det =    26

%Teorema Laplasa

>>Laplas=A(1,1)*(A(2,2)*A(3,3)-A(2,3)*A(3,2))-A(1,2)*(A(2,1)*A(3,3)->>A(2,3)*A(3,1))+A(1,3)*(A(2,1)*A(3,2)-A(2,2)*A(3,1))

Laplas =    26

%DetA

>>DetA=det(A)

DetA =    26

ЗАДАНИЕ 2. Вычислить определитель 4-го порядка двумя способами: по

теореме Лапласа раскладывая по элементам 1-ой строки и сравнить результат со

значением, вычисленным по команде det(A).

4)

%%Opredelitel 4*4

>>A=[[1 -2 8 2];[2 2 3 -5];[-1 -1 4 1];[4 6 -2 7]];

>>A11=A(2:4,2:4)

>>A12=[A(2:4,1),A(2:4,3:4)]

>>A13=[A(2:4,1:2),A(2:4,4)]

>>A14=A(2:4,1:3)

>>d1=Deter(A11);d2=Deter(A12);d3=Deter(A13);d4=Deter(A14);

>>DEL=d1*A(1,1)-d2*A(1,2)+d3*A(1,3)-d4*A(1,4)

A11 =

    2     3    -5

   -1     4     1

    6    -2     7

A12 =

    2     3    -5

   -1     4     1

    4    -2     7

A13 =

    2     2    -5

   -1    -1     1

    4     6     7

A14 =

    2     2     3

   -1    -1     4

    4     6    -2

delta A11 =  209

delta A12 =   163

delta A13 =  6

delta A14 = -22

DEL =   627

>>Det=det(A)

Det =  627

ЗАДАНИЕ 3. Решить СЛАУ 3-го порядка по формулам Крамера.

4)

>>Glavniy=[[2 -1 5];[5 2 13];[3 -1 5]];

>>DOP=[14;-15;-4];

>>Dopx=[DOP, Glavniy(:,2:3)]

>>Dopy=[Glavniy(:,1), DOP,Glavniy(:,3)]

>>Dopz=[Glavniy(:,1),Glavniy(:,2),DOP]

>>DET=det(Glavniy);

>>DETx=det(Dopx);

>>DETy=det(Dopy);

>>DETz=det(Dopz);

>>x=DETx/DET

>>y=DETy/DET

>>z=DETz/DET

Dopx =

   14    -1     5

  -15     2    13

   -4    -1     5

x =   -18

y =  -11.9565

Dopy =

    2    14     5

    5   -15    13

    3    -4     5

Dopz =

    2    -1    14

    5     2   -15

    3    -1    -4

z =    7.6087

%%Proverka

>>Glavniy*X

ans =

  14.0000

 -15.0000

  -4.0000

ЗАДАНИЕ 4. Проверить результат предыдущего задания, применив встроенную

операцию левого деления матриц.

>>Glavniy=[[2 -1 5];[5 2 13];[3 -1 5]];

>>DOP=[14;-15;-4];

>>X=Glavniy\DOP

X =

 -18.0000

 -11.9565

   7.6087

ЗАДАНИЕ 5. Найти наилучшее решение переопределенной системы уравнений.

Сделайте проверку. Для этого полученное решение X подставьте в исходную

систему линейных уравнений A* X = B и найдите значение невязок   для каждой

строчки. Решение правильно, если произведение  равно нулю.

>>A=[[2;5;3;2],[-1;2;-1;-3]];

>>B=[10;31;8;6];

>>P=A\B

>>EX=A*P-B

>>ATEX=A'*EX

P =

   5.0731

   2.0715

ATEX =

 1.0e-014 *

        0

   0.5329

EX =

  -1.9253

  -1.4913

   5.1479

  -2.0684


 

А также другие работы, которые могут Вас заинтересовать

20717. ПОЛНЫЕ МЕТРИЧЕСКИЕ ПРОСТРАНСТВА 57 KB
  Чтобы разобраться в этом вопросе рассмотрим понятие фундаментальной последовательности на R. Определение: последовательность {xn} называется фундаментальной если выполняется Пример. ТЕОРЕМАпринцип сходимости Коши Для сходимости последовательности необходимо и достаточно чтобы она была фундаментальной. Понятие фундаментальной последовательности переносится на метрические пространства.
20718. Формула и ряд Тейлора. Биномиальный ряд 130.5 KB
  Формула и ряд Тейлора. Биномиальный ряд. Теорема о разложении функции в ряд Тейлора: пусть функция имеет в некотором интервале производные до порядка включительно а точка находится внутри этого интервала. Используя эту теорему можно сделать следующий вывод: если функция имеет на некотором отрезке производные всех порядков раз они имеются все то каждая из них будет дифференцируемой и поэтому непрерывной то можно написать формулу Тейлора для любого значения .
20720. Линейные дифференциальные уравнения второго порядка с постоянными коэффициентами 72.5 KB
  Вопрос о том является ли это решение общим приводит к понятию линейной независимости системы частных решений линейно независимых функций 1 и фундаментальной системы решений 2. Совокупность всех линейнонезависимых частных решений уравнения называется фундаментальной системой решений этого уравнения тогда есть общее решение для уравнения . Таким образом для решения нужно: найти частные решения; выяснить их линейную независимость ; найти общее решение согласно .
20721. Мощность множества. Арифметика счетной мощности 59.5 KB
  Пусть A некоторое счетное мнво тогда по определению A N.Из всякого бесконечного мнва можно выделить счетное подмново.Сумма конечного числа счетных мнв есть счетное мнво. Сумма счетного числа конечных мнв есть счетное мнво.
20722. Предел и непрерывность функции в точке. Основные свойства функции непрерывной на отрезке 29.5 KB
  Иногда говорят что предел функции в точке а : fx=b      х: ха ха и fxb Данное определение называется определением предела функции на языке .3 Если fx=fa то функция назся непрерывной в точке а.4 Если использовать предел функции в точке то определение функции в точке можно оформить в виде:    : ха х[ аb] и fxb Опред.
20723. Предел числовой последовательности. Необходимый и достаточный признак сходимости числовой последовательности 62 KB
  Определение: Если каждому по определённому закону можно поставить в соответствие то числа получающиеся при каждом конкретном n образуют числовую последовательность. Если такое имеет место то пишут что последовательность расходится. Теорема Необходимое условие сходимости числовой последовательности: если последовательность {Xn} сходится то она ограничена. Определение 2: Если предел сходящейся последовательности равен 0 то она называется бесконечно малой последовательностью.
20725. Замечательные пределы 40.5 KB
  Замечательные пределы Существует 4 замечательных предела: I. Покажем доказательство первого предела. ; ; ; ; ; ; ; по свойству функции имеющей предел имеем предел зажатой последовательности ч.