Добавил:
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

книги из ГПНТБ / Доценко С.В. Теоретические основы измерения физических полей океана

.pdf
Скачиваний:
7
Добавлен:
23.10.2023
Размер:
5.68 Mб
Скачать

только от свойств датчика и показывает, какую долю в общем зна ­ чении выходной величины Y занимает измеряемая величина X (рь р 2 , Рз; т ) , находившаяся в момент т в данном элементе объема

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

* ( P i ,

Ра, Рз;

^) =

5 ( Р 1 ) Ч Р а ) Ч Р з ) Ч - 0 .

 

то

 

 

 

 

 

 

 

 

Г = / / 0 ( г ь

го, г3 ;

t),

 

 

 

т. е. аппаратная функция есть реакция датчика

на

четырехмерный

импульс.

 

 

 

 

 

 

 

Несмотря на то что вывод соотношения

(1.2)

проведен

не строго,

точные в ы р а ж е н и я преобразования, рассчитанные

для

различных

линейных датчиков,

имеют

именно такой

вид,

а

достигнутая при

этом наглядность позволяет быстрее прийти к цели в более слож ­

ных случаях

(гл. V I I ) .

 

 

 

 

 

 

Отметим,

что аналогично может быть введено понятие

аппарат ­

ной функции

всего прибора в целом,

учитывающее

свойства не

только

датчика, но

и инерционность

остальной

части

прибора.

В дальнейшем, если

не

оговорено противное, под Н0 понимается

именно т а к а я

функция.

 

 

 

 

 

 

З а п и ш е м

в ы р а ж е н и е

(1.8) в более компактном

виде. Обозначим

радиус-вектор центра

датчика

 

 

 

 

 

 

 

 

 

r = r 1 i 0 - f - r 2 j 0 + r 3 k o ,

 

 

 

где i 0 l jo, k 0 единичные орты, г радиус-вектор точки

измерения

 

 

 

 

Р = Р 1 io -f- Р2 Jo+Рз^о,

 

 

 

элемент

объема

 

 

 

 

 

 

 

 

 

 

 

dp=dpi

dp2

dp3.

 

 

 

Тогда

 

 

 

t

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

К = j dp j X (P ;

*) tf0

(r - Pi t ~ -0 d-s

 

(1 -9)

 

 

 

— 00

 

 

 

 

 

где интегрирование по пространству распространяется на всю об­ ласть совместного существования функций Х(р; т) и Но —р; t — т ) .

Измерители и регистраторы, входящие в измерительные при­ боры, достаточно полно могут быть охарактеризованы коэффициен­ тами передачи и частотными характеристиками . Эти параметры, несомненно, д о л ж н ы быть согласованы с параметрами датчика и выбраны таким образом, чтобы обеспечить наилучшие метрологи­ ческие показатели всего измерительного прибора в целом. По -

20

скольку эти узлы прибора представляют собой электронные устрой­

ства,

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

очень

гибки

и могут изменяться

в

широ­

ких пределах в отличие от датчиков,

имеющих, как правило,

жест­

кие

параметры,

определяемые

конструктивными

требованиями.

Вследствие

этого

часто

оказывается,

что

самым

неуправляемым

элементом прибора является датчик,

а измеритель

и

регистратор

д о л ж н ы рассчитываться на основании его параметров .

 

 

 

Б о л ь ш а я

часть

авторов, занимающихся

гидрофизическими ис­

следованиями, не принимает во внимание свойство датчиков

осред-

нять

измеряемую

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

молчаливо

предполагая, что датчики являются точечными. П р е д е л ы

примени­

мости этого

предположения, изложенные выше, требуют детального

анализа работы конкретных датчиков. В некоторых работах (напри­ мер, в [41]) говорится о возможности исследования неточечных дат ­ чиков, однако само это исследование не проводится.

§4. Математическое описание аппаратных функций

имногомерных полей

При анализе функций одной переменной (например, времени) оказывается полезным спектральное р а з л о ж е н и е этой функции (ряд

или

интеграл Фурье) .

П о л я

измеряемой

физической

величины

Х(р;

т)

и аппаратные функции

 

# о ( р ; т) в

общем

случае

зависят

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

Согласно

теории

кратных

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

Фурье,

функция

четырех переменных

(трех

координат

и времени)

f (pi, рг, рз', х), удовлетворяющая усло­

виям

Д и р и х л е по всем переменным,

может

быть

представлена

в виде

[59, 74]

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

со

 

 

 

 

 

 

 

 

/ ( P i .

Р2> Рз!

*)=~(ЩГ

Щ

d

a i d

a 2 d

a i X

 

 

 

 

 

 

 

 

 

 

 

 

— СО

 

 

 

 

 

 

 

 

 

со

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

X

J

/ ( « . ,

« 2 ,

« 3

;

 

 

 

w)eHa'9'+°*H+a3fMda,

 

 

 

 

— со

 

 

 

 

 

 

 

 

 

 

 

 

где

четырехмерный спектр

этой

функции

дается

выражением

 

 

 

 

 

 

со

 

 

 

 

со

 

 

 

 

 

 

 

7(«ь а 2 , а 3 ;

)=I

I J

d P i dpodp3

j

/ ( Р ь Р 2 , р 3 ;

"ОХ

 

ООСО

e~j (a,p,+a-iP:+a3P3 + »-c)

Здесь величина

со есть круговая

частота,

а величины « ь аг, а з —

волновые числа,

связанные

соответственно

с координатами pi, рг

и р3 . Обозначив:

 

 

. в

 

 

do. = dax

doL2

daz,

 

21

м о ж но записать эту пару

преобразовании в векторной форме:

/ ( р ;

-•) =

, r , l v l

( / ( « ;

(о) exp[y(ap-4-tot)]afarfi.

 

 

 

(2r0-i

 

 

 

 

 

 

(1.10)

 

 

 

 

 

 

 

 

 

/ ( « ;

< • > ) = / ( р ; х ) е х р

[ — У (ap +

шт)] ofp fif-c.

 

Предположим,

что функция

f (pi, р2 , рз; т)

не зависит

от р3 , но

записана в четырехмерной форме. Тогда

 

 

 

 

 

7 ( а ъ а 2 . а з ; w ) = J / ( P i - Р_>; Т ) Х

 

 

 

 

 

 

 

 

со

 

 

 

 

X в " у

rfPl

rfP2

 

J <T; 'a *' tfp3

 

 

 

 

 

 

 

— 00

 

 

 

 

 

 

 

 

00

 

 

 

 

и, учитывая

известное соотношение

Je- Ja P^p = 2n6 (a), где б (a) —

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

найдем

—со

 

 

 

 

 

 

 

 

 

 

/ ( « ) ,

an, a 3 ; ш ) = / ( а ь

а 2 ;

w )

(<?-3).

 

Это соотношение дает

возможность совмещать

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

Фурь е функций разного числа измерений.

 

 

 

Отметим, что аппаратные функции приборов всегда

удовлетво­

ряют условиям Д и р и х л е

(в силу ограниченности объема

измерения

и инерционности),

и для них всегда

справедливо

представление

(1.10). Этим ж е условиям удовлетворяют

и многие

детерминирован­

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

Рассмотрим однородные случайные «замороженные» поля, т. е. такие не зависящие от времени поля Х(р), для которых плотность

вероятности р ( | )

{N)(X\,...,

XN) не меняется при любых парал ­

лельных переносах точек р ( 1 ) , ... , p(jV>

[44, 51]. В силу однородности

поля

Х(р)

его среднее значение

Х(р)

д о л ж н о быть

постоянным.

При

этом

без ограничения

общности

можно

считать, что Л"(р) = 0 ,

заменяя,

если надо,

исходное

поле

Х(?)

полем

Х'(р)=Х(р)

Однородное случайное поле может быть представлено, в отличие •от (1.10), в виде спектрального разложения (стохастического инте­ грала Фурье — Стильтьеса)

X ( p ) = j e x p ( ; a p ) r f Z ( « ) ,

(1.11)

где интегрирование распространяется на все пространство волно­ вых векторов, а функция Z(a) является комплексной случайной

'22

функцией. Если поле Х(р) вещественно, а пространственная корре ­ ляционная функция

 

Я ( г ) = * ( Р ) * ( Р + г )

 

удовлетворяет условию

J" | В (г) | d r < o o ,

то случайные

амплитуды

разложения (1.11) обладают следующими свойствами:

 

 

d Z ( - a ) =

rfZ*(a),

(1.12)

 

d Z ( a ) = 0 ,

 

 

dZ* (a) dZ ( о , ) = 8 (a - 04) (7 (a) da da,,

(1.13)

где 8 (a) трехмерная

дельта - функция .

В х о д я щ а я в

в ы р а ж е н и е

(1.13) неслучайная функция G(a) носит название трехмерной спект­

ральной

плотности (или

трехмерного

спектра)

поля Х(р).

Эта

функция д о л ж н а

быть всюду неотрицательна.

 

 

 

 

Из перечисленных соотношений вытекает, что корреляционная

функция

поля связана

с трехмерным

спектром

выражением

 

 

 

В (г) =

j exp (Jar) G (a) da,

 

 

(1.14)

и, в свою очередь,

 

 

 

 

 

 

 

 

 

 

О ( а ) = - щ г J exp (-jar)

В (г) dr.

 

 

 

Таким образом, задание трехмерного спектра G(a)

равносильно

заданию корреляционной функции В(г).

 

 

 

 

 

Наиболее простым

типом однородных

полей

являются однород­

ные изотропные

поля

Х(р),

дл я которых

плотности

вероятности

/7? ( 1 )

р(лг)№> • • •' ^n)

не меняются

при любых

параллельных пе­

реносах,

вращениях

и зеркальных

отражениях

системы

точек

р( 1 ) , ... , p(JV>. Такие

поля дл я краткости

называются

просто изотроп­

ными. И х корреляционные функции зависят только

от одного

пере­

менного г = | г |

 

 

 

 

 

 

 

 

 

Д ( г ) = * ( р ) * ( р + г ) .

При этом трехмерный спектр поля G (а) зависит лишь от а= | а | . Рассмотрим значения поля Х(р) в точках какой-либо прямой, например, в точках прямой рг = Рз = 0. Они представляют собой од­ нородное поле на этой прямой (т. е. стационарный процесс от пере ­

менного pi) . Одномерное преобразование Фурье функции В (г)

со

Oi Ы=-ъг

j

e-*rB(r)dr

 

 

—со

 

 

должно быть всюду неотрицательным и носит название

одномерной

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

одномерного спектра)

поля Х{р).

 

 

 

О т с ю да корреляционная функция поля

со

 

5 ( r ) = j е>""0, («,) da,

(1.14а)

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

со

В (г) =

j j J е > <«•'•+•«">> С? (l/"a?+al+as ) da, da2 da3 .

 

 

 

— CO

 

 

 

 

 

 

 

 

 

 

Корреляционная функция поля на рассматриваемой прямой да­

ется частным случаем этой

формулы

при /*2 = гз = 0, т. ё.

 

 

 

 

со

 

 

 

 

 

 

 

 

 

£

(r) = j

J j eJa'r

 

G ( j / a l + a g + a j j ) da, da2

da3 .

 

 

 

 

CO

 

 

 

 

 

 

 

 

 

Д л я упрощения

этого

в ы р а ж е н и я

перейдем

в

плоскости

(а», а3)

к полярной системе

координат

 

 

 

 

 

 

 

 

 

 

 

a 2 = 7 . cos

ср,'

 

 

 

 

 

 

 

 

 

a3 =v. sin «р.

 

 

 

 

 

Учитывая, что при этом

элемент

площади

есть

xdxdcp,

и

инте­

грируя по ф, получим

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

da,

 

 

Производя

во внутреннем

интеграле замену

переменной

а =

= Уа^ + и 2 , окончательно

найдем

 

 

 

 

 

 

 

 

 

со Г

со

 

 

 

 

 

 

 

5 ( r ) = 2 i c

j

e h

' r \ О (a) a. da.

da,.

 

 

 

Сравнение полученного в ы р а ж е н и я с формулой (1.14а) дает иско­ мое в ы р а ж е н и е одномерного спектра трехмерного изотропного поля через трехмерный спектр

G, ( а , ) = 2 т с J (у (a) a da.

Дифференциру я эту формулу по ai, получим

в ы р а ж е н и е трех­

мерного спектра через

одномерный

 

0

(

а ) =

*

(1.15)

 

4

'

2яа da

 

24

Обобщением однородных случайных полей являются локальн о

однородные случайные поля

Х(р),

все распределения

вероятностей

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

пар то­

чек не меняются при любом

параллельном

переносе всех

рассмат ­

риваемых точек. Среднее значение

приращений

 

 

 

 

такого

поля является

линейной

функцией

вектора

г

 

 

 

 

 

 

т (r) =

[ ^ ( p + r ) - ^ ( p ) j = c , r >

 

 

 

 

где ci постоянный

вектор,

а общий

второй

момент

приращений

АхХ(д)

поля Х(?)

 

 

 

 

 

 

 

 

 

 

 

 

 

£ > ( Р 2 - Р , , г ь

г2 ) = [X(р,

+

г,) - X ( Р 1 ) |

(?о +

Го)

- X ( Р 2

) ]

в ы р а ж а е т с я через структурную функцию этого поля

 

 

 

 

 

 

D(r)

=

 

 

 

[X(P+r)-X(P)}*

 

 

 

 

при помощи равенства

 

 

 

 

 

 

 

 

 

 

 

D(p,

г,,

ro) = 4 - [ D (P - r . ) +

^ ( P + r 2

) - D ( p - r , + r o ) - £ > ( p ) ] .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(1.16)

Структурная

функция

однозначно

определяется

своей

спект­

ральной плотностью

(трехмерным

спектром)

G(a)

 

 

 

 

 

 

 

D(r)=2

j ( l - c o s a r ) G(a)rfa .

 

 

(1.17)

Локальн о однородное

поле,

зависящее

только

от г = | г | ,

назы ­

вается локально

изотропным.

поле У (t)

 

 

 

 

 

 

Если

локально однородное

одномерно,

то оно

носит

название

случайного

процесса

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

приращениями .

Среднее значение его приращений

 

 

 

 

 

 

 

 

 

 

туЦ{)=^7Щ1^¥Ш=а^,

 

 

 

 

 

 

(1.18)

где di — константа, а структурная

функция

 

 

 

 

 

 

£Ж ) = [ П Н - ' . ) - П 0 ] 2

вы р а ж а е т с я через свой энергетический спектр 5 (со) следующим об­ разом:

со

Dy(({)=2

j (1 — coscoz',) 6* (ш) й?ш.

(1.19)

 

— оэ

 

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

25

ГЛАВА Л

УР А В Н Е Н ИЕ СВЯЗИ ВЫХОДНОГО СИГНАЛА

ДВ И Ж У Щ Е Г О С Я ПРИБОРА С ВЕЛИЧИНОЙ ИЗМЕРЯЕМОГО ПОЛЯ

§1. Вывод уравнения

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

Поэтому следует признать наиболее

типичным измерение

поля

д в и ж у щ и м с я прибором и рассматривать

неподвижный прибор

как

частный случай движущегося . Р а з н и ц а

между

движущимся и не­

подвижным прибором существенна только для

датчика прибора,

так как его выходной сигнал зависит и от его положения в простран­

стве и от времени. Остальные ж е части измерительного

прибора

не­

чувствительны

к положению относительно

исследуемой

среды,

а,

следовательно, и к характеру своего движения в ней.

 

 

 

 

 

Н а й д е м связь между величиной измеряемого поля и величиной

сигнала на выходе прибора при условии

произвольного

поступа­

тельного (без вращений) движения прибора в море.

 

 

 

 

 

Если датчик прибора неподвижен относительно начала

коорди­

нат, в

которых

рассматривается поле

Х(р;

т), то

нормированный

сигнал на выходе прибора дается выражением (1.9), причем радиус-

вектор

центра

датчика

г есть

постоянная

величина.

При

этом

У (г; t)

является явной

функцией координат

центра

датчика

и вре­

мени.

 

 

 

 

 

 

 

 

 

 

 

Если ж е датчик с течением

времени

перемещается в

рассматри­

ваемой

системе

координат, то

радиус-вектор его центра

является

функцией времени r(*f), и выходной сигнал становится явной функ­

цией только времени. Выясним, как изменится при этом

выражение

(1.9). Д л я этого представим его в иной форме. Считая

г = const,

26

произведем замену переменных г' = г — р и t' = t — т:

 

со

 

 

 

 

 

Y=\dr'

j

x { r - r '

;

t-t')H0(r';

 

t')df.

 

 

о

 

 

 

 

 

 

Найдем вид этого в ы р а ж е н и я

в

том случае, когда датчик пере­

мещается в пространстве так, что

радиус-вектор

его

центра меня­

ется по закону г (/).

 

 

 

 

 

 

 

П р е д п о л о ж и м вначале, что поле и датчик

одномерны. При не­

подвижном датчике

 

 

 

 

 

 

 

 

СО

СО

 

 

 

 

 

У (г,; 0 = j

dr\\x{rx-r\;

 

t-t')H0(r[;

 

i)

dt .

 

—со

0

 

 

 

 

 

Обозначая

 

 

 

 

 

 

 

 

 

x = t —

f,

 

 

(2.1)

 

 

P i = r i — r [ .

 

 

 

 

 

 

 

получим

 

 

 

 

 

 

 

К (г,;

0 =

] dr[]x(Pl;

i)H0(r\;

t)

dt

(2.2)

Пусть центр датчика

движется

по закону n = r i ( 0 -

Координата

pi и время х связаны с неподвижным полем, координата г' ряется относительно центра

иосителы-ю момента регист­

рации сигнала

на выходе

прибора (рис. 6).

 

 

Н а й д е м

связь

м е ж д у

пе­

ременными,

входящими

в

выражение (2.2), в некото­ рый момент времени т. В этот момент распределение

Рис. 6. Движение датчика в одно­ мерном поле.

Но

Pi

изме­

Pi.rt

поля

имеет вид X (pij т ) , местоположение центра датчика i\ =

= г\

(т), и соотношения (2.1) принимают

форму:

 

Pi — '"i (*)—>"i=r, (t —

t)—r[.

Подстановка этих величин в формулу (2.2) и возвращение к преж ­ ним переменным дает значение сигнала на выходе одномерного

27

прибора, движущегося в одномерном поле

со со

V(t)=

J

dPl\X[ri(t-^)-Pu

t--.\M0(Pl;

 

— со

0

 

Обобщим найденное выражение на случай трехмерного поля и

поступательного движения прибора

 

 

со

 

 

y ( 0 = J r f P J > [ r ( / — с ) - р ; t—z\H0(p-

-с) Л .

(2.3)

о

 

 

Полученное соотношение является исходным для всех последую­ щих построений. Оно определяет взаимосвязь между выходной ве­ личиной поступательно движущегося прибора (обычно электриче­ ским сигналом) и измеряемым физическим полем.

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

цифика ж е датчиков

заключается

в том, что они

осуществляют

трансформацию функций четырех

переменных

Х{р;

т) в

функцию

одной переменной У (г). Поскольку

измерительный прибор

создается

с целью определения

структуры именно поля

А'(р; т ) , необходимо

выяснение отличия величины выходного сигнала прибора от вели­ чины измеряемого поля на траектории движения прибора при изве­ стных параметрах прибора и законе его движения г (t), и оценка пределов применимости данного прибора при данной методике из­ мерения.

П р е ж д е чем переходить к дальнейшим расчетам, введем одно обозначение, а именно, определим аппаратную функцию прибора следующим образом:

Это дает возможность в выражении (2.3) изменить нижний пре­ дел интегрирования по времени и получить окончательный вид уравнения связи выходного сигнала прибора с измеряемым полем

V(t) = ^X [г ((-'.)-р;

t-i\H{?;

z)dpdt,

(2.4)

где интегрирование по всем переменным производится в бесконечных пределах. Сделаем одно замечание относительно поля X, входящего в это соотношение. При выводе формулы (2.4) не делалось никаких предположений о характере поля X, только считалось, что осуще­ ствима операция свертки его с аппаратной функцией прибора Н. Это условие выполняется как для детерминированных полей, удов­ летворяющих условиям Дирихле, так и для широкого класса слу­ чайных полей. Поэтому на соотношении (2.4) будет основан анализ полей первого и второго типов, т. е. будут решаться задачи, возни-

28

к а ю щ ие как при вертикальном зондировании, так и при горизон­ тальном движении прибора относительно среды.

Найдем спектр выходного сигнала прибора Y (t) при условии из­ мерения им детерминированного поля.

Пользуясь первой формулой (1.10), выразим функцию простран­ ственно-временного распределения поля Х(р; т) через его спект­ ральную плотность X (а; со)

^ ( р ; х ) = (2я)4 § х ( а >

° > 1 ) е х р [ / ( о р + м 1 т ) ] ^ 0 ^ ш 1

и подставим это в ы р а ж е н и е в формулу (2.4)

YW=~J5K)*~§X(a'

ш 0 е х р { / [ а г ^ - ^ - а р + ш ^ - ю , ^ } X

X

И (р; t ) dp dt da е Ц .

Поскольку одномерный спектр выходного сигнала есть

 

 

 

K ( w ) = J

y(t)e~Jal dt.

то отсюда

получим

 

У(ш)=

 

1

J -А'(а; ш,) е х р | у

[ar( * — <с) — ap-j-co,* ш^ —и/!]] X

( 2

i l ) 4

X Я (р; "0 dp d- da с Ц dt,

где, как указывалось выше, интегрирование по всем переменным производится в бесконечных пределах. Это выражение можно запи­ сать как

^ с ° ) = - ( 2 ^ г !

^ ^ ( P i ' ^ ' e x p l - y f a p + H J r f p X

X Л

j

exp {J [at (t -

x) + (со, - M ) (t - *)]} dt da

dw{.

Отсюда,

учитывая, что

 

 

 

H{u;

a)) = J / / ( p ; x ) e x p [ - y ( a p + ^)Jrfprft

(2.5)

есть спектральная характеристика аппаратной функции прибора, найдем спектр выходного сигнала в виде

У (to) 9 ^ 4 J-У («; u>,)//(a; со) С [a; (со, to)] da dco,, (2.6)

где функция

со

С (а; to) = j exp [у [ar (0 + ( u ^l) dt

—оэ

целиком определяется характером движения датчика r{t). В даль ­ нейшем будем называть ее функцией движения прибора.

29

Соседние файлы в папке книги из ГПНТБ