Лабораторные работы 1-7
.pdfЛАБОРАТОРНАЯ РАБОТА №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× 10− 15 |
|
|
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