Файл: Крылов, В. И. Методы приближенного преобразования Фурье и обращения преобразования Лапласа справочная книга.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 (р),

д. е. задача обращения преобразования Лапласа. Имеются