книги из ГПНТБ / Исакович, М. А. Общая акустика учеб. пособие
.pdf§ 23. Разложение Фурье волны с произвольной зависимостью от времени
Покажем, что при соблюдении известных условий, налагае мых на временную зависимость волны р (t, г), которые будем считать выполненными, можно представить волну в виду супер позиций гармонических волн различных частот путем разложе ния по Фурье функции р. Эти условия таковы: если функция
периодична по времени, то она разлагается в ряд Фурье; |
если |
|||
функция не |
периодична, но достаточно |
быстро убывает |
при |
|
t —>—оо и t |
—>+ оо (например, является |
ограниченным по |
|
вре |
мени импульсом), то она разлагается в интеграл Фурье. Если спа дание на бесконечности недостаточно быстрое, то разложение в интеграл Фурье неосуществимо. Будем пользоваться комплекс ным представлением волн.
Периодическая функция |
с |
периодом Т = |
2я/ш0 |
разлагается |
|
в ряд |
|
|
|
|
|
P(t, г) = Ё Рп(г)е~1па°*, |
|
где р„ (г) = |
_£»_ Гр п |
г)е*пш'* dt. |
|
П= —со |
|
|
Z n |
JО х |
1 |
Коэффициенты рп (г), меняющиеся от точки |
к точке, — ампли |
||||
туды спектра волны в каждой точке. |
в интеграл |
||||
Непериодическая функция |
разлагается |
||||
|
СО |
|
|
|
|
p(t, Г) = |
J |
рш(г) е~ш |
dm, |
|
|
где
СО
Ав(г) = -2 )Г J р Ѵ’ Г) e'ai(it
Коэффициенты рш (г) называют спектральной плотностью ампли туды разложения. Элементы интеграла в области от —оо до нуля —
волны с отрицательными частотами. В этой области ра = р1ш. Покажем, что каждое из гармонических слагаемых, т. е. член ряда р„ (г) или элемент интеграла da рш(г) е~ш , яв ляется волной, способной распространяться в данной среде. Математически это значит, что гармонические слагаемые должны
каждое в отдельности удовлетворять уравнению (2 2 .2 ).
Волна р (t, г) удовлетворяет по условию волновому уравнению
Ар |
1 |
д*р |
= 0 . |
с» діг |
|||
Гармоническую компоненту ряда или интеграла можно записать,
70
опуская постоянный множитель, в виде
|
ь |
Л |
Рсо = } реш dt, |
|
а |
считая, что для |
периодической функции b — а = Т, а частота |
для компоненты номера п есть /гео 0; для непериодической функции
интеграл берется в бесконечных |
пределах (а = —оо, Ь = +оо), |
а частота принимает все значения от —оо до +оо. |
|
Умножим волновое уравнение |
на еш и проинтегрируем по |
времени в пределах от а до Ь. В первом члене, меняя порядок
интегрирования и |
дифференцирования, найдем |
ь |
ь |
IАр еш dt = А J р еш dt = Арш.
Во втором члене дважды произведем интегрирование по частям:
|
о |
|
|
|
|
|
|
С&Р |
dt = - ^ е ш — г© j |
|
dt = |
|
|
|
|
J дР |
|
|
|
|
|
|
|
|
ь |
dp |
|
ш |
и |
ш |
|
|
ре |
со ре |
dt. |
||||
|
а |
dt |
|
|
|||
|
|
|
|
|
|
Но для периодической функции значения функции и ее производ ных на концах интервала длиной в один период равны между собой, поэтому первые два члена исчезают. Для непериодической функции эти члены исчезают потому, что по условию сама функ ция р исчезает на бесконечности. Следовательно,
ь |
ь |
j | f e^ |
= _C 0 2 | р Л г = : - С 0 2 рШ. |
Отсюда получим, подставляя в проинтегрированное уравнение, Арео + ^Рш = О,
т. е. компоненты Фурье действительно удовлетворяют уравнению Гельмгольца, а значит, члены разложения данной волны действи тельно являются волнами, каждая из которых может распростра няться независимо от других.
Доказанная теорема имеет важнейшее значение: она придает физический смысл разложению Фурье по времени. Эта математи ческая операция имеет смысл замены волны с произвольной вре менной зависимостью суперпозицией волн со стандартной зави симостью от времени — гармонической зависимостью.
71
Разложение Фурье имеет физический смысл не только для волн, удовлетворяющих волновому уравнению, но, как можно показать (см. § 26), и для волн в более сложных средах. Необхо димо только, чтобы уравнение, которому удовлетворяет давление (или какая-либо иная характеристика волны), было линейным и однородным.
Помимо рассмотренных типов волн, возможны еще волны, которые неразложимы ни в ряд, ни в интеграл Фурье, но все же могут быть представлены в виде суперпозиции некоторого дискрет ного набора гармонических волн: это суммы гармонических волн несоизмеримых частот. Такие волны называют почта периоди ческими, потому что, как можно показать, любой отрезок такой волны повторяется со сколь угодно большой точностью через достаточно большой промежуток времени. Простейший пример почти периодической волны — биения между двумя волнами близких, но несоизмеримых частот. Аналитически почти перио дическую волну можно записать в виде
Pit, r )= È Рп(г)е~
П
где частоты соЛ несоизмеримы.
§ 24. Спектральные разложения волн
Разложение волн с любой зависимостью от времени на гармо нические волны разных частот — это пример так называемого спектрального разложения: представления данной функции в виде линейной суперпозиции (ряда или интеграла) стандартного на бора функций с более простыми свойствами. Если эти вспомога тельные функции изучены, то исследование других функций сводится к определению коэффициентов в спектральном разло-. жении. В акустике (и в других волновых науках) в качестве такого стандартного набора удобно пользоваться гармоническими функ циями времени, представляя заданную волну в виде интерферен ционной картины гармонических волн разных частот. Спектраль ный подход освобождает нас от необходимости исследовать каж дую волну со своей зависимостью от времени в отдельности: каждая звуковая волна оказывается представленной в виде су перпозиции гармонических функций, и рассмотрение временной зависимости оказывается упрощенным до предела.
Но поле гармонической волны зависит вообще от трех коорди нат-, и при одной и той же частоте зависимость от координат может быть самой разной. Возникает вопрос о возможности даль нейшего упрощения изучения волн: возможности представления произвольных гармонических по времени функций от координат также в виде суперпозиции некоторого набора гармонических волн (конечно, той же частоты), стандартно зависящих от коорди нат, — вопрос о пространственном спектре гармонической волны.
Ответ на этот вопрос зависит от акустической ситуации. Если известно поле данной гармонической волны на плоскости, то в качестве стандартного набора можно взять плоские гармони ческие волны (мы увидим в § 33, что это возможно только при некотором обобщении понятия плоских волн), если известно поле на сфере, то удобно производить разложение в спектр по набору так называемых сферических волн, и т. п. В данной главе рассмот рим разложение поля, по плоским волнам. Для этого теперь изучим подробно плоские волны и их обобщения.
§ 25. Плоские гармонические волны
Произвольную плоскую волну можно разложить в спектр, т. е. ее можно представить в виде суперпозиции плоских же гармонических волн (оговорки — те же, что и выше для спек трального разложения любой волны). Напишем в комплексной форме бегущую плоскую гармоническую волну. Как указано в § 17, выражение для плоской волны (в векторной записи) полу чается из выражения для временной зависимости в точке путем
замены времени t |
на бином t— Sr, где 5 — вектор медленности |
||
волны. |
В гармонической волне временная зависимость дается |
||
множителем е~ш . |
Значит, |
бегущую плоскую гармоническую |
|
волну |
можно записать в виде |
||
|
|
р = р о ехр |
[—іа (t— Sr) ], |
где ро — постоянная (вообще— комплексная). Эту постоянную будем называть комплексной амплитудой плоской волны. Для гармонических волн удобно пользоваться вместо вектора медлен ности S пропорциональным ему волновым вектором к = aS. Плоская гармоническая волна записывается тогда в виде
р = Ро ехр (— iat 4- ikr), |
(25.1) |
что можно считать комплексной записью |
формулы (18.4) при |
Ро = |Р о Іе‘е- Фронты волны совпадают с |
плоскостями kr = |
= const, перпендикулярными к к. Комплексная амплитуда коле бания в какой-либо точке есть р 0 ехр (ikr) и равна амплитуде волны р о (общей для всех точек), умноженной на фазовый множи тель ехр (ikr). Если в той или иной задаче амплитуда волны не существенна (например, при нахождении коэффициента отражения волны от препятствия), то амплитуду, как и временной множитель, опускают:
р = ехр (ikr). |
(25.2) |
Для волны, бегущей, например, вдоль оси х вправо |
или влево, |
р = ехр (±ikx). |
(25.3) |
Для того чтобы вернуться к вещественной записи плоской волны, необходимо предварительно восстановить оба комплекс-
73
ных множителя, и р 0 и е~ш , |
и лишь тогда брать вещественную |
часть от получившегося выражения: |
|
р = Re [ра>ехр (— Ш + |
ikr)] = I Po I cos (соt — kr — е), (25.4) |
В плоской гармонической волне зависимость от координат в дан ный момент времени также является синусоидальной, как и вре менная зависимость в каждой данной точке. Фаза е комплексной амплитуды волны окажется существенной, только если придется иметь дело одновременно с несколькими волнами. Имея дело только с одной волной, всегда можно выбрать начало отсчета времени, например, так, чіюбы амплитуда волны была веще ственна (е = 0 ).
Скорость частиц |
в волне р |
= |
exp {ikr), согласно |
(22.5), вы |
ражается формулой |
|
|
|
|
ѵ = - ^ ех^ |
^ |
= - т ^ - |
<25-5) |
|
Для плоской волны можно написать разложение Фурье, при нимая за аргумент вместо времени линейную комбинацию t — Sr. Тогда спектр разложения не будет зависеть от координат. Разло жение для волны, бегущей вдоль оси х, имеет вид
( |
/ ! = - { - СО |
|
|
|
£ рпехр (— іпщі + |
ink0x) |
|||
t ---- -) = |
||||
для периодической волны с периодом Т = 2я/со0 |
(при k 0 = а»„/с) |
|||
и вид |
|
|
|
|
р(і---- J-) = |
1 А»ехР |
Ш + ikx)d& |
||
|
— CD |
|
|
|
(при k = м/с) для непериодической волны, разложимой в интеграл Фурье по времени. Отсюда видно, что волна в каждый момент времени оказывается разложенной в спектр по координатам, т. е. представлена в виде суперпозиции синусоидальных простран ственных распределений. Амплитуды спектров разложения р„ (или Po) не зависят ни от координат, ни от времени.
Кроме плоских волн, разложимых в ряд или в интеграл Фурье, возможны еще волны, хотя в таком виде и не представимые, но которые все же можно выразить в виде суперпозиции гармониче- ' ских плоских волн. Такими волнами будут, в частности, волны вида
|
Р ( t — -7 |
-) = S Рпехр (— m nt + iknx), |
|
где |
р„ — постоянные, |
величины м„ несоизмеримы, |
а м„/£„ = с |
для |
всех и. |
|
|
|
Переходя к декартовой системе координат, в которой направ |
||
ляющие косинусы волнового вектора равны cos а, |
cos ß, cos у, |
||
74
получим из (25.2) координатное представление плоской гармони ческой волны:
р = exp (ikcosa-x -(- cosß- г/ -f- ikcosy-z) = |
|
|
= exp (ikxx-\- ikyy |
ikzz), |
(25.6) |
где проекции волнового вектора на оси координат |
kx = |
k cos а, |
ky — k cos ß, kz = k cos у равны волновым числам следов волны на координатных осях. Так как эти проекции меньше модуля волнового вектора, то медленность следов меньше, чем медлен ность волны, а скорость следов больше скорости волны.
На координатных плоскостях следы волны представляют собой
двухмерные |
гармонические волны. |
Например, на плоскости |
z = 0 бежит |
волна р = exp (ikxx + |
ikyy) с волновым числом |
kx -(- kl — ^ sin у, равным проекции волнового вектора k на пло
скость 2 = 0. Скорость этого следа также больше скорости волны. Из уравнения (22.2) следует
&* + % + £ = fe2 = co2/c2- |
(25.7) |
Каждая комбинация трех вещественных чисел kx, ky, kz, удовлетворяющих уравнению (25.7), соответствует плоской гар монической волне данной частоты, бегущей в направлении, опре деляемом направляющими косинусами
cos а = kjk, cos ß = kylk, cos у = kjk.
Часто располагают какую-либо координатную плоскость (на пример, плоскость xz) параллельно волновому вектору данной
плоской |
волны. Тогда движение |
не зависит |
от третьей коорди |
|||
наты |
(у) |
и волну |
можно записать в виде р |
= exp (ikxx + |
ikzz), |
|
где |
kx = |
k cos Ѳ, |
kz = k sin Ѳ. |
Угол Ѳ между волновым |
век |
|
тором и осью X называют углом скольжения данной волны |
отно |
|||||
сительно оси X или относительно плоскости ху.
Комплексная форма записи удобна не только для звуковых, но и для любых гармонических волн. Так, температурную волну
(19.2) можно |
записать, опуская временной множитель, |
в виде |
|
Т = Т0ехр(Щхх — \ %х). |
(25.8)- |
Вязкая волна |
(19.4) запишется в виде |
|
|
V = v0 exp (i^yX — %ух). |
(25.9) |
§ 26. Сохранение формы бегущих гармонических плоских волщ Дисперсионное уравнение
В средах, подчиняющихся волновому уравнению, плоская волна любой формы распространяется без искажения. В других средах этим свойством обладают только гармонические плоские волны. Единственное условие, налагаемое, при этом на среду, —
75
это ее линейность: в среде должен быть справедлив принцип су перпозиции, и, кроме того, если в среде может распространяться какая-либо волна, то волна, отличающаяся только множителем при давлении,- скорости частиц и т. п., также может распростра няться в данной среде. Уравнение, которому подчиняется, напри мер, давление р в такой среде, можно записать в виде
5 ’ (Р )= 0, |
(26.1) |
где S ’ — линейный однородный оператор |
(например, для сред, |
подчиняющихся волновому уравнению, |
і? = А---- |
Докажем свойство сохранения формы гармоническими волнами произвольной частоты или произвольной длины волны для такой среды. Для этого нужно доказать, что в числе решений уравне
ния (26.1) есть волны вида |
|
р = ехр(— Ш -\-ікг) |
(26.2) |
для любой длины волны (и соответственной частоты) или для любой частоты (и соответственной длины волны).
Пусть оператор 9? переводит некоторую функцию р в некото рую функцию q:
£ (Р) = q. |
(26.3) |
Если q = 0, то соответственная функция р представляет собой свободную волну в среде. Возьмем р в виде (26.2) и, подставив в (26.3), продифференцируем обе части этого равенства по времени.
Всилу линейности и однородности оператора дифференцирование
идействие оператора можно переставлять между собой, так что
= |
д%\р) |
= |
tCöp)= |
— |
— |
Таким |
образом, |
dq/dt = —шр, |
откуда |
следует |
|
|
|
q = q e |
\ |
|
|
где q от времени уже не зависит, но вообще зависит от со. Под ставляя это выражение в (26.3), получим уравнение, в котором временные зависимости отсутствуют и содержащее частоту какпараметр. Например, для волнового уравнения получится урав
нение Гельмгольца относительно амплитуд колебаний q в разных точках. Далее, дифференцируя (26.3) по координатам, найдем снова в силу тех же свойств оператора
Ѵр = 99£ (р) = 3? (Ѵр) = 3 (ikp) = ikSS (р) = ikq,
т. е. Ѵр = ikq. Отсюда следует, что q можно представить в виде
q = f (о, k) exp (— m t + ikr),
где / может зависеть, помимо со и k, только от коэффициентов опе ратора.
76
При произвольных а» и k функция ехр (— m t + ikr), конечно, не представит волну, которая могла бы распространяться в дан ной среде, так как вообще эта функция не является решением уравнения (26.1). Но если выбрать со и А; так^чтобы выполнилось условие
/((о, Ä) = 0,_ |
(26.4) |
то будет и q = 0 , и для такой комбинации со и к уравнение (26.1) будет удовлетворено: соответственная монохроматическая волна в среде сможет распространяться как свободная. Вообще каждому значению со будет соответствовать решение уравнения (26.4) относительно к и каждому значению k будет соответствовать решение уравнения (26.4) относительно частоты со. Для изотроп ной среды это уравнение может содержать только модуль волно вого вектора, так что уравнение можно привести к виду
’ F (со, k 2) = 0. |
(26.5) |
Уравнение (26.5) называют дисперсионным уравнением для данной среды. Например, дисперсионное уравнение, соответст вующее волновому уравнению (16.1), есть
k 2 — со Ус2 = 0 ,
где с — постоянная. Дисперсия в средах, описываемых волновым уравнением, таким образом, отсутствует: фазовая скорость гармо нической волны любой длины есть с.
Пример дисперсионного распространения дают изгибные волны на стержне. Как известно из курса сопротивления материалов, уравнение для поперечного смещения £ стержня при малых коле баниях можно записать в виде
д% , _P _É^ _n
дх* ~r G dt2 ~ ’
где G — коэффициент изгибной жесткости. Отсюда, подставляя решение вида ехр (— m t + ikx), найдем дисперсионное уравне ние изгибных волн в виде.
Ä*---- g-- со2 = 0 |
(26.6) |
(в § 7 мы применили другой способ нахождения этого дисперси онного уравнения). Зависимость фазовой скорости от частоты или волнового числа имеет вид
4 |
|
с = а/k = kV~Gfp = і/ю К С /р. |
(26.7) |
Еще пример — морские волны на поверхности воды («грави тационные волны»). Как известно из гидродинамики, на поверх ности несжимаемой жидкости, находящейся в поле силы тяжести,
77
могут распространяться поверхностные волны. Потенциал ско рости таких волн удовлетворяет уравнениям
Ä<p=0. 3 + T - S * “ 0-
Первое из этих уравнений — условие несжимаемости среды, вто рое — уравнение движения поверхностной волны. Поверхност ную волну можно найти в виде ф = ехр (— m t + ikx -f- kz). Подставляя это выражение, автоматически удовлетворяющее тре- ' бованию несжимаемости среды, в уравнение движения, найдем
дисперсионное уравнение в виде
k ---- 1 - со2 = 0 . |
(26.8) |
Скорость поверхностных волн зависит от частоты или волнового
числа по закону |
|
с = g/<o = V gW , |
(26.,9) |
что, в отличие от изгибных волн на стержне, соответствует нор мальной дисперсии.
Спектральный подход к решению задач акустики требует нахождения всех монохроматических волн, способных распро страняться в данной среде. В принципе это можно сделать, решая дисперсионное уравнение относительно со или относительно k (фактически такое решение может оказаться очень трудным).
Если решение получено, т. е. известно <и = |
со (k) или k |
= k |
(со), |
|||
то |
фазовая скорость получается в виде |
с = соIk (со) |
или |
с = |
||
= |
со (k )/k — как функция либо |
частоты, либо |
волнового числа. |
|||
Вообще для любого линейного |
уравнения |
для |
волн 3? (р) — О |
|||
всегда имеет место дисперсия: случай волнового уравнения, когда дисперсия отсутствует, — исключение, но исключение очень важ ное: это, как мы видели, уравнение волн в свободной неограни ченной среде.
В акустике встречаются два принципиально различных типа дисперсии. Один тип обусловлен физическими свойствами среды: зависимостью упругих напряжений не только от деформаций, но и от скорости изменения деформации. В плоской звуковцй волне в неограниченной среде возможен только этот тип диспер- •сии. Он всегда сопровождается поглощением звуковой энергии. Классические примеры таких сред-— лед, вар. При малой ско рости деформирования возникающие упругие силы малы, и за достаточно долгое время эти тела могут растекаться подобно жид костям под действием собственного веса. Но при резком ударе возникающие силы — такие же, как в обычных твердых телах: кусок льда или вара разбивается при таком ударе на осколки. Поэтому в таких телах при разной частоте колебаний скорость волн различна: с ростом частоты всегда растут'и упругие силы,
78
скорость звука увеличивается (аномальная дисперсия). Опыт показал, что такой дисперсией обладает также ряд жидкостей и многоатомные газы. «Область дисперсии», т. е. частотный диапа зон, в котором зависимость упругости от частоты заметна, ме няется в различных веществах от долей герца до тысяч мегагерц. Более подробно этот тип дисперсии рассмотрен в гл. XII.
Другой тип дисперсии обусловлен границами среды, в которой распространяется волна, и не зависит от свойств среды. Этот тип с поглощением звука не связан и целиком определяется кинема тикой волнового движения в ограниченной среде. Такова, напри мер, рассчитанная выше дисперсия скорости изгибных волн в стержне. Физическая картина дисперсии для изгибных волн заключается в том, что коэффициент упругости стержня растет при уменьшении длины изгибаемого участка; поэтому с умень шением длины волны, т. е. с увеличением частоты, скорость волн
растет. Дисперсия наблюдается и при |
распространении волн |
в жидких средах, заключенных в трубах, |
и т. д. Более подробно |
эти вопросы рассмотрены в гл. VIII. |
|
§ 27. Групповая скорость. Распространение узкополосного сигнала
Монохроматическая волна не может передать никакой инфор мации, никакого сигнала: в такой волне в каждой точке проис ходили, происходят и всегда будут неизменно происходить гармо нические колебания. Чтобы передать сигнал при помощи волны, необходимо, чтобы в ней что-нибудь менялось, чтобы волна была модулирована тем или иным способом, например, чтобы она дли лась ограниченный промежуток времени. Это уже не будет моно хроматическая волна; такой сигнал можно рассматривать как интерференционную картину, образованную суперпозицией гар монических волн разных частот. Информацию передаст именно эта интерференционная картина.
Но в диспергирующей среде сама интерференционная картина меняется, так как компоненты разных длин волн распростра няются с разной скоростью. Таким образом, в диспергирующей среде передаваемая информация оказывается искаженной.
Какова глубина этого искажения и в какой мере все-таки можно передавать сигнал в диспергирующей среде, можно выяс нить при помощи фурье-представления волны.
Выясним раньше всего, как найти изменение данного профиля волны при ее распространении в среде с заданным законом диспер сии. Для этого достаточно выполнить следующие действия: раз ложим по Фурье данный профиль на сумму синусоид различных длин волн и припишем каждой синусоиде временной множитель соответственно дисперсионному уравнению среды, как сказано в предыдущем параграфе. Каждая из полученных таким образом компонент — свободная гармоническая волна, фазовая скорость
79
