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

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

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

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

Добавлен: 09.04.2024

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

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

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

 

 

 

 

 

Стратегия

эффективного

 

 

 

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

 

 

 

 

635

 

Brooks

S. H . , Mickey M . R., Optimum Estimation of Gradient Direction

in

Steepest

 

Ascent

Experiments,

Biometrics,

17,

48

(1962).

 

 

 

 

 

 

H i l l

W . .T., Hunter

W . G., A

 

Review

of Response Surface

Methodology;

A

Literature

Survey.

Technometrics,

8,

571

(1966).

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Эволюционное

планирование

 

 

 

 

 

 

 

 

 

 

Box G. E. P., A Simple System of

Evolutionary

Operation

Subject

to

Empirical

Feedback,

Technometrics,

 

8,

 

19

(1966).

 

 

 

 

 

 

 

 

 

 

 

Box G. E. P., Hunter J. S., Condensed Calculations for

Evolutionary

Operation

Programs,

Technometrics.

 

1,

 

77

(1959).

 

 

 

 

 

 

 

 

 

 

 

Delver A . , Organon N . V . , Gradual Optimization of Production Processes,

Sigma,

7 (3). 45

(1961).

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Techno­

 

Hunter

W . G.,

K i t t r e l l J. R., Evolutionary Operation,

A

Review,

metrics,

8,

389

(1966).

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Симплексный

метод

 

 

 

 

 

 

 

 

 

 

 

Carpenter В . H . , Sweeny H . С , Process Improvement

w i t h

Simplex Self-

Directing

Evolutionary Operation,

Chem.

Eng.,

1 1 7 (July

5,

1965).

 

 

 

Lowe C. W . , Trans.

Inst.

Chem.

Eng.

(London),

42,

T334

(1964).

 

 

Spendley W . , Next

G. R., Himsworth

F. R., Sequential

 

Application

of Simplex

Designs

in

Optimization

and

Evolutionary

Operation,

 

Technomet­

rics,

4,

441

(1962).

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Оптимизация

с

адаптацией

 

 

 

 

 

 

 

 

 

 

Rees N . W . , Self-Adaptive Control Systems, British

Chem.

Eng.,

5,

106

(Feb. 1960).

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Последовательные

планы для увеличения точности

 

 

 

 

 

 

 

 

 

 

 

 

 

 

оценивания

параметров

 

 

 

 

 

 

 

 

 

 

 

Atkinson А. С , Hunter W . G., The Design of Experiments

for

Parameter

Estimation.

 

Technometrics,

10,

271

 

(1968).

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Box G. E . P., Lucas H . L . , Design

of

Experiments

i n

Nonlinear

Situ­

ations,

Biometrics,

 

46,

77 (1959).

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Draper

N . R., Hunter W . G., The Use

of Prior

Distribution i n the

Design

of Experiments

for Parameter Estimation

in Nonlinear Situations,

 

Biometrika,

54.

147 (1967).

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Ann.

 

Elfving

 

G., Optimum Allocation in

Linear

Regression

Theory,

Math.

Stat.,

 

23,

255

(1952).

 

 

 

 

 

 

 

 

/ . Royal

 

Stat.

Soc,

 

 

 

Kiefer

J. С ,

Optimum Experimental Designs,

 

B21,

272

(1959).

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Kiefer

J. C , Optimum Experimental Designs V ,

w i t h Applications

to

Systematic

and

Rotatable

Designs,

Proc. Fourth

Berkeley

Symposium,

1,

381

(1961).

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

K i t t r e l l

J. R.,

Hunter W . G., Watson

С. C ,

Obtaining Precise

Parameter

Estimates for Nonlinear Catalytic Rate

Models,

AIChE

J.,

 

12,

5

(1966).

 

Stone M . , Application of a Measure of Information to the

Design and Com­

parison

of

Regression

Experiments,

 

Ann.

 

Math.

Stat..

30,

55

(1959).

Ann.

 

Wald A . , On the Efficient Design

of

Statistical

Investigations,

Math.

Stat.,

 

14.

134

(1943).

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 


63(5

 

 

 

 

 

 

 

 

 

 

Глава

8

 

 

 

 

 

 

 

 

 

 

 

 

 

Последовательные планы для различения моделей

 

 

 

Behnken

D . W . ,

Estimation of Copolymer

Reactivity

Ratios,

J.

Polymer

Sei., A2 ,

645

(1964).

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Box

G.

E.

P.,

H i l l

J.

W . ,

Discrimination

among

Mechanistic

Models,

Technometrics,

 

9,

57

(1967).

 

 

 

 

 

 

 

 

 

 

 

 

 

Box

G. E . P., Lucas H . L . , Design of Experiments in Nonlinear

Situations,

Biometrika,

 

 

46,

77

(1959).

 

 

 

 

 

 

 

 

 

 

 

 

 

Hunter

 

W .

G.,

Mezaki

R., A n

Experimental

Design

Strategy

for

Dis­

tinguishing

Among Rival

Mechanistic

Models,

Can. J.

Chem.

Eng.,

45, 247

(1967).

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Hunter

W .

G., Reiner A . M . , Designs

for Discriminating

between

Two

R i v a l Models,

Technometrics,

1, 307

(1965).

 

 

 

 

 

 

 

Ann.

Keifer J., Wolfowitz J.,

Optimal

Designs

in

Regression Problems,

Math.

Stat.,

 

30,

271

(1959).

 

 

 

 

 

 

 

 

 

 

 

 

 

Nelson

 

A . C , Some Optimum Regression Designs for Discriminating

between Two Models, Research Triangle

Institute,

Techn. Rept. 4,

Durham,

N . C . , Jan. 27, 1967.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Peterson

 

D .

W . ,

Discriminant

Functions

(AD465281),

Stanford

Elect­

ronics

Lab.

Stanford,

Calif.,

A p r i l

1965.

 

 

 

 

 

 

 

 

 

Roth

P.

 

M . , Design

of

Experiments

for

Discrimination

Among

Riva l

Models,

 

Ph.

 

D .

Thesis,

Princeton

U n i v . ,

1966.

 

 

 

 

 

 

 


Часть III

ОЦЕНИВАНИЕ МОДЕЛЕЙ, ОСНОВАННЫХ НА ПРИНЦИПАХ ЯВЛЕНИЙ ПЕРЕНОСА

В разд. 1.1 уже отмечалось, что модели процессов могут осно­ вываться на эмпирических соотношениях, балансах популяций и принципах явлений переноса или на их комбинациях. В части I I I описывается, каким образом можно оценивать параметры основан­ ных на принципах явлений переноса моделей, включающих обыкновенные дифференциальные уравнения, уравнения в част­ ных производных и передаточные функции (частотную характери­ стику). Д л я исследователя могут быть полезными четыре различ­ ных вида информации: 1) сама форма модели, 2) параметры модели (коэффициенты), 3) входы модели и 4) выходы модели. Задачи, которые здесь рассматриваются, можно классифицировать по сле­ дующей схеме (известные факторы отмечены символом Ѵ ) :

Вид информации

И д е н т и ф и к а ц и я м о д е л и О ц е н и в а н и е п а р а м е т р о в П р е д с к а з а н и е О б р а т н а я з а д а ч а

Модель,

 

 

 

включаю­

Парамет­

Вход

Выход

щая

ры

модели

модели

граничные

модели

 

 

условия

 

 

 

V

 

V

V

V

V

V

V

V

V

V

V

 

Весьма общей и чрезвычайно сложной задачей является

описа­

ние процесса с помощью подходящей модели, или проблема

иденти­

фикации.

Если задан определенный класс моделей процессов

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


638 Глава 9

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

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

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

Глава 9

ОЦЕНИВАНИЕ В МОДЕЛЯХ ПРОЦЕССОВ,

ПР Е Д С Т А В Л Е Н Н Ы Х О Б Ы К Н О В Е Н Н Ы М И

ДИ Ф Ф Е Р Е Н Ц И А Л Ь Н Ы М И У Р А В Н Е Н И Я М И

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

1. Анализ баллистических траекторий.

2.Анализ работы химических реакторов.

3.Обработка сигналов.

4.Кибернетика.

5.Явления переноса в человеческом обществе.

Здесь будет рассмотрено оценивание параметров и несколько менее подробно — оценивание состояния. Оценивание «состояния» процесса означает оценивание зависимых переменных этого про­ цесса. Можно различить три типа оценивания состояния. Если


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

обыкновенных дифференциальных

уравнений

639

проведены наблюдения

в интервале от t0 до tf, то оценивание

век­

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

как:

 

1.

Интерполяция,

или

сглаживание, если t

< tf.

 

 

2.

Фильтрация, если t

= tf.

t >

tf.

 

3.

Предсказание,

или

экстраполяция, если

 

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

Здесь нет возможности достаточно широко осветить ряд прак­ тических расчетных аспектов оценивания, таких, как:

1.Время, требуемое для проведения измерений и вычислений.

2.Выбор используемых пробных входных функций.

3.Эффективность установки для сбора данных.

4.

Скорость сходимости

процедуры оценивания.

5.

Тип и расположение

приборов для обработки данных.

Все эти факторы не очень просто поддаются обобщению, и луч­ ше всего их иллюстрировать специальными примерами, имеющи­ мися в литературе. Требования к экспериментированию и оцени­ ванию моделей, изложенные в предыдущих главах, предъявляются также и к моделям, основанным на явлениях переноса. Экспери­ менты следует проводить в соответствии с хорошим планом (гл. 8), который должен быть составлен таким образом, чтобы ошибки при измерениях были независимыми (гл. 5), а модель, коэффици­ енты которой требуется оценить, должна быть адекватной (гл. 7).

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

Естественно начать изложение с обзора типов моделей, которые будут рассматриваться в дальнейшем, и описания способа, с по­ мощью которого предполагается ввести ненаблюдаемые ошибки в детерминированную модель процесса. Затем будут применены два метода оценивания — метод наименьших квадратов и метод максимального правдоподобия, описаны и проиллюстрированы три несколько различающихся метода оценивания: последователь­ ное оценивание, квазилинеаризация и метод «ошибки уравнения». Будут рассмотрены случаи как дискретных, так и непрерывных