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

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

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

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

Добавлен: 09.04.2024

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

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

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

718 Глава 10

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

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

10.1. В Х О Д Н Ы Е С И Г Н А Л Ы

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

следующие три входных

сигнала:

ступенчатый,

импульсный

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

и синусоидальный

(фиг. 10.1.1). Из стохастиче­

ских входных

сигналов

наиболее

часто используется случай­

ная последовательность прямоугольных импульсов (случайная

бинарная последовательность);

иногда может использоваться

сам входной сигнал. Случайные

сигналы рассматриваются в

гл. 12.

 

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


x(t)-U(t-a)

Вреліенная

Пространство изод-

ражений по Лапласу

 

область

 

 

U(t-

,

(0,t<a

 

 

 

 

a

 

 

 

a

 

 

 

ô(t-ot)=0,

t*a

x(t)=acosù)t

Ф и г . 10.1.1. Типичные входные

сигналы .

а единичный ступенчатый сигнал; б импульсный

сигнал; в — синусоидальный1

сигнал; г — случайная последовательность прямоугольных импульсов (случайная бинар­ ная последовательность).

720

Глава

10

области

любой из трех входных

сигналов позволяет эффективно

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

ных систем, показана

в работе

Найквиста

и

др.

[1].

На данном

этапе

необходимо

остановиться

на

некоторых

преимуществах

и недостатках использования

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

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

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

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

рений,

чем при

ступенчатом входном сигнале.

И

наконец,

в случае реальных процессов для генерирования

или подачи синусоидального входного сигнала требуется более сложное оборудование, чем в случае ступенчатого или импульс­ ного сигнала, и большие затраты времени, поскольку необходимо рассматривать несколько частот и на каждой частоте должно быть достигнуто установившееся состояние (для чего может потребовать­ ся несколько часов). При синусоидальном входном сигнале допу­ скается внесение в процесс лишь небольших возмущений. Однако, как указывалось ранее, в случае нелинейного процесса частота выходного сигнала смещается относительно частоты входного сигнала. Хоуген [2] приводит дополнительные сведения о некото­ рых практических вопросах, связанных с выбором входного сигнала.



Оценивание параметров уравнений

в

частных производных

721

10.2. Р Е А К Ц И Я С И С Т Е М Ы

Н А

В Х О Д Н О Й С И Г Н А Л

 

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

Таблица 10.2.1

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

сраспределенными параметрами г )

Процесс

Передача тепла к трубе с по­ током ж и д к о ­ сти

Модель процесса

dT(t,z)

 

,

dT(t,z)

dt

1

1

dz

= h[Ta-T(t,

z)\,

T{0,z)

=

T a

[ \ - e - ^ z ] ,

 

T(t, 0)=T0U

(t)

Адсорбция

в Жидкость:

 

 

слое насадки

dcL

 

 

dcL

 

 

- d f + V L

^ z -

=

 

= —k(cL

— c L i e )

Решение

T(t,z) = Ta[l-e-WVz]

+

+ T0u ( < - e - ^ z

z*

 

= 1 - (

e-(t*+7*)x

0

 

X I 0

(2]Л*2*) dz*,

Твердое вещество:

z* =

 

 

W Z

с = 0, * < — , z < 0 ,

 

 

c(t,0) = cLOU(t)

!) U — единичная ступенчатая функция, T — температура, с концентрация, г — координата вдоль оси, ѵ — скорость, t — время, h — межфазный коэффициент теплопере­ дачи, h — межфазный коэффициент массообмена, /о функция Бесселя, знаком ~ обо­

значена переменная интегрирования.

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

Импульсный входной сигнал можно рассматривать либо как граничное условие для модели, либо как член, характеризующий источник в самом дифференциальном уравнении. Например, импульсный входной сигнал при z = 0 (начало оси z) и t = 0


722

 

 

 

 

 

Глава 10

 

 

 

 

можно

ввести

в модель

как

источник:

 

 

 

 

 

 

- g - + I , . g - ==aô(*)ô(s),

y(0-,z)

= 0,

(10.2.1)

где

член aô (t) ô (z)

описывает импульс

в

точке z = 0

в момент

£ =

0 (в соответствии со свойствами дельта-функции).

Входной

сигнал

также

можно

рассматривать

как

граничное

условие:

 

 

 

л + у

&

- и '

y = (t,

0) = ßo(0-

 

 

В уравнении (10.2.1) ѵ обозначает скорость потока

жидкости,

а у

— зависимую переменную. Поскольку

по определению ô (х)

^ ô (x) dx = 1, то ô (t) имеет размерность (время)"1 , а ô (z) —

оо

(длина)- 1 ; размерности a и ß должны быть такими, чтобы диф­ ференциальное уравнение было согласованным с точки зрения размерностей. Допустим, что некоторое фиксированное количе­ ство индикаторного вещества (трасера) т, равномерно распреде­ ленное по поперечному сечению А трубы или канала (по которому движется поток), вводится в трубу в некоторой фиксированной точке (z = 0) в момент t = 0. В этом случае коэффициент в члене, характеризующем источник, принимает значение a = ml А.

Уравнение (10.2.1) можно решить, выполняя последовательно преобразования Лапласа для обеих частей уравнения по t и z, группируя подобные члены и переходя затем к обратным пре­ образованиям:

y(t,z)=-%.ô(t-±). (10.2.2)

Соотношение (10.2.2) можно интерпретировать следующим образом: в момент t = z/v входной импульс появляется в точке z, все еще сохраняя форму импульса. Поскольку в этой модели не учитывается влияние диффузии или продольного перемешивания, такая интерпретация вполне допустима. Импульс распространяет­

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

и

его можно

измерить

на

выходе

трубы.

 

 

 

 

 

 

 

В более сложной модели учитывается осевая дисперсия (связан­

ная с

диффузией

или

продольным

перемешиванием):

 

 

^

L +

V ^ L =

D

^

+ Источник,

(10.2.3

где D — эффективный коэффициент

дисперсии, а член, характе­

ризующий источник, такой же, как и ранее, если вещество равно­

мерно

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

сечении

z =

z0 трубы или канала:

 

 

 

Источник = - ^ ô ( z — z 0 (t),

(10.2.4)