УДК 519.876.2:658.5
ОПТИМАЛЬНОЕ УПРАВЛЕНИЕ ПРИ ДОСТИЖЕНИИ ЦЕЛЕЙ В ОБЛАСТИ КАЧЕСТВА ПРОМЫШЛЕННОГО ПРЕДПРИЯТИЯ
© 2021 А.Г. Ивахненко, О.В. Аникеева
Юго-Западный государственный университет, г. Курск, Россия
Статья поступила в редакцию 24.07.2021
Назначение обоснованных значений целей в области качества отнесено к одной из основных функций на тактическом уровне управления. Рассмотрена математическая модель при целевом управлении в области качества, являющаяся системой обыкновенных дифференциальных уравнений с постоянными коэффициентами. Дана характеристика параметров классического квадратичного функционала, соответствующего функции потерь качества Тагути. Эти параметры определяют затраты, связанные с отклонениями каждой цели в области качества на конец планового периода и в течение него, а также затраты связанные с отклонениями величин управляющих воздействий от статических значений, для достижения поставленных целей в переходных процессах. Приведены решения этой системы для наиболее распространенных на практике ступенчатого и линейного законов управления. Наличие такого решения позволило применить стандартный метод поиска экстремумов функции многих переменных при оптимизации без использования вариационного исчисления. Рассмотрен пример поиска минимума составляющих и самого классического квадратичного функционала, с использованием фактических данных по результатам деятельности промышленного предприятия при различных значениях параметров.
Ключевые слова: оптимальное управление, цели в области качества, математическая модель, функционал.
Б01: 10.37313/1990-5378-2021-23-4-18-26
Исследование выполнено при финансовой поддержке РФФИ в рамках научного проекта № 19-01-00015.
ВВЕДЕНИЕ И ПОСТАНОВКА ЗАДАЧИ
Назначение обоснованных значений целей в области качества является одной из основных функций на тактическом уровне управления качеством [1]. Наименования целей и их значений определяют содержание документа «Цели в области качества», как правило, со сроком выполнения на один год. Вопросы, связанные с достижимостью поставленных целей при наиболее распространенных ступенчатом и линейном законах управления, рассмотрены в работах [2, 3]. В данной работе рассмотрен характер поведения составляющих функционала, предложенного при постановке общей задачи оптимизации для линейной модели системы менеджмента качества (СМК) [4] с акцентами на безусловную достижимость поставленных целей и/или минимизацией затрат.
Задачей данной работы является исследование классического квадратичного функционала и
Ивахненко Александр Геннадьевич, доктор технических наук, профессор, профессор кафедры машиностроительных технологий и оборудования. E-mail: [email protected]
Аникеева Олеся Владимировна, кандидат технических наук, доцент, доцент кафедры дизайна и индустрии моды. E-mail: [email protected]
его составляющих при ступенчатом и линейном законах управления для линейной модели СМК.
ПОСТАНОВКА ОБЩЕЙ ЗАДАЧИ ОПТИМИЗАЦИИ ДЛЯ ЛИНЕЙНОЙ МОДЕЛИ СМК
Сама математическая модель при целевом управлении СМК имеет вид:
х (?) = ах (?) + ви (/% (1)
где составляющие вектора (переменные состояния) X = (Х(1), Х(2))т являются текущими значениями целей в области качества Х(1) и скоростями их изменения Х(2), каждый из которых имеет размерность п; п - количество целей; А - системная матрица (2п*2п); В - матрица параметров управления; иф - вектор управляющих воздействий.
Система уравнений (1) совместно с заданными начальными условиями для переменных состояния Х(0) представляет собой задачу Коши. Решение этой системы для произвольного управляющего воздействия имеет следующий вид:
t
x^) = ехр^ X(0) +1 ехрА ви(т)ёт, (2)
о
где выражение ехрА - матричная экспонента, вычисление которой связано с преобразованием матрицы А, для представления ее в виде произведения жордановых матриц, или же разложением в ряд по степеням данной матрицы [6].
Рассмотрим второе слагаемое в правой части системы (1). В практике управления для достижения значений поставленных целей в области качества наибольшее распространение получили ступенчатый и линейный законы управления. Ступенчатый закон управления отражается в указанном выше документе только требуемыми значениями целей на конец планируемого периода [Х(1)], при этом явно не учитываются скорости их изменения. Линейный закон управления отражается в документе не только конечными значениями [Х(1)], но и их промежуточными значениями в течение планируемого периода, то есть неявно задает скорости изменения значений поставленных целей. Составляющие матрицы В для этих законов управления определяются в два этапа. На первом этапе рассматривается решение уравнения (1) после завершения переходных процессов в СМК, то есть при выполнении условия dX/dt = 0. Определенные таким образом компоненты матрицы В обеспечивают потенциальную достижимость целей в области качества. На втором этапе в матрицу В вводятся коэффициенты усиления для обеспечения фактической достижимости требуемых значений целей в течение планового периода времени, то есть задаются более высокие значения целей в области качества по сравнению со статической задачей.
Тогда математическая модель целевого управления СМК из (1) преобразуется к виду:
x (г) = ах {г) + в *и(г), (3)
где В - матрица, включающая элементы матрицы В и коэффициенты усиления к, г = 1...п.
При ступенчатом законе иф = 1(Е) (где 1(Е) = 0 при Ь < 0, 1(Е) = 1 при t > 0) правую часть выражения (3) представим так:
В'и (г) = (0, В )Т1(г), (4)
где 0 - нулевой вектор размерностью п, В1 - вектор размерностью п, компоненты которого содержат введенные коэффициенты усиления кг.
При линейном законе управления и^) = (1(Е), t)т, правую часть (3) представим таким образом: В 'и(г) = (0, В1 )Т1(г) + (0, В2 )Т г, (5)
где В2 - вектор размерностью п, компоненты которого также содержат введенные коэффициенты усиления кг.
Необходимость введения коэффициентов усиления обоснована тем, что действующая СМК не находится в стабильном состоянии, наоборот, ей присущи постоянные переходные процессы обусловленные изменением целей в
области качества. В связи с этим, для фактического достижения требуемых значений целевых показателей необходимо задавать в планах их повышенные значения. Отметим, что действующие СМК являются устойчивыми, то есть собственные значения матрицы A, входящей в выражения (1) и (3), имеют отрицательные действительные части.
Хотя векторы Bj различаются в выражениях (4) и (5), однако схожая структура этих выражений позволяет представить систему (2) в следующем виде:
X (t) = exp at( x (0) + ллв! + лгвг) -
- ллв1 - л'2вг - A"B2t . (6)
Решения (6) системы уравнений (3) для каждой составляющей вектора X являются функциями не только времени t, но и введенных коэффициентов усиления X = X(t, k1, k2, ..., kn). Это позволяет использовать стандартный метод поиска экстремумов функции многих переменных для нахождения оптимальных решений.
В работе [4] при постановке задачи оптимизации с использованием функции потерь качества Тагути [5] было предложено определять минимум функционала I (I ^ min) I = Vt + L1 + L2 следующего вида:
1 = £P^D^T) -[x(1)l])2 + i=1
+JZ ßiOWe) - [X(1);I])2d(e)+, (7) 0 i=1
l 2n
+JX KI(u(e)l - ustl)2d(e)
о i=n+1
где р. и р. - параметры, определяющие затраты, связанные с отклонениями каждой цели в области качества на конец планового периода, и в течение него, соответственно (г = 1 ... п); К - параметры, определяющие затраты, связанные с отклонениями величин управляющих воздействий от статических управляющих воздействий, для достижения поставленных целей в переходных процессах (г = п + 1 ... 2п); им - вектор управляющих воздействий, определенный при решении задач статики и квазистатики качества из (1) из условия достижения значений поставленных целей [Х(1)] при dX(t)/dt = 0.
ИССЛЕДОВАНИЕ КЛАССИЧЕСКОГО КВАДРАТИЧНОГО ФУНКЦИОНАЛА И ЕГО СОСТАВЛЯЮЩИХ ПРИ СТУПЕНЧАТОМ ЗАКОНЕ УПРАВЛЕНИЯ
Рассмотрим на примере характер поведения этого функционала, являющегося суммой составляющих, определенных выражением (7). При анализе фактических результатов деятельности
предприятия ЗАО «Салют» [7] в работе [2] были выделены две цели в области качества (п = 2).
Начальные значения в 2017 году составили: для целей в области качества предприятия х1(0) = 0.206, х2(0) = 0.90; для скоростей их достижения х3(0) = 0.008, х4(0) = 0.02. Фактические значения целей в области качества, которые были достигнуты к концу 2017 года: [х1] = 0.233; [х2] = 0.93. Плановый период времени Т = 1 год. Подход к определению составляющих матрицы А представлен в [3], а сама матрица имеет вид: ( 0 0 1 0 ^ 0 0 0 1 - 0.815 0.090 - 0.889 0.094 0.750 - 0.896 0.818 - 0.938
(11)
A =
(8)
При ступенчатом законе управления коэффициенты матрицы управляющих воздействий В, обеспечивающие фактическую достижимость поставленных целей х1 и х2 для выражения (4), определяются так:
b131 = 0.190^ - 0.084k2, b141 , = -0.175^ + 0.883k2, ,
(9)
где к и к2 - коэффициенты усиления управляющих воздействий для целей х1 и х2, соответственно.
Часть системы уравнений (6) для поставленных целей представляет собой выражения: х^) = [0.530F1•sin(0.873t + 1.021) -- 0.208F2•sin(0.873t + 0.977)+ +0.374•10-3•sin(0.873t +1.103)]е-а596г + (10) + [0.619^т(0.70к +1.088)+ ( ) +0.189£^т(0.70к +1.152) + +0.011*т(0.70к - 0.029)]е-03Ш - ^ + 0.206;
x2(t) = [-1.737F1•sin(0.873t + 0.979) +
+ 0.682F2•sin(0.873t + 0.935) -- 0.123•10-2•sin(0.873t + 1.061)] е-0-596 + +[1.575^т(0.70к + 1.155) + +0.481Е^т(0.70и + 1.219) +
+ 0.029*т(0.70к + 0.037)] е-03Ш - F2 + 0.900, где 0.206 - 0.233к1, F2= 0.9 - 0.930к2.
На рис. 1 представлено изменение целей х1 и х2 во времени при коэффициентах усиления управляющих воздействий к1 = к2 = 1 для ступенчатого закона управления.
Как видно из рис. 1, обе цели в течение планового периода времени не достигаются. Первый момент достижения плановых значений показателей соответствует t * 2.2 года для первой цели, и t * 1.9 года - для второй цели. Минимальные значения коэффициентов к1 = 1.25 и к2 = 1.1, обеспечивают достижение требуемых значений целей в области качества при t = 1 год, что связано не с выходом на установившийся режим, а с первым достижением значений поставленных целей при затухающих колебаниях, происходящих от начальных значений переменных х1 и х2 к их плановым значениям.
На рис. 2 представлены графики зависимостей всех частей функционала I = 1(кр к2) - при принятых значениях весовых коэффициентов затрат: р1 = р2 = Р1 = Р2 = К = К2 = 1, т.е. при условии равных удельных затрат на потери качества и управление, вызванные их отклонениями от номинальных значений.
Составляющие функционала I при приведенном выше равенстве единице значений весовых
Рис. 1. Результаты моделирования целей в области качества при ступенчатом законе управления (k1 = k2 = 1):
а) цель x1; б) цель x2 Fig. 1. The results of modeling quality goals with a step-by-step control law (kl = k2 = 1): a) goal x1; b) goal x x2
коэффициентов удельных затрат на управление имеют вид:
V = 0.613-Ю-2^2 +0.081к22 -- (0.027к2 - 0.013)к1 - 0.142к2 + 0.068, (12) Ь1 = 0.148-10-2к12 + 0.019к22 -- (0.663к2 - 0.337)10-2к1 - 0.034к2 + 0.018, (13) Ь2 = 0.054(1 - к1)2 + 0.865(1 - к2)2. (14) Минимальные значения терминальной части V функционала определялись из решения системы уравнений дУ1/дк1 = 0 и дУ/дк2 = 0, при соответствующих условиях определяемых вторыми производными, в области действительных чисел и условия к1 > 0 и к2 > 0 для р1, р2 е (0..2]. Полученные минимальные значения УШп, которым соответствуют значения коэффициентов к1 = 1.246 и к2 = 1.079 при достигнутых требуемых значениях целей [х1] = 0.233 и [х2] = 0.93, представлены в табл. 1.
Минимальные значения составляющей функционала I определялись из решения системы уравнений д!1/дк1 = 0 и д!1/дк2 = 0 в области действительных чисел и условия к1 > 0 и к2 > 0 для Р1, Р2 е (0..2]. Полученные минимальные зна-
Таблица 1. Минимальные значения составляющей Ушп функционала I при ступенчатом законе управления
P2=0.1 P2=0.5 P2=1 P2=1.5 P2=2
P1=0.1 1.96-10-21 9.61-10-21 1.121-10"20 7.764-10-20 1.289-10"20
Pi=0.5 10-21 2.04-10-20 1.55-10-20 3.36-10-20 2.36-10-20
Pi=1 4.4-10-21 2.04-10-20 1.96-10"20 2-10-20 1.49-10-20
Pi=1.5 10-21 8-10-22 10-20 2-10-20 2-10-20
pi=2 10-21 2-10-20 1.84-10-20 2-10-20 3.84-10-20
Таблица 2. Минимальные значения составляющей функционала I при ступенчатом законе управления
p2=0.1 p2=0.5 p2=1 p2=1.5 p2=2
Px=0.1 0.580-10-4 при k1 = 1.486 k2 = 1.171 1.807-10-4 при k = 1.422 k2 = 1.161 3.329-10-4 при k = 1.343 k2 = 1.148 4.841-10-4 при k1 = 1.266 k2 = 1.136 6.342-10-4 при k1 = 1.190 k2 = 1.124
pi=0.5 1.669-10-4 при k1 = 1.499 k2 = 1.171 2.9-10-4 при k1 = 1.486 k2 = 1.171 4.437-10-4 при k1 = 1.470 k2 = 1.168 5.971-10-4 при k = 1.454 k2 = 1.166 7.503-10-4 при k = 1.438 k2 = 1.163
Pi=1 3.030-10-4 при k = 1.500 k2 = 1.169 4.262-10-4 при k = 1.494 k2 = 1.172 5.8-10-4 при k1 = 1.486 k2 = 1.171 7.338-10-4 при k1 = 1.478 k2 = 1.170 8.874-10-4 при k1 = 1.470 k2 = 1.168
Pi=1-5 4.391-10-4 при k1 = 1.500 k2 = 1.167 5.623-10-4 при k1 = 1.497 k2 = 1.172 7.162-10-4 при k = 1.492 k2 = 1.171 8.701-10-4 при k = 1.486 k2 = 1.171 10.238-10-4 при k1 = 1.481 k2 = 1.170
pi=2 5.752-10-4 при k = 1.499 k2 = 1.165 6.984-10-4 при k1 = 1.498 k2 = 1.171 8.524-10-4 при k1 = 1.494 k2 = 1.172 10.063-10-4 при k1 = 1.490 k2 = 1.171 11.601-10-4 при k = 1.486 k2 = 1.171
Рис. 2. Графики зависимостей составляющих
функционала I при ступенчатом законе управления и условии равных удельных затрат на потери качества и управление:
Pi = Р2 = Pi = Р2 = К = К2 = 1
Fig. 2. Graphs of the dependencies of the components of functional I with a stepwise control law and the condition of equal unit costs for quality loss and management:
Pi = P2 = Pi = P2 = Ki = K2 = 1
чения которым соответствуют значения
коэффициентов к1 и к2 , представлены в табл. 2. При этом только при Р1 = 0.1 и Р2 = 2 на конец года первая цель не достигается: [х1] = 0.228. В остальных случаях 0.233 < [х1] < 0.248. Вторая цель достижима при всех рассматриваемых значениях р1 и Р2: 0.944 < [х2] < 0.946.
Минимальные значения составляющей Ь2 функционала I определялись из решения системы уравнений д!2/дк1 = 0 и д!2/дк2 = 0 в области действительных чисел и условия к1 > 0 и к2 > 0 для К1, К2 е (0..2]. Полученные минимальные значения Ь^^ совпадают при всех К1, К2 е (0..2], и равны нулю. При этом им соответствуют одинаковые значения коэффициентов к1 = к2 = 1, т.е. обе цели на конец года недостижимы: для всех случаев [х1] = 0.218, [х2] = 0.918.
Сам функционал 1(к1,к2) при принятых значениях весовых коэффициентов удельных затрат: р1 = р2 = Р1 = Р2 = К1 = К2 = 1, имеет вид: 1 = 0.062кк12 +0.965к22 -- (0.092 + 0.033к2)к1 - 1.906к2 + 1.006. (15)
Минимальное значение функционала 1тЬ = 0.126-10-2 соответствует значениям коэффициентов усиления к1 = 1.011, к2 = 1.005, при которых обе цели в области качества не достижимы на конец года.
В определенной ранее области достижимости целей (ОДЦ) для значений коэффициентов к1 = 1.25 и к2 = 1.1, минимальное значение функционала составляет I . = 0.127-10-1.
^ тт
ИССЛЕДОВАНИЕ КЛАССИЧЕСКОГО КВАДРАТИЧНОГО ФУНКЦИОНАЛА И ЕГО СОСТАВЛЯЮЩИХ ПРИ ЛИНЕЙНОМ ЗАКОНЕ УПРАВЛЕНИЯ
При линейном законе управления коэффициенты матрицы управляющих воздействий В, обеспечивающие фактическую достижимость поставленных целей х1 и х2 для выражения (5) определяются так:
Ь2 31 = 0.190к1 - 0.084к2 - 0.087, Ь2^ = -0.175^ + 0.883к2 - 0.652, (16) а система (3) при линейном законе управления имеет следующий вид:
Г 1 ( 0 0 1 0 ^ Г 211 ( 0 1 Г о ^
х2 0 0 0 1 х2 + 0 1(0 + 0
хъ — - 0.815 0.090 - 0.889 0.094 х3 0.087 Ь2,31
V 24 0.750 - 0.896 0.818 - 0.938 V Х4 0.652 1Д41)
(17)
Соответствующая часть системы (6) для рассматриваемых целей представляет собой: х^) = [0.530•Я2•sin(0.873t + 1.021) -- 0.208•Я2•sin(0.873t + 0.977) + + 0.509•Я3•sin(0.873t + 0.070) -0.195•Я4•sin(0.873t - 0.013)] е-0-596 + + [0.619-Я^т(0.70и + 1.088) + + 0.189•Я2•sin(0.701t + 1.152) + (18)
+ 0.808-Я^т(0.70и - 0.044) + +0.243-Я4^т(0.70и - 0.010)] е-0318' -
- '-(Я3 - 0.008) - Я1 + 0.206;
х2(') = [-1.737Я^т(0.873' + 0.979) + +0.682-Я2^т(0.873' + 0.935) -- 1.668Я3^т(0.873' + 0.028) + +0.638•Я4•sin(0.873t - 0.055)] е-0596' +
+ [1.575-Я^т(0.70и + 1.155) + (19) + 0.481•Я2•sin(0.701t + 1.219) + + 2.058Я^т(0.70П + 0.022) + + 0.620Я^т(0.70П + 0.056)] е-0318' -
- '-(Я4 - 0.02) - Я2 + 0.900,
где Я1 = -0.225 + 0.254к1 + 0.0003к2, Я2 = -0.942 + 0.00003к1 + 0.974к2, Я3 = 0.214 - 0.233к1, Я4 = 0.920 - 0.930к2.
На рис. 3 представлено изменение целей х1 и х2 во времени при коэффициентах усиления управляющих воздействий к1 = к2 = 1 при линейном законе управления.
Из рис. 3 видно, что обе цели в течение планового периода также не достигаются. Оба плановых значения показателей достигаются лишь при ' > 2.2 года.
На рис. 4 представлены графики зависимостей всех частей функционала Г = Г(кр к2) при принятых значениях весовых коэффициентов удельных затрат: р1 = р2 = Р1 = Р2 = К1 = К2 = 1 для линейного закона управления.
В аналитическом виде для линейного закона управления составляющие функционала Г при приведенном выше равенстве единице значений весовых коэффициентов удельных затрат на управление имеют вид: У; = 0.892-Ю-^2 +0.011 к22 - (0.406к2 - 0.187)10-2к1 -- 0.021к2 + 0.118, (20) = 0.141•10-(k22 + 0.172-10-2к22 --(0.648к2 - 0.271)10-(k2 -- 0.36740-2к2 + 0.323-10-2,
¿2*= ^
Минимальные значения терминальной части У* функционала также определялись из решения системы уравнений дУ*/дк1 = 0 и дУ*/дк2 = 0 в области действительных чисел и условия к1 > 0 и к2 > 0 для р1, р2 е (0..2]. Полученные минимальные значения У*Шп, которым соответствуют значения коэффициентов к1 = 1.887 и к2 = 1.290, а также достигнутые требуемые значения целей [х1] = 0.233 и [х2] = 0.93, представлены в табл. 3.
Минимальные значения составляющей функционала Г определялись из решения системы уравнений дL1'/дk2 = 0 и дЬ1'/дк2 = 0 в области действительных чисел и условия к1 > 0 и к2 > 0 для Р1, Р2 е (0..2]. Полученные минимальные значения Ь1'т2п, которым соответствуют значения коэффициентов к1 и к2, представлены в табл. 4. Обе цели достижимы при любых Р1, Р2 из указанного интервала: 0.237 ^ [х1] ^ 0.250, 0.945 ^ [х2] ^ 0.946. Однако следует отметить, что здесь, по сравне-
(21) (22)
X2(t) 0.94 Н
0.93
0.92-
0.91 ■
0.90-
о
3 t
а) б)
Рис. 3. Результаты моделирования целей в области качества при линейном законе управления (k1 = k2 = 1):
а) цель x1; б) цель x2
Fig. 3. The results of modeling quality goals with a linear control law (k1 = k2
a) goal x1; b) goal x2
1):
Рис. 4. Графики зависимостей составляющих функционала Г при линейном законе управления Fig. 4. Graphs of the dependencies of the components of the functional Г with a linear control law
нию со ступенчатым законом управления, значения коэффициентов усиления намного выше: так, если при ступенчатом законе максимальное значение коэффициента к1 = 1.5 и к2 = 1.172, то для линейного закона и приведенных Р1 и Р2 они составляют: к1 = 1.683 и к2 = 1.573.
Минимальные значения составляющей Ь2* функционала Г определялись из решения системы уравнений дЬ2/дк1 = 0 и дЬ2/дк2 = 0 в области действительных чисел и условия к1 > 0 и к2 > 0 для К1,К2 е (0..2]. Полученные минимальные значения Ъ2тП, также, как и при ступенчатом законе управления, совпадают при всех К1, К2 е (0..2], и равны нулю. При этом им соответствуют одинаковые значения коэффициентов к1 = к2 = 1, т.е. обе цели на конец года недостижимы: для всех случаев [х1] = 0.214, [х2] = 0.915.
Общий функционал /*(к1,к2) при принятых значениях весовых коэффициентов удельных затрат на управление, равных единице, имеет вид: Г = 0.055к12 + 0.878к22 - (0.106 + + 0.471-10-2к2)к1 - 1.755к2 + 0.934. (23) Минимальное значение функционала Гт1п = 0.163-10-2 достигается при значениях коэффициентов усиления к1 = 1.005, к2 = 1.002, при которых обе цели в области качества не достижимы на конец года. В определенной ранее ОДЦ для значе-
Таблица 3. Минимальные значения составляющей У^ функционала Г* при линейном законе управления
Р2=0.1 р2=0.5 Р2=1 Р2=1.5 Р2=2
Р1=0.1 1.16-10-21 2.26-10-20 3.97-10-20 7.29-10-21 10-23
pi=0.5 8.6-10-22 5.97-10-21 1.16-10-20 3.33-10-21 6.8-10-22
Pi=1 10-22 10-22 1.16-10-20 9-10-22 9-10-22
Pi=1.5 0 6.1-10-21 1.28-10-20 0 0
pi=2 1.2-10-21 0 8-10-22 6-10-22 1.5-10-21
Таблица 4. Минимальные значения составляющей функционала Г при линейном законе управления
p2=0.1 p2=0.5 p2=i p2=i.5 p2=2
Pi=0.1 7.049-i0-5 при ki = 2.66i k2 = i.572 2Л4Ы0-4 при ki = 2.549 k2 = i.553 3.933-i0-4 при ki = 2.4i0 k2 = i.530 5.722-i0-4 при ki = 2.272 k2 = i.506 7.508-i0-4 при ki = 2.i35 k2 = i.483
pi=0.5 2.087-i0-4 при ki = 2.682 k2 = i.572 3.524-i0-4 при ki = 2.66i k2 = i.572 5.320-i0-4 при ki = 2.633 k2 = i.567 7.ii5-i0-4 при ki = 2.605 k2 = i.563 8.9i0-i0-4 при ki = 2.577 k2 = i.558
Pi=1 3.8i5-i0-4 при ki = 2.683 k2 = i.569 5.253-i0-4 при ki = 2.674 k2 = i.573 7.049-i0-4 при ki = 2.66i k2 = i.572 8.845-i0-4 при ki = 2.647 k2 = i.570 i0.640-i0-4 при ki = 2.633 k2 = i.567
pi=1.5 5.543-i0-4 при ki = 2.683 k2 = i.565 6.98Ы0-4 при ki = 2.679 k2 = i.573 8.777-i0-4 при ki = 2.670 k2 = i.573 i0.573-i0-4 при ki = 2.66i k2 = i.572 i2.369-i0-4 при ki = 2.65i k2 = i.570
pi=2 7.27Ы0-4 при ki = 2.682 k2 = i.56i 8.709-i0-4 при ki = 2.68i k2 = i.573 i0.505-i0-4 при ki = 2.674 k2 = i.573 i2.30i-i0-4 при ki = 2.668 k2 = i.573 i4.098-i0-4 при ki = 2.66i k2 = i.572
ний коэффициентов к1 = 1.9 и к2 = 1.3, минимальное значение функционала составляет Г^ = 0.123.
АНАЛИЗ РЕЗУЛЬТАТОВ
На рис. 5 представлены графики зависимостей составляющих: У1(к1,к2),^1(к1,к2), Ь1(к1,к2) и 11*(к1,к2) функционалов I и Г при двух законах управления.
Графики зависимостей Ь2(к1,к2) и Ь2'(к1,к2) на рис. 5 не представлены, т.к. полностью совпадают друг с другом.
Итак, для рассмотренной математической модели при целевом управлении в области качества являющейся системой обыкновенных дифференциальных уравнений с постоянными коэффициентами, представлены решения для ступенчатого и линейного законов управления. Это позволило использовать стандартный метод поиска экстремумов функции многих переменных для нахождения оптимальных решений без использования вариационного исчисления. Рассмотрен пример применения классического квадратичного функционала, соответствующего
а) б)
Рис. 5. Графики зависимостей составляющих функционалов I и Г при двух законах управления:
а) Vt(kpk2) и Vt(k1,k2) при Pl = р2 = 1; б) ¿1(k1,k2) и ¿1*(k1,k2) при = Р2 = 1 Fig. 5. Graphs of the dependencies of the components of the functionals I and Г for two control laws: a) V^k) and V^k) for ^ = P2 = 1; b) ^(k^) and L/^^) for P1 = P2 = 1
функции потерь качества Тагути, с использованием фактических данных по результатам деятельности промышленного предприятия.
ЗАКЛЮЧЕНИЕ
Результаты моделирования подтвердили, что для достижения поставленных целей в течение планового периода времени необходимо задавать их повышенные значения. При этом для линейного закона управления значения коэффициентов усиления (k1 = 1.683 и k2 = 1.573) выше, чем для ступенчатого закона управления (k1 = 1.5 и k2 = 1.172), хотя линейный закон предполагает более «спокойную» работу по выполнению плановых показателей. Применение этих законов в практике СМК безусловно связано со стилем управления и корпоративной культурой. При ступенчатом законе управления руководство предприятия имеет некоторые основания требовать выполнения плановых показателей уже на следующий день после их утверждения, то есть в течение всего планового периода времени. Использование функции потерь качества Тагути позволило раскрыть различие между двумя подходами, используемыми в практике планирования. При первом походе необходимо безусловное достижение утвержденных показателей, а второй подход допускает некоторое их недовыполнение при снижении затрат. Существенный интерес представляет отдельное исследование составляющей функционала, определяемой способами реализации управляющих воздействий в социально-экономических системах и соответствующими им затратами на управление, в частности взаимосвязи между коэффициентами усиления в матрице параметров управления и элементами системной матрицы, отражающими сопротивление персонала.
СПИСОК ЛИТЕРАТУРЫ
1. Анцев В.Ю., Иноземцев А.Н. Всеобщее управление качеством: учеб. пособие для студентов вузов, обучающихся по специальности 220501 Управление качеством. Тула: ТулГУ, 2005. 243 с.
2. Аникеева О.В., Ивахненко А.Г., Сторублев М.Л. Моделирование влияния значений параметров взаимодействия потенциала и организационного сопротивления на достижимость целей в области качества при ступенчатом виде управления // Известия ТулГУ. Технические науки. 2020. Вып. 10. С. 3-9.
3. Аникеева О.В., Ивахненко А.Г. Обеспечение достижимости целей в области качества с помощью целенаправленного изменения значений показателей подсистем промышленных предприятий при линейном законе управления // Фундаментальные и прикладные проблемы техники и технологии. 2020. № 6 (344) С. 156-165.
4. Ивахненко А.Г., Аникеева О.В. Постановка задачи
оптимизации при управлении качеством в пространстве состояний // Избранные научные труды двадцатой Международной научно-практической конференции «Управление качеством». М.: Пробел-2000, 2021. С. 154-159.
5. Taguchi Genichi, Subir Chowdhury, Yuin Wu. Taguchi's Quality Engineering Handbook. Wiley-Interscience. 2004. 1696 p.
6. Федорюк М.В. Обыкновенные дифференциальные уравнения. 3-е изд. СПб: Лань, 2003. 447 с.
7. Максимова Н.А. Разработка методов и моделей принятия оптимальных управленческих решений для обеспечения организационной устойчивости предприятий текстильной и легкой промышленности на базе совершенствования организации складского хозяйства: дисс. ... канд. техн. наук. Санкт-Петербург, 2019. 154 с.
REFERENCES
1. Antsev V.Yu., InozemtsevA.N. Vseobshchee upravlenie kachestvom: uchebnoe posobie [Total Quality Management]. Tula, TSU, 243 p. (in Russ.).
2. Anikeeva O.V., Ivakhnenko A.G., Storublev M.L. [Modeling the influence of the values of the interaction parameters of potential and organizational resistance on achievement of goals in the area of quality with step-by-step control] // Izvestiya Tula State University. Technical sciences,
2020, no. 10, pp. 3-9. (in Russ.).
3. Anikeeva O.V., Ivakhnenko A.G. [Ensuring the achievability of quality goals by purposefully changing the values of indicators of industrial enterprises subsystems under the linear management law] // Fundamental and applied problems of engineering and technology, no. 6 (344), 2020, pp. 156-165. (in Russ.).
4. Ivakhnenko A.G., Anikeeva O.V. [Setting optimization problem in quality management in the state space] // Izbrannye nauchnye trudy dvadcatoj Mezhdunarodnoj nauchno-prakticheskoj konferencii «Upravlenie kachestvom». Moscow: Probel-2000,
2021, pp. 154-159. (in Russ.).
5. Taguchi Genichi, Subir Chowdhury, Yuin Wu. Taguchi's Quality Engineering Handbook. Wiley-Interscience. 2004. 1696 p.
6. Fedoryuk M.V. Obyknovennye differencial'nye uravneniya [Ordinary differential equations]. St. Petersburg: Lan', 2003, 447 p. (in Russ.).
7. Maksimova N.A. Razrabotka metodov i modelej prinyatiya optimal'nyh upravlencheskih reshenij dlya obespecheniya organizacionnoj ustojchivosti predpriyatij tekstil'noj i legkoj promyshlennosti na baze sovershenstvovaniya organizacii skladskogo hozyajstva [Development of methods and models for making optimal management decisions to ensure the organizational stability of textile and light industry enterprises on the basis of improving the organization of warehouse management]: PhD thesis. St. Petersburg, 2019, 154 p. (in Russ.).
OPTIMUM MANAGEMENT IN ACHIEVING THE INDUSTRIAL QUALITY GOALS
© 2021 A.G. Ivakhnenko, O.V. Anikeeva
Southwest State University, Kursk, Russia
A mathematical model for targeted control in the field of quality, which is a system of ordinary differential equations with constant coefficients, is considered. The solutions of this system are given for the most common stepwise and linear control laws in practice. An example of optimization based on the classical quadratic functional corresponding to the Taguti quality loss function is considered, using actual data based on the results of an industrial enterprise.
Keywords: optimal management, quality goals, mathematical model, functional. DOI: 10.37313/1990-5378-2021-23-4-18-26
Alexander Ivakhnenko, Doctor of Technical Sciences, Professor at the Department of Mechanical Engineering Technologies and Equipment. E-mail: [email protected]
Olesya Anikeeva, Candidate of Technical Sciences, Associate Professor at the Department of Design and Fashion Industry. E-mail: [email protected]