ЧИСЛЕННОЕ ДИФФЕРЕНЦИРОВАНИЕ

Для численного дифференцирования MathCAD применяет довольно сложный алгоритм, вычисляющий производную с колоссальной точностью до 7-8-го знака после запятой. Этот алгоритм (метод Риддера) описан во встроенной справочной системе MathCAD, доступной через меню Help (Справка). Погрешность дифференцирования не зависит от констант TOL или CTOL, в противоположность большинству остальных численных методов, а определяется непосредственно алгоритмом. Исключение составляют функции, которые дифференцируются в окрестности сингулярной точки; например функции f(x)=l/x это будут точки вблизи х=о. При попытке найти ее производную при х=о будет выдано сообщение об одной из ошибок деления на ноль "Can't divide by zero" (Деление на ноль невозможно) или "Found a singularity while evaluating this expression. You may be dividing by zero" (Найдена сингулярность при вычислении этого выражения. Возможно, вы делите на ноль) Если попробовать численно определить производную очень близко к нулю, например, при х=10100, то может появиться сообщение об ошибке "Can't converge to a solution" (Невозможно найти решение). Встретившись с одной из упомянутых ошибок, присмотритесь повнимательнее к дифференцируемой функции и убедитесь, что вы не имеете дело с точкой сингулярности.

В заключении этой главы хочу привести пример ТИПИЧНОЙ ОШИБКИ:

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

Его первые две строки вычисляют производную cos(x) в точке х=0.5. Последняя строка демонстрирует неправильное применение оператора дифференцирования. Вместо вычисления производной cos(x) в той же точке, как этого можно было ожидать, получено нулевое значение. Это случилось потому, что аргумент функции cos(x) введен не в виде переменной х, а в виде числа. Поэтому MathCAD воспринимает последнюю строку как вычисление сначала значения косинуса в точке х=0. 5, а затем дифференцирование этого значения (т. е. константы) также в точке х=0.5, в соответствии с требованием первой строки. Поэтому ответ, на самом деле, неудивителен - в какой точке ни дифференцируй константу, результатом будет ноль.

ГЛАВА III

ИССЛЕДОВАНИЕ ФУНКЦИЙ ОДНОЙ ПЕРЕМЕННОЙ. ПОСТРОЕНИЕ ГРАФИКОВ.

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

В чём же состоит основная задача исследования функции? Как я считаю, что эта задача является выполненной, когда по найденным аналитическим путём данным (максимумов, минимумов, промежутков возрастания, и т. д.), мы знаем поведение функции на всех промежутках интервала, на котором она определена, т. е., можем построить её график. Поэтому, прежде чем пытаться исследовать функцию на MathCADе, советую на первое время проанализировав её вручную и построив график, сравнить полученные результаты с графиком, построенным MathCAD. Если видно несоответствие, следует искать ошибку (прежде всего у себя, т. к. вероятность ошибки MathCAD мала – он строит графики, разлагая функции в ряд и аппроксимируя их).

В принципе, в MathCAD, для того, чтобы построить график, необязательно исследовать функцию. Можно всё сделать буквально нажатиями трёх кнопок (описано ниже). Но это не интересно. Гораздо интереснее, сначала проанализировав функцию (возможностями программы), прикинуть её график, а потом построить его, и посмотреть, всё ли в порядке.

ПЛАН ИССЛЕДОВАНИЯ ФУНКЦИИ

Сразу отмечу, что некоторые пункты плана можно опускать. Например, зачем у степенной функции у=а^х искать период?

I.

1) область определения функции;

2) область значения функции;

3) чётность – нечётность;

4) периодичность;

5) асимптоты;

6) точки пересечения графика с осями координат;

II. Исследование функции по 1-ой производной:

1) участки монотонности (возрастания, убывания);

2) максимумы, минимумы;

III. Исследование функции с помощью второй производной:

выпуклости, вогнутости, точки перегиба.

Поясню некоторые пункты плана. С областью определения, значения, чётностью, периодичностью, точками пересечения графика с осями координат, надеюсь, всё понятно (кому нет – обратитесь за помощью к учебным пособиям по мат. анализу, благо таковых куча).

Затруднения, на мой взгляд, из I части исследования могут вызвать асимптоты.

Асимптотой кривойназывается прямая, обладающая тем свойством, что расстояние между переменной точкой кривой и этой прямой стремиться к нулю при неограниченном удалении точки на кривой от начала координат.

Асимптоты разделяются на вертикальные и наклонные (горизонтальные асимптоты есть частный случай наклонных асимптот).

Откуда приходим к ПРАВИЛУ 1: для отыскания вертикальных асимптот, нужно найти те значения х, при которых функция обращается в бесконечность (lim y(x) à ). Их может быть бесконечно много.

ПРАВИЛО 2. Уравнение наклонной асимптоты: у=kх+b, где , .

Поиск участков монотонности, максимумов и минимумов тоже не должен вызывать затруднений (всё это должно быть известно со школы). - Ищем критические точки (точки в которых производная равна нулю или не существует) – они будут у нас «подозреваемыми» в том, что являются максимами или минимами. Затем смотрим, меняет ли знак производная в окрестности этих точек.

Наибольшее затруднение вызовет, наверное, исследование функции по второй производной – оно позволяет нам найти вогнутости, выпуклости, точки перегиба кривой и подтвердить, найденные с помощью первой производной максимумы и минимумы.

Кривая называется выпуклой на интервале, если она расположена ниже любой своей касательной.

Кривая называется вогнутой на интервале, если она расположена выше любой своей касательной.

Необходимое условие точки перегиба (точки, отделяющей выпуклую часть от вогнутой):Если х=с - точка перегиба, то f ’’(с) либо равно нулю, либо не существует.

Достаточное условие точки перегиба:Если вторая производная при переходе через точку меняет знак – то эта точка есть точка перегиба.

Если f ’’ (х)>0, то кривая вогнутая на интервале, если f ’’ (х)<0 – выпуклая.

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

Пусть нам требуется найти наибольшее и наименьшее значение функции на отрезке [1, 5] (пример 3.15) . Для этого необходимо найти максимумы и минимумы функции (или хотя бы просто критические точки, не проверяя их на максимум или минимум) на этом промежутке (обращаю внимание, что максимум не означает максимальное возможное значение, как минимум и не означает минимально возможное значение), и значения функции на концах промежутка. На языке MathCAD это будет выглядеть так:

1) Определяем функцию: .

2) Ищем производную (подробно описано в главе дифференцирование). Для этого можно просто нажать на значок оператора производная из панели CALCULUS: . В метку, обведённую кружочком можно вводить только символ (т. к. непосредственно сам вид функции мы определили заранее. Нажимаем CTRL + ю ( в лат регистре) и получаем ответ: . Упрощать не будем – это в данном случае не имеет смысла

3) Скопируем в буфер обмена полученный результат (чтобы не вводить его вручную с клавиатуры), для чего нажимаем левой кнопкой мыши на конце выражения (подчёркнуты места нажатий для копирования – можно нажимать в любом месте) и, перемещая её с нажатой кнопкой (либо вправо либо влево), выделим кусок, необходимый для копирования. Затем выбираем из главного меню пункт «редактирование», а из него пункт «копировать».

4) Набираем служебное слово Given (используется для решения уравнений и т. п.) и нажимаем ENTER.

5) Если метка, отмеченная красным крестиком, оказалась выше служебного слова, перетащим её под слово Given, для чего следует просто щелкнуть мышью под ним.

6) Выбираем пункт «редактирование» - «вставка». Наша производная должна появиться в месте, где была красная метка.

7) «Цепляем» синей меткой всё выражение, нажимая пробел (подробно описано в примерах выше).

8) Из панели инструментов Булевая ( BOOLEAN) выбираем кнопку с символом = или нажимаем CTRL + =. Равно должно появиться жирное.

ВНИМАНИЕ. Не вводите численное равно с панели Evaluation (или просто = с клавиатуры) и равно - определитель ( ) – будет сообщение об ошибке!

9) Вводим 0.

Математическая справка: по существу мы решаем уравнение y’=0. Т.е. ищем критические точки функции.

10) Набираем на клавиатуре служебное слово Find(x) (найти что-то), где х – это та переменная, значения которой нам нужно найти. Ставим à.

11) Нажимаем ENTER и получаем ответ: . Переведём этот результат на обычный математический язык. Это значит, что в точке х=4 производная равна нулю, следовательно это и есть критическая точка. Но существует и ещё одна критическая точка – х=2, т.к., вспомнив определение критических точек, найдём, что производная при х=2 не существует – знаменатель обращается в ноль при х=2. Следовательно, х=2 – тоже критическая точка.

12) Ищем значение функции на концах промежутка и в критических точках: поочерёдно вводим y( x), где вместо х подставляем абсциссу, от которой хотим найти значение функции и нажимаем численное равно (просто нажимаем равно на клавиатуре =) Т. е.: .

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

Для убедительности привожу график этой функции (подробно построение графиков будет разобрано в следующем примере):

Из графика сразу же становится всё понятно. Например, производная в х=2 не существует, т. к там острый пик у функции, и, следовательно, нельзя провести касательную в данной точке.

Рассмотрим теперь подробно построение графиков функции и асимптот к ним на примере 6.15.

1. Определим функцию

2. Ищем вертикальные асимптоты ( те х, при которых знаменатель обращается в ноль и функция не существует: вводим . (см. выше).

3. Вводим и сразу получаем ответ. Следовательно, уравнения вертикальных (их две) будут выглядеть так: .

4. Ищем наклонные асимптоты. Коэффициент k есть , где вместо f(x) программа автоматически подставляет определенную в пункте 1 нами функцию и считает предел (о вычислении пределов см. гл.I). Коэффициент b есть .

5. Теперь определим функцию наклонной асимптоты: (у=kx+b)- она потребуется нам для построения графиков асимптот. Вместо а можно ввести любую букву лат. алфавита. Не забывайте в скобках вводить независимую переменную.

6. Строим график:

a) нажимаем на панели Graf (графики) кнопку «декартов график»:

b) в появившемся окне: начинаем заполнять метки: В метку, обведённую кружочком, вводим аргумент первой функции ( ), т. е. х. В метку, обведённую квадратом вводим f(x). (Если бы мы не определили функцию заранее, то пришлось бы вводить в метку сам вид функции сам вид функции, т.е. , что громоздко и неудобно). Перемещение между метками осуществляется нажатием TAB.

c) Ставим запятую после f(x) (запятая не видна) нажимая на б .

d) В появившуюся ниже метку вводим a(x).

e) Переместившись на нижнюю метку (надо нажать на х левой кнопкой мыши), ставим там запятую, и вводим ещё один х.

f) Возвращаемся на вертикальную метку (нажимаем на a(x)) и ставим запятую.

g) Вводим х, но уже на вертикальной оси

h) Возвращаемся на горизонтальную ось, нажимая на х левой кнопкой мыши.

i) Ставим запятую и вводим .

j) Возвращаемся на вертикальную ось и вводим там ещё один х

k) На горизонтальной оси вводим .

Пояснения: лучше всего( чтобы избежать возможных ошибок) вводить столько аргументов на оси Ох, сколько функций будет у вас на оси Оу. В данном случае, т. к. помимо обычных линий необходимо построить вертикальные линии х=а, нужно чтобы на оси Ох стояла а (в нашем случае в роли а и ), а по оси у стояла х. Честно признаться, я потратил часа два, на то, чтобы догадаться как построить вертикальную линию – в Центре Ресурсов и в Помощи это нигде не описано. Здесь видна явная недоработка разработчиков MathCAD.

7) Итоговый график будет выглядеть так: Надо признаться, совсем неказисто. Дело в том, что MathCAD. Изменим их:

a) дважды щелкаем левой кнопкой мыши в любом месте на графике (в прямоугольном окне или щелкаем в окне правой кнопкой мыши и в появившемся всплывающем меню выбираем пункт Формат (FORMAT))

b) Появится одно из следующих окон (в зависимости от того, является ли ваша версия русской или нет). Ниже приведено соответствие между двумя окнами:

Далее по мере необходимости я буду в скобках дублировать английскую команду.

c) из вкладки х-у оси (x-y axes) в блоке «стиль осей графика» (Axes Style) выбираем пункт «пересечение» (crossed) и нажимаем кнопку «применить» (наблюдаем изменения, не закрывая диалог);

d) теперь закроем окно и изменим масштаб по оси Оу. Для этого нажмём мышью в окне где построен график. Появится рамочка с синей меткой:

e) в метках, обведённых кружочком, вводим удобный нам в данном случае диапазон изменения величины y (в нашем примере путём подгонки я нашёл, что это будет промежуток [-10; 10]). Обычно MathCAD автоматически расставляет пределы, позволяющие «понять» функцию, но в данном случае это у него не получилось

f) если хотите, можете изменить пределы по оси x, но в данном примере программа выбрала их оптимально

g) изменим трассировки (представление, цвет и жирность) линий. Для чего дважды щёлкаем левой кнопкой мыши по окну, где построен график, - для открытия окна форматирования графика (п. b). Выбираем вкладку «трассировка» (Traces). В первой колонке (под названием Legend Label) по существу указаны порядковые номера (trace 1, 2, 3 и т.д.) графиков (сверху вниз по оси Оу). Т.е. сначала идёт трассировку графика f(x), trace 2 – метка графика а(х) и т. д.. (Кстати в окне, где построен график, по оси Оу слева можно увидеть соответствующую ему трассировку). Во второй колонке указан символ, который соответствует данной трассировки. Можете поэкспериментировать с его изменением, но в данном случае оптимальным будет выбор параметра NONE(нет). В колонке LINE (линия) представлена текущее отображение типа линии графика. Вы можете, нажимая мышью по выделенным овалом на рисунке стрелочкам выбирать из всплывающего меню различные пункты (можно поэкспериментировать и с этим параметром, выбирая соответствующий из всплывающего меню, и нажимая кнопку «применить» (чтобы не выходить каждый раз из окна). В столбце COLORS (цвета) можно тем же способом выбрать соответствующий цвет для каждой кривой. В столбце TYPE(тип) можно выбрать тип представления графика (точками, диаграммой и т. д.). В нашем примере, конечно оптимальным будет выбор LINES (линия). Итак, для каждой функции в окне, где построены их графики «сконструируем» из предложенных MathCADом параметров соответствующий вид.

h) В итоге, после всех преобразований должно получиться что-то похожее на представленный ниже рисунок:

На первый взгляд кажется, что построение графиков в MathCAD довольно сложная задача. Но выше я описал весь процесс очень подробно, и, стоит только «приноровиться» к возможностям программы всё пойдёт очень быстро, фактически нажатием трёх-четырёх кнопок.

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

Пусть дана функция (пример 7.15). Требуется повести полное её исследование.

Вспомнив наш план, начнём с области определения. Можно сразу догадаться. Что это любое действительное число (знаменатель x²+4 не равен нулю ни при каком действительном х), а можно проверить это с помощью вычислительного блока и операторов GIVEN-FIND).

Определим функцию: .

Решим уравнение . Вводим служебное слово начала вычислительного блока (Given), выражение с булевым равно, служебное слово FIND(х), нажимаем CTRL+ю (подробно этот процесс описан в пункте, где мы искали наибольшее и наименьшее значение функции на отрезке) и получаем ответ: . Т.е. уравнение не имеет действительных корней.

МАТЕМАТИЧЕСКАЯ СПРАВКА i – мнимая единица, где i²=-1 (не равносильно из-за неоднозначности определения операции извлечения корня у мнимых чисел).

Следующий шаг – проверка функции на чётность. Набираем с клавиатуры у(-х), нажимаем CTRL+ю (à), и программа в символьной форме выводит результат: . Мы видим, что , откуда следует, что данная функция является нечётной (её график должен быть симметричен относительно начала координат).

Найдём асимптоты кривой (подробно описано на задаче 6.15). Ясно, что вертикальных асимптот данная кривая не имеет ( функция не «вылетает» в бесконечность ни в одной точке). Ищем наклонные асимптоты. Коэффициент k= , коэффициент b= . Откуда вроде бы следует, что асимптотой данной кривой является ось Ох. Да, действительно, как мы потом увидим из графика, кривая асимптотически приближается к оси Ох, но найдя точки пересечения с осями координат (решая уравнение у(х)=0 через GIVEN-FIND операторы: ), получим, что график проходит через начало координат (х=0, у=0). Так что с точки зрения математики, ось Ох вряд ли можно назвать асимптотой кривой.

Найдем производную в символьном виде с помощью оператора дифференцирования (гл.II) (здесь и далее, чтобы не повторяться, в скобках указываю, где данная тема была подробно разобрана): .

Т.к. производная непрерывна во всех точках, то критические точки следует приравнивая её нулю (если бы производная где-либо была разрывной точку её разрыва нужно было бы включить в критические). Повторяем уже знакомые нам действия и получаем ответ: . Обращаю внимание, что в вычислительный блок в данном случае я не стал копировать и затем вставлять производную (как мы делали в примере 3.15), а просто написал оператор дифференцирования. MathCAD автоматически вычислил производную, а затем приравнял её к нулю.

ВНИМАНИЕ Не всегда это получается, иногда выводятся сообщения об ошибки.

Итак, теперь мы имеем две критические точки, которые могут быть максимами или минимами, а могут, вообще говоря, не являться ни тем ни другим. Поэтому, проверим знак производной при переходе через эти точки. Я сделал это следующим образом:

1. MathCAD отказался почему-то отказывается решать неравенства с помощью блока GIVEN-FIND, поэтому приходиться воспользоваться панелью Символических операторов (Symbolic). Если она не видна, следует на панели Математика (Math) нажать на кнопку с изображением шапочки: .

2. Далее, в появившейся панели символических операторов следует нажать на кнопку с надписью Solve (решить): . В появившемся на рабочем листе служебном слове: начинаем заполнять метки. В левую метку вводим наше неравенство, которое требуется решить, т.е. можно ввести это таким образом: ( т.е., решая такое неравенство, мы по существу ищем промежутки возрастания и убывания функции).

ВНИМАНИЕ Между оператором символьного дифференцирования и нулём необходимо поставить БУЛЕВОЕ «больше или равно» ( на панели BOOLEAN необходимо нажать соответствующую клавишу: ). В правую метку заносим переменную, которую необходимо «отыскать», т.е., в нашем случае х.

3. Нажимаем ENTER и получаем ответ: . Осмыслим его. По существу, это означает, что наша функция возрастает на промежутке [- ; -2] и [2; ]. Значит, она непременно должна убывать на интервале [2; 2]. Проверим это (можно, конечно, не делать этого, но для тренировки лучше повторить решение неравенств, в тем более, что я покажу далее как можно решать их без помощи блока SOLVE). Итак, решим неравенство . (решить неравенство в виде данным методом не получается).

4. Скопируем в буфер обмена выражение, стоящее справа от оператора дифференцирования (à): (подробно о копировании выражений см. гл. II пример 3.15 и в примере из этой главы № 3.15). Вставляем выражение в свободном месте на рабочем листе, для чего из пункта редактировать в главном меню выбираем пункт «вставка». Должно появиться только выражение, выделенное прямоугольником. Если опять появился оператор дифференцирования – сотрите его, нажимая клавишу Backspace на клавиатуре. Далее, мышью поставим метку-заполнитель так, чтобы она «цепляла» только ту переменную, которую необходимо решить, т. е. х в любом месте в выражении (подчёркнуто) . БУДЬТЕ ВНИМАТЕЛЬНЫ! НЕ ЗАЦЕПИТЕ МЕТКОЙ КВАДРАТ!

5. Из пункта «символика» главного меню выбираем пункт «переменная», а далее из выплывающего подменю пункт «решение» («разрешить»).

6. Получаем ответ: . Как и должно было получиться, функция убывает на промежутке [-2; 2]

СОВЕТ Если вам непонятно как перевести ответ из языка MathCAD в обычный язык поясняю – начертите координатную линию и расставьте точки –2, 2 и т.п.. А дальше смотрим так: первая часть решения х≥-2 – отмечаем всё то, что больше –2, вторая часть решения х≤2 – отмечаем всё то, что меньше 2, затем смотрим, где пересекаются наши отметки )

-2 2

Теперь мы уже имеем в наличии участки монотонности, максимумы и минимумы функции ( -2 – точка максимума т. к. производная меняет знак с + на -, точка 2 – точка минимума – производная меняет знак с – на +). Но всё же, этого не достаточно, чтобы построить приличный график. Найдём вогнутости и выпуклости кривой, для чего посчитаем вторую производную (см. математическую справку).

Используем уже знакомый нам по главе Дифференцирование оператор производной N-ого порядка. Щелкаем его из панели CALCULUS и в появившемся на листе выражении заполним метки соответствующим образом (подробнее см. главу II) à . Напомню, что вместо выражения я вставил символ у(х), т. к. функция была определена нами заранее. Нажимаем CTRL+ю (лат.), ENTER и получаем результат: .

Приравняем нулю вторую производную (необходимое условие точки перегиба) и составим вычислительный блок. Получим ответ: . Посмотрим на выполнение достаточного условия точки перегиба, т. е. решим неравенство: любым из предложенных выше способов. Должно получиться . Переведу это на обычный математический язык: у’’>0 при х€ (-∞;- )U[0; ) Т. е. у наша кривая будет вогнутой на указанных интервалах. Найдём теперь участки, где кривая будет выпуклой, т. е. решим неравенство . Сразу привожу ответ. который должен получиться: . Переводя это на понятный язык, получим. Что кривая будет выпуклой на интервале (- ; 0}U( ; ∞).

Ну что. Всё готово. Чтобы построить график на бумаге и, соответственно на MathCADе тоже. Подробнейшим образом процесс построения графика был описан на примере из данной главы № 6.15, поэтому далее я привожу готовый график функции (MathCAD практически всё сделал автоматически – выбрал удобно масштабы и т.п., я лишь изменил представление осей на CROSSED – пересечение из меню форматирования графиков):

Вспомним то, что я говорил насчёт асимптоты (оси Ох). Да, действительно, график приближается неограниченно к этой оси на бесконечности, но он её пересекает в начале координат, поэтому говорить об оси Ох как об асимптоте, неверно.

Повторим исследование для разрывной функции . ( пример 8.15)

Найдем её область определения, решая неравенство (логарифмическая функция определена только тогда, когда аргумент её строго больше нуля) - . Т.е. ( - ; -1)U(0; )

Теперь исследуем функцию на чётность – нечётность, для чего вводим с клавиатуры f(-x), нажимаем CTRL+ю, ENTER, и получаем результат: - функция не является ни чётной ни нечетной.

Найдём теперь асимптоты кривой. Решим такое уравнение: (мы ищем такие значения х, которые «выкидывают» функцию в бесконечность, а, как известно, ). Мы, проводя вычисления в SOLVE-блоке, должны получить следующий результат:

.

Значит, точка х=0 является вертикальной асимптотой данной функции. Но это ещё не всё. Внимательно посмотрев на уравнение. Задающее функцию, найдем, что она уходит в бесконечность при х=-1. (По-существу мы должны решить такое уравнение: , где х₀ его корни, т.е. найти те значения х_0, при которых функция обращается в бесконечность – но в таком виде MathCAD решить его не может, поэтому приходиться додумываться самому, глядя на уравнение искать «необычные» точки). Имеем две вертикальные асимптоты: . Наклонные асимптоты находим, используя оператор предела (подробно см. пример 6.15): , .

Точки пересечения с осями координат вычислим следующим образом:

. Обратите внимание на строку 3. В ней я использовал оператор присвоения (определения). Таким образом. Присвоив в строке 3 найденному корню букву у, в строке 4 программа автоматически посчитала значение функции в этой точке. И мы сразу получили точку пересечения с осью Ох: ( , 0 ). Для того, чтобы найти точку пересечения с ось Оу вычислим . Правильно, MathCAD выдаёт ошибку т.к. функция в нуле неопределенна.

Вычисляя производную и решая следующее уравнение (критические точки):

получим сообщение об ошибки. Да, действительно, производная не равна нулю нигде, но это не означает, что не существует критических точек – они есть: x = 0 x=-1 – в них выражение, определяющее производную не существует.

Посмотрим, где производная меньше, а где больше нуля, решая неравенства:

1. (ответ )

2. (ответ )

С учётом области определения имеем: возрастание (- ;-1]U(0; ). Промежуток убывания [-1; 0] не входит в область определения.

Вычисляя вторую производную, а затем. Приравнивая её к нулю, получим:

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

Из неравенства

находим вогнутости и выпуклости (ответ на неравенство т. е. (- ; -1)U(-1; -1.2)). Но промежуток (- ; -1)U(-1; -1.2) не входит в область определения и, окончательно имеем вогнутость кривой на промежутке (- ; -1).

Из неравенства:

находим, что кривая выпукла на промежутке (0; ) (учётом области определения).

Строим график. Выбирая соответствующие трассировки и пункт «пересечение осей» (подробнейшим образом описано на примерах выше) (crossed) получаем:

.ПРИМЕЧАНИЕ Коэффициенты k и b наклонной (горизонтальной) асимптоты были определены нами заранее.

На этом считаю необходимым завершить данную главу, т.к. думаю, что более чем подробно описал исследование функций и построение их графиков. На мой взгляд MathCAD обладает достаточно мощными возможностями, чтобы провести полное исследование функции, но всё таки кое-что необходимо додумывать самому, подключая знания из математики (критические точки, область определения и т. п.).

Настало время перейти к не менее важному разделу математики, к действию «обратному дифференцированию», как говорят математики (хотя интегрирование первоначально не имело к дифференцированию никакого отношения, и лишь благодаря Лейбницу и Ньютону, мы теперь смотрим на интегрирование и дифференцирование как на неразрывные вещи, а то бы и сейчас писали все интегралы через интегральные суммы).

ГЛАВА IV.

ПРОСТЕЙШЕЕ ИНТЕГРИРОВАНИЕ ФУНКЦИЙ ОДНОЙ ПЕРЕМЕННОЙ.

Рассмотрим теперь простейшие приёмы интегрирования в среде MathCAD, а вместе с этим и приложения определенных интегралов, которые лично для меня, как физика особенно интересны.

Вот простейший пример (2.15) интегрирования в среде MathCAD. Для его расчёта используем (т.к. интеграл определённый) оператор «определённый интеграл» из уже знакомой панели CALCULUS:

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

Если вы хотите получить символьный результат (т.е. не в виде числа в десятичной форме, а в виде ln, экспонент, sin, и т.д.) нужно «зацепить» меткой всё выражение, нажать CTRL+ю и ENTER. В результате получиться следующее выражение: . Чтобы получить ответ в виде числа достаточно нажать равно на клавиатуре после символьного выражения. Т.е. .

Следующий пример (4.15) выполняется аналогично: Внимание следует обратить лишь на вставку функции asin(x) ( подробно на примерах 11.15 12.15 в главе I). Здесь мы также получили символьный результат - 4/3, переведём теперь его в десятичную дробь: . Можно было получить этот результат сразу, вводя = сразу после интеграла (а не ставить -->)

В примере 5.15 следует вместо оператора определённого интеграла использовать оператор неопределенного интеграла:

Остальные же действия выполняются аналогично, поэтому далее я привожу сразу ответ: .

В примере 7.15 обратить внимание следует на правильность расставления скобок, а в остальном он не должен вызвать затруднений: .

Ответ на следующий пример 8.15 .

Если вы во время вычисления следующего примера 9.15 используете символьный вывод результатов (à), то ничего хорошего из этого не получиться. MathCAD выведет в качестве ответа что-то большое и страшное. Я здесь даже не буду приводить ответ на этот пример в символьном виде. Хотя, если вы задумаете просчитать этот пример вручную (с помощью универсальной тригонометрической подстановки), то вы получите вполне нормальный ответ. Всё дело как раз и состоит в том, что MathCAD не знает подстановок – программа пытается всегда всё разложить в ряд. Это и приводит к такому результату. В чиленном же виде всё получается довольно просто: .

Примеры 11.15 и 12.15 вычисляются стандартным способом: (11.15)

(12.15)

ПРИЛОЖЕНИЯ ОПРЕДЕЛЁННОГО ИНТЕГРАЛА.

В этом подпункте, специально мною выделенном, мы поговорим о некоторых прикладных задачах, решающихся с помощью определённого интеграла (вычисление площадей, длин дуг кривых и т.д.)

Определённыq интеграл (геометрическое истолкование) – это площадь, ограниченная графиком функции, соответствующей осью (не обязательно Ох – может быть уравнение в виде х=у²) и пределами интегрирования.

Пусть требуется вычислить площадь фигуры, ограниченной линией х= , у=0 и графиком функции у= . Для наглядности построим график этой функции и линий, которые её ограничивают (см. Главу III):

Искомая площадь фигуры есть площадь треугольника, обведённого овалом. По существу здесь мы должны проинтегрировать функцию у= в пределах от 0 до ( где - абсцисса точки пересечения графика функции и прямой): .

В примере 15.15 кривая задана в параметрической форме: ; . – эта кривая называется циклоида. Её ограничивает линия у=6 и по условию у≥6, а 0<x<12π (х, но не t!!!). Для того чтобы построить график параметрически заданной функции необходимо заранее определить функцию (задать как зависит х и у от параметра t), по горизонтальной оси поставить вместо привычной метки х метку х(t), по вертикальной вместо у(х) y(t). Затем на вертикальной оси ставим запятую (подробно см. гл.III о просмотре графиков нескольких функций на одном поле) и вводим число 6 (т.к. по условию у=6 – ограничивающая линия). В результате этих действий должно получиться следующее:

Ищем площадь, выделенную на рисунке эллипсом (из площади циклоиды вычитаем площадь прямоугольника). Но для начала определим пределы интегрирования, для чего решим следующее уравнение: Учитывая период функции cos (т.к. MathCAD выдаёт только главный корень), найдём, что t= ; .

ПРИМЕЧАНИЕ Вообще-то период cos , но как мог бы заметить внимательный читатель, следующий корень отличается от предыдущего не на , а на , это связано с тем, что решая уравнение вида cos a=0 следует помнить, что следующий корень будет отличаться от предыдущего на , а не на 2π – частный случай тригонометрического уравнения, вообще-то, решая такие частные случаи, для наглядности рисуют единичную окружность, где х-ось cos, у-ось sin, и смотрят, когда повториться тот или иной корень)

Найдём площадь по формуле: S= .: . (я полагаю, что мы определили заранее все функции). Это мы нашли площадь, ограниченную осью Ох, циклоидой и вертикальными линиями х=х(t₁) и х=х(t₂). Вычитаем из неё площадь прямоугольника и получаем искомую площадь: .

ВНИМАНИЕ В выражении, обведенном овалом я нашёл площадь этого самого прямоугольника, ограниченного линиями у=6, х=х(t₁) и х=х(t₂). Всего лишь длину (у=6) умножил на ширину: . Нужно помнить, что вычитать надо не параметры t₂ - t₁ , а значение функции х от этих параметров!!!

Поработаем теперь с полярными координатами (16.15). Здесь нужно найти площадь, ограниченную следующими функциями, заданными в полярных координатах: .

Для начала построим их графики:

Искомая площадь заштрихована на рисунке.

МАТЕМАТИЧЕСКАЯ СПРАВКА В полярных координатах S= .

Как видно из рисунка (это можно легко доказать и с помощью математики: поу словию мы считаем радиус r всегда положительным, т.е. и . Значит, решая тригонометрическое уравнение и, находя так называемое главное решение, имеем: .) пределы интегрирования есть и . Искомая площадь получается вычитанием из большой окружности маленькой, т.е. S= .

В задании 17.15 нужно вычислить длину дуги, заданной в декартовом ортонормированном базисе.

МАТЕМАТИЧЕСКАЯ СПРАВКА L_{в декарт. коорд.}= .

Действуем стандартным методом:

1. Определяем функцию , где по условию. Здесь даже нет необходимости строить график – пример только лишь на применение формулы.

2. Вводим соответствующую формулу:

3. Нажимаем CTRL+à, ввод и получаем ответ: . Для удобства я посчитал его в численной форме.

ПРИМЕЧАНИЕ У меня почему-то в этом примере, когда я после определения функции, стал считать отдельно символьно y’, программа в качестве ответа писала следующее: . Наверное, это связано с тем. что она брала какие-то данные из других номеров (я делал их все в одном файле), в то время, как создав новый документ. Расчёт был произведён правильно.

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

МАТЕМАТИЧЕСКАЯ СПРАВКА L_{пар. ф.}= .

Определяем зависимость у и х от t, т.е функции , . По условию - это и будут наши пределы интегрирования. Далее вводим соответствующую формулу, нажимаем à и получаем символьный результат: . Опять же, считаю, что нет необходимости строить график – пример только на применение формулы.

Разберём последний пример в этой главе – 19.15. Здесь мы имеем дело с функцией, заданной в полярных координатах: .

МАТЕМАТИЧЕСКАЯ СПРАВКА L_{пол. к.}= .

По условию имеем, что . Для наглядности привожу график этой функции:

Синим и зелёным пунктиром построены пунктиром построены графики прямых соответственно φ= и φ=0, являющихся, по-существу, пределами интегрирования. Сама кривая называется кардиоида.

Осталось ввести формулу и получить результат: .

ЗАКЛЮЧИТЕЛЬНЫЕ ПРИМЕЧАНИЯ

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

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

Можно вычислять интегралы с одним или обоими бесконечными пределами. Для этого на месте соответствующего предела введите символ бесконечности, воспользовавшись, например, той же самой панелью Calculus (Вычисления). Чтобы ввести -°° (минус бесконечность), добавьте знак минус к символу бесконечности, как к обычному числу.

Чтобы получить результат интегрирования, следует ввести знак равенства или символьного равенства. В первом случае интегрирование будет проведено численным методом, во втором - в случае успеха, будет найдено точное значение интеграла с помощью символьного процессора MathCAD. Конечно, символьное интегрирование возможно только для небольшого круга несложных подынтегральных функций. Но далее мы углубляться не будем. Ибо это уже не простейшие вычисления.

ГЛАВА V

НЕМНОГО ЛИНЕЙНОЙ АЛГЕБЫ.

И в заключении хочу рассмотреть один пример из линейной алгебры (9.15). Одной из основных задач линейной алгебры является нахождение собственных значений и собственных векторов матрицы. Нам дана квадратная матрица: . Требуется найти её собственные вектора и соответствующие им собственные значения.

Очевидно, первое, что нам нужно сделать – это ввести матрицу. Набираем на клавиатуре любую, желательно, латинскую букву (этой букве мы присвоим (определим) нашу матрицу – также бы мы сделали на бумаге для уменьшения громоздкости выкладок), оператор присваивания ,нажимаем на кнопку с изображением матрицы (если данная панель не видна, щёлкаем левой кнопкой мыши на похожем значке в панели математика). В появившемся диалоговом окне «вставить матрицу» в русской версии: или в соответствующем окне в английской версии: задаем количество строк и количество столбцов в матрице и нажимаем OK или INSERT (вставить) и матрица появляется на рабочем листе в таком виде: . Щёлкаем мышью по первой метки и вводим первый элемент матрицы с клавиатуры. Перемещаясь клавишей TAB по меткам вводим последовательно остальные её элементы: . Далее используем спец.функции MathCAD.

МАТЕМАТИЧЕСКАЯ СПРАВКА: По существу мы должны решить следующее уравнение в матричной форме: . Чтобы столбец векторов имел не тривиальное (нулевое) решение, необходимо и достаточно, чтобы определитель матрицы A был равен 0.

Нам потребуются всего лишь две функции: - eigenvals (A) – находит вектор собственных чисел матрицы А, и eigenvec(A, λ) – находит собственные вектора. Соответствующие данному числу λ.

Из меню вставка выбираем пункт «функция» (или щелкаем на значок функции в панели STANDART (стандартная)). В появившемся диалоговом окне из меню «категория функций» выбираем пункт «VECTOR and MATRIX», а из меню «названия функции сначала функцию с названием eigenvals и нажимаем ENTER:

На рабочем листе появится следующее: . В метку вводим ту букву, какую присвоили в начале матрице. Нажимаем = с клавиатуры и получаем ответ: .

Далее вставляем функцию с названием eigenvec по той же схеме. В первую метку вводим букву, соответствующую матрице, во вторую первое собственное значение, т.е. 7. Нажимаем = и получаем ответ: Эти действия необходимо повторить для двух других собственных значений: , .

ЗАКЛЮЧЕНИЕ

Описанные мною действия, как многим может показаться, демонстрируют использование MathCAD в качестве обычного калькулятора с расширенным набором функций. Но это не так!

MathCAD может использоваться для оцифровки звука из звуковых файлов, для вставки и редактирования рисунков и т.д. К сожалению, остались не рассмотрены эти и многие другие возможности MathCAD - линейная и не линейная регрессия, возможности по работе с трехмерной графикой, функции нескольких переменных, работа с циклами… Это связано с ограниченностью по времени и размеру данной работы, ибо вообще-то это – курсовой проект. Эти и многие другие темы подробно рассмотрены в Центре ресурсов (к сожалению, он на английском языке).

Конечно, те возможности, которые я попытался описать выше (удачно или нет – судить вам) составляют вряд ли 10 % от всего того, что можно делать с помощью программы. Ведь, по существу, считая даже простенький интеграл на MATHCAD, мы составляем именно свою программу (жалко, что она не может быть откомпилирована, т.е. стать независимым приложением).

Может быть, опытным пользователям покажется, что данное пособие написано уж через чур подробно, но я рассчитывал на неподготовленных читателей, и вообще это пособие может послужить началом к более крупной и серьёзной работе, в которой будут описываться все, в том числе и расширенные возможности MATHCAD. А пока на это нет времени…

Теперь о том, что я думаю о 2001 версии MathCAD. Быть может это моё субъективное мнение, но MathCAD 2000 мне показался лучше. Я работал с полурусифицированной английской версией MathCAD 2000, затем с 2001-й версией, но она хоть и русифицирована больше и качественней, каких-то сверхновых функций не имеет (добавлены спецфункции по работе со звуком, графикой, WEB), а установка MATHCAD 2001 стала намного сложнее. Насчёт возможностей: какие-то номера я делал в 2001, какие-то в 2000-й и разницы не почувствовал. Плохо то, что документ, сделанный в 2001 версии никогда не открывается в 2000-й, а документ, сделанный в 2000-й не всегда откроется в 2001. Кроме всего этого, MathCAD-2001 думает почему-то намного дольше 2000-го (возможно. Это связано с реализацией новых алгоритмов)

Можно сделать вывод, который уже как бы неявно фигурировал с самого вступления - для профессионального математика, возможно, что функций и возможностей MathCADа покажется недостаточным, а для студента, не проводящего ни каких суперрасчётов - чтобы проверить тот же самый Типовой Расчёт, вполне хватит. Более серьезным вычислителям можно порекомендовать использовать такие пакеты как MATLAB, MAPLE и др.

На этом считаю необходимым завершить этот кратенький очерк о MathCAD. Спасибо, до свидания…