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

книги из ГПНТБ / Авдеев, Н. Я. Аналитико-статистические исследования кинетики некоторых физико-химических процессов учеб. пособие

.pdf
Скачиваний:
22
Добавлен:
21.10.2023
Размер:
10.28 Mб
Скачать

(4)

Параметры распределения а, р, п функции (1) в каждом кон­

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

по одному из

следующих способов.

 

через х моду

 

Способ характеристических

точек. Обозначая

и

х', х" абсциссы точек перегиба функции (1) и

приравнивая

к

нулю производные f'(x),

f"(x) в этих точках, получим си­

стему уравнений, которая для вычисления параметров распреде­ ления а, р, п по характеристическим точкам (л;, х', х") функции

(1) приводится к численному решению относительно параметра р "уравнения

{аЬ)р + ^ ± ^ = Р ± 1

(б)

 

р — 1

Р — 1

 

и вычислению параметров а,

п по формулам:

 

п = 1 +

1

а = П— 1

(6)

1 — (ab)P ’

р х Р

 

где введены обозначения: а =

х' : X-, b = х"

: х.

Способ пропорциональных точек. Если выбор точек интер­ полирования кривой плотности распределения (1) подчинить условию пропорциональности

хг : х2 = х2 : хз = хз : х4 = D0

(7)

и задать значения функции / (х) в этих точках, то. после соот­

ветствующих преобразований получим формулы (D2 • D 2 >

0):

р =

ln ( D i : Da) : ln D 0, а

= — D l

(Da — Dt) * x ~

( 8)

n -

l ==[~b f - D~ + D3)

:ln D °’

D3 =

ln (h ■/а).

 

Di =

\n(fifz :fl),

 

О, =

ln (fj, : fl)

 

Способ моментов. Полагая в (3) т — 1, 2, 3, получим систему трех уравнений, которая может быть использована для вычисле­

10

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

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

 

 

 

 

 

X

 

 

 

 

 

 

 

Ф

(X) = А

I" е ~ а*Р х п ~1сіх

 

(9)

 

 

 

 

 

о

 

 

 

 

рассмотрены в работах [1, 2].

 

 

 

 

Способ кратных

параметров.

Если в

(9) положить п =

кр,

где /г — натуральное число,

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

по­

лучим:

 

 

 

 

 

 

 

 

 

Ф (х) — 1

 

 

 

,к -\

( 10)

1 — 2Г-|------1------ )—... —J—

 

 

 

 

 

 

2!

31

(й— 1)1

 

где введена подстановка

 

 

 

 

 

 

 

 

 

г — ахр

 

 

(11)

Применяя

к

уравнению

(11)

метод

наименьших квадратов

[12], получим

[13]:

 

 

 

 

 

 

 

X

г‘

хі — “ lg Пг«lg Пхі

 

 

 

 

 

 

т

 

 

 

 

( 12)

 

 

 

 

 

 

 

 

 

 

 

 

Z

ig“ xl — lg'-Пх/

 

 

 

 

 

 

 

Ml

 

 

 

 

 

где m >

2 — число точек интерполирования; zt— находятся

по

опытным

данным f

{x) в точках

хь х2, хз,...,

хт и табличным

значениям функции

(10).

 

 

 

 

 

При т =

2 формулы (12)

принимают более

простой вид:

 

=

lg(га : Zi)

 

 

 

 

 

 

 

lg(*.:*i)'

 

 

 

 

 

При вычислении параметров а, р, п = kp по трем и большему числу точек интерполирования кумулятивной кривой распре­ деления по формулам (12) можно пользоваться, например, сле­ дующей схемой записи вычислений:

11

Лэ

Фэ(-Ѵ)

Ч

п/п

 

 

 

1

2

3

4

1

Л1

Ф(Ч)

Ч

2

Л*о

Ф(х„)

Ч

3

А‘з

Ф(Ч)

Ч

тхт Ф{хт) гт

Суммирование по столбцам

Вспомогательные

вычисления

Вычисления' по формулам

lg л-

ig«

lg* lg*

lg2 *

гР

фр (*)

\ф р -

 

 

 

 

 

 

-ф ,\

5

6

7

8

9

10

11

lg Л‘і

ig ч

lg*i lg«,

lg“ gj

Ч

Ф(хі)

ДФ,

lg *'о

lg г.,

lg*a lgzä

lg2 *2

Ч

Ф(х„)

АФ,

lg *3

lg 4

lg*3 lgz3

lg“ *3

ч

Ф(*з)

АФІ

lg*m

lg 4n

 

lg2 *m

г,и

Ф(хт)

АФт

а

b

А

В

Среднеабсолютная

погрешность

 

 

 

 

Ä = — ---- ",

%

 

 

 

 

 

rn

 

 

 

 

 

Среднеквадратичес­

 

ab

а1

 

кая

погрешность

 

C = —

D = —

 

-V- m

 

 

rn

m

 

 

 

 

 

 

 

 

А — С

I

Ф—ра)

Критерий согласия

Р =

m

— —

а = 10

 

по К. Пирсону

 

B — D

 

 

 

 

 

Графы

1, 2, 3 заполняются по данным

наблюдений; графа

4 — по табличным ізначениям функции

(10),

графы 5

и

6 — по

табличным

значениям логарифмов xt

и г,-,

графы 7

и

8 — по

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

Графы 9 и 10 заполняются для получения оценки (графа 11) согласованности опытных (графа 3) и расчетных (графа 10) ре­ зультатов.

При вычислении параметров распределения по формулам (13) по двум точкам интерполирования узловые точки следует выбирать в наиболее характерных местах кумулятивной кривой. Например, если кумулятивная кривая распределения задана графически, то узловые точки следует брать по разные стороны от точки перегиба в местах наибольшей кривизны кривой. Если кумулятивная кривая не имеет точки перегиба, то узловые точки выбираются по разные стороны от наибольшей кривизны кривой так, чтобы выполнялось условие Ф (х2) Ф (хх) « 5 0 % .

12

Таким образом, мы видим, что если случайная величина распределена по закону (9), то при я = кр параметры (а, р, я) могут быть вычислены по способу наименьших квадратов по двум или большему числу точек интерполирования кумулятив­ ной кривой распределения с использованием табличных значе­ ний функции (10). При к = 1, т. е. при я = р из (1) и (10) как частный случай получаются известные в дисперсионном анализе уравнения Розина-Раммлера [14] или закон Вейбулла [9]

/ (л) =

аре~ахІ>хР -\ Ф(х) = 1— е~ а*Р.

(14)

Параметры а и р

функций распределения

(14) вычисляются

по формулам (12, 13), в которых значения zt

находятся в явном

виде из соотношения

 

 

 

zi = — ln [1 — Ф (Xi) ] = — \nq (Xj).

(15)

Статистические характеристики случайной величины, рас­

пределенной по закону (14), мода (х),

медиана (х\ ), математи-

_

2

ческое ожидание (х), среднеквадратическое отклонение (а) и коэффициент вариации (ул.) определяются по формулам:

 

 

 

,,т( +р}

(16)

 

х ~ а

р

і

І \

 

 

\ 1

1

 

а = а р

 

 

1

ъ =

о

(17)

7 ) - гТ +

 

>

 

У Т +

7 )' ъ - т

 

Практическая ценность применения уравнений (14) и их модификаций для выравнивания опытных кривых распределе­ ния большого диапазона случайных величин в настоящее время хорошо известна [1, 2, 9, 10, 11, 13, 14, 15, 16, 17, 18]. Примеры применения уравнений (1) и (9) в их более общем виде пока не­ многочисленны [1, 2]. Объясняется это, по-видимому, тем, что пока нет достаточно удобного в практическом отношении ал­ горитма вычисления параметров распределения (а, р, я) по ми­ нимальному числу опытных данных. Кроме того, если аналити­ ческое выражение уравнений (14) в настоящее время находит теоретическое обоснование в статистике экстремальных значе­ ний [9, 19 ], то уравнения (J) и (9) такого обоснования пока не имеют.

13

§ 2. Кривые Столетова намагничивания ферромагнетиков

По существующей теории процесса намагничивания ферро­ магнетиков [20] и статистической обработки большого экспери­ ментального материала [21, 22, 23] следует, что характер за­ висимости величины намагничивания / от напряженности на­ магничивающего поля Н может быть выражен уравнением вида (14):

 

 

• / = / Д

і - е~анР),

(18)

где Is

— магнитное насыщение; а и р — параметры, определяе­

мые из опыта по формулам (12,

13).

 

По известным функциональным соотношениям

[20, 24]

 

 

В = ptf = # + 4 я /,

(19)

где В — магнитная индукция;_р— магнитная

проницаемость,

и по уравнению

(18), находим:

 

 

 

 

В = Н + 4 я / Д і - е - анР),

(20)

 

 

|А = 1 + і ^

( і _ е-аЯ';)<

(21)

Из

уравнения

(21) и

равенства [20, 24]

 

 

 

р

= 1 +

4л%

(22)

получается явное аналитическое выражение магнитной воспри­ имчивости ферромагнитного вещества:

1 = Ь ( \ - е - * нР).

(23)

Исследуя уравнения (21) и (23) на экстремум, находим, что магнитнаяпроницаемость р и магнитная восприимчивость % достигают максимума для значений Н, определяемых из уравне­ ния

е~анР(\ — арНР) = 1.

(24)

В случае нулевого решения уравнения (24) из (21) и (23) пре­ дельным переходом имеем:

1

1 при р >

1

 

I

0

при р >

1

(25)

1 4- 4л/. при р =

1

X = { a/s

при

Р =

1

1-j-oo при

р <

1

^

оо

при

р <

1.

 

14

Разлагая экспоненциальный множитель (24) в ряд и

ограни­

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

приближенное

не нулевое решение уравнения (24)

 

 

 

Н* = и ха р ^ /24Р~ ‘5 ~ 3-

 

(26)

Геометрическая интерпретация уравнёний (18, 20,\ 21, 23)

показывает, что

кривая намагничивания (18) при р >

1 имеет

S -образную форму с одной точкой перегиба,

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

приближается к предельному горизонтальному положению I —

= I/, при р <

1 кривая намагничивания не имеет точек пере­

гиба, но для больших значений Н также принимает горизонталь­

ное положение. Кривая магнитной индукции (20)

при р > 1 и

р <

1 имеет аналогичную форму с кривой намагничивания, но

не

имеет

горизонтальной части, при больших

Я

вырождается,

в прямую

В = Я + 4 лIs. Кривые

магнитной

проницаемости

(21)

и магнитной восприимчивости

(23) при р >

1 характери­

зуются тем, что (.1 и %сначала сильно возрастают с увеличением напряженности поля Н, а затем, достигнув максимума при Я = == Н\а — Н% (26), начинают падать. При этом для больших зна­ чений Н величина р, стремится к единице, а величина %— к ну­

лю. При р = 1 имеем случай, когда р и %при Я =

0 принима­

ют максимальные значения р = 1 + 4nls, %= als и

сразу же

начинают падать, монотонно приближаясь к своим предельным значениям 1 и 0. При р < 1 величины р и %обладают бесконеч­ ным скачком, но с увеличением Я начинают падать, как и в пре­ дыдущих случаях. Важно заметить, что приведеннаія геометриче­ ская характеристика уравнений (18,- 20, 21, 23) находится в соответствии с характеристикой экспериментальных кривых Столетова для ряда образцов ферромагнетиков [21 ], представ­ ленных табл. 1, 2 и рис. 1—5.

Из табл. 1 видна хорошая согласованность расчетных (Ір) и опытных (/0) определений намагничивания сплава Гейслера [21 ]. Среднеабсолютная погрешность взаимного отклонения расчетных и опытных определений составляют около 1 %; мак­ симальная погрешность по отдельным измерениям менее 3%. Данные табл. 2 показывают на принципиальную зависимость характеристики намагничивания сплава Гейснера от времени старения и других факторов.

В табл. 3 и 4 приведены результаты расчета по формулам (14—17 при X = і) и опытных [22, 23] измерений остаточной на­ магниченности активных масс железного электрода щелочных

15

Рис. 1. Кривые намагничивания сплава Гейслера, построенные по уравнению (18)

6 9 Н-ЮЛктеі

С и старения

 

600°

 

от

 

закалки

 

после

 

1,0% А1)

2100 часов

Си,

и

М„6%

2, 6, 50

Гейслера (1,4%

209° С в течение

сплава

при

намагничивания

 

Кинетика

 

7

о

O O r f O c O ( N - - H O O

см —

•'^‘ СМ—' O C M C O l - П С М О

—-OOOOOOOO

О О О О О О О О Ю

со " " " О С О С О О Ю С М

“ O * - - * * - *

CD^OTFCOCOCMCMCM

o o o o o o o o o coowN-HTfiocom

C O O W C O l O t D S O i l N

 

 

 

 

X

 

 

 

 

0)

 

 

 

 

о .

и

 

 

 

са

 

 

 

н

 

o o o o o o o o o

и

 

coocsN ^’tmcoin

 

 

М ^ О О - « - н - н М

 

 

o o o o o o o o o

 

 

C O O O N O C O L O С О Ю

 

 

N O

l O O " «

i CM

 

о

O^NKjnW -OO

 

 

Tt* —<

 

 

 

CO

 

 

 

 

inNcocoooiMom

 

 

^ ( M - . - H - O O O O

 

 

o o o o o o o o o

 

 

O O i O O O O C O W O C O

 

 

CD

О О

О CD

 

 

со

со см см см - - — —

 

оа

o o o o o o o o o

 

C O ^ ( D N 0 5 0 -

COCD

 

 

СО "Т СО 00 —• т г Ю Ю СО

 

Н

 

 

 

 

и

o o o o o o o o o

 

 

о о т ^ с о с о — ^ ю с о ю

 

 

WCO'ifVlOiniOlßlD

 

 

o o o o o o o o

 

 

S 'Ч'ЮСО—'ГЮОЮ

 

 

СМСОтГтГЮЮЮСОСО

 

 

o o o o o o o o o

 

 

L O

O O O O O O O O

 

 

^

(N СО'СГ Ю И О СО О

 

o o o o o o o o o

O O O O O ^ C S C O l O С О О О О - н ^ О С О С Ч О )

СМСОСО^трт^сОСОСМ

o o o o o o o o o

C O O C O C O N W O I I M I O

<- i C O N ( N 0 O ( N ( O O 5

~« CM CM CM CM

o o o o o o o o o

C O O inOMNiOCOOOO)D C O M N n t N l

W

o o o o o o o o o

O O O i ß i n N O N N l O CS I C 05 W Ю CD 00 05

CMCOOOCMOOCOT PCM —

CO CM —1—■<

Ю 00 CM OO тр —«CDrJ*

ютгсососмсмсм-ч—< o o o o o o o o o

О О Ю П Ю О С О О Ю

TTNSCOTTOIO-N

N-CDtOlO^TPCOCOCM

o o o o o o o o o

NNlO^OOOOONlO COCO—lßNO —

— CM CM CM CM

o o o o o o o o o

WNin-COOOONin

C O l O O C M C O l O l O C O N -

o o o o o o o o o

cCOlßOJWCO^WCDNo o o - " - « t o o o i n i m o

o o o o o o o o o

Ю О О О О О О О О

ен (N со т г 1Л (D со О

Обозначения: / г — опытные данные намагничивания (гаусс); / р — расчетные данные намагничивания; W — «интенсивность» намагничивания определяется по уравнению W = f (14).

18

Время

старения,

час

2

6

50

2100

Т а б л и ц а 2

Статистическая

характеристика

намагничивания

 

 

 

сплава Гейслера

(табл.

1)

 

 

 

Параметры

Характеристика намагничивания

 

Р

а • ІО-6

Нэ

" .

 

Яор

О

Ун

/

 

 

 

2

 

 

 

 

 

0,48

836

0

85

 

378

894

2,40

23

0,39

21100

0

23

 

191

618

3,23

50

1,00

390

0

177

 

266

256

1,00

4,5

1,61

6,74

211

307

 

347

221

0,64

2,5

Обозначения: Нэ — мода кривых распределения;

 

 

 

 

— медиана

«

«

;

 

 

 

Г

Нср — математическое ожидание;

а— среднеквадратическое отклоненію;

у„

— коэффициент вариации;

/

— показатель неоднородности.

аккумуляторов в зависимости от силы тока (і) в обмотке электро­ магнита. Среднеабсолютная погрешность расчетных и опытных определений (табл. 3) по всем испытанным образцам меньше 2%; -максимальная погрешность по отдельным замерам не превыша­ ет 5%. Из табл. 4 видно, что статистическая характеристика остаточной намагниченности активных электродных масс зави­ сит не только от природы компонентов исходного вещества, но и от технологии их получения [22—28].

§ 3. Коалесценция капель эмульсий на границе раздела двух жидких фаз

Кинетику коалесценции капель эмульсий на границе разде­ ла двух жидких фаз [29, 30] можно представить аналитическим выражением вида

^

N = Wo (1 -

(27)

где Na — число всех капель; N — число капель, коалесценировавших за время т > т0; а и р — параметры, зависящие от условий коалесценции и от природы вещества.

Уравнение (27). получается из статистического закона тео­ рии «наиболее слабого звена» [19], примененного к распреде­ лению «времени жизни» капель у плоской поверхности раздела двух жидких фаз. Схематически это можно представить следую-

19

Соседние файлы в папке книги из ГПНТБ