Файл: Доценко С.В. Теоретические основы измерения физических полей океана.pdf

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

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

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

Добавлен: 15.07.2024

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

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

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

только от свойств датчика и показывает, какую долю в общем зна ­ чении выходной величины 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