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

книги / Математические методы и планирование эксперимента в грунтоведении и инженерной геологии

..pdf
Скачиваний:
3
Добавлен:
12.11.2023
Размер:
16.89 Mб
Скачать

при 81Oil

О

_

 

X

г

J f tt?

g U iU l-X -u ^

Ai%

 

 

Sx~

. „ . ^ г

 

 

 

 

 

JL _ 1.

 

 

/Т\1J

X k

 

 

 

21а:«.

of(x-r

 

 

 

 

 

«Як

 

 

 

 

 

 

- экспериментальный

коэффициент регрессии x

на

, остальные

обозначения

те

же,

что

в

4*1.

Экспериментальная прямая регрес­

сии ос на ц.

 

 

 

 

 

 

 

 

 

 

 

 

х -

X

=:

*>х

 

 

 

 

 

 

 

 

 

 

 

 

где

 

 

 

_5ос _

£ u ^ c - ^ u i f

 

 

 

 

 

 

 

A t

 

 

4 * / з = z

g ,

 

 

 

%

 

 

 

 

 

 

 

 

 

 

 

u«< *-

4 *

 

 

 

 

 

 

 

 

 

 

 

 

- экспериментальный

коэффициент регрессии х

в а ^

. Коэффици­

енты регрессии

связаны

соотношением

 

 

 

 

 

4 .3 .

Множественная линейная корреляция

 

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

ходимые

для решения формулы для случая трех в е л и ч и н и £ .

В целях

сокращения записи

приводятся только формулы, относя­

щиеся к

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

регрессии % на х и

.

Плоскость

регрессии.

Плоскость регрессии

Н на (зс, ^ ) ха­

рактеризуется

уравнением

,

. д

,

ч

 

% -% =

 

V *

(* - ? )»

причем, коэффициенты регрессии -^я/х и -бг/а- определяются через коэффициенты корреляции между парами величин х и <&, х и , $ и з следующим образом

о

1-2*г) • . а _ ( z ai,-2 * s'2 * % )- З г

V " ..^

(7 “^

^

Коэффициенты корреляции между парами величия

3s * ( K - l ) S j S z

где JO - общее количество результатов опыта, а

.

Zx3 = C«-L) SX %

'

Сводный и частный коэффициент корреляции. Мерой зависи­ мости между величинами 3 и х , ^ служит сводный коэффициент корреляции

R =

V

1'

 

 

'

Сводный коэффициент корреляции Я

всегда заключён между 0 и ! •

Если R=Q

,

то между величиной 1

и величинами X , у,

нет линей­

ной корреляциОииой зависимости. При R - i

существует

линейная

функциональная зависимость % от

х и

.-

 

Для изучения, влияния только одного И8 факторов,

например

на величину^ , вводят частный коэффициент корреляций

величин

х % по отношению в величине ^

 

 

 

г ,* А '= . g x * - Z « y Z 4 g__________ _

Аналогично вычисляются частные коэффициенты корреляции

* Ц г/Х

Р А З Д Е Л П. ПЛАНИРОВАНИЕ ЭКСПЕРИМЕНТА

Глава I . ОСНОВНЫЕ ПОНЯТИЯ МАТЕМАТИЧЕСКОГО ПЛАНИРОВАНИЯ ЭКСПЕРИМЕНТА

1.1» Факторы,. параметры оптимизации и модели

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

В общем случае планирование эксперимента мояно подразде­ лить наг

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

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

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

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

факторами могут быть гранулометрический и

петрографический

составы,

влажность,

основность

(кислотность), содержание ор­

ганических примесей,

прочность

крупного заполнителя, возраст

и д р ., а

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

- объемная масса, плотность,

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

ния и сцепления, модули упругости и деформации, прочность

и

Др.

 

 

 

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

схе­

мой (рис. 18)• Переменные Х{ ,34, • . Хи. называют входными

 

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

ми (факторами),

переменные

 

Ъ л -

входными не­

контролируемыми (случайны­

ми переменными или помеха­

ми), переменные

,ц г , . . .

* * > * ^п. “ переменными

со­

стояния

(выходными параме­

трами),

которые характери-

вуют состояние

объекта

в

зависимости от

изменения

входных переменных.

 

Чаще всего

требуется

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

Каждый фактор может принимать в опыте одно из нескольких значений, называемых уровнями, фиксированный набор уровней всех изучаемых факторов в пределах одного опыта представляет условия проведения этого опыта (точка % \ рио» 19, а , б) . Про­ странство, ограниченное фиксированным набором уровней всех опытов плана эксперимента, называется факторным пространством (рис. 19, а ) . Оно задается координатными ооями, на которых от­ кладываются кодированные значения факторов. Для двух факторов факторное пространство ограничивается плоскостью, для трех и выше - более сложными многомерными геометрическими фигурами.

Каждому набору уровней изучаемых факторов в пределах од­ ного опыта соответствует свое определенное значение параметра

оптимизации» Связь между уровнями всех изучаемых факторов и соответствующими им значениями параметра оптимизации характе­

ризуется математической

моделью

у,

о с * , . - , x*v).

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

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

У-= 4 . * £ « 1 * 1 ,

(59)

для которой поверхностью отклика служит пдоокость. Ее исполь­

зуют при исследовании объектов с четко выраженной

линейной

зависимостью между факторами и параметром оптимизации.

 

Однако исследователю чаото приходится изучать объект

в

некоторой локальной области факторного пространства,

для

ко-

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

в котором 1 =1 ,

2, . . . , и. -количество управляемых факторов;

4 0 ,

i i A q tin . -

коэффициенты модели.

 

Для ^писания поверхности отклика полиномами первой

сте­

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

опи­

сания криволинейной поверхности отклика полиномами второй степени - центральные некомпозиционные планы на шестиугольни­ ке, центральные композиционные планы на кубе типа 63 «Ьц »&5 »

Нау и центральные композиционные ротатабельные планы.

При

изучении свойств (качества) материалов, зависящих только

от

соотношений компонентов, используют симплексные или симплексреиетчатые планы.

1 .2 . Последовательность планирования и обработки результатов эксперимента

Общая последовательность планирования эксперимента пока­

зана на рис. 20. Для построения модели исследователь

долхен

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

количестве

опытов, варьируя значения факторов по специально

разработан­

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

определить

с требуемой надежностью коэффициенты уравнения. При этом

каж­

дый фактор может иметь три уровня - минимальный ( - 1) ,

нулевой

( 0) и максимальный (+1) , которые называются кодированными

зна-г

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

max

X, Ч г Г - г \) /^ (zТ -гТ )!г-

 

 

 

П о с т а н о в к а

зад ачи

 

67

 

 

----- 7Г~ ~ Ь !

о

р

Основны х .

 

 

]

<2

ф сиаяоро&

—►2*

 

 

Ф а к т о р о в

ы л ж к ж г I ы i i m s s t s s s s U * - i> t\

CL

 

CL

 

 

 

 

Ъ

ЬыЪор формымодели

 

 

 

 

 

 

 

 

l

 

 

 

 

 

 

 

 

 

&ы&ор плаща.

 

 

 

5 Jрандомизоиеизсолыщ

Уетыно4лжна« числа

 

 

6 ftctp#AJtShHUX 0**TV4

 

 

7 j П о с т а н о в к а

эксперим ент а

 

 

 

<?

I

расчет

построчных

Ъиеяерсий

 

 

 

°

1

 

пяряжлеЛнщ олы/яоо

 

 

 

?<Х

 

 

П роверки

-

Ует

 

 

числя пяряллелмнх опЫтоо

 

 

 

 

 

 

 

Д а .У

Vtr =g Vr

Нет

 

 

SS

 

ftp»Бернес значимости чиеленкьа

 

 

значения параллельных ольто4

 

 

 

 

 

 

 

А # Л

г*-*? 5гт

Н*т

 

 

 

 

 

 

 

Ц

 

 

 

 

 

 

оЪнорt>btfости е££&(* Эисаереий

 

 

р т

[ИС4*Ни*Яи?/ГШРС4си

\^я

Гп Г

~~pcfC4€tn

 

з :'*СЛР0Ц1*+Ъимсс*м/

 

\коэрфициенто4 регрессии

 

и

р а с ч е т

Д и с п е р с и й

 

 

п р • < * / > & * ^ н т и и м е е г ч

 

К о о ф с р и ц м е н т о 4

 

Н а к о у р ф и ч и с / г г ** р е гр е с с и и

 

 

j/Й | Расист дисперсий адекбгтнисгпи. ]

л

 

 

 

 

 

 

Т

 

 

 

 

flpofcptier cietc/amHoctrtUмоЯели

Ну»

 

-

w

* 1

 

 

2 p a rp u < fc c K u jr ш -

\ S -

Тщ р н р л гя ч и я & t»% t4 u \

К Ч

 

 

______2Sk

Рио. 20.

П о с т р о е н и е м од е- л и с Н Я Т У Р Ы Л * -

Н Ш M p tM tH H H X

где

0CL - кодированные значения I - г о фактора;

натураль­

ные

значения L - ro фактора; OL “ интервал варьирования

1 -го

фактора.

 

 

 

После установления основных факторов и параметра

оптими­

зации исследователь выбирает форму модели и соответствующий

ей план эксперимента, реализация которого позволяет

осущест­

вить комплексное воздействие факторов на параметр

оптимиза­

ции.

 

План эксперимента записывается в виде таблицы - матрицы планирования. В табл. 2 приведена в качестве примера матрица планирования двухфакторного эксперимента с минимальным числом опытов, необходимых для построения модели первого порядка.

Факторное пространство этого плана представляет собой квадрат (рис. 21).

Цифрами I , 2 , 8, 4 на рис. 21 показана номера опытов,

Т а б л и ц а 2

Л

 

Кодированные

Пара-

Кодовое

 

опы­

значения

метр

обозна­

 

та

-факторов•

оптими­

чение

 

 

 

 

зации

строк

 

 

х<

ха.

ч

 

 

 

 

 

 

 

I

-X

- I

 

( I )

>3

2

+1

- I

Чг

а

Рис. 21,

8

- I

+1

Чъ

б

4

+1

+1

Чь

аб

 

 

 

 

 

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

мер,

опыт №I проводится при кодированных значениях факторов

Э£|

, Ij, —I •

 

 

 

В последнем столбце матрицы приведены кодовые обозначе­

ния

ее строк. Буква "а" обозначает, что соответствующий

ей

фактор x L находится на верхнем уровне

(х ^ + 1) ; паб" показыва­

ет,

что на верхних уровнях находятся

оба фактора; символ

(I )

обозначает, что оба фактора - на нижних уровнях. Кодовые обо­

значения значительно сокращают запись матриц.

Так, матрица,

приведенная в табл. 2, записывается следующим

образом:

( I ) ,а ,б ,а б .

 

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

При проведении эксперимента исследователю необходимо для каждого опыта при фиксированных, уровнях факторов (в соответст­ вии с планом) получить среднее значение параметра оптимизации с требуемой надежностью. При этом можно воспользоваться мето­ дом [ 8] , основанном на предварительном (из рекогносцировочной серии измерений параметра оптимизации в одном из .опытов пла­ на) определении меры изменчивости (\г) как отношения, установ­ ленного в опыте среднеквадратического отклонения (<э) к стати­ стическому среднему, т . е . :

V = ( ^ /

(62)

По полученному значению меры изменчивости и требуемой доверительной вероятности (£) по табл. XI Приложения £8] , ус­ танавливается необходимое число измерений параметра оптими­ зации в данном опыте.

Очень важно исключить из полученных данных грубые ошибки. Для оценки сомнительных максимальных или минимальных резуль­ татов можно воспользоваться Z -критерием [ I ] :

где $

- среднеквадратическая ошибка измерения.

 

 

Полученное значение *4 -критерия сравнивают с табличным,

принимаемым по табл, ХП Приложения в зависимости от

уровня

значимооти и количества измерений в опыте. Необходимое

усло­

вие

- Z p ^ Z r . Если

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

ся грубым ( "промахом11) и в раочет не включается. После

под­

счета

дисперсий и определения средних значений параметра оп-

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

K - i

Затем из всех дисперсий находится наибольшая, которая делится на сумму всех дисперсий плана

Дисперсии однородны, если расчетное значение критерия Кохрена

(Кр)

не превышает табличное (K j),

которое устанавливают

по

табл.

УШ Приложения. Бели КР> ^т•

необходимо

выявить и устра­

нить

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

использовать

бо­

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

Для оценки значимости коэффициентов модели и проверки ее

адекватности

нужно подсчитать дисперсию среднего

но формуле

 

 

 

 

(63)

в которой X -

количество опытов, К - количество

измерений

в

опыте.

 

 

 

 

Методы расчета коэффициентов модели и проверки их значи­

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

пост­

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

 

каж­

дого плана. v

 

 

 

 

Соседние файлы в папке книги