Файл: Григоркина Р.Г. Прикладные методы корреляционного и спектрального анализа крупномасштабных океанологических процессов.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