УДК: 637.1:65.011.56
Профессор А.А. Хвостов, доцент Д.И. Ребриков, аспирант В.Е. Мерзликин
(Воронеж. гос. ун-т. инж. технол.) кафедра информационных и управляющих систем.
тел. (473) 255-38-75
E-mail: rebrikov-vgta@vgta.vrn.ru
Professor A.A. Khvostov, associate Professor D.I. Rebrikov, graduate V.E. Merzlikin
(Voronezh state university of engineering technologies) Department of information and control systems. phone (473) 255-38-75 E-mail: rebrikov-vgta @vgta.vrn.ru
Идентификация массовых распределений жировой фазы в молоке с помощью универсальных распределений Пирсона
Identification of the mass distributions of the milk fat phase using universal Pearson distributions
Реферат. В работе рассмотрена задача по аппроксимации экспериментальных значений коэффициента затухания ультразвуковых колебаний и массовых распределений жировых шариков в молоке и молочных продуктах. Произведен анализ экспериментальных данных с точки зрения выбора метода аппроксимации. Предложена аппроксимирующая зависимость, которая основана на решении дифференциального уравнения Пирсона. Рассмотрены преимущества предложенного метода аппроксимации с учетом вида получаемых экспериментальных данных. Разработан алгоритм построения математической модели, описывающей спектр времен релаксации и массовых распределений жировых шариков в молоке и молочных продуктах. В результате аппроксимации семейством кривых Пирсона экспериментальных данных показана возможность качественно верно описать изменение распределения жировой фазы в ходе процесса гомогенизации. Оценена погрешность аппроксимирующей зависимости, которая составила 18 %. Показано, что с ходом процесса гомогенизации молочных продуктов изменяется вид кривой, описывающей распределение жировых шариков, в виду того, что появляется локальный экстремум, вызванный наличием части негомогенизированных жировых шариков. При этом точность выбранной математической модели резко снижается. Вместе с этим, теряется физический смысл и её параметров. Для устранения выявленных отклонений в работе предложено массовое распределение жировых шариков рассматривать как функцию с двумя модами. Обосновано, что усложнение модели не только увеличивает в два раза количество её параметров, но и усложняет интерпретацию результатов измерений в системе управления и делает сложным анализ получаемых параметров аппроксимации лицом, принимающим решение. По результатам аппроксимации экспериментальных данных и по их виду предложено использовать статистические моменты распределения для решения поставленной задачи.
Summary. In this paper we consider the problem of approximating the experimental values of the coefficient of attenuation of ultrasonic oscillations and the mass distribution of the fat globules in the milk and milk products. The analysis of experimental data in terms of the choice of the method of approximation was done. A approximating dependence is based on the solution of Pearson differential equations. The advantages of the proposed method for the type of approximation of the experimental data obtained. An algorithm for constructing a mathematical model describing the relaxation spectrum and mass distribution of the fat globules in the milk and milk products was implemented. As a result, a family of Pearson approximation curves of the experimental data shows the ability to qualitatively correctly describe the change in the distribution of the fat phase in the process of homogenization. It estimates the error of approximating dependence, which amounted to 18 %. It is shown that during of the process of homogenization of dairy products changes shape of the curve describing the distribution of the fat globules, in view of the fact that there is a local extremum, caused by the presence of the non-homogenized fat globules. The accuracy of the selected mathematical model is significantly reduced. At the same time, it loses its physical meaning and its parameters. To address the identified deviations in the proposed mass distribution of fat globules as a function with two modes. It is proved that the complexity of the model is not only doubles the number of its parameters, but also complicates the interpretation of measurement results in a control system, and makes it difficult to analyze the obtained parameters of approximation by decision-maker. As a result of approximation of experimental data suggested to use statistical moments of the distribution for problem decision.
Ключевые слова: молочные продукты, акустические свойства, массовое распределение, аппроксимация, распределение Пирсона.
Keywords: dairy products, acoustic properties, distribution of mass, approximation, Pearson distributions.
© Хвостов А.А., Ребриков Д.И., Мерзликин В.Е., 2015
Использование связи значений масс распределения жировых шариков по массовым или объемным фракциям [1, 2] со спектром времен релаксации молока от акустических свойств раствора [3] при параметрической идентификации представляет трудности по следующим причинам:
- число точек Н(т) ограничено количеством их измерений;
- разная разрешающая способность по массе и частоте усложняет сопоставление точек, полученных экспериментально, на кривых релаксационного спектра и распределения масс жировых шариков по массовым или объемным фракциям с релаксационным спектром молока;
- вероятное увеличение случайной ошибки без фильтрации исходных данных.
То есть необходима аппроксимация экспериментальных значений Н(т) и распределения жировой фазы. Рассмотрим синтез аппроксимирующей зависимости на примере массового (объемного) распределения жировой фазы в молоке. Поскольку полная функция распределения начинается с нулевого значения и заканчивается нулевым значением, проходит через один или несколько максимумов и при этом вид спектра заранее неизвестен, то для его аппроксимации можно использовать решения дифференциального уравнения Пирсона [4, 5]:
dH (т) т + а
dт
т•Н(т), (1) Ь + Ьт + Ьт т
-'о 1 'Л^ 1 2 где а, Ьо, Ьи, Ь2, — постоянные.
Общее решение этого уравнения может быть представлено в виде:
-I
(г+а) dт Ьо+Ьт+Ыт2
Н (т)= Н0е Ьо+Ьт+Ь2т (2)
Параметры а, Ьо, Ьи, Ь2, определяются по методу моментов:
а = ¡и3 (ц4 + З^2) / А,
Ь0 = -Н (4 Н2Н4 - 3н3, ) / A, < Ь = (ц4 + 3$. ) / А,
Ь2 = - (- 3н3 - 6М2 ) / A,
А = 10н2М4 -18^2 -12 Н
(3)
где Ц2, Из, Н-4 - первые четыре центральных момента распределения измеряемой величины:
н = Е (т -н)к • н т):
(4)
1=1
где к - порядок рассматриваемого момента; п - количество экспериментальных точек.
Структура решения дифференциального уравнения (1) зависит от полученных корней ¥\, ¥2 уравнения [5]:
Ь0 + Ьтт + Ь2т2 = 0 (5)
Для проверки корней данного уравнения вводится параметр:
ь2
ж= -
4ЬоЬ2
(6)
Величина ж называется критерием Пирсона (каппа Пирсона) и ее значения определяют следующие свойства корней уравнения:
- Если ж < 0, уравнение (2) имеет действительные корни разных знаков.
- Если 0< ж<1, уравнение (2) имеет комплексные корни.
- Если ж>1, уравнение (2) имеет действительные корни одного знака.
Всем этим вариантам Пирсон ставит в соответствие три основных типа своих кривых, которые он назвал кривыми I, IV и VI типа. Также ж может принимать значения 0,1 ±да, что дает переходные типы кривых. Всего семейство кривых Пирсона включает 12 типов, плюс нормальную кривую.
Имея рассчитанный параметр ж, можно определить вид кривой, аппроксимирующей экспериментальные данные.
В случае нескольких основных времен релаксации согласно принципу суперпозиции релаксационных процессов в молоке суммарная функция ^ (Т,с) будет аппроксимироваться взвешенной суммой функций плотности в двумерном виде:
IV
Г (Т ,) = Е£/1 (Т с),
(7)
где 1 - номер рассматриваемого релаксационного процесса, N - количество рассматриваемых релаксационных процессов, ^ - весовой коэффициент для 1-го релаксационного процесса в общие механические потери.
Таким образом, алгоритм, с помощью которого осуществляется построение математической модели, которая описывает спектр времен релаксации в молоке, будет выглядеть следующим образом (рисунок 1):
1. Разделение множества всех экспериментально точек в плоскости (Т,с) на подмножества, соответствующие отдельным выборкам и относящиеся к конкретному релаксационному механизму.
Начало
/ гг. /
/ ^ /
Разбиение множества всех наблюдаемых экспериментально точек в плоскости на подмножества исходя из принципа унимодальности (т)
/ ' \ j=l
j<m
Ввод исходных данных £эксп (экспериментальные значения распределения жировой фазы или релаксационного спектра) при заданных аргументах Раг<1 (частота, время релаксации или масса жирового шарика)
Для каждого релаксационного механизма
Определение матрицы моментов Мш, Мт j-го распределения , расчет по ним значений параметров уравнения Сш, Сти нахождение корней и
Определение типа распределения по значению корней ^ и Ш2
1*^2<=С Ш1 и ш *ш >0 11 1комплексные
Распределение Пирсона I типа Распределение Пирсона VI типа Распределение Пирсона IV типа
* 1
j=j+1
I
fl, (Par), в Л Par) YiPar), Г}(Раг\
^аал 1 ^у ^¡яап.жп
Уточнение параметров суммарного распределения методом сопряженных градиентов по критерию МНК
Аппроксимация параметров распределения. Построение модели f(Par)
Определение точности и адекватности модели
Представление результатов аппроксимации
Рисунок 1. Алгоритм построения математической модели
2. Создание из множества экспериментальных точек, соответствующих каждому механизму релаксации, подмножеств экспериментальных точек, связанных с фиксированным частотам или температурам, для создания выборок условных распределений по частотам и температурам.
3. Приведение к нормированным значениям для расчета выборочных моментов условных распределений.
4. Расчет выборочных моментов
Мш, МТ условных распределений.
5. Расчет коэффициентов квадратного уравнения в знаменателе (1) Сш, СТ и вычисление его корней \ус, хуТ
6. Выбор типа распределения из семейства универсальных распределений Пирсона.
7. Расчет параметров распределения
Г с, Ут VI .
8. Дополнительное уточнение коэффициентов суммарного многомодального распределения с помощью метода покоординатного спуска по среднеквадратичному критерию и начальными приближениями, полученными на этапе 7.
9. Оценка адекватности и точности модели.
По проведенным экспериментальным
исследованиям установлено, что полученные корни уравнения (2) - вещественные числа разных знаков. Поэтому для моделирования спектра времен релаксации необходимо использовать бета-распределение первого рода:
Хт/1 • (1 - у(т)У2
Н (т) = :
(8)
к • В^,
где s1, s2 - параметры рассматриваемого распределения; у(т) - аргумент; B(s1,s2) - значение полученной бета-функции.
Параметры для распределения вычисляются по формулам:
я, =
1 В2 (С1 + С2 )'
(9)
В2 (С1 + С2 )
Аргумент распределения можно определить по формуле:
т-01
у(т) = ■
(10)
где вх = н - с1, д2 = с1 + с2 .
Для того чтобы идентифицировать параметры модели релаксационного спектра был использован среднеквадратичный критерий:
£ = 1 (Н(т, )эксп - Н(т)расч )2
(11)
При проведении экспериментальных исследований использовано 42 образца молока с различным массовым распределением жировой фазы вследствие проведения процедуры гомогенизации при разных давлениях.
Необходимо заметить, что в процессе параметрической идентификации использованы предварительно нормированные данные эксперимента. Переход от абсолютных значений Н(т) к нормированным вызван тем, что площадь под кривой распределения жировой фазы должна
тггах
быть равна 1, однако условие | Н(т^т = 1 не
ттт
выполняется, поэтому для дальнейшего расчета массового распределения жировой фазы по спектру времен релаксации необходимо осуществить нормировку исходных данных.
Для чего необходимо соблюдение следующего условия:
1\
I Н (т, ) = 1,
(12)
"норм N
где N - количество точек измерения.
Для чего необходимо осуществить введение дополнительного нормирующего коэффициента "норм, который обеспечивает выполнение равенства (4):
-— (13)
IН (т)
1=1
При дальнейших расчетах используемые экспериментальные значения спектра времен релаксации подразумеваются как нормированные значения.
На рисунке 2 представлены экспериментальные значения распределений жировой фазы молока до процесса гомогенизации и при разных давлениях проведения процесса. Приведенные данные показывают увеличение содержания низкоразмерных фракций с увеличением давления гомогенизации.
Аппроксимация семейством кривых Пирсона с использованием зависимости структуры (3) позволила качественно верно описать изменение распределения жировой фазы в ходе процесса гомогенизации (рисунки 3, 4, 5). Из представленных графиков видно, что в начале процесса гомогенизации вид распределения имеет ярко выраженный экстремум, благодаря чему распределение можно считать унимодальным. В этом случае математическая модель на основе кривых Пирсона с приемлемой погрешностью (порядка 18%) описывает распределение жировой фазы.
,52 ,к ,в1,02
->тт
^ =
i 1)
1x10
В 3 п 1x103 о
hn4
h90
;i) I
m
hi 801
н о
10 0
10
rt- 13
1x10 12 1.5x10 12
0 5x10 Площадь жировых ш ариков, мкм2
--негомогенизированное молоко,
- . - . - - гомогенизированное при давлении 90 и
180 МПа, соответственно
Рисунок 2. Сравнительные распределения жировых
шариков в объеме, пробы по диаметру,, ,, ^ апоГепетоеу еоеаю 1е6ппа
Л
н о о o
Ё
0 а
1
0.2
0.15
0
д
х\ < -V v v * v .. v w .
0 10 20 Диаметр жировых шариков, мкм
30
UVUIUUW. IV1
вая аппроксимации,
nerteiaoey
Экспериментальные данные
л туетаоеию
Рисунок 3. Аппроксимация экспериментальных данных распределения жировых шариков негомо-
генизированного молока ДУ Прсо®— -
аио1епе1аоеу еоеаге 1еоппа
0 10 20 Диаметр жировых шариков, мкм кривая аппроксимации,
..ЙЖршентальные данные п л дУегаабёшб
Рисунок 4 Аппроксимация экспериментальных данных распределения жировых шариков гомогенизированного молока при давлении 90 МПа ДУ Пирсона.
aii6ieneiaoey eoeaie iesniia
10
20
Диаметр жировых шариков, мкм aeaiaoo, iei
aiMray1™^™™,
экспериментальные данные _ „ yeniaOeiaiO
Рисунок 5/Аппроксимация экспериментальных данных распределения жировых шариков гомогенизированного молока при давлении 180 МПа ДУ Пирсона
Однако в ходе процесса гомогенизации вид кривой распределения изменяется, и экстремум смещается к левой границе области определения функции распределения фазы по размерам. Кроме этого, ввиду наличия как гомогенизированной, так и негомогенизирован-ной части в общем объеме начинает проявляться остаток негомогенизированного молока в виде небольшого локального экстремума. В этом случае для того, чтобы достигнуть приемлемого уровня отклонений модельных данных от экспериментальных в математическую модель приходится вводить вторую моду. Такое усложнение модели не только увеличивает в два раза количество её параметров, но и усложняет интерпретацию результатов измерений в системе управления и ЛПР.
Проведенные исследования показали возможность качественного описания изменения объемного (массового) распределения жировой фазы в молочных продуктах с использованием семейства распределений Пирсона. Для более точного описания необходимо вводить в математическую модель описание второй моды или использовать описание функции распределения статистическими моментами.
30
1
0
ВестникВГУИТ, №2, 205
ЛИТЕРАТУРА
1 Крусь Г.Н. и др. Технология молока и молочных продуктов. М.: Колосс, 2003. 315 с.
2 Kress-Rogers E., Brimelow C. J. B. Instrumentation and sensors for the food industry. Second edition. Abington: Woodhead Publishing Limited, 2000.
3 Honerkamp J. et al. Determination of the relaxation time spectrum from dynamic module using an edge preserving regularization method // Rheologica Acta. 2000. V. 39. № 2. Р. 163-173.
4 Крамер Г.Математические методы статистики. М.: Мир, 1975. 658 c.
5 Кендалл М., Стьюарт А. Теория распределений. М.: Наука, 1966. 588 c.
REFERENCES
1 Krus' G.N. et al. Tekhnologiya moloka i molochnykh produktov [Technology of milk and dairy products]. Moscow, Koloss Publ, 2003, 315 p. (In Russ.).
2 Kress-Rogers E., Brimelow C. J. B. Instrumentation and sensors for the food industry. Second edition. Abington: Woodhead Publishing Limited, 2000.
3 Honerkamp J. et al. Determination of the relaxation time spectrum from dynamic module using an edge preserving regularization method. Rheologica Acta, 2000, vol. 39, no. 2, pp. 163-173.
4 Kramer G.Matematicheskie metody statistiki [Mathematical methods of statistics]. Moscow, Mir, 1975, 658 p. (In Russ.).
5 Kendall M., Stuart A. Teoriya raspredele-nii [Theory of distributions]. Moscow, Nauka, 1966, 588 p. (In Russ.).