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


 

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

24698. 34 ЗАЩИТА АД 110 KB
  Наиболее просто токовая отсечка выполняется с реле прямого действия встроенными в привод выключателя. С реле косвенного действия отсечка выполняется с независимыми токовыми реле по схемам на рис.7; Iпуск пусковой ток электродвигателя; k0TC коэффициент отстройки Токовую РЗ электродвигателей мощностью до 2000 кВт следует выполнять как правило по наиболее простой и дешевой однорелейной схеме рис. На электродвигателях мощностью 20005000 кВт токовая отсечка выполняется двухрелейной.
24699. Основные особенности выполнения РЗ на блоках 88 KB
  2 отсутствие электрической связи между генератором и сетью имеющее место в блочных схемах облегчает решение вопросов селективности РЗ генератора от замыканий на землю вследствие высокой стоимости мощных генераторов и трансформаторов повышенные требования в части чувствительности быстродействия и надежности на блоках без поперечных связей необходимость действия на останов блока в целом; На блоках малой мощности до 30 МВт включительно в качестве РЗ от внешних КЗ применяется МТЗ с комбинированным пуском по напряжению. На блоках...
24700. ЗАЩИТА РОТОРА от замыкания на корпус 63 KB
  Для периодического контроля за состоянием изоляции цепей возбуждения используется вольтметр один зажим которого соединен с землей а второй поочередно подключается к полюсам ротора. Если изоляция ротора достаточно высока замеры вольтметра в обоих случаях будут близки к нулю. Второй конец обмотки токового реле заземляется через специальную щетку имеющую электрический контакт с валом ротора.
24701. Защита ротора от перегрузки 38 KB
  Для предотвращения повреждения ротора при перегрузке предусматривается специальная РЗ а также выполняется ограничение длительности форсировки возбуждения. Наиболее полноценную РЗ ротора от перегрузки можно осуществить с помощью реле имеющего характеристику соответствующую перегрузочной характеристике ротора. Выдержка времени первой ступени при одних и тех же значениях тока ротора примерно на 20 меньше выдержки времени второй ступени.
24702. ПОВРЕЖДЕНИЯ И НЕНОРМАЛЬНЫЕ РЕЖИМЫ РАБОТЫ ГЕН-В, ОСНОВНЫЕ ТРЕБОВАНИЯ К ЗАЩИТЕ ГЕНЕРАТОРОВ 41.5 KB
  Обмотка ротора гена находится под сравнительно невысоким напряжением и поэтому ее изоляция имеет значительно больший запас элой прочности чем изоляция статорной обмотки. Однако изза значительных механических усилий обусловленных большой частотой вращения роторов турбогенов относительно часто наблюдаются случаи повреждения изоляции и замя обмотки ротора на корпус т. Замыкание на корпус в одной точке обмотки ротора неопасно так как ток в месте замыкания очень мал и нормальная работа генератора не нарушается. При двойных...
24703. Общие принципы работы реле. Работа реле на переменном токе 91.5 KB
  Общие принципы работы реле. Работа реле на переменном токе. В устройствах РЗ и электрической автоматики применяются реле на базе электромеханических конструкций полупроводниковых приборах из отдельных диодов транзисторов и др. Электромеханические реле обладают большими габаритами значительным потреблением мощности требуют тщательного ухода имеют ограниченное быстродействие и чувствительность.
24704. ИНДУКЦИОННЫЕ РЕЛЕ 220 KB
  ИНДУКЦИОННЫЕ РЕЛЕ Работа индукционных реле основана на взаимодействии переменных магнитных потоков с токами индуктированными ими в подвижной системе реле. Основными элементами реле являются два электромагнита 1 и 2 и подвижная система 3 расположенная в магнитном поле электромагнитов рис. С осью 4 жестко связан подвижный контакт реле 5 замыкающий при повороте неподвижные контакты 6. Момент Мэ приводит в движение подвижную систему 3 которая в зависимости от знака направления Мэ действует в сторону замыкания или размыкания контактов...
24705. МТЗ. Структурная и принципиальная схема 154.5 KB
  МТЗ. Селективность действия МТЗ достигается с помощью выдержки времени. МТЗ являются основным видом РЗ для сетей с односторонним питанием. Соответственно при КЗ в точке К2 быстрее всех сработает МТЗ 3.
24706. Погрешности ТН. Повреждения в цепях ТН 124.5 KB
  Повреждения в цепях ТН. Во вторичных цепях ТН могут возникать повреждения КЗ и обрывы. Для предупреждения ложных действий РЗ предусматриваются блокирующие устройства которые реагируют на появление U0 и I0 при повреждениях в цепях напряжения во вторичных цепях ТН и подают сигнал. Недостатком таких устройств блокировки является то что при КЗ в цепях фазных напряжений они не действуют.