Помощничек
Главная | Обратная связь


Археология
Архитектура
Астрономия
Аудит
Биология
Ботаника
Бухгалтерский учёт
Войное дело
Генетика
География
Геология
Дизайн
Искусство
История
Кино
Кулинария
Культура
Литература
Математика
Медицина
Металлургия
Мифология
Музыка
Психология
Религия
Спорт
Строительство
Техника
Транспорт
Туризм
Усадьба
Физика
Фотография
Химия
Экология
Электричество
Электроника
Энергетика

Вычисление определённого интеграла с помощью пакета Maxima



Определенный интеграл в символьном режиме (аналитически) вычисляется с помощью команды

integrate(f, x, a, b),

где f – подынтегральная функция, x – переменная интегрирования, a и b соответственно верхний и нижний пределы интегрирования.

Примеры. Вычислить интеграл:

1)

В вводной ячейке задаём команду

(%i4) integrate((1+cos(x))^2,x,0,%pi);

получаем результат

(%o4) (3*%pi)/2.

2)

(%i14) integrate(cos(2*x)*exp(-x),x,0,3);

(%o14) (%e^(-3)*(2*sin(6)-cos(6)))/5+1/5

(%i15) %o14, numer;

(%o15) 0.184874675854

Численное интегрирование выполняется функцией romberg или при помощи функций пакета quadpack.

Пример. Вычислить интеграл

В стоке меню выбираем кнопку Анализ → Integrate…(Рис. 7.2), в результате возникает вспомогательное окно для ввода подынтегральной функции и пределов интегрирования Интегрировать(Рис. 7.3), здесь же указывается режим интегрирования (численное), а также метод интегрирования romberg или quadpack.

Рис. 7.2

Рис. 7.3

Нажав на клавишу Okв рабочем окнепоявитсявводная ячейка и результат интегрирования

(%i2) quad_qags(tan(x)/(sin(x)^2-5*cos(x)^2+4), x, %pi/4, acos(1/sqrt(3)));

(%o2) [0.081093021621633,9.0031339740762459*10^-16,21,0]

В ячейке вывода (%о2) массив результата вычисления содержит:

0.081093021621633 – приближённое значение интеграла;

9.0031339740762459*10^-16 – относительная погрешность вычислений;

21 – число интервалов разбиения;

0 – признак корректности вычислений (0 – без проблем).

Примеры для самостоятельного решения

Вычислить определённый интеграл с помощью пакета Maxima

1) ( ) dx . 2) ( ) dx . 3) ( ) dx .

4) ( ) dx . 5) ( ) dx . 6) ( ) dx .

7) ( ) dx . 8) ( ) dx . 9) ( ) dx . 10) ( ) dx . 11) ( ) dx . 12) ( ) dx . 13) ( ) dx . 14) ( ) dx . 15) ( ) dx .

Применение определённого интеграла

7.4.1. Вычисление площадей плоских фигур

Прямоугольные координаты

Известно, что определенный интеграл на отрезке представляет собой площадь криволинейной трапеции, ограниченной графиком функции f(x). Если график расположен выше оси Ох (см. рис. 7.4), т.е. f(x) > 0,

Рис. 7.4

то площадь имеет знак “+“ и в этом случае искомая площадь определяется формулой

. (7.5)

Если график функции f(x) расположен ниже оси Ох (см. рис. 7.3), т.е. f(x) < 0,

Рис. 7.5

то площадь имеет знак “ - ” и .

Площадь фигуры, ограниченной кривыми y = f1 (x) и y = f2 (x) (при условии f2 (x)f1 (x) см. рис. 7.6), прямыми x = a и x = b может быть найдена с помощью определенных интегралов

. (7.6)

Рис. 7.6.

Примеры.

1. Вычислить площадь фигуры, ограниченной графиком функции y = sin x и осью Ох при . Построим данную фигуру (рис. 7.7)

Рис. 7.7.

Используя формулу (7.5), находим искомую площадь фигуры

(кв. ед.)

2. Вычислить площадь фигуры, ограниченную линиями у = х, , х = 2.

Изобразим данную фигуру (рис. 7.8)

Рис. 7.8.

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

(кв. ед.)

3. Вычислить площадь фигуры, ограниченную линями и .

Построим фигуру, площадь которой необходимо вычислить (рис. 7.9).

Рис. 7.9.

Находим точку пересечения линий и . Решаем уравнение

, имеем или , откуда . Из рисунка видно, что пределами интегрирования являются . Определяем площадь фигуры, используя формулу (7.6)

(кв. ед.)

Параметрические координаты

Если криволинейная трапеция ограничена кривой, заданной параметри-чески

,

прямыми и и осью Ох, то её площадь определяется по формуле

, (7.7)

где α и β определяются из равенств и .

Примеры.

1. Найти площадь фигуры ограниченной эллипсом, который задан в параметри-ческой форме

.

Изображаем фигуру, площадь которой необходимо определить (рис. 7.10)

Рис. 7.10.

Найдём четвёртую часть площади S, расположенной в первой четверти ко-ординатной плоскости (на рисунке она окрашена серым цветом). Здесь х изменяется от 0 до а и, следовательно, t изменяется от π / 2 до 0. По формуле (7.7) находим

=

= . Таким образом, . Значит .

2. Вычислить площадь фигуры, ограниченной одно аркой циклоиды и осью Ох.

Циклоида есть траектория точки, расположенной на ободе колеса радиуса а , при равномерном качении колеса по оси Ох. При одном обороте колеса центр колеса переместится на расстояние (рис. 7.11).

Рис. 7.11.

Уравнение циклоиды в параметрической форме имеет вид

,

При изменении параметр t изменяется в пределах . По формуле (7.7) находим искомую площадь

=

= =

= .

Полярные координаты

В некоторых случаях вычисление площадей криволинейных фигур удобно проводить в полярных координатах.

Полярная система координат определяется заданием точки О (полюс)луча Ор, исходящего из точки О (полярной оси) и масштаба для измерения длины. Положение точки М на плоскости определяется в полярной системе координат двумя числами: полярным радиусом (рис. 7.12), выражающим длину отрезка ОМ в выбранном масштабе, и полярным углом φ = .

Рис. 7.12.

Из рис. 7.10 видно, что независимо от расположения точки М на плоскости имеют место следующие формулы перехода от полярных координат ( ) к декартовым (х, у) : , ;

и от декартовых к полярным: , .

Найдём площадь криволинейного сектора.Пусть кривая АВ за­дана в полярных координатах уравнением , , причем функция непрерывна и неотрицательна на отрезке [α, β]. Плоскую фигуру, ограниченную кривой АВ и двумя лучами, составляющими с полярной осью углы α и β, будем называть кри­волинейным сектором (рис. 7.13). Площадь S криволинейного сек­тора определяется формулой

. (7.8)

Рис. 7.13.

Доказательство. Разобьем произвольно отрезок [α, β] иа п частей точками ол , выберем на каждом частичном отрезке [ ] произвольно точку ( ) и построим круговые сек­торы с радиусами . В результате получим веерообразную фигуру, площадь которой приближенно равна площади S криволи­нейного сектора:

,

где . С другой стороны, площадь веерообразной фигуры является интегральной суммой для интеграла (7.8). Так как функция непрерывна на отрезке [α, β ], то предел этой суммы при существует и равен интегралу (7.8).

Следовательно, и площадь криволинейного сектора численно равна этому определенному интегралу:

.

Пример 1.Вычислить площадь фигуры, ограниченной поляр­ной осью и первым витком спирали Архимеда: , где а — положительное число (рис. 7.12).

Рис. 7.12.

Решение. При изменении от 0 до 2π полярный радиус описывает кривую, ограничивающую криволинейный сектор ОАВС. Поэтому по формуле (7.8) имеем

.

Расстояние от точки С до полюса равно . Поэтому круг радиуса ОС имеет площадь πOC2 = 4 π3 a2 = 3 ∙ , т. е. площадь фигуры, ограниченной полярной осью и первым витком спирали Архимеда, равна 1/3площади круга с радиусом, равным наибольшему из полярных радиусов витка. К этому вы­воду пришел еще Архимед.

Пример 2.Второй закон Кеплера (закон площадей) о движении планет солнечной системы гласит: площадь, описываемая радиусом-вектором планеты, проведенном из центра Солнца, возрастает пропорционально времени.

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

Рис. 7.15.

Рассмотрим перемещение планеты в окрестностях точек А (афелий) и П (перигелий), по закону Кеплера площади секторов и равны между собой, т.е.

,

где - площадь сектора, опирающегося на дугу , длина этой дуги равна , аналогично - площадь сектора, опирающегося на дугу , длина этой дуги равна .

Из формулы для площади криволинейного сектора (7.8) следует, что

или , т.е. . Здесь и перемещение планеты за один и тот же промежуток времени в окрестности точек А и П орбиты. Разделим предыдущее равенство на промежуток времени :

.

Отношение перемещения планеты ко времени есть скорость планеты в точке А, аналогично отношение перемещения ко времени есть скорость планеты в точке П, т.е.

, .

В результате или , откуда следует, что

.

 




Поиск по сайту:

©2015-2020 studopedya.ru Все права принадлежат авторам размещенных материалов.