Файл: Григоркина Р.Г. Прикладные методы корреляционного и спектрального анализа крупномасштабных океанологических процессов.pdf

ВУЗ: Не указан

Категория: Не указан

Дисциплина: Не указана

Добавлен: 29.06.2024

Просмотров: 211

Скачиваний: 1

ВНИМАНИЕ! Если данный файл нарушает Ваши авторские права, то обязательно сообщите нам.

состояния, диффузии и теплопроводности.

'Воздействие атмо­

сферных процессов на поверхность океана,

а

т а к ж е массовых

сил различного происхождения (сила тяжести,

приливообразу-

ющие силы и т. д.) является входом этой системы. Выходными

процессами системы являются течения, колебания

уровня моря,

изменение температуры воды и т. д. П а р а м е т р ы ,

определяющие

характер воздействия входных процессов на океан как динами ­ ческую систему, зависят главным образом от его морфометрических особенностей, географического положения бассейна, типа невозмущенной стратификации водных масс и их устойчивости.

Идентификация

гидродинамических

и

термодинамических

моделей океанологических

процессов с

учетом вероятностных

свойств входных и выходных процессов

представляется в а ж н о й

задачей . П е р в ы м этапом в

решении этой задачи д о л ж н о быть

исследование связи

м е ж д у

статистическими

характеристиками

входных и выходных процессов. Некоторые успехи в этом от­ ношении, как отмечалось в предисловии, у ж е имеются при ис­ следовании морского волнения и колебаний уровня моря.

Наиболее простым классом динамических систем являются линейные системы, о б л а д а ю щ и е свойствами аддитивности и од­

нородности

(принцип

суперпозиции).

Если

через

f(t)

и

g(t)

обозначить два любых

внешних

 

воздействия

(входные сигналы),

а через а и

b — некоторые

постоянные

параметры,

то

свойство

аддитивности

 

можно записать

(Бендат,

1965)

в

виде

 

 

 

 

 

H[f(t)+g(t)]=H[f(t)]+H[g(t)],

 

 

 

 

 

 

(5.1)

т. е. если внешнее воздействие

 

f(t)

производит

эффект

 

H[f(t)]r

а воздействие

g(t)—H[g(t)],

то

совместное их

воздействие

даст

эффект H[f(t)-{-g(t)].

Свойство

однородности

 

 

 

 

 

 

 

 

H[af(t)]=aH[f(t)]

 

 

 

 

 

 

(5.2)

означает, что

 

если воздействие

f(t)

производит

эффект

 

H[f(t)],

то при

любом

действительном

числе

а

воздействие

af (t)

произ ­

водит

эффект

aH[f(t)].

 

 

 

 

 

 

 

 

 

 

Б л а г о д а р я

свойствам аддитивности и однородности линейных

динамических

систем,

для

таких

систем могут быть

относитель­

но легко получены аналитические и численные решения диффе ­ ренциальных уравнений, которыми описываются эти системы, что является значительным преимуществом линейных систем с математической точки зрения.

Д и н а м и ч е с к а я

система, результат воздействия на которую не

зависит от начала

этого воздействия, а зависит лишь от интер­

вала времени м е ж д у его началом

и данным моментом, относит­

ся к классу

стационарных. Линейные

стационарные

системы

представляют

наиболее изученный

класс

динамических

систем,

и могут быть полностью описаны

линейными дифференциаль ­

ными уравнениями с постоянными

коэффициентами.

 

7*

 

 

 

99-



Линейные модели динамических процессов в океане часто яв­ л я ю т с я достаточным приближением к реальным условиям и хо­

рошо описывают

широкий класс крупномасштабных явлений.

Именно в

таком

приближении, например, получены основные

результаты

теории

течений, приливов. Многие океанологические

з а д а ч и могут быть сведены к исследованию одномерных стати­

стических связей м е ж д у внешним воздействием

(входом)

и его

результатом (выходом) .

 

 

Рассмотрим вопрос о том, к а к связаны м е ж д у

собой статисти­

ческие характеристики входных и выходных процессов.

Пусть

произвольное внешнее возмущение представляет собой единич­

ную импульсную

функцию

6(т)

(дельта-функция Д и р а к а ) ,

мате­

матически определяемую

как

 

 

 

 

 

 

6 ( т ) = 0

для

хфО,

 

 

е

 

 

 

 

 

 

 

j б ( т ) с ? т = 1

для

всех

е > 0 ,

(5.3)

 

—е

 

 

 

 

 

 

о

 

е

 

 

 

 

 

§&(x)dx=

jb(x)dx=

д л я

всех е > 0 .

 

-s

 

о-

 

 

 

 

 

Дельта - функция с физической точки зрения может быть ин­ терпретирована как мгновенный импульс, приложенный к входу системы. Если на линейную систему с постоянными параметра ­ ми действует не единичный, а произвольный импульс, то реакция системы на этот импульс может быть охарактеризована при по­ мощи весовой функции h(t) (импульсная переходная ф у н к ц и я ) , которая определяется как реакция системы на единичную им­ пульсную функцию по истечении времени т. Таким образом, весо­ вая функция есть реакция системы в момент х на единичную им­ пульсную функцию, приложенную к системе на время т раньше

(Бендат,

1965).

 

 

 

 

Пусть

x(t) есть

сигнал

на входе системы (внешнее

возмуще ­

ние на входе), a y(t)

— на выходе. Через весовую функцию вход­

ной и выходной сигнал

могут быть связаны следующим

образом:

 

 

 

 

со

 

 

 

y(t)

=

\ h{x)x{t-x)dx.

(5.4)

 

 

 

 

о

 

П р и этом

учитывается,

что

поскольку реальная физическая си­

стема не может реагировать в данный момент на будущие воз­

мущения, весовая функция

таких систем д о л ж н а удовлетворять

условию

 

ft(T)=0

при т < 0

100


(условие физической возможности системы) . И з (5.4) следует, что весовая функция полностью определяет реакцию линейной системы на внешние возмущения .

В качестве характеристики линейных динамических

систем

часто

 

используется т а к ж е передаточная функция

Н(Р),

которая

м о ж е т

быть

определена

как

преобразование

Л а п л а с а

весовой

функции

 

 

 

 

со

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

H(P)=\h(x)-er**dx.

 

 

 

 

 

 

 

 

(5.5)

 

 

 

 

 

 

 

о

 

 

 

 

 

 

 

 

 

где р— комплексная частота,

p — b-\-i®.

 

Действительная

часть

комплексной

частоты

б

есть

коэффициент

затухания,

мнимая

часть

 

соответствует частоте

гармонического

сомножителя

зату­

х а ю щ е г о колебания . Д л я

гармонического

колебания

действитель­

н а я

часть

соответствующей

ему

комплексной

частоты

равна

нулю .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

П р и б-*-0, т. е. при

замене

р на

/со, преобразование

Л а п л а с а

совпадает с

преобразованием

Фурье, и

передаточная

функция

м о ж е т

быть записана в виде обращения

Фурье

весовой

функции

 

 

 

 

 

 

со

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Я(со) =

j h ( x ) e - i ° " d x = y ( i a ) .

 

 

 

 

 

(5.6)

 

 

 

 

 

 

о

 

 

 

 

 

 

 

 

 

 

 

В форме (5.6) передаточную функцию называют

т а к ж е

час­

тотной

характеристикой

 

y(ia)

 

системы. Таким

образом,

в

част­

ном

случае

(6 =

0) понятия

передаточной

функции

и частотной

характеристики

совпадают.

Передаточная

функция,

к а к

следует

из (5.5), является более широким понятием, чем частотная ха­ рактеристика, и представляет собой аналитическое продолжение частотной характеристики на всю комплексную плоскость.

•В дальнейшем термины «передаточная

функция»

и «частот­

н а я характеристика» мы будем употреблять

в смысле

(5.6).

Если внешнее возмущение задать в виде гармонической функции

р е а к ц и я системы на возмущение такого вида запишется:

 

 

со

 

 

 

 

u{t)

=

§ h(x)ei<*lt-Vdx=y(ia)eia>t

.

(5.7)

 

 

о

 

 

 

 

Ч а с т о т н а я характеристика

может

быть представлена

в виде

У(«о) =

М с о )

I

•ei^=Pia)-\-iQ(a),

 

где | * / ( с о ) | — м о д у л ь ,

определяющий амплитуду на выходе си­

стемы при частоте

со

(амплитудная частотная

х а р а к т е р и с т и к а ) ,

101


qp (со) — сдвиг

ф а з ы реакции

на выходе системы

по отношению-

к фазе внешнего

возмущения

((разовая частотная

характеристи ­

к а ) . Р(а)

называют вещественной частотной характеристикой,.

Q(co) — м н и м о й

частотной характеристикой

 

 

 

 

Щи)

=

\у(<о)\

созф(со),

(5.9)

 

 

 

Q(co) =

ly(to) [ зт-ф(со),

(5.10)

 

 

 

|y(co)| =

[ ^ J c o ) + Q 3 ( t o ) ] 2 ,

(5.11)

 

 

 

c p ( c o ) = a r c t g - | g | - .

(5.12)

Амплитудная

и

вещественная

частотная характеристики я в л я ­

ются четными функциями со, т. е.

 

 

 

 

 

Ы с о ) | =

[у( - <со)|,

(5.13)

 

 

 

Р ( с о ) = Р ( - с о ) ,

(5.14)

а ф а з о в а я

и

мнимая

частотные

характеристики — нечетными

функциями

со,

 

ф(со)=—<р(—со),

(5.15)

 

 

 

 

 

 

Q ( c o ) = - Q ( - c o ) .

(5.16)

Весовую и

передаточную

функции,

а т а к ж е частотную х а р а к т е ­

ристику называют динамическими характеристиками линейной

стационарной системы, с помощью которых эта

система м о ж е т

быть полностью описана.

 

Весьма в а ж н ы м является то обстоятельство,

что динамиче ­

ские характеристики одномерной линейной системы вполне опре­ деленным образом связаны сравнительно простыми соотношени­ ями со статистическими характеристиками процессов на входе-

выходе

системы

 

(с дисперсиями,

функциями

корреляции и

спектральной

плотности) .

 

 

 

 

 

Если известны

динамические характеристики

системы, а т а к ­

ж е

статистические

характеристики

процессов

на

входе

системы,

то

могут

быть

найдены характеристики на

выходе

системы,

С другой стороны, по известным статистическим

характеристи ­

кам

входного

и

выходного процессов определяются

динамиче ­

ские характеристики системы. Динамические характеристики си­ стемы могут быть найдены либо теоретически по уравнениям, описывающим процесс, либо эмпирически по статистическим ха­ рактеристикам входного и выходного процессов. Сравнение ре­ зультатов этих двух подходов может быть основой для опреде­ ления неизвестных параметров системы.

Удачная попытка связать трехмерные спектры скорости и ка ­ сательного н а п р я ж е н и я ветра и описать трансформацию энергии

102