Добавил:
Upload Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

Лабораторные работы 1-7

.pdf
Скачиваний:
16
Добавлен:
27.02.2016
Размер:
610.74 Кб
Скачать

ЛАБОРАТОРНАЯ РАБОТА №2 Тема: Работа с массивами данных

Цели работы:

1.Научиться создавать массивы и производить с ними вычисления.

2.Научиться использовать векторные и матричные операторы и функции.

2.1. Создание и изменение массива

Матрицей А называется таблица вида

a11

a21

А= ...

ai1

...

am1

a12

...

a1 j

...

a1n

 

a

...

a

2 j

...

a

 

 

 

22

 

 

 

2n

 

...

... ...

...

...

 

,

a

...

a

...

a

 

 

i2

 

 

ij

 

in

 

...

... ...

...

...

 

 

a

 

...

a

 

...

a

 

 

m2

mj

 

 

 

 

 

 

mn

 

содержащая т строк и п столбцов. В этом случае говорят, что матрица А размерности m ×n .

аij - это элемент матрицы А, стоящий в строке с номером i и столбце с номером j.

Элементами матрицы могут быть действительные числа или другие математические объекты.

Матрица А называется квадратной, если в ней число строк равно числу столбцов (m=n).

п- мерным вектором V называется упорядоченная совокупность п действительных чисел V (x1, x2 ,..., xn ), где x1, x2 ,..., xn - координаты вектора.

Общее название для матрицы и вектора – массив.

В MathCAD для создания массивов и выполнения действий с ними удобно использовать палитру «Матрицы».

14

Для создания или изменения массива необходимо вызвать окно создания массива. На палитре «Матрицы» ему отвечает крайняя левая кнопка в верхнем ряду. Кроме того, окно создания массива можно вызвать через главное меню: Insert –Matrix.

В появившемся окне: Rows – число строк; Columns – число столбцов; Insert – вставить;

Delete – удалить; Close – закрыть.

Для создания массива в выбранном месте рабочего документа достаточно в окне создания массива указать число строк и число столбцов и нажать ОК или Insert. При заполнении появившегося шаблона для перемещения между элементами используются клавиши «вверх», «вниз», «вправо», «влево» или клавиша табуляции.

Имеющийся массив можно изменить, добавив или удалив строки и столбцы.

Для изменения массива необходимо:

-выбрать один из элементов массива, щелкнув по нему мышью;

-вызвать окно создания массива;

-указать число строк и число столбцов, которое необходимо добавить или удалить;

-нажать Insert - для добавления, Delete – для удаления строк и столбцов. Обратите внимание, строки и столбцы добавляются и удаляются ниже и

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

15

Пример 1.

1) Создадим массив размерности 3×3 и заполним его произвольными числами.

1

2

3

4

5

6

 

 

 

7

8

9

2) Встанем на элемент a21 , в окне создания массива укажем 2 строки и 1

столбец и нажмем Insert. Получим

1

 

2

3

 

 

 

4

 

5

6

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

7

 

8

9

 

 

Заполним свободные места произвольными числами:

1

1

2

3

 

 

4

2

5

6

 

 

 

3

4

5

6

 

7

8

9

10

 

 

 

11

8

9

 

7

 

3) Встанем на элемент a43 , в окне создания массива укажем 1 строку и 2

столбца и нажмем Delete. Получим

1

1

 

 

4

2

 

 

 

 

3

4

 

 

 

7

11

Массивам можно задавать имена:

A :=

.....

 

, и в дальнейшем использовать

 

 

.....

 

 

имя массива, а не вводить каждый раз массив заново.

Пример 2. Зададим три массива: А – размерности 4×5 ; В и С – размерности 4×4 .

 

 

1

2

3

4

5

 

 

1

2

1

2

 

 

5

6

5

6

 

 

2

1

5

8

7

 

 

 

3

4

3

4

 

 

 

2

4

6

8

 

A :=

 

 

B :=

 

 

C :=

 

 

 

9

4

6

8

2

 

 

5

6

7

8

 

 

0

5

1

4

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

2

8

9

7

 

1

3

5

7

 

9

3

1

0

16

2.2.Действия с массивами

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

Знаки операций с массивами можно набирать вручную, используя соответствующие клавиши или их сочетания на клавиатуре (см. таблицу 2), либо при помощи палитр «Матрицы» и «Калькулятор» .

Вприведенной таблице: А, В – произвольные матрицы, V, U п-мерные векторы.

 

 

 

 

 

 

 

 

 

 

 

Таблица 2

 

 

 

 

 

 

 

 

 

 

 

 

Операция

Клавиши

 

 

 

 

 

Описание

 

 

 

 

 

 

 

 

 

 

A+B; V+U

+

Сумма двух матриц или двух векторов.

 

 

 

 

 

 

 

 

 

А-В;

V-U

Разность двух матриц или двух векторов.

 

 

 

 

 

 

 

 

 

 

 

α A ;

α V

*

Произведение

матрицы или

вектора на

число.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Произведение двух матриц. (Число

A B

*

столбцов матрицы A должно быть равно

 

 

 

 

 

числу строк матрицы B.)

 

 

 

 

 

 

 

 

 

 

 

 

 

A

 

Вычисление

определителя

квадратной

 

 

матрицы.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Возведение квадратной матрицы в степень

 

An

^

п (n- целое число). При n=-1 получим

 

 

 

 

 

обратную матрицу.

 

 

 

 

 

 

 

 

 

 

 

 

AT

Ctrl+1

Транспонирование матрицы.

 

 

 

 

 

 

 

 

 

V

 

 

Вычисление длины вектора:

 

 

 

 

 

V

 

= x 2

+ x 2 +...+ x 2 .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

2

n

 

 

 

 

 

 

 

 

 

 

 

 

 

17

 

 

 

V

Ctrl+4

Суммирование элементов вектора.

 

 

 

U V

*

Скалярное произведение векторов.

 

 

 

U ×V

Ctrl+8

Векторное произведение двух трехмерных

векторов.

 

 

 

A<k >

Ctrl+6

Извлечение из матрицы А столбца с

номером k.

 

 

 

Ai, j

[

Извлечение из матрицы А элемента aij .

Vi

[

Извлечение из вектора V элемента xi .

 

 

 

Следует помнить, что нумерация элементов массива начинается с нуля. Чтобы перейти к привычной нумерации, необходимо изменить значение встроенной константы ORIGIN, присвоив ей значение 1.

Пример 3.

1) Выполним действия с заданными ранее массивами А, В, С.

 

 

 

 

 

5

10

15

20

25

 

 

 

 

 

13

14

13

14

 

 

27

25

20

26

5 A =

10

5

25

40

35

 

2 B 3 C =

 

0

4

12

16

B C =

59

61

46

62

 

20

30

40

10

 

 

 

 

3

11

4

 

 

113

76

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

45

 

 

 

 

 

 

10

 

 

 

109

106

 

 

 

 

 

0

10

40

45

35

 

 

 

 

 

25

3

7

14

 

 

74

64

35

50

 

 

 

 

14

22

24

32

 

 

 

 

 

1.531×

3

 

 

3

1.207×

3

 

3

 

 

 

 

 

 

 

 

 

 

 

10

1.588× 10

10

1.61× 10

 

 

2

=

34

52

56

74

 

3

 

1.226× 103

1.312× 103

1.058× 103

1.388× 103

 

B

 

66

100

112

 

 

 

 

C

=

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

700

 

 

731

 

 

531

 

712

 

 

 

 

 

 

 

146

 

 

 

 

 

 

 

 

 

 

 

 

 

 

42

65

80

103

 

 

 

 

 

1.135×

3

 

 

3

 

827

 

 

3

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

10

1.156× 10

 

 

1.13× 10

 

 

 

B

 

=

1.776× 1015

 

 

C

 

= −276

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2) Поскольку для матрицы С определитель отличен от нуля, найдем для нее обратную матрицу и сделаем проверку.

 

 

0.261

0.174

0.043

0.217

 

 

 

1

0

0

0

 

 

1

0

0

0

1

=

0.493

0.384

0.029

0.188

1

=

0

1

0

0

1

C =

0

1

0

0

C

 

0.413

0.478

 

 

C C

 

0

1

 

C

 

0

1

 

 

 

0.87

0.391

 

 

0

0

 

 

0

0

 

 

0.833

0.583

0.333

0.333

 

 

 

0

0

0

1

 

 

0

0

0

1

18

3)Извлечем из матрицы А столбец с номером 4:

5

A 4 = 72

7

Поскольку константа ORIGIN не была переопределена, то извлекся последний столбец. Присвоим ORIGIN значение 1 и снова извлечем из матрицы А столбец с номером 4:

ORIGIN:= 1

4

A 4 = 88

9

4) Зададим два вектора: V- равный третьему столбцу матрицы А, U- равный второй строке матрицы В; и выполним с ними действия.

V := A 3

U := (BT)2

 

 

3

 

 

 

3

 

 

 

5

 

 

 

4

 

V =

 

 

U =

 

 

 

6

 

 

3

 

 

 

 

 

 

 

 

 

 

8

 

 

4

 

 

 

 

 

 

9

 

 

 

 

 

6

 

 

 

15

 

 

 

 

 

 

 

15

 

 

 

 

 

 

8

 

 

 

21

 

 

3 V =

 

 

 

 

 

2 U =

 

 

V + 4 U =

 

 

 

 

18

 

 

 

 

 

6

 

 

18

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

24

 

 

 

 

 

8

 

 

24

 

 

V

 

= 11.576

 

U

 

= 7.071

 

 

V = 22

U = 14

V U = 79

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

5) Извлечем из заданных матриц и векторов некоторые элементы.

B1,2 = 2 A3,4 = 8 B2,2 = 4 U4 = 4 V2 = 5

19

2.3. Векторные и матричные функции

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

Напомним, что список встроенных функций можно вызвать через главное меню (Insert - Function) или используя кнопку с надписью f(x) на панели инструментов.

Рассмотрим основные матричные и векторные функции: rows (A) - определяет количество строк в массиве А; cols (A) - определяет количество столбцов в массиве А;

length (V) - выдает размер вектора V (число элементов в векторе); max (A) - выдает максимальный элемент массива;

min (A) - выдает минимальный элемент массива; rank (A) - определяет ранг матрицы А;

augment (A,B) – создает новый массив, располагая массив А слева от массива В (А и В должны содержать одинаковое число строк);

stack (A,B) – создает новый массив, располагая массив А над массивом В (А и В должны содержать одинаковое число столбцов);

submatrix (A,ir,jr,ic,jc) – создает подматрицу матрицы А, состоящую из элементов, расположенных в строках с номерами от ir до jr и столбцах с номерами от ic до jc.

Пример 4. Применим рассмотренные функции к заданным ранее матрицам и векторам.

rows(A) = 4

cols(A) = 5 rows(B) = 4

cols(B) = 4

length (U) = 4

length (V) = 4

min(A) = 0

max(A) = 9

 

min(V) = 3

max(V) = 8

rank(A) = 4

rank(B) = 3

rank(C) = 4

 

 

20

 

1 2

3

4

5

1

2

1

2

 

1 2 1

2

3

augment(A ,B) =

2 1

5

8

7

3

4

3

4

augment(B,U) =

3 4 3

4

4

 

4

6

8

2

5

6

7

 

 

6

7

8

 

 

 

 

9

8

 

5

3

 

0

2

8

9

7

1

3

5

7

 

1

3

5

7

4

 

 

 

(

 

 

 

 

 

 

)

1

1

5

3

3

 

 

 

 

 

 

augment

A

1

 

1

1

2 3 2 5 4

 

 

 

 

 

 

 

 

,B

 

,C

,V,U

=

9

5

0

6

3

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

1

9

8

4

 

 

 

 

 

 

1

2

1

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

3

4

3

4

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

5

 

6

5

6

 

 

5

6

7

8

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

4

6

8

 

 

 

 

1

3

5

7

 

 

 

stack(C,VT)

 

 

 

 

 

stack(B,C) =

 

 

 

=

0

 

5

1

4

 

 

 

5

6

5

6

 

 

 

 

 

 

 

 

 

9

 

3

1

0

 

 

 

2

4

6

8

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0 5 1

4

 

 

 

 

 

 

 

 

 

3 5 6 8

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

9

3

1

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

1

5

 

 

 

 

 

 

 

 

 

2

1

2

submatrix(A ,2,4,1,3) =

 

9

4

6

submatrix(B,1,3,2,4) = 4

3

4

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

2

8

 

 

 

 

 

 

 

 

 

6

7

8

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

1

5

2

1

2

 

augment(submatrix(A ,2,4,1,3) ,submatrix(B,1,3,2,4))

= 9

4

6

4

3

4

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

2

8

6

7

8

 

2.4.Задания для самостоятельного решения

1.Задайте три массива: А – размерности 4×5 ; В и С – размерности 4×4 .

2. Вычислите: 6A , 4C 5B , AT C , B A , C B , B C , 2C3 5B2 , B , C .

3.Для тех матриц, для которых возможно, найдите обратные матрицы и сделайте проверку.

4.Задайте два вектора: V- равный второму столбцу матрицы А, U-равный третьей строки матрицы С.

5.Вычислите: 3V , 5U , 5V +3U , V , U , V , U , U V , V U .

21

6. Найдите значения выражений:

а) (A4,5 + B4,4 +C3,3 )2 (V1 )2 (U3 )2 ; б)

B2,3 A1,2 +U4

.

 

 

(V3 )2 +(C3,4 )

Сделайте проверку.

 

 

7.Найдите максимальные и минимальные элементы во всех массивах, вычислите ранги матриц А, В, С и длины векторов V и U.

8.Создайте новые массивы:

-объединив матрицы С и А, матрицы В и С, матрицу А с векторами V и U, векторы V , U и последние столбцы матриц А, В, С;

-расположив матрицу С над матрицей В;

-добавив к матрице В вектор U в качестве первой строки;

-выделив из матрицы С подматрицу, состоящую из второй и третей строк, и второго, третьего и четвертого столбцов;

-выделив из матрицы А произвольную подматрицу размерности 3×2 .

22

ЛАБОРАТОРНАЯ РАБОТА №3

Тема: Решение систем алгебраических уравнений

Цели работы:

1.Рассмотреть способы решения систем уравнений.

2.Рассмотреть способы решения систем линейных уравнений.

3.1. Основные теоретические положения

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

F1 (x1 , х2 ,........, xn )= 0,

F

(x

, х

2

,........, x

n

)= 0,

 

2

1

 

 

 

 

 

....................................,

 

 

(x

 

 

 

 

 

 

 

 

)= 0.

F

 

, х

2

,........, x

n

 

т

1

 

 

 

 

 

Решением системы уравнений называется любой набор значений x1 , x2 ,....., xn , при подстановке которых уравнения системы обращаются в верные равенства.

Системой линейных уравнений называется система вида

a11 x1 + a12 x2

+..... + a1n xn

= b1,

 

 

 

 

 

 

 

 

 

a21 x1 + a22 x2

 

+..... + a2n xn

 

= b2,

 

 

 

 

 

 

 

 

 

...................................................

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

a

x + a

m2

x

2

+..... + a

mn

x

n

= b

 

 

 

 

 

 

 

 

 

 

m1 1

 

 

 

 

 

m,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

a ...

 

a

b

 

b

 

 

x

 

 

 

 

 

 

 

11

 

 

 

 

1n

1

 

 

или в матричной форме A X = B , где A = a12 ...

 

a2n

b2

,

 

1

 

,

 

1

 

 

B =

...

X =

... .

 

 

 

 

 

... ...

 

...

...

 

 

 

 

 

 

 

 

 

 

 

 

 

a

m1

...

 

a

mn

b

 

bn

 

 

xn

 

 

 

 

 

 

 

 

 

 

 

n

 

 

 

 

 

 

 

 

Система, не имеющая ни одного решения, называется несовместной. Система, обладающая хотя бы одним решением, называется совместной.

23