Matlab  |  Mathematica  |  Mathcad  |  Maple  |  Statistica  |  Другие пакеты Поиск по сайту
Internet-класс  |  Примеры  |  Методики  |  Банк задач  |  Форум  |  Download  |  Ссылки  |  Конкурсы


10 ноутбуков
В.П.Дьяконова
 
Интегралы и Mathematica 4.0

Введение

Вычисление интегралов - привлекательная возможность системы Mathematica 4. Однако при этом важно понимать, как выполняется это тонкое действие. Начнем с вычисления определенных интегралов вида

[Graphics:Images/index_gr_1.gif]

Здесь f(x) - подынтегральная функция, a- нижний предел интегрирования, b- верхний предел интегрирования и x - переменная, по которой производится интегрирование. Зададим простую функцию f(x)=x^2.

[Graphics:Images/index_gr_2.gif]

Исторически под интегралом понималась площадь m прямоугольников, высоты левых граней у которых равны значениям функции. Это соответствует затемненной площади на следующем рисунке (m=4, a=1, b=3).

[Graphics:Images/index_gr_3.gif]
[Graphics:Images/index_gr_4.gif]

[Graphics:Images/index_gr_5.gif]

[Graphics:Images/index_gr_6.gif]

Такой способ приближенного вычисления интеграла сейчас известен как способ интегрирования методом прямоугольников. При числе прямоугольников m=10 (можно выбрать любое m) m он реализуется следующим образом:

[Graphics:Images/index_gr_7.gif]
[Graphics:Images/index_gr_8.gif]
[Graphics:Images/index_gr_9.gif]
[Graphics:Images/index_gr_10.gif]
[Graphics:Images/index_gr_11.gif]

Заметим, что точное значение интеграла равно 8.6667, таким образом простейший метод прямоугольников дает большую погрешность (примерно 1/m). Можно усовершенствовать метод, беря для высот прямоугольников отсчеты функции в середине каждого прямоугольника. Графически это иллюстрирует следующий рисунок:

[Graphics:Images/index_gr_12.gif]

[Graphics:Images/index_gr_13.gif]

[Graphics:Images/index_gr_14.gif]

Такой метод интегрирования называется модифицированным методом прямоугольников. Ниже дан пример его реализации:

[Graphics:Images/index_gr_15.gif]
[Graphics:Images/index_gr_16.gif]
[Graphics:Images/index_gr_17.gif]
[Graphics:Images/index_gr_18.gif]
[Graphics:Images/index_gr_19.gif]

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

[Graphics:Images/index_gr_20.gif]
[Graphics:Images/index_gr_21.gif]
[Graphics:Images/index_gr_22.gif]
[Graphics:Images/index_gr_23.gif]
[Graphics:Images/index_gr_24.gif]

Мы можем продолжать улучшение методов интегрирования. К примеру можно вычислять интегралы методом трапеций, используя кусочно-линейную аппроксимацию f(x), методом Симсона (параболическая аппроксимация f(x) на каждых смежных трех точках) и т. д. Описание этих методов можно найти в литературе по численным методам. Чаще всего на практике используется метод Симпсона. Он не очень сложен и дает вполне приемлемую точность. Однако, при описании возможностей интегрирования в системе Mathematica все эти методы имеют исключительно иллюстративное значение. Ибо встроенные функции интегрирования системы Mathematica реализуют лучшие из адаптивных методов интегрирования, для которых указанные выше методы являются частными случаями.

Геометрическая трактовка интеграла (общий случай)

Геометрически значение определенного интеграла представляет собой площадь, ограниченую кривой f(x), осью x и вертикалями c x=a и y=b. Это иллюстрирует приведенный ниже рисунок, в котором функция f(x) задает параболу,смещенную вниз, а пределы a=1.5 и b=3.5.

[Graphics:Images/index_gr_25.gif]
[Graphics:Images/index_gr_26.gif]
[Graphics:Images/index_gr_27.gif]

[Graphics:Images/index_gr_28.gif]

[Graphics:Images/index_gr_29.gif]

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

Вычисление определенных интегралов в Mathematica 4

Теперь вычислим определенный интеграл для функции f(x) в общем случае и при заданных пределах интегрирования a и b:

[Graphics:Images/index_gr_30.gif]
[Graphics:Images/index_gr_31.gif]
[Graphics:Images/index_gr_32.gif]

Мы получили аналитическое выражение для определенного интеграла. Задав конкретные a и b можно вычислить интеграл как число:

[Graphics:Images/index_gr_33.gif]
[Graphics:Images/index_gr_34.gif]
[Graphics:Images/index_gr_35.gif]
[Graphics:Images/index_gr_36.gif]
[Graphics:Images/index_gr_37.gif]

Разумеется после задания конкретных пределов интеграл также вычисляется как число

[Graphics:Images/index_gr_38.gif]
[Graphics:Images/index_gr_39.gif]

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

[Graphics:Images/index_gr_40.gif]
[Graphics:Images/index_gr_41.gif]

Обратите внимание и на такой пример:

[Graphics:Images/index_gr_42.gif]
[Graphics:Images/index_gr_43.gif]

Функции Integrate и NIntegrate

Для вычисления данного интеграла (да и других) можно использовать функции:

[Graphics:Images/index_gr_44.gif]
[Graphics:Images/index_gr_45.gif]
[Graphics:Images/index_gr_46.gif]
[Graphics:Images/index_gr_47.gif]
[Graphics:Images/index_gr_48.gif]
[Graphics:Images/index_gr_49.gif]

Функция Integrate дает аналитическое значение интеграла, если оно существует, тогда как функция NIntegrate, эквивалентная N[Integrate[f[x],{x,a,b}]], дает численное значение интеграла.

Интегралы с переменным пределом

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

[Graphics:Images/index_gr_50.gif]
[Graphics:Images/index_gr_51.gif]

Здесь интересно и задание функции пользователя fun[y], которая выражается через интеграл. Может быть построен график этой функции.

[Graphics:Images/index_gr_52.gif]

[Graphics:Images/index_gr_53.gif]

[Graphics:Images/index_gr_54.gif]

Приведем еще один пример, в котором оба предела интегрирования являются функциями переменной y:

[Graphics:Images/index_gr_55.gif]
[Graphics:Images/index_gr_56.gif]
[Graphics:Images/index_gr_57.gif]

[Graphics:Images/index_gr_58.gif]

[Graphics:Images/index_gr_59.gif]

А следующий пример показывает случай, когда верхний предел, как и функция f[x], является функцией переменной x:

[Graphics:Images/index_gr_60.gif]
[Graphics:Images/index_gr_61.gif]

Определенный интеграл от функции синуса

Итак, площадь при вычислении определенных интегралов со знакопеременной f(x) образуется алгебраическим суммированием сегментов над и под осью абсцисс x. С этим связаны некоторые тонкие особенности интегрирования. Рассмотрим к примеру f(x) в виде двух периодов синусоиды. Очевидно, что общая площадь в алгебраическом смысле в этом случае равна 0.

[Graphics:Images/index_gr_62.gif]

[Graphics:Images/index_gr_63.gif]

[Graphics:Images/index_gr_64.gif]

Вычисление определенного интеграла с этой функцией и с пределами от -2p до 2p и впрямь дает нуль:

[Graphics:Images/index_gr_65.gif]
[Graphics:Images/index_gr_66.gif]

Mathematica неплохо справляется с вычислением интеграла от функции синуса даже при пределах от -1000000p до 1000000p, поскольку очевидно, что при любом конечном числе периодов синусоиды их общая площадь равна 0.

[Graphics:Images/index_gr_68.gif]
[Graphics:Images/index_gr_69.gif]

Проблема 1. Многие интегралы могут не иметь значений при пределах, стремящихся к бесконечности. С первого взгляда может показаться, что и при бесконечном "числе" пределов мы должны получать знание интеграла, равное 0. Однако на деле такой интеграл просто не существует, что и видно из следующих примеров:

[Graphics:Images/index_gr_70.gif]
[Graphics:Images/index_gr_71.gif]
[Graphics:Images/index_gr_72.gif]

Замена одного предела конечным (нулем) не меняет дела:

[Graphics:Images/index_gr_73.gif]
[Graphics:Images/index_gr_74.gif]
[Graphics:Images/index_gr_75.gif]
[Graphics:Images/index_gr_76.gif]
[Graphics:Images/index_gr_77.gif]
[Graphics:Images/index_gr_78.gif]

Интеграл от функции x[Graphics:Images/index_gr_79.gif]

Теперь рассмотрим функцию, площадь которой равна 1:

[Graphics:Images/index_gr_80.gif]
[Graphics:Images/index_gr_81.gif]

[Graphics:Images/index_gr_82.gif]

[Graphics:Images/index_gr_83.gif]

Специфика этой функции заключается в том, что ее значение быстро уменьшается с увеличением x. Например, при интервале изменения x от 0 до 10000 кривой функции уже просто не видно на графике (точнее она совпадает с осью x):

[Graphics:Images/index_gr_84.gif]

[Graphics:Images/index_gr_85.gif]

[Graphics:Images/index_gr_86.gif]

Это обстоятельство довольно коварно. Многие методы интегрирования основаны на обработке конечного числа отсчетов функции f(x). Однако в подобной ситуации все отсчеты могут быть "нулевыми" с позиций ограниченной точности вычисления очень малых чисел. Например, в системе Mathcad 7 вычисление интеграла [Graphics:Images/index_gr_87.gif]  давало грубейшую ошибку в виде значения 0.

Чтобы лучше понять возникновение данной ошибки рассмотрим вычисление определенного интеграла описанным выше методом прямоугольников:

[Graphics:Images/index_gr_88.gif]
[Graphics:Images/index_gr_89.gif]
[Graphics:Images/index_gr_90.gif]

Результат не так уж плох. Но теперь, сохранив число прямоугольников m=100, изменим предел интегрирования b на 100000:

[Graphics:Images/index_gr_91.gif]
[Graphics:Images/index_gr_92.gif]
[Graphics:Images/index_gr_93.gif]

Вместо ожидаемого приближения к 1 мы получаем "неожиданный" результат - близкое к нулю значение интеграла. Причина этого уже указывалась - отсчеты f2(a+i*dx) попадают на участки функции f2(t), где ее значения очень малы. А первый отссчет попадает на нуль функции. Сумма почти (мли равных) нулю отсчетов дает близкое к 0 значение.

Mathematica 4 имеет усложненный алгоритм вычисления интегралов, который, по мере возможности, анализирует поведение f(x) и старается исключить описанную ситуацию. Нетрудно убедиться,чтоMathematica 4 знает особенность данной функции (f2[x]) и вычисляет интеграл верно:

[Graphics:Images/index_gr_94.gif]
[Graphics:Images/index_gr_95.gif]

Вычисление интегралов с f(x) вида y(x)*sin(x)

А что будет, если интегрируемая функция синусоида, амплитуда которой медленно уменьшается или нарастает? Может показаться, что ситуация при этом остается прежней - как при интегрировании синусоиды. Но это далеко не так! Изменение амплитуды синусоиды превращает ее в новую функцию с принципиально новыми свойствами. Рассмотрим, например, функцию в виде синусоиды с экспоненциально уменьшающеся до 0 амплитудой:

[Graphics:Images/index_gr_96.gif]
[Graphics:Images/index_gr_97.gif]

[Graphics:Images/index_gr_98.gif]

[Graphics:Images/index_gr_99.gif]

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

[Graphics:Images/index_gr_101.gif]
[Graphics:Images/index_gr_102.gif]
[Graphics:Images/index_gr_103.gif]
[Graphics:Images/index_gr_104.gif]
[Graphics:Images/index_gr_105.gif]

Поскольку функция f2(x) быстро затухает, то можно ожидать, что даже при конечном b, например b=100, будет получено то-же значение интеграла, что и при точном расчете. И впрямь:

[Graphics:Images/index_gr_106.gif]
[Graphics:Images/index_gr_107.gif]

Теперь рассмотрим функцию в виде синусоиды, амплитуда которой нарастает по экспоненциальному закону, приближаясь к значению 1:

[Graphics:Images/index_gr_108.gif]
[Graphics:Images/index_gr_109.gif]

[Graphics:Images/index_gr_110.gif]

Нетрудно заметить, что здесь площадь сегментов под осью x всегда больше, чем сегментов над осью x. И в самом деле, вычислим интеграл в пределах от 0 до 6p - три полных периода синусоиды. Получим отрицательное число.

[Graphics:Images/index_gr_112.gif]
[Graphics:Images/index_gr_113.gif]

Нетрудно сообразить, что это число по модулю равно значению интеграла для функции f2. И далее:

[Graphics:Images/index_gr_114.gif]
[Graphics:Images/index_gr_115.gif]
[Graphics:Images/index_gr_116.gif]
[Graphics:Images/index_gr_117.gif]
[Graphics:Images/index_gr_118.gif]
[Graphics:Images/index_gr_119.gif]

А теперь вычислим интеграл с этой функцией в пределах от 0 до Ґ. Результат оказывается неожиданным - положительным числом:

[Graphics:Images/index_gr_120.gif]
[Graphics:Images/index_gr_121.gif]

Разум подсказывает, что в данном случае Mathematica 4 дает явную ошибку. Система Mathcad 2000, кстати, в данном случае честно сообщила, что не в состоянии вычислить данный интеграл.

Возьмем еще один пример:

[Graphics:Images/index_gr_122.gif]
[Graphics:Images/index_gr_123.gif]

[Graphics:Images/index_gr_124.gif]

[Graphics:Images/index_gr_125.gif]

С вычислением интеграла от этой функции Mathematica справляется легко:

[Graphics:Images/index_gr_126.gif]
[Graphics:Images/index_gr_127.gif]

Первообразная и неопределенный интеграл

Пусть функция f(x) есть производная от F(x), то есть f(x)dx есть дифференциал функции F(x). Следовательно

f(x)dx=dF(x). Функция F(x) называется первообразной для функции f(x). Нахождение наиболее общей первообразной для заданной функции f(x) называют вычислением неопределенного интеграла. Последний обозначается как [Graphics:Images/index_gr_128.gif]. Есть множество способов аналитического вычисления неопределенного интеграла. Mathematica 4 позволяет найти значения неопределенного интеграла для многих тысяч видов таких интегралов. Ниже дан пример этого для функции f(x)=x^n:

[Graphics:Images/index_gr_129.gif]
[Graphics:Images/index_gr_130.gif]
[Graphics:Images/index_gr_131.gif]
[Graphics:Images/index_gr_132.gif]

Пусть теперь конкретно n=2:

[Graphics:Images/index_gr_133.gif]
[Graphics:Images/index_gr_134.gif]

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

[Graphics:Images/index_gr_135.gif]

Или в нашем случае:

[Graphics:Images/index_gr_136.gif]
[Graphics:Images/index_gr_137.gif]

Наличие константы С говорит о существовании бесконечно большого числа первообразных, графики которых просто смещены по вертикали на величину C (от -Ґ до +Ґ). Ниже представлены графики первообразных функции x^2 для трех значений C = -1,0 и1:

[Graphics:Images/index_gr_138.gif]

[Graphics:Images/index_gr_139.gif]

[Graphics:Images/index_gr_140.gif]

Mathematica может вполне заметить самые объемные таблицы неопределенных интегралов. Ниже дано несколько примеров, поясняющих простоту вычисления интегралов:

[Graphics:Images/index_gr_141.gif]
[Graphics:Images/index_gr_142.gif]
[Graphics:Images/index_gr_143.gif]
[Graphics:Images/index_gr_144.gif]
[Graphics:Images/index_gr_145.gif]
[Graphics:Images/index_gr_146.gif]
[Graphics:Images/index_gr_147.gif]
[Graphics:Images/index_gr_148.gif]
[Graphics:Images/index_gr_149.gif]
[Graphics:Images/index_gr_150.gif]

Интегралы от функций с условиями сравнения

На практике многие функции приходится задавать, используя те или иные условия сравнения. Вначале рассмотрим пример при задании подынтегральной функции с помощью функции условных выражений If. Функция f5[x]задает треугольную площадку с площадью 1/8.

[Graphics:Images/index_gr_151.gif]
[Graphics:Images/index_gr_152.gif]
[Graphics:Images/index_gr_153.gif]

Посмотрим как справляется Mathematica с вычислением интегралов (определенного и неопределенного) от этой функции:

[Graphics:Images/index_gr_154.gif]
[Graphics:Images/index_gr_155.gif]
[Graphics:Images/index_gr_156.gif]
[Graphics:Images/index_gr_157.gif]

Теперь используя функцию знака Sign(дает 1 при положительном аргументе и -1 при отрицательном)зададим периодичекую функцию, интегрирование которой по периода должно дать площадь 2*p.

[Graphics:Images/index_gr_158.gif]
[Graphics:Images/index_gr_159.gif]
[Graphics:Images/index_gr_160.gif]

Вот как Mathematica вычисляет интегралы от этой функции:

[Graphics:Images/index_gr_161.gif]
[Graphics:Images/index_gr_162.gif]
[Graphics:Images/index_gr_163.gif]
[Graphics:Images/index_gr_164.gif]

Еще один пример показывает задание функции "двухполупериодного выпрямления" синусоиды. Такая функция встречается повсеместно в электротехнике, описывая выпрямление переменного тока.

[Graphics:Images/index_gr_165.gif]
[Graphics:Images/index_gr_166.gif]
[Graphics:Images/index_gr_167.gif]

Вычисление интегралов от нее показано ниже:

[Graphics:Images/index_gr_168.gif]
[Graphics:Images/index_gr_169.gif]
[Graphics:Images/index_gr_170.gif]
[Graphics:Images/index_gr_171.gif]

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

Вычисление кратных интегралов

Средства системы Mathematica 4 позволяют вычислять кратные, например двойные и тройные интегралы. Для этого операторы интегрирования используются один в другом. Например, если подынтегральная функция f(x,y) есть функция двух переменных, то двойной интеграл вычисляется как в следующем примере:

[Graphics:Images/index_gr_172.gif]
[Graphics:Images/index_gr_173.gif]

Аналогично для функций трех переменных можно вычислять тройные интегралы:

[Graphics:Images/index_gr_174.gif]
[Graphics:Images/index_gr_175.gif]
[Graphics:Images/index_gr_176.gif]

В кратных интегралах пределы могут быть функциями, например

[Graphics:Images/index_gr_177.gif]
[Graphics:Images/index_gr_178.gif]

Приведенными примерами возможности системы Mathematica 4 далеко не ограничиваются.

В начало страницы

Карта сайта | На первую страницу | Поиск | О проекте | Сотрудничество | e-mail
Корпоративная почта | ActiveCloud | Антивирус Касперского | Matlab | Подписка на MSDN для вузов | ИТ-ПРОРЫВ

Исправляем ошибки: Нашли опечатку? Выделите ее мышкой и нажмите Ctrl+Enter

Наши баннеры


Copyright © 1993-2017. Компания Softline. Все права защищены.

Дата последнего обновления информации на сайте: 04.03.17
Сайт начал работу 01.09.00

Softline – программное обеспечение, IT-консалтинг, лицензирование, обучение

подарки – подарочные сертификаты

 

            Rambler's Top100