Главная Информатика
Вычисление интегралов в Mathcad
|
|
|||||
ИнтегрированиеИнтегрирование в MathCAD реализовано в виде вычислительного оператора. Допускается вычислять интегралы от скалярных функций в пределах интегрирования, которые также должны быть скалярами. Несмотря на то, что пределы интегрирования обязаны быть действительными, подынтегральная функция может иметь и комплексные значения, поэтому и значение интеграла может быть комплексным. Если пределы интегрирования имеют размерность, то она должна быть одной и той же для обоих пределов. Оператор интегрированияИнтегрирование, дифференцирование, как и множество других математических действий, устроено в MathCAD по принципу "как пишется, так и вводится". Чтобы вычислить определенный интеграл, следует напечатать его обычную математическую форму в документе. Делается это с помощью панели Calculus (Вычисления) нажатием кнопки со значком интеграла или вводом с клавиатуры сочетания клавиш <Shift>+<7> (или символа "&", что то же самое). Появится символ интеграла с несколькими местозаполнителями (рис. 1), в которые нужно ввести нижний и верхний интервалы интегрирования, подынтегральную функцию и переменную интегрирования. Можно вычислять интегралы с одним или обоими бесконечными пределами. Для этого на месте соответствующего предела введите символ бесконечности, воспользовавшись, например, той же самой панелью Calculus (Вычисления). Чтобы ввести - (минус бесконечность), добавьте знак минус к символу бесконечности, как к обычному числу. ![]() Рис. 1. Оператор интегрирования Чтобы получить результат интегрирования, следует ввести знак равенствили символьного равенства. В первом случае интегрирование будет проведено численным методом, во втором - в случае успеха, будет найдено точное значение интеграла с помощью символьного процессора MathCAD. Конечно, символьное интегрирование возможно только для небольшого круга несложных подынтегральных функций. Результат численного интегрирования - это не точное, а приближенное значение интеграла, определенное с погрешностью, которая зависит от встроенной константы TOL. Чем она меньше, тем с лучшей точностью будет найден интеграл, но и тем больше времени будет затрачено на расчеты. Поумолчанию TOL=O.OOI. Для того чтобы ускорить вычисления, можно установить меньшее значение TOL. Отдавайте себе отчет в том, что при вводе в редакторе MathCAD оператора численного интегрирования, вы, фактически, создаете самую настоящую программу. Пользователь имеет возможность выбирать сам алгоритм численного интегрирования. Для этого:
Обратите внимание, что перед тем как один из алгоритмов выбран впервые, как показано на рис. 2, флажок проверки в контекстном меню установлен возле пункта AutoSelect (Автоматический выбор). Это означает, что алгоритм определяется MathCAD, исходя из анализа пределов интегрирования и особенностей подынтегральной функции. Как только один из алгоритмов выбран, этот флажок сбрасывается, а избранный алгоритм отмечается точкой. ![]() Рис. 2. Выбор алгоритма численного интегрирования Разработчиками MathCAD 2001 запрограммированы четыре численных метода интегрирования:
Старайтесь все-таки оставить выбор численного метода за MathCAD, установив флажок AutoSelect (Автоматический выбор) в контекстном меню. Попробовать другой метод можно, например, чтобы сравнить результаты расчетов в специфических случаях, когда у вас закрадываются сомнения в их правильности. Если подынтегральная функция "хорошая", т. е. не меняется на интервале интегрирования слишком быстро и не обращается на нем в бесконечность, то численное решение интеграла не принесет никаких неприятных сюрпризов. Приведем основные идеи итерационного алгоритма Ромберга, который применяется для большинства таких функций.
Если интеграл расходится (равен бесконечности), то вычислительный процессор MathCAD может выдать сообщение об ошибке, выделив при этом оператор интегрирования, как обычно, красным цветом. Чаще всего ошибка будет иметь тип "Found a number with a magnitude greater than 10Л307" (Найдено число, превышающее значение 10307) или "Can't converge to a solution" (Не сходится к решению), как, например, при попытке вычислить интеграл [- dx Тем не менее, символьный процессор справляется с этим Jo Vx интегралом, совершенно правильно находя его бесконечное значение. |
<< | СОДЕРЖАНИЕ | >> |
---|