Файл: Крылов, В. И. Методы приближенного преобразования Фурье и обращения преобразования Лапласа справочная книга.pdf
ВУЗ: Не указан
Категория: Не указан
Дисциплина: Не указана
Добавлен: 22.10.2024
Просмотров: 108
Скачиваний: 0
212 |
В Ы Д Е Л Е Н И Е О С О Б Е Н Н О С Т Е Й И З О Б Р А Ж Е Н И Я F (р) |
[ Г Л , И |
|
F ( P ) |
|
|
|
— \ - ё - е ~ а Р , |
|
а > 0 |
|
|
Р + Р |
|
|
|
|
^-Р + 11 |
е - а р |
* |
а > |
0 |
р 2 + р 2 |
е |
° > |
и |
|
L P + R 2 < ra p > |
а > 0 |
|
||
Р2 — Р 2 |
|
|
|
|
Продолжение
|
|
|
|
|
Ц х ) |
|
( |
о, |
|
0 |
< |
х С |
а |
1 |
е~ Р |
|
а |
< |
ж |
|
1" |
0 , |
0 |
< ж |
< |
а |
|
< |
A, co s |
[Р |
(дс— а )] |
+ |
ц р - 1 sin [Р (ж — а ) ] , |
|
V- |
а |
< х |
/■ 0 , |
0 < х < а |
|
| A c h [р (л: — a )] + n p - i Sh [ р ( ж — а ) ] , |
||
' |
а < |
ж |
Г Л А В А 12
ВЫДЕЛЕНИЕ ОСОБЕННОСТЕЙ ФУНКЦИИ ПРИ ПРЕОБРАЗОВАНИИ ФУРЬЕ
Вычисление интегралов Фурье по значениям преобра зуемой функции / основано, как и вычисление интеграла Меллина, преимущественно на интерполировании функции f при помощи многочленов или рациональных функций, при этом интерполирование выполняется либо во всей области интегрирования, либо в ее частях. Точность вычисления интегралов зависит как от избранного правила, так и от свойств функции /, но обычно, чем более гладкой является функция f и чем быстрее она стремится к нулю при уда лении аргумента ее в бесконечность, тем более точный результат может быть получен при вычислениях.
Предварительную подготовку функции f к вычислению интегралов Фурье можно проводить в следующих трех направлениях:
1) Улучшение дифференциальных свойств функции /,
вчастности повышение порядка ее дифференцируемости.
2)Увеличение плавности изменения функции /. Чтобы пояснить содержание этой краткой фразы, достаточно, повидимому, привести примеры. Хорошо известна возмож ность сколь угодно точного и равномерного приближения на замкнутом отрезке непрерывной функции многочленом.
Внекоторых случаях для достижения заданной точности приближения необходимо строить многочлен высокой сте пени с большими и сильно изменяющимися значениями производных от него. С таким затруднением можно встре титься даже в случае аналитической функции /, когда ее особые точки лежат вблизи отрезка, на котором строится аппроксимирующий многочлен.
В случаях такого рода естественно пытаться облегчить
задачу приближения путем предварительного выделения
214 |
ВЫДЕЛЕНИЕ ОСОБЕННОСТЕЙ ФУНКЦИИ |
[ГЛ. 12 |
|
из / «особой, негладкой части» с тем, чтобы вычисления |
|||
нужно |
было производить с более плавно |
изменяющимся |
|
«остатком». С задачей такого увеличения плавности |
изме |
||
нения мы встречались выше в § 11.2. |
|
|
|
3) |
Ускорение стремления к нулю |
функции |
/, когда |
|х | неограниченно возрастает.
Мы остановимся только на двух вопросах: устранении разрывов первого рода и ускорении стремления к нулю преобразуемой функции.
§ 12.1. Устранение разрывов первого рода
Пусть функция / задана на отрезке [а, Ь] и х есть внутренняя точка этого отрезка. Предположим, что f имеет в точке х правостороннее предельное значение
/(х + 0) = lim |
f(x') |
х' -*■ X, X' > X
илевостороннее предельное значение
fix —0) = |
lim |
fix'). |
|
х' -* X, X' < X |
|
Если значения f (л; 4- 0) и / (х — 0) существуют, но хотя |
||
бы одно из них отлично от f (х), то |
говорят, что f имеет в |
|
точке х разрыв непрерывности первого рода. |
||
Если точка х является одним из |
концов отрезка [а, Ь\у |
|
то при определении разрыва рассматривается лишь одно из предельных значений: /(х + О) или fix —0).
Почти аналогично определяются разрывы первого рода у производных функции /. Например, рассмотрим первую производную f и будем считать, что она существует во всех точках некоторой окрестности значения х, исключая, может быть, саму точку х. Предположим также, что существуют предельные значения f' (х + 0) и /' (х — 0). Если окажется, что /' (х —0) ^ / '( х + 0), то говорят, что первая производная f имеет в точке х разрыв первого рода со скач ком f ( х + 0) — /' (х —0).
Для определенности записи будем иметь в виду косинуспреобразование. Предположим, что / и ее производные до
некоторого порядка т непрерывны на |
полуоси [0, |
оо) |
|
всюду, |
за исключением q точек х}, х2, |
... , хч (х0 = |
0-< |
< хх < |
х2 < . . . < xq), где они имеют разрывы первого рода. |
||
§ 12.1] |
УСТРАНЕНИЕ РАЗРЫВОВ ПЕРВОГО РОДА |
215 |
Величину скачка f(‘>(х) в точке х. обозначим
ttp = P ( x !+ 0 ) - P ( x l - 0 )
(/ = 0, j= l,2 , ... ,q ) .
Когда необходимо устранить разрывы только у функ ции /, достаточно ввести кусочно-постоянную функцию F0(x), абсолютно интегрируемую на [0, со) и имеющую такие же скачки, как и /. За такую функцию можно, очевидно, взять
|
|
2 [ E (x - X j) - \] k 7 , |
|
|
(12.1.1) |
||
|
|
/=1 |
|
|
|
|
|
где Е (х) есть «гасящая функция» (см. стр. |
143). |
разрыва |
|||||
На отрезках полуоси |
[0, со) между точками |
||||||
F0 имеет указанные |
ниже значения: |
|
|
|
|
||
|
(^0,+ ^ 0,+ ---+ ^0,)= а '0), |
0 < * < |
хг |
|
|
||
|
- W ' + ... + k'*‘)=a»\ |
Xl< x < х2, |
|
( 12. 1. 2) |
|||
Fо М = ' |
|
|
kq— ajj\ |
< х < Xq, |
|
||
|
|
|
|
|
|||
|
|
|
0, |
Xq< X. |
|
|
|
Косинус-преобразование функции F0 находится просто: |
|||||||
оо |
q |
|
xj |
|
|
|
|
jj F0 (х) cos рх clx= 2 |
аТ |
] cos Рх dx — |
|
|
|
||
о |
/=i |
*/_! |
^ |
q |
|
|
|
|
q |
|
|
|
|
||
= |
p- 2 a'i ‘ (sin PXJ ~ sin PXJ - J |
= " T |
2 |
^ |
sin pXF |
||
|
И /= l |
|
|
|
/= i |
|
|
Разность |
Ф0 (x)=f(x)—F0(x) непрерывна всюду, |
кроме то |
|||||
чек Xj, |
а в этих точках для нее ф0 (х,-+ |
0) = |
ф0 (х,-— 0). |
||||
Если ее доопределить в точках Xj, |
положив там ф0 (xj) = |
||||||
=фо (* ;+ 0), она будет непрерывна везде на полуоси [0, оо). Так как вне точек х<функция F0 (х) кусочно-постоян
на, |
то всюду, кроме точек Xj, производные всех поряд |
ков |
от фо (х) и f (х) совпадают, если производные от / су |
ществуют.
Если мы хотим устранить у f только разрывы первой
производной |
оставив неизменными разрывы самой функ |
ции / и ее |
производных выше первого порядка, можно |
215 |
ВЫДЕЛЕНИЕ ОСОБЕННОСТЕЙ ФУНКЦИИ |
[ГЛ. 12 |
воспользоваться функцией *) |
|
|
|
F1(x)= £ [ E ( x - x ; ) - \ ] k ' y ( x - x j ) . |
(12.1.3) |
|
/= i |
|
Она кусочно-линейная, непрерывная на [0,оо) и обращается
в |
тождественный нуль |
при x > x q. Первая производная |
от |
нее |
|
|
F[(x)= j ] |
[ E ( x - x j ) ~ 1]А“' |
; = i
кусочно-постоянна вне точек Xj и имеет в этих точках ска чки, равные соответственно величинам к)\
Разность |
<Pi {х) = f (х) — |
(х) будет |
иметь производную |
|||||||
ф !(х) везде |
при х ф х Л ) |
= |
1,2...... q). Кроме |
того, |
в точ |
|||||
ках х, |
будет cpi (ху + 0) = |
(pj (Xj —0), |
и |
если |
доопределить |
|||||
ф!(х) в точках Xj, положив там ф\(xj) |
= |
{xj-\- 0), функ |
||||||||
ция ф1 (х) будет непрерывной на полуоси [0,оо). |
|
|||||||||
Косинус-преобразование |
функции |
|
Fl {x) находится |
|||||||
с помощью несложных вычислений. Так как |
|
|
||||||||
СО |
|
|
|
|
|
|
|
|
|
|
^ [Е {х — Xj) — 1] (х — Xj) cos рх dx = |
|
|
|
|
|
|||||
о |
|
|
х! |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
= ^ {Xj — х) cos pxdx == р2 (1 — cos pxj), |
|||||||
то |
|
|
о |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
оо |
|
|
|
|
|
|
|
|
|
|
^ Fx(х) cos рх dx = |
|
|
|
|
|
|
|
|
||
о |
q |
со |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
= |
^ |
^ ^ |
('r ~ */) — 11 (* — Xj) cos рх dx = |
|
||||||
|
/ = 1 |
о |
|
|
ч |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
= ? |
2 |
- z ° s PXj). |
(12.1.4) |
||||
|
|
|
|
|
/=1 |
|
|
|
|
|
Отметим также, |
что разность |
|
|
|
|
|
||||
Фи (x) = f( x ) - F 0(x)~F1(x)
*) Ввиду простоты мы не поясняем соображения, на основании Которых строится /д (х).
§ |
12.2] |
УВЕЛИЧЕНИЕ СКОРОСТИ СТРЕМЛЕНИЯ К НУЛЮ |
217 |
|||||||||
будет |
обладать |
свойствами: |
Ф01 |
+ 0) = |
Ф01 (-^у — 0), |
|||||||
tpo'i (х/ + |
0) = |
<po'i (xj — 0), и если ср01 доопределить в точках |
||||||||||
Xj |
равенствами cp01 (xj) = <p0i (xj + 0), то <р01 |
|
будет |
непре |
||||||||
рывно дифференцируемой при х ^ О . |
|
|
При |
всяком |
||||||||
|
Эти |
рассуждения |
можно |
продолжить. |
|
|||||||
i = 0, 1 |
,т положим |
|
|
|
|
|
|
|
|
|||
|
|
Fi (x)= 2 |
[ £ |
( * - |
* |
, ) - |
( 1 |
2 |
. |
1 |
. 5 ) |
|
|
|
|
/= 1 |
|
|
|
|
|
|
|
|
|
|
Такая функция Ft непрерывно дифференцируема по л: |
|||||||||||
до порядка |
/ — 1 |
включительно, |
производная |
же порядка |
||||||||
i |
при х ф х , |
( / = 1 , 2 , . . . , (?) |
имеет значение |
|
|
|
||||||
|
|
|
Ff]М = |
2 [ E ( x - Xj) - |
l]kjC) |
|
|
|
|
|||
|
|
|
|
|
i= i |
|
|
|
|
|
|
|
и |
является |
кусочно-постоянной |
функцией |
со |
скачками |
|||||||
в точках Xj, величины которых равны соответственно kSp.
Для х > х д верно равенство Fff>(x) = 0. |
|
|||||||
|
С помощью Fi (х) у функции / могут быть удалены |
|||||||
скачки производной |
порядка i. |
|
|
|||||
|
Если |
же |
мы положим |
|
|
|
||
|
|
Ф01 |
m(x)=f (х) — [F0(x) -\-F1(x) |
. .-jFт(х)] |
||||
и |
затем |
в |
точках |
Xj |
доопределим |
ф01... т, полагая |
||
Фм ... т(лу)= |
Ф01... т (*/ + 0)> получим функцию, т раз не- |
|||||||
прерывно |
дифференцируемую на [0, оо). Этим будут устра |
|||||||
нены |
разрывы у / и ее производных до порядка т. |
|||||||
|
§ |
12.2. Увеличение скорости |
стремления к нулю |
|||||
|
|
|
|
преобразуемой функции |
|
|||
|
Пусть |
f(x) стремится |
к нулю по степенному закону |
|||||
так, что для некоторого s будет (х + |
a)s f (х) -»■ С Ф 0 (а > 0). |
|||||||
Во |
многих |
случаях |
параметр а, |
зависящий от поведения |
||||
f вблизи начала координат, можно принять равным еди нице. Функция f представима в форме
/ (х) — + /1 (х)>