Файл: Крылов, В. И. Методы приближенного преобразования Фурье и обращения преобразования Лапласа справочная книга.pdf
ВУЗ: Не указан
Категория: Не указан
Дисциплина: Не указана
Добавлен: 22.10.2024
Просмотров: 103
Скачиваний: 0
В. И. КРЫЛОВ, Н. С. СКОБЛЯ
МЕТОДЫ
ПРИБЛИЖЕННОГО ПРЕОБРАЗОВАНИЯ ФУРЬЕ И ОБРАЩЕНИЯ ПРЕОБРАЗОВАНИЯ ЛАПЛАСА
СПРАВОЧНАЯ КНИГА
ИЗДАТЕЛЬСТВО «НАУКА»
ГЛАВНАЯ РЕДАКЦИЯ ФИЗИКО-МАТЕМАТИЧЕСКОЙ ЛИТЕРАТУРЫ
М о с к в а 1974
517.2 К 85
УДК 517.5
A t
I ъа, г, С’Лк
»-'?Гчне.тл<иИ.(
Ы&!НГ1 .. :<f,
4.V
у/сГ/ S ¥ y -3 S ¥ ¥■ £ |
||
4 л |
' |
**'“^»»** |
Методы приближенного преобразования Фурье и об |
||
ращения преобразования |
Лапласа (справочная книга). |
|
В. И. К р ы л о в , |
Н. С. |
С к о б л я . Главная редакция |
физико-математической литературы издательства «Нау ка», М., 1974.
Гармонический анализ и преобразование Лапласа очень часто применяются для решения многих теорети ческих и прикладных вопросов.
В книге содержится изложение большинства извест ных методов приближенного обращения преобразования Лапласа и вычисления интегралов Фурье.
Книга предназначена для научных и инженернотехнических работников, которым в их деятельности приходится иметь дело с теорией или приложениями преобразования Лапласа и интегралов Фурье. Она будет полезным справочником для работников вычислительных центров и конструкторских бюро.
В книге 2 рис., библ. 11 названий.
©Издательство «Наука», 1974.
20204—091 К 053(02)-74 83-74
ОГЛАВЛЕНИЕ
Предисловие.............................. ...................................................... 5
ЧАСТЬ ПЕРВАЯ
ОБРАЩЕНИЕ ПРЕОБРАЗОВАНИЯ ЛАПЛАСА
Глава |
1. Введение |
................................................................................... |
понятия |
|
|
|
|
|
9 |
||||
§ |
1.1. |
Основные |
теории преобразования Лапласа |
||||||||||
(9). § 1.2. Комплексный интеграл, осуществляющий обраще |
|||||||||||||
ние преобразования Лапласа (15). |
§ 1.3. |
Представление функ |
|||||||||||
ций интегралом Лапласа (18). § 1.4. Некорректность задачи |
|||||||||||||
обращения преобразования Лапласа (22). |
|
|
|
|
|||||||||
Глава 2. Некоторые аналитические методы обращения преобра |
|||||||||||||
|
зования Лапласа.................................................................. |
|
|
|
|
|
24 |
||||||
§ 2.1. |
Нахождение оригинала |
с помощью формулы обра |
|||||||||||
щения (24). § 2.2. Разложение оригинала в степенные |
ряды |
||||||||||||
(27). |
§ 2.3. |
Разложение |
оригинала в обобщенные степенные |
||||||||||
ряды (29). |
|
|
|
|
|
|
|
|
|
|
|
|
|
Глава 3. |
Методы численного обращения преобразования Лап |
||||||||||||
|
|
ласа, |
основанные на использовании специальных раз |
||||||||||
§ 3.1. |
ложений ................................................................................... |
|
|
|
|
|
|
|
|
31 |
|||
Обращение преобразования Лапласа с помощью |
|||||||||||||
многочленов, |
ортогональных |
на |
конечном |
промежутке |
(31). |
||||||||
§ 3.2. |
Обращение |
преобразования |
Лапласа |
с помощью |
ряда |
||||||||
Фурье |
по |
синусам (52). § 3.3. Обращение преобразования Лап |
|||||||||||
ласа с помощью рядов по обобщенным многочленам Чебышева — |
|||||||||||||
Лагерра (54). |
|
|
|
|
|
|
|
|
|
|
|
||
Глава |
4. |
Методы |
вычисления интеграла Меллина при помощи |
||||||||||
§ 4.1. |
интерполяционных ....................квадратурных формул |
59 |
|||||||||||
Общая теория |
интерполяционных |
методов (59) § 4.2. |
|||||||||||
Интерполяционный |
метод с равноотстоящими узлами (62) . § 4.3. |
||||||||||||
Интерполяционный |
метод |
с |
неравноотстоящими |
узлами |
(63). |
||||||||
§ 4.4. Замечания о других интерполяционных методах. Приме |
|||||||||||||
нение |
отрезка |
ряда |
Тейлора |
(72). |
§ 4.5. |
Некоторые теоремы |
|||||||
о сходимости интерполирования (73). § 4.6. |
Теоремы о сходимо |
||||||||||||
сти интерполяционных методов обращения (85). |
|
|
|||||||||||
Глава 5. Методы численного обращения преобразования Лапласа |
|||||||||||||
|
при |
помощи квадратурных формул, |
имеющих наивыс |
||||||||||
|
шую степень ......................................................... |
точности |
|
(90). |
|
90 |
|||||||
§ 5.1. |
Теория |
квадратурных формул |
§ 5.2. Ортого |
||||||||||
нальные многочлены, |
связанные с квадратурной формулой наи- |
1*
4 |
ОГЛАВЛЕНИЕ |
высшей степени точности (97). § 5.3. Методы вычисления коэф фициентов и узлов квадратурной формулы (112).
Глава 6. Методы обращения преобразования Лапласа с помощью квадратурных формул с равными коэффициентами . . . 121
§ 6.1. Построение вычислительной формулы (121). § 6.2. Замечание о расположении узлов (124).
ч а с т ь в т о р а я
ПРЕОБРАЗОВАНИЯ ФУРЬЕ И ИХ ПРИМЕНЕНИЕ К ОБРАЩЕНИЮ ПРЕОБРАЗОВАНИЯ ЛАПЛАСА
Глава 7. |
Введение |
................................................................................... |
Фурье |
(126). |
§ 7.2. |
Приведение |
125 |
|
§ 7.1. |
Преобразования |
|
||||||
интеграла типа Меллина к преобразованию Фурье (131). |
|
|||||||
Глава 8. |
Обращение преобразования Лапласа с помощью ряда |
133 |
||||||
|
Фурье......................................................................... |
|
убывающего оригинала |
f (х) (133). |
||||
§8.1. Случай быстро |
|
|||||||
§ 8.2.Случай быстрого убывания модуля изображения F (р) (135). |
|
|||||||
Глава 9. |
Интерполяционные формулы для вычисления интегра |
137 |
||||||
§ 9.1. |
лов Ф урье............................................................................. |
|
|
|
|
(137). § 9.2. |
||
Несколько предварительных замечаний |
|
|||||||
Вычислительные |
формулы , |
основанные |
на |
алгебраическом |
|
|||
интерполировании |
функции |
f (х) (138). § 9.3. |
Вычислительные |
|
||||
формулы, основанные на интерполировании рациональными |
|
|||||||
функциями (166). |
|
|
|
|
|
|
|
|
Глава 10. |
Формулы для вычислений, |
имеющие наивысшую сте |
192 |
|||||
|
пень точности . . |
. . ......................................................... |
||||||
§ 10.1. Введение (192). |
§ 10.2. |
Построение формулы наи |
|
|||||
высшей степени точности (194). |
|
|
|
|
|
ч а с т ь т р е т ь я
ВЫДЕЛЕНИЕ ОСОБЕННОСТЕЙ ФУНКЦИИ ПРИ ВЫЧИСЛЕНИЯХ
Глава 11. Выделение особенностей изображения F ( р ) ................ |
|
202 |
||
§ 11.1. Введение (202). |
§ 11.2. Устранение |
и ослабление |
|
|
особенностей изображения F (р) (204). § 11.3. Замечание об уве |
|
|||
личении скорости стремления |
к нулю изображения F (р) |
(208). |
|
|
§ 11.4. Таблица изображений F (р) и соответствующих оригина |
|
|||
лов f (х) для построения особой части изображения Ft (р) (209). |
|
|||
Глава 12. Выделение особенностей функции при преобразовании |
213 |
|||
Фурье...................................................................................... |
|
(214). § |
12.2. |
|
§ 12.1. Устранение разрывов первого рода |
|
|||
Увеличение скорости стремления к нулю преобразуемой функ |
|
|||
ции (217). |
|
|
|
220 |
Л итература............................................................................................. |
|
|
|
|
Список обозначений ............................................................................ |
|
|
|
221 |
Предметный указатель . ...................................................................... |
222 |
ПРЕДИСЛОВИЕ
Проблема приближенного обращения преобразования Лапласа и, в особенности, численного его обращения воз никла из потребности довести решение до числа в том случае, когда существующие таблицы функций и их изо бражений не дают возможности по изображению найти оригинал или требуют очень больших вычислений.
Было построено, особенно за два последних десятиле тия, много способов обращения. Они опубликованы раз розненно в специальных журналах и книгах и известны, как правило, неширокому кругу лиц. Насколько осведом лены авторы, в научной литературе нет книг, содержащих систематическое изложение всех этих методов *).
Авторы хотели написать книгу, где было бы дано достаточно полное описание современного состояния про блемы обращения, и хотели сделать книгу полезной тем, кто применяет преобразование Лапласа к решению тех или иных задач. Первое из этих намерений можно было осу ществить просто, так как литература по проблеме обраще
ния |
еще невелика и сравнительно |
легко |
обозрима. |
Что |
же |
касается вопроса о полезности |
книги, |
то самый |
бес |
пристрастный ответ здесь даст время, но авторы считают, что в этом отношении в адрес книги может быть сделано несколько замечаний. Наиболее важным из них, несом ненно, будет справедливое замечание о недостаточной исследованности проблемы и о еще малом числе получен ных результатов. Чтобы быть лучше понятыми читателя ми, не знакомыми с проблемой обращения преобразования
*) Некоторым исключением здесь является книга [8], но ее целью было объединение всех известных вспомогательных численных таблиц, предназначенных для вычисления интеграла Меллина. Математическая же теория обращения преобразования Лапласа в ней дама в очень сокращенном виде в форме краткого справочника, поясняющего таблицы.
6 ПРЕДИСЛОВИЕ
Лапласа, авторы должны здесь сделать некоторые пояс нения. Для этого авторы прибегли к нестрогим, но на глядным соображениям.
Проблема обращения есть не что иное, как задача отыскания решения f(x) интегрального уравнения первого рода
|
|
СО |
|
|
|
\f{x)e-pxdx=F{p), |
(*) |
|
|
о |
|
где F (р) считается известной функцией комплексного аргу |
|||
мента р, |
аналитической в некоторой полуплоскости вида |
||
R e p > у (у < о о ). |
|
|
|
Равенство (*) есть преобразование Лапласа функции |
|||
f(x) в F (р). Ядро |
интеграла егр,с является целой |
анали |
|
тической |
функцией |
аргументов х и р с плавным |
измене |
нием, и операция интегрирования, которая выполняет усреднение / с весом е~рх, может значительно сгладить особенности в поведении преобразуемой функции /.
В задаче обращения нужно по изображению F (р), изменяющемуся, ввиду его аналитичности, очень плавно, восстановить все возможные неровности поведения ориги нала f{x). Поэтому можно ожидать, что аппарат, при помощи которого решается задача обращения, не может быть простым и грубым, а должен быть достаточно слож ным и чувствительным даже к малым оттенкам в поведе нии F (р), чтобы выполнить «тонкую работу» восстановле ния оригинала f(x).
Обратим внимание еще на одно свойство проблемы обращения — на неустойчивость оригинала / относительно малых изменений изображения F. Такая неустойчивость
становится |
видной при |
первом же |
взгляде |
на |
преобразо |
|
вание (*). |
В самом деле, |
пусть |
/ — заданный |
оригинал |
||
и F —соответствующее |
ему |
изображение. |
Подвергнем / |
любому изменению, быть может даже сильному, но на очень малом отрезке. Новый оригинал назовем fv Такое изменение незначительно повлияет на интеграл и мало изменит изображение F (р), так что новое изображение Рг(р) будет близким к F (р).
Мы могли бы построить новый оригинал flt изменяя / не на одном участке полуоси 0= ^х -< оо , а на многих
.ее участках, но так, чтобы общая сумма длин этих участков
ПРЕДИСЛОВИЕ |
7 |
была бы достаточно малой величиной, все равно при этом новое изображение Ft (р) можно было бы сделать сколь угодно близким к F (р), так что близким изображениям F (р) и F1(p) могут отвечать оригиналы f (х) и Д (л:), сильно отличающиеся между собой на многих участках полуоси О г ^ х с о о .
Теперь можно более полно объяснить, в каком отноше нии содержание книги следует считать недостаточным, и попытаться указать причины, по которым авторы нахо дят, что эти недостатки вряд ли могут быть устранены
вближайшие годы.
Внастоящее время мы находимся в самой начальной стадии построения теории приближенного обращения пре образования Лапласа. В большинстве публикуемых работ либо указываются видоизменения уже известных методов вычислений, либо создаются новые вычислительные схемы, не имеющие аналогов в прошлом. Примерами методов первого вида могут служить методы вычисления комплекс ного интеграла Меллина, основанные на идее интерполи рования или на. идее Гаусса достижения наивысшей сте пени точности, и др. Использование же ортогональных многочленов для обращения преобразования возникло только в связи с преобразованием Лапласа и не имело близких аналогов в прошлом.
Построение вычислительного метода или указание на возможность построения есть лишь первый шаг в созда нии теории метода. За ним должно следовать выяснение условий сходимости правила, определение скорости схо димости, нахождение оценки погрешности априорной и апостериорной, выработка способов улучшения сходимости, если последняя окажется недостаточно быстрой, и т. д. Исследования проблемы обращения в этом направлении особенно важны ввиду ее неустойчивости.
Как читатель увидит из содержания книги, результаты
такого рода получены в очень небольшом числе и лишь в самых простых случаях.
Главными причинами, препятствующими быстрому реше нию стоящих здесь задач, являются указанные выше два свойства проблемы обращения: неизбежная сложность математического аппарата обращения преобразования Лап ласа и неустойчивость проблемы обращения относительно изменений функций F, которая должна вызвать ту или
8 ПРЕДИСЛОВИЕ
иную форму неустойчивости в любом вычислительном про цессе решения этой проблемы.
Недостаточная изученность правил приближенного обращения преобразования делает справочную книгу непол ной во многих отношениях. По-видимому, трудно ожи дать быстрого улучшения наших знаний в этих вопросах, и в ближайшие годы книгу, вероятно, невозможно будет освободить от такой неполноты.
До сих пор говорилось только об обращении преобра зования Лапласа, о гармоническом же анализе не упомина лось вообще. Это связано с тем, что в книге гармониче ский анализ рассматривается только в той мере, в какой он связан с проблемой обращения, а именно, из гармони ческого анализа взяты только задачи вычисления интегра лов Фурье, тесно связанные с обращением преобразова ния Лапласа, и вычисление интегралов Фурье есть для нас один из возможных методов решения проблемы обра щения. Поэтому все сказанное об обращении относится и к вычислению интегралов Фурье.
Книга предназначена широкому кругу читателей, зани мающихся теорией преобразования Лапласа или его науч ными и техническими приложениями. Поэтому авторы не стремились к лаконичности изложения и старались сделать его доступным не только для специалистов-математиков. Книга предполагает знакомство читателя с основами анализа и теории функций комплексного переменного в объеме втузовского курса математики с расширенной
программой.
Авторы
Ч А С Т Ь П Е Р В А Я
ОБРАЩЕНИЕ ПРЕОБРАЗОВАНИЯ ЛАПЛАСА
Г Л А В А 1
ВВЕДЕНИЕ
§1.1. Основные понятия теории преобразования Лапласа
Втечение нескольких последних десятилетий в мате матике, механике и технике для решения многих задач стали особенно часто и успешно применяться операцион ные методы на основе преобразования Лапласа. Эти методы
нашли |
широкое |
приложение |
в теории теплопроводности, |
|
в электро- и радиотехнике, |
при |
исследовании нестацио |
||
нарных |
явлений |
в электрических |
цепях, в задачах дина |
мики систем автоматического регулирования, в теории линейных дифференциальных, интегральных и разностных уравнений и во многих других задачах.
Операционный метод решения задач можно подразде лить на 4 этапа:
1)от искомой функции-оригинала f(t) переходят к функ ции-изображению F (р);
2)над F (р) производят операции, соответствующие
операциям над / (/), после чего получают уравнение отно сительно F (р), которое часто бывает значительно проще уравнения для оригиналов, например, обыкновенное диф ференциальное уравнение заменяется алгебраическим, ур в- нение в частных производных — обыкновенным дифферен циальным и т. д.;
3)полученное уравнение для изображений решают относительно F (р);
4)от найденного изображения F (р) переходят к ориги налу f{t), который и является искомой функцией.
Во многих случаях самым |
трудным является 4-й |
этап — нахождение оригинала f (t) |
по изображению F (р), |
д. е. задача обращения преобразования Лапласа. Имеются