- •Федеральное агентство по образованию
- •«Национальный исследовательский томский политехнический университет»
- •Практикум на эвм
- •«Векторные и матричные функции в пакете maple»
- •Лабораторная работа №2 векторные и матричные функции в пакете maple
- •Массивы, векторы и матрицы
- •Функции для работы с векторами и матрицами
- •Операции с векторами
- •Символьные операции с матрицами
- •Ход работы
- •Контрольные вопросы
- •Индивидуальные задания
Ход работы
Предварительно изучив функции, выполнить следующие преобрзования:
1.1.Задать
вектора с использованием шаблона.![]()
1.2. Вычислить их скалярное произведение.
1.3.
Вычислить элементы нового вектора
.
Заметим, что понятие «вектор» в пакете
аналогично определению «массив». Для
выполнения этой операции воспользуемся
оператором циклаfor.
> A:=<cb^2|az^2|2*z|2*az>;

> B:=<4|5|4|8>;

> X:=array(1..4);

> for i to 4 do X[i]:=(A[i])^3*cos(B[i])^2 od;


> evalm(X);

Вектора A и B заданы. Прежде, чем вычислить по заданному правилу новый вектор, опишем его как массив заданной размерности с не определёнными значениями элементов. В операторе цикла выражение между do и od выполняется, пока переменная i изменяется (в данном примере по умолчанию) от1 до4 с шагом 1.
Загрузить пакет расширения linalg. Ввести в 2-х режимах (используя командную строку ввода и указанную пиктограмму) матрицыCиD.
2.2Построить транспонированные матрицы.
2.3. Вычислить матрицу, обратную C для 2-х случаев:
а) с использованием в качестве элементов матрицы целых чисел;
б) с заданным хотя бы одним элемента числом с десятичной точкой (2.).
2.4. Умножить матрицу C на обратную для случаев (а) и (б) п.
2.5. Построить характеристический полином для матриц C и D.
2.6. Найти характеристические числа.
2.7. Используя функцию цикла for, построить новую матрицу, где Y[i,j]=F(C[i,j]).
Пример 1.
> restart;
> with(linalg):
> C:=array(1..3,1..3,[[5,-2,2],[1,5,1],[0,2,3]]);

> transpose(C);

> det(C);
75
> Hp:=charpoly(C,lambda);

> solve(Hp,lambda);
3,5,5
C1:=evalm(C^(-1));

C1:=evalf(evalm(C^(-1)));

> (v,e):=Eigenvectors(C);

> with(LinearAlgebra):
> W := <<-1|-3|-1>,<3|5|6>,<-3|-3|-4>>;

(v,e):=evalf(Eigenvectors(W));

> (v,e):=Eigenvectors(W);

Пример2. Формирование новой матрицы по заданному правилу, с использованием вложенных циклов.
> Y:=Matrix(1..3,1..3);

> for i to 3 do for j to 3 do Y[i,j]:=i^2*sin(C[i,j]+j) od od;
> evalm(Y);

Контрольные вопросы
1. Чем отличаются разные способы задания векторов и матриц?
2. Назначение пакетов расширения?
3. Чем объяснить различие результатов п.2.4?
4. Как организовать с помощью изученных функций вложенные циклы?
Индивидуальные задания
1.
.
;
;
.
2.
.
;
;
.
3.
.
;
;

4.
.
;
;
.
5.
.
;
;
.
6.
.
;
;

7.
.
;
;

8.
.
;
;

9.
.
;
;

10.
.
;
;

11.
.
;
;
.
12.
.
;
;
.
13.
.
;
;
.
14.
.
;
;

15.
.
;
;

16.
.
;
;
.
17.
.
;
;
.
18.
.
;
;
.
19.
.
;
;
.
20.
.
;
;
.
21
.
;
;
.
22.
.
;
;
.
23.
.
;
;
.
24.
.
;
;
.
25.
.
;
;
Векторные и матричные функции в пакете MAPLE
Методические указания
по выполнению лабораторной работы
Составил:
Козловских Александр Владимирович
Подписано к печати
Формат 60х84.16. Бумага ксероксная.
Плоская печать. Усл. печ.л.. Уч.-изд. л
Тираж зкз. Заказ . Цена свободная.
ИПФ ТПУ. Лицензия ЛТ №1 от 18.07.2011.
Типография ТПУ. 634034, Томск, пр. Ленина, 30.
