Файл: ЭО Лаб работа 2.doc

Целью лабораторной работы №2 является освоить технологию вычиcления конечных и бесконечных сумм и произведений с помощью встроенных функций sum() и product(), пределов функций и числовых последовательностей с помощью встроенной функции limit, дифференцирования и интегрирования функций, работы с матрицами и векторами, комплексными числами, дифференциальными уравнениями в программе Maple.

Вычисление сумм и произведений

Для выполнения операции суммирования в Maple используется функция sum(f,k=m..n), в которой f – выражение, зависящее от переменной суммирования k, m и n – пределы суммирования. Пределы суммирования могут быть как конечными, так и бесконечными, т.е. функция sum() может быть использована для суммирования рядов.

Например, нам необходимо найти сумму

>sum(2^n*(n+1)/n!, n=0..20)

>evalf(%)

22.16716830

Вычислим сумму ряда

> sum(n/(n+1)!, n=0..infinity)

Вычислим сумму ряда

c переменной x.

> sum(n*(n+2)*x^n, n=1..infinity)

Функция product(), выполняющая операцию умножения, аналогична функции sum(). Аргументы функции product(f, k=m..n) следующие: f – выражение, зависящее от индекса k и являющееся k-ым членом умножения, m и n – интервал изменения k.

Например, докажем что

> product(1-1/n^2, n=2..infinity)

Разложение в ряд Тейлора

Разложение функции f(x) в степенной ряд в окрестности точки а

осуществляется командой series(f(x), x=a, n), где а – точка, в окрестности которой производится разложение, n – число членов ряда.

Аналогичного действия команда taylor(f(x), x=a, n) раскладывает функции f(x) в окрестности точки x=a до порядка n-1 по формуле Тейлора.

Команды series и taylor выдают результат, имеющий тип series. Для того, чтобы иметь возможность дальнейшей работы с полученным разложением, его следует преобразовать в полином с помощью команды convert(%,polynom).

Функцию многих переменных f(x₁,…,x_n) можно разложить в ряд Тейлора по набору переменных (x₁,…,x_n) в окрестности точки (a₁,…,a_n) до порядка n с помощью команды mtaylor(f(x), [x1,…,xn], n). Эта команда находится в стандартной библиотеке, поэтому перед использованием должна быть вызвана readlib(mtaylor).

Пример. Разложить в степенной ряд в окрестностих₀=0, удерживая 5 первых членов.

> f(x)=series(exp(-x)*sqrt(x+1), x=0, 5);

Разложить в ряд Тейлора в окрестности точки (0, 0) до 6-ого порядка.

> readlib(mtaylor):

> f=mtaylor(sin(x^2+y^2), [x=0,y=0], 7);

Вычисление пределов

В Maple для некоторых математических операций существует по две команды: одна прямого, а другая – отложенного исполнения. Имена команд состоят из одинаковых букв за исключением первой: команды прямого исполнения начинаются со строчной буквы, а команды отложенного исполнения – с заглавной. После обращения к команде отложенного действия математические операции (интеграл, предел, производная и т.д.) выводятся на экран в виде стандартной аналитической записи этой операции. Вычисление в этом случае сразу не производится. Команда прямого исполнения выдает результат сразу.

Для вычисления пределов имеются две команды:

прямого исполнения – limit(expr,x=a,par), где expr – выражение, предел которого следует найти, a – значение точки, для которой вычисляется предел, par – необязательный параметр для поиска односторонних пределов (left – слева, right – справа) или указание типа переменной (real – действительная, complex – комплексная).
отложенного исполнения – Limit(expr,x=a,par), где параметры команды такие же, как и в предыдущем случае. Пример действий этих команд:

> Limit(sin(2*x)/x,x=0);

> limit(sin(2*x)/x,x=0);

С помощью этих двух команд принято записывать математические выкладки в стандартном аналитическом виде, например:

> Limit(x*(Pi/2+arctan(x)),x=-infinity)= limit(x*(Pi/2+arctan(x)), x=-infinity);

Односторонние пределы вычисляются с указанием параметров: left – для нахождения предела слева и righ – справа. Например:

> Limit(1/(1+exp(1/x)),x=0,left)= limit(1/(1+exp(1/x)),x=0,left);

> Limit(1/(1+exp(1/x)),x=0,right)=limit(1/(1+exp(1/x)), x=0,right);

Вычисление производных

Для вычисления производных в Maple имеются две команды:

прямого исполнения – diff(f,x), где f – функция, которую следует продифференцировать, x – имя переменной, по которой производится дифференцирование.
отложенного исполнения – Diff(f,x), где параметры команды такие же, как и в предыдущей. Действие этой команды сводится к аналитической записи производной в виде . После выполнения дифференцирования, полученное выражение желательно упростить. Для этого следует использовать команды simplify factor или expand, в зависимости от того, в каком виде вам нужен результат.

Пример:

> Diff(sin(x^2),x)=diff(sin(x^2),x);

Для вычисления производных старших порядков следует указать в параметрах x$n, где n – порядок производной; например:

> Diff(cos(2*x)^2,x$4)=diff(cos(2*x)^2,x$4);

Полученное выражение можно упростить двумя способами:

> simplify(%);

> combine(%);

Большинство задач дифференциального исчисления функций многих переменных решается в Maple теми же командами, что и для функций одной переменной, только с указанием дополнительных параметров.

Для вычисления частных производных функции f(x₁,…, x_m) используется уже хорошо известная вам команда diff. В этом случае эта команда имеет такой формат: diff(f,x1$n1,x2$n2,…, xm$nm), где x1,…, xm – переменные, по которым производится дифференцирование, а после знака $ указаны соответствующие порядки дифференцирования. Например, частная производная записывается в виде:diff(f,x,y).

Найти все частные производные 2-го порядка функции .

> restart; f:=(x-y)/(x+y):

> Diff(f,x$2)=simplify(diff(f,x$2));

> Diff(f,y$2)=simplify(diff(f,y$2));

> Diff(f,x,y)=diff(f,x,y);

Исследование функций

Исследование функции необходимо начинать с нахождения ее области определения, но, к сожалению, это трудно автоматизируемая операция. Поэтому при рассмотрении этого вопроса приходится решать неравенства (см. тему II). Однако, ответить на вопрос, определена ли функция на всей числовой оси, или нет, можно исследовав ее на непрерывность.

Проверить непрерывность функции f(x) на заданном промежутке [x₁,x₂] можно с помощью команды iscont(f,x=x1..x2). Если функция f непрерывна на этом интервале, то в поле вывода появится ответ true – (истина); если функция f не является непрерывной на этом интервале, то в поле вывода появится ответ false – (ложь). В частности, если задать интервал x=-infinity..+infinity, то функция f будет проверяться на всей числовой оси. В этом случае, если будет получен ответ true, то можно сказать, что функция определена и непрерывна на всей числовой оси. В противном случае следует искать точки разрыва. Это можно сделать двумя способами:

с помощью команды discont(f,x), где f – функция, исследуемая на непрерывность, x – переменная. Эта команда пригодна для нахождения точки разрыва первого и второго родов.
с помощью команды singular(f,x), где f – функция, x – переменная. Эта команда годится для нахождения точек разрыва второго рода как для вещественных значений переменной, так и для комплексных.

Обе эти команды выдают результаты в виде перечисления точек разрыва в фигурных скобках. Тип такой записи называется set. Для того, чтобы в дальнейшем можно было использовать полученные значения точек разрыва, следует из типа set с помощью команды convert перевести их в обычный числовой тип.