 |
Дифференцирование выражений
|
|
|
|
Вычисление производных функций n-го порядка -одна из самых распространенных задач математического анализа. Для ее реализации Maple V имеет следующие основные функции:
| diff(F,x1,x2,...,xn) | Diff(F,x1,x2,...,xn) |
| diff(F,[x1,x2,...,xn]) | Diff(F,[x1,x2,...,xn]) |
Здесь F -дифференцируемое алгебраическое выражение, в частности, функция F(x1,x2,...,xn) ряда переменных, по которым производится дифференцирование.
Функция Diff является инертной формой вычисляемой функции diff и может использоваться для естественного воспроизведения производных в документах.
Первая из этих функций (в вычисляемой и в инертной форме) вычисляет частные производные для выражения F по переменным x1,x2,...,xn. В простейшем случае diff(F(x),x) вычисляет первую производную функции F(x) по переменной x. При n большем 1 вычисления производных выполняются рекурсивно, например diff(f(x,y),x,y) эквивалентно diff(diff(f(x,y),x),y).
| Оператор $ можно использовать для вычисления производных высокого порядка. Для этого после имени соответствующей переменной ставится этот оператор и указывается порядок производной, например diff(f(x),x$4). |
Дифференциальный оператор
Для создания функции с производными может также использоваться дифференциальный оператор D. Порою он позволяет создавать более компактные выражения, чем функции diff и Diff.
| Дифференциальный оператор можно записывать в следующих формах: D(f) или D[i](f), где параметр f-выражение или имя функции. i-положительное целое число, выражение или последовательность. |
| Оператор D(f) просто вычисляет имя производной от f. В форме D(f)(x) этот оператор подобен diff(f(x),x). |
Поиск экстремумов, анализ функций не непрерывность
Фуекция extrema() позволяет найти экстремумы выражения или функции, а также вывести координаты точек экстремумов. Её синтаксис
|
|
|
Extrema(expr, условия, переменная, 'nv');
expr - выражение, экстремум которого следует найти,
условия - ограничения (по области и др.) задаются в виде множества уравнений или неравенств,
переменные - переменные, по которым ищется экстремум, задаются в виде множества (необязательный параметр)
'nv' - имя переменной, ккоторой будут присваиваться координаты точек экстремума (необязательный параметр).
Для поиска максимального и минимального значений функции или выражения в Maple предназначены функции maximize () и minimize(). Их синтаксис:
minimize (expr, vars, опция);
maximize (expr, vars, опция);
expr - выражение, экстремум которого следует найти,
vars - список или множество переиенных
опция - параметр, указывающий область изменения независимых переменных выражения.
-
Параметр опция может принимать значения: 'infinity' - поиск максимума или минимума выполняется по всей числовой прямой, или представляет собой множество уравнений вида x=a..b, определяющих область изменения независимых переменных.
Функция iscont() позволяет исследовать выражение (функцию) на непрерывность. Синтаксис этой команды:
iscont(expr, x=a..b)
iscont(expr, x=a..b,'closed')
expr - выражение, исследуемое на непрерывность,
x - независимая переменная
a..b - промежуток изменения независимой переменной
'closed' - параметр, указывающий, что конечные точки также должны проверяться.
Если выражение непрерывно, возвращается логическоезначение true, иначе - false.
Другая функция discont() позволяет определить точки, в которых нарушается непрерывность функции f(x). Она вычисляет все точки в пределах изменения x от минус бесконечности до плюс бесконечности. Обращение к этой команде:
Discont(f(x),x)
Многие операции, например, интегрирования и дифференцирования, чувствительны к особенностям функций, в частности, к ее разрывам и особым точкам. Команда singular(expr, vars) позволяет найти особые (сингулярные) точки выражения expr, в которых оно испытывает разрывы.
|
|
|
Вычисление интегралов
Для вычисления неопределенных и определеных интегралов Maple V предоставляет следующие функции:
| int(f,x) | int(f,x=a..b) | int(f,x=a..b,опция) |
| Int(f,x) | Int(f,x=a..b) | Int(f,x=a..b,опция) |
Здесь: f -подынтегральная функция,
x -переменная, по которой выполняются вычисления,
a и b - нижний и верхний пределы интегрирования (эти пределы могут быть и аналитическими выражениями, и бесконечностью, обозначаемой как infinity)
опция -опционально заданное условие.
По умолчанию команда интегрирования проверяет выражение на непрерывность в области интегрирования и вычисляет интеграл как сумму отдельных определенных интегралов на промежутках непрерывности функции. Опция continuous отключает этот режим и вычисляет интеграл как разность значений первообразной подынтегральной функции в точке начала и конца промежутка интегрирования.
Еще одна опция CauchyPrincipalValue вычисляет несобственные интегралы первого и второго рада в смысле главного значения Коши.
Если Maple не может найти замкнутую форму выражения для определенного интеграла, то команда интегрирования в области вывода повторяет саму себя. В подобных случаях можно вычислить значение определенного интеграла численным способом с помощью команды evalf(). Синтаксис подобной конструкции следующий:
evalf(int(f,x=a..b));
evalf(Int(f,x=a..b));
evalf(Int(f,x=a..b, digits, flag));
где digits задает число значащих цифр при вычислении приближенного значения интеграла,
flag - параметр, позволяющий явно задать метод численного интегрирования.
| При численном интегрировании по умолчанию используется квадратурная формула Кленшо-Куртиса (Clenshaw-Curtis). Если в подынтегральном выражении встречается сингулярность, то применяется специальная методика символьного анализа для ее разрешения. Для задач с неустранимыми сингулярностями используется адаптивный метод двойных экспоненциальных квадратур. Параметр flag, явно задающий метод численного интегрирования, может принимать значения, представленные в таблице: |
| Значения параметра flag. | |
| Значение | Смысл |
| _Ccquad | Применяется только квадратура Кленшо-Куртиса без вызова процедуры обработки сингулярности. |
| _Dexp | Применяется адаптивный метод двойных экспоненциальных квадратур. |
| _Ncrule | Применяется метод квадратурной формулы Ньютона-Кортеса, являющийся методом фиксированного порядка, и не эффективен для высоких точностей. |
|
|
|
| Функции int и Int могут использоваться для вычисления кратных интегралов. Для этого они используются неоднократно. В пакете student имеются дополнительные функции интегрирования, которые дополняют уже описанные возможности. |
Разложение функции в ряд
Система Maple предлагает множество функций для разложений функций во всевозможные ряды. Познакомимся с некоторыми из них.
Команда series() дает разложение функции в форме степенного многочлена в окрестности некоторой точки. Её синтаксис:
series(выражение, условие)
series(выражение, условие, n)
| выражение - разлагаемое выражение, условие - уравнение вида x=a или имя переменной (например, x) n - неотрицательное число, задающее число членов ряда (по умолчанию = 6, можно переустановить системной переменной Order) |
Если в качестве параметра условие задано имя перемнной, то будет построено разложение функции в рад по этой переменной в области точки с ее нулевым значением. Задав условие ввиде равенства x=a, будем получать разложение функции по переменной x в окрестности точки x=a.
Остаточная поргешность при разложении функции в степенной ряд задается членом вида O(xn). При точном разложении этот член отсутствует. В общем случае для его удаления можно использовать функцию convert(выражение, polynom).
Для разложения функции в широко распространенный ряд Тейлора используется команда taylor(). Ее синтаксис:
|
|
|
