Файл: Химмельблау Д. Анализ процессов статистическими методами.pdf

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

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

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

Добавлен: 09.04.2024

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

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

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

Оценивание

параметров

уравнений

в частных

производных

745

где С — наблюдаемый отклик, функция времени; — оценка момента т^; М# — оценка момента т # , то необходимо исследовать

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

свойства

интегралов І0, It и І2. Интегралы І0,

It и / 2 статистически

независимы.

Рассмотрим

интеграл / :

ь

I = j ф (t) X {t) dt,

а

где ф (t) — детерминированная функция, X (t) — нестационарная случайная функция. При заданной функции ф (t) и заданных пределах интегрирования интеграл / является случайной величи­ ной. Распределение случайной величины / зависит от выбранной функции ф (г), функции распределения случайной функции X (t) и пределов интегрирования. В частном случае, когда X (t) — нормально распределенная случайная функция, интеграл / явля­ ется нормально распределенной случайной величиной и

 

ь

 

g { / } =

]<p(t)${X(t)}dt,

(10.4.8)

 

а

 

 

b

 

Var {/} =

j ф (t) Var {X (t)} dt.

(10.4.9)

 

а

 

Если заданы математическое ожидание и дисперсия случайной функции X (t), то математическое ожидание и дисперсию инте­ грала / можно вычислить по формулам (10.4.8) и (10.4.9). Д а ж е в том случае, когда X (t) не является нормально распределенной случайной функцией, но интеграл аппроксимируется суммой и от­ дельные наблюдения Xt (t) независимы, согласно центральной предельной теореме, случайная величина / имеет приближенно нормальное распределение.

Поскольку функция X (t) является нестационарной,

то такие

ее статистические характеристики, как математическое

ожидание

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

как

описывалось в разд. 2.1. Следовательно,

математическое ожи­

дание случайной

функции

X (t) является

функцией

времени,

и

несмещенную

оценку

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

ожидания

% {X (t))

можно вычислить по п результатам, заоегистоиоованным в любой


746 Глава 10

заданный момент времени t, по

формуле

 

п

 

 

£ * ( 0 = 4-2Xt<*)'

*' = ! . • • • . * .

(Ю.4.10)

г=1

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

Дисперсия оценки х (t) в любой момент времени t имеет вид

пп

Var &х (t)} = ±

2 2 C o v

WХ * (<)) =

 

 

 

 

 

 

i=l ft=l

 

 

 

 

 

 

2

 

 

n

n

 

 

 

 

 

 

= ^

 

+ 7^2 2 8 { № W - | » i W ] f t W - | 4 W ] } .

(10.4.11)

Если

Xt

(t)

и

X f e

(t) независимы,

то

равенство (10.4.11) прини­

мает

вид

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Var

(*)} = ±

й .

(Ю.4.12)

Если

Аі->-0,

т. e,

n =

(tf/At)

—>- о о , то выборочное

 

среднее

•сходится

к

среднему

непрерывной

случайной величины.

Если

ХІ (t)

и Xj, (t)

коррелированы, то необходимо вычислить

двойную

сумму, входящую в формулу (10.4.11), и при увеличении

объема

выборки п

дисперсия

оценки \ix (t) может не снизиться до доста­

точно малой величины, что зависит от результата двойного сум­

мирования. Если при большом

объеме выборки п

с увеличением

к происходит

уменьшение г х х

(к, t),

то формула

(10.4.11) при­

нимает вид

 

 

 

 

 

ѵ„&(і»-^+42 ( І _ І ) Х

 

 

 

 

ft=i

 

 

 

 

X [ r „ ( f t , 0 - n i ( < ) ] « ^ ) + -|-S

* " ( M ) - i & ( 0 ,

(10.4.13)

 

 

fc=i

 

 

 

где R x x (к, i)

— оценка автокорреляционной функции в момент t.

Применяя

формулу (10.4.8) к интегралам / 0 , І1

и / 2 ,

получаем

g { / 0 } = Ç g { C } Ä = j с A ,

оо

oooo

Ô0

oo

oo

g { / 2 } = j *»g{C}d* = j *c a A .


Оценивание параметров уравнений в частных производных 747

Если Мі

вычисляется как (q2imV)

Ilf то Mi

является

несмещенной

оценкой:

 

 

 

оо

 

 

 

 

 

 

 

 

 

Ш

<Л / Л = -&

%

= -& \

tc dt = ті.

 

 

 

 

 

о

 

 

Если же

Ali

вычисляется

как

отношение IJIJIQ),

ТО

»<•*->--?-«(-&)•

и для определения величины смещения оценки момента піі необ­ ходимо использовать формулу (10.3.9).

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

ных может

оказаться, что

вместо вычисления u c (t)

и о2-

(t)

и определения с

помощью

этих величин оценки % {Ih}

проще

вычислить Ik

для

нескольких

экспериментов

и найти выборочное

среднее Ih:

 

 

 

 

 

 

 

 

п

 

 

 

 

 

 

і = 1

 

 

 

 

Трудно что-либо сказать

о возможности

использования

М^

для оценивания коэффициентов дифференциального уравнения.

Поскольку

соотношение

между ЛІ#

и интегралами / 0 ,

It и І2

является

нелинейным,

то величину

смещения оценки

g

определить трудно. Дисперсию для Л/# можно оценить путем проведения повторных экспериментов.

Какими бы ни были распределения С (t), Ik и AIh, для опреде­ ления дисперсий величин р с (t), Ik и AIh соответственно можно использовать неравенство Чебышева [формула (3.3.8)].

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

ші =1 + -р--

После вычисления

t=i

по результатам независимых экспериментов имеем

Ѵ а і Г { Л / і } = і £^ _ (Mii-MiY ,


748

Глава 10

и, следовательно,

V a r { / > } « ( - Ç ) 2 V a r > f i } -

Как работа Левеншпиля и Смита, так и работа Ван-дер-Лаана основаны на предположении о том, что ввод индикаторного веще­ ства описывается дельта-функцией. Дельта-функция представляет собой математическую идеализацию, так как требует введения конечного количества индикаторного вещества за нулевое время,, и в действительности ее можно лишь аппроксимировать. В работах [11—13] описан метод, не требующий применения идеальной дельта-функции. Метод предусматривает измерение концентра­

ции в двух точках,

находящихся на

контрольном участке, а н&

в одной точке, как

предполагалось

выше. Детерминированные

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

находится

разность

между ними. Эту разность можно

связать

с параметром Р и,

следовательно, с коэффициентом дисперсии.

Не имеет

значения,

где именно индикаторное вещество

вводится

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

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

через

z4 координату первой точки, а через z2 — координату вто­

рой.

Пусть

 

 

 

Ami = (mi)Z2

 

(10.4.14)

 

Am# = (m#)Z2~(m#)zl.

(10.4.15)

Математическое ожидание и дисперсия АМ± (или AM#)

имеют вид

 

Ш {AMJ = % { ( i ^ U -

% {{МЛ*},

(10.4.16)

 

Var {АМІ} = Var {(Mt)Z2}

- j - Var « M , ) Z l } .

(10.4.17)

Бишофф и Левеншпиль [13] представили также решения, свя­ зывающие уравнения (10.4.14) и (10.4.15) с безразмерным пара­ метром Р для проб, взятых как в пределах, так и за пределами контрольного участка. Два таких соотношения приведены в таб­ лице 10.4.1. При проведении экспериментов крайне желательно выполнять измерения в точках, достаточно удаленных от концов сосуда, с тем чтобы краевые эффекты были незначительны; в этом случае могут применяться исключительно простые выражения для бесконечной трубы. Другой, еще более важной причиной использования выражений для бесконечной трубы является то,


Оценивание

параметров

уравнений

в частных

производных

749

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

Пример 10.4.1. Оценивание коэффициента осевой

дисперсии по экспериментальным данным

Результаты, полученные в двух контрольных точках, для экспериментальной схемы, соответствующей конфигурации 4 в табл. П.10.4.1, представлены на фиг. П.10.4.1; предполагается, что регистрировалась непосредственно зависимая переменная С.

/ \

а

 

А

s,.

 

\ i

 

 

 

 

 

 

 

 

 

 

t,-a

^

1

I

t

 

/

 

 

 

 

 

 

 

 

 

5

А

 

іV

 

 

 

 

 

 

 

 

 

 

t

Ф и г . П.10.4.1.

данные для первой контрольной точки, Л — начало регистрации данных; б — дан­ ные для второй контрольной точки.

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

h = к + (ô2 - Ôt ), .