CC BY
53
Введём в рассмотрение бинарное отношение 7 на Р3 как
Y = {(1,0), (1,X), (1,1'), (X,0), (0', 0)}
и распространим его покомпонентно на векторы с компонентами в Р3, в том числе на элементы полу-решётки состояний. Пусть далее для элементного базиса B через хв обозначается вектор значений функций всех элементов из B на их аргументе х.
Теорема 2. Пусть множество переходов T совместимо для функции f и RS-базис B содержит элемент, не сохраняющий отношения 7. Тогда f реализуется функционально устойчивой к состязаниям на всех переходах из множества T схемой в базисе B, если и только если для любого у є М = UxeUfT m(xB) существует нижняя грань множества Ty = {fT(х) : х є UfT & у є m(xB)}.
Доказательство. По теореме 1 доказываемое утверждение равносильно реализуемости в базисе B расширения fr. В свою очередь, в условиях теоремы необходимым и достаточным условием реализуемости fr в базисе B является квазимонотонность отношения Гft,B, которое содержит пары (хв, fT(х)) для всех х є UfT [8, теорема 3]. Пусть отношение Гft,в квазимонотонно. Тогда поскольку для любого у є М существует нижняя грань множества {хв : у є т(хв)} (например, у), то по тесту квазимонотонности [1] существует и нижняя грань множества Ty. Необходимость доказана.
Достаточность. Построим функцию h : М ^ L, положив Н(у) = inf Ty для каждого у є М, и ее расширение g по правилу точечного продолжения [1]: g(a) = £tem(a) g(t). Согласно [1, теорема 3.2], функция g является монотонной. Рассмотрим произвольное значение х є UfT. Имеем:
д(хв) = £yem(xB) д(у) = £y£m(xB) h(^) < fT(х) что означает g < Гft ,в . СлеДовательно, отношение f,в квазимонотонно, и fr реализуется в базисе B. Теорема доказана.
Библиографический список
1. Агибалов Г.П. Дискретные автоматы на полурешет-ках. Томск: Изд-во Томск. ун-та, 1993.
2. Eichelberger E.B. Hazard Detection in Combinational and Sequential Switching Circuits // IBM Journal of Research and Development. 1965. V. 9, № 2. P. 90-99.
3. Миллер Р. Теория переключательных схем. М.: Наука, 1971. Т. 2.
4. Яблонский С.В. Основные понятия кибернетики//
Проблемы кибернетики. М.: Наука, 1959. Вып.2.
С. 7-38.
5. Рогинский В.Н. Основы дискретной автоматики. М.: Связь, 1975.
6. Чеботарев А.Н. Риск в асинхронных логических схемах // Кибернетика. 1976. № 4. С. 8-11.
7. Агибалов Т.П., Комаров Ю.М., Липский В.Б. Синтез комбинационных схем, свободных от статических состязаний // Автоматика и вычислительная техника. 1979. № 1. С. 1-6.
8. Панкратова И.А. Условия реализуемости функций на полурешётке в реальных базисах переключательных элементов // Дискретный анализ и исследование операций. Сер. 1. 2006. Т. 13, № 3. С. 40-61.
УДК 681.5:51-74
МЕТОДЫ И АЛГОРИТМЫ ИДЕНТИФИКАЦИИ ПАРАМЕТРОВ ВНЕШНЕГО ТЕПЛООБМЕНА ПРОЦЕССОВ ТЕПЛОМАССОПЕРЕНОСА
В.Н. Ткаченко, А.А. Иванова, О.В. Шуба
Институт прикладной математики и механики НАНУ, Украина E-mail: [email protected], [email protected]
Предложены методы и алгоритмы идентификации параметров внешнего теплообмена на основе метода наименьших квадратов. Проведён сравнительный анализ методов идентификации распределенных параметров.
Methods and Algorithms of Identification of External Heat Exchange Parameters for Heat and Mass Tranfer Processes
V.N. Tkachenko, A.A. Ivanova, O.V. Shuba
The methods and algorithms of identification of external heat exchange parameters based on least square method are proposed. The comparative analysis of distributerd parameters identification methods is performed.
© В.Н. Ткаченко, АА. Иванова, О.В. Шуба, 2008
В.Н. Ткаченко п др. Методы п алгоритмы идентификации параметров внешнего теплообмена ВВЕДЕНИЕ
Трудность практического использования математических моделей связана с отсутствием необходимых значений (физических констант или зависимостей) для коэффициентов уравнений модели. В связи с этим возникает задача идентификации параметров модели. Высокие требования к точности моделей вынуждают рассматривать неизвестные величины параметров распределенными в пространстве либо во времени. Одним из примеров может служить процесс тепломассопереноса в слитке, движущемся внутри машины непрерывного литья заготовок (МНЛЗ), математическая модель которого разработана в работе [1].
1. ИДЕНТИФИКАЦИЯ ПАРАМЕТРОВ ТЕПЛООБМЕНА МНЛЗ
Рассмотрим некоторый участок непрерывно движущегося стального слитка в системе координат привязанной к конструкции машины непрерывного литья (рис. 1).
Рис. 1. Участок непрерывного слитка: 1 — слиток, 2 — направление движения слитка, 3 — форсунки, распыляющие охлаждающую водо-воздушную смесь, 4 — опорные и приводные ролики
Уравнение тепломассопереноса для двумерной модели в прямоугольнике (0,1) х (0, т) выглядит следующим образом:
дТ(т,х,у) + ) дТ(т,х,у)
1
дт
дх
д
А (Т)
дТ
дх
+
д_
ду
А(Т)
дТ ду _
с(Т)р(Т) [дх
где у(т) — скорость движения среды, Т(т, х, у) — температура, с(Т) — удельная теплоёмкость, г(Т) — плотность, и 1(Т) — теплопроводность сплошной среды. Заданы начальные условия:
Т(0, х, у) = То (х, у)
и граничные условия:
дТ
А(Т) дх
= 0,
х=0
дТ
А(Т) дх
= 0,
дТ
А(Т) ду
= 0,
У=0
А(Т)
дТ
ду
У=1
= а(х) • (Гое Т 1У = 1^
(1)
Здесь а(х) — коэффициент конвективной теплоотдачи, Т — температура внутри слитка, ТУ= — температура на поверхности слитка, Тов — температура окружающей среды.
Требуется определить коэффициент теплоотдачи а(х). В качестве дополнительной информации выступают данные измерений температуры на поверхности слитка.
Такие задачи называются граничными обратными задачами [2,3]. Они являются некорректными в классическом смысле. Чтобы получить представление о степени неустойчивости решения обратной
х=т
задачи, воспользуемся методом прямого обращения. Для этого из граничного условия (1) получим выражение для вычисления а(х):
Для применения метода наименьших квадратов примем во внимание следующую информацию. Коэффициент а(х) имеет специальное распределение вдоль поверхности слитка. Известно, что на участке, накрываемом факелом форсунки, его можно приблизить параболической функцией, которая приобретает максимальное значение в точке, соответствующей координате сопла форсунки, а на остальных участках — константой. Будем искать а(х) в виде
Величина h определяется полушириной захвата факела форсунки. Следовательно, неизвестными являются всего два параметра — А и ас.
Рассмотрим сначала участки, на которых а(х) = ас = const. Обозначим K — множество узлов x, в которых мы считаем КТ постоянным. Множество остальных узлов, где КТ распределяется согласно параболическому закону, обозначим В.
После проведения конечно-разностной аппроксимации получаем формулу для невязки тепловых потоков на границе:
У
Требуется найти такое ас, чтобы сумма квадратов невязок была минимальной, т.е. чтобы выполнялось условие
Каждому узлу Хі из множества В поставим в соответствие точку уі на отрезке [-к, к] таким образом, чтобы |уі| равнялось расстоянию от соответствующего хі до координаты ближайшей форсунки. Аналогичным способом находим
Необходимо также отметить, что для определённых при помощи МНК значений ас и А выполняется достаточное условие существования минимума функции Б. Легко проверить, что частные производные второго порядка Б по каждому их этих параметров строго больше нуля.
2. РЕЗУЛЬТАТЫ ВЫЧИСЛЕНИЙ
Методом наименьших квадратов и методом прямого обращения были проведены численные расчёты. В качестве теплофизических параметров для модели были выбраны данные процесса непрерывной разливки стали (для марки ст 40), ширины сляба 1м, полутолщины сляба I = 0,1м и скорости движения слитка V = 1м/мин. Данные вычислений представлены на рис. 2. Здесь хорошо видно, что решение, полученное методом прямого обращения, является неустойчивым и непригодным для практического использования. Вторая кривая представляет сплайн аппроксимацию, которая является результатом решения той же задачи методом наименьших квадратов.
60 Научный отдел
Ó P1 і — acP2i ,
S — (P1i — acP2i )2 ^ min, Vi : x Є K.
Из необходимого условия существования экстремума S(а) находим ас
A
2
Рис. 2. Результаты определения коэффициента теплоотдачи: 1 мого обращения, 2 — методом наименьших квадратов
методом пря-
Таким образом, мы нашли сплайн-аппроксимацию распределённого в пространстве коэффициента теплоотдачи на поверхности движущегося слитка, которая даёт нам минимальное среднеквадратичное отклонение между температурой поверхности измеренной и вычисленной по модели в результате решения прямой задачи.
3. МЕТОД РЕГУЛЯРИЗАЦИИ ТИХОНОВА
Метод регуляризации Тихонова позволяет получить устойчивое решение без жестких ограничений на вид функций, описывающих изменение коэффициентов теплообмена по длине зоны теплообмена. Распределение коэффициентов а(х) отыскивается путем минимизации следующего регуляризирую-щего функционала:
т
I(а ) = I [Трасч(х ) - Гизм(ж )] 2 йх + в О [а(х )],
0
где Трасч (х) — температура, получаемая в результате решения прямой задачи теплопроводности, Тизм (х) — измеряемая температура; О [а(х)] — стабилизирующий функционал; в — параметр регуляризации, в > 0.
Структурно регуляризирующий функционал состоит из двух функционалов, первый из которых оценивает величину невязки между измеренной и расчетной температурами, второй предназначен для стабилизации получаемого решения. Величина вклада стабилизирующего функционала определяется параметром регуляризации в, значение которого предлагается определять в соответствии с «принципом невязки». Величина в должна быть такова, чтобы выполнялось следующее равенство:
т
J [Трасч (х ) - Тизм (х)] 2йх = 82,
0
где 82 — интегральная характеристика среднеквадратической погрешности измерений температур.
Предполагая р-кратную дифференцируемость искомых функций, можно предложить в качестве стабилизирующего функционала следующий функционал достаточно общего вида:
О [а(х )] =
дг а(х) дхг
йх,
где рг — весовые коэффициенты, г = 1,... ,р; р — порядок стабилизации функционала.
Задача определения параметров путем минимизации регуляризирующего функционала со стабилизирующим функционалом становится «условно корректно поставленной».
Недостатком метода регуляризации является большой объем вычислений, связанный с процедурой поиска нужного значения параметра регуляризации, а также с овражностью регуляризирующего функционала. Другой существенный недостаток метода связан с самой идеей регуляризации: сглаживания решения в пределах погрешности измерений. Чем больше погрешность, тем можно получить более гладкую кривую, но при этом возрастает опасность получения хотя и более плавной кривой, но все в большей степени отклоняющейся от истинной. В этом смысле метод регуляризации принципиально отличается от ряда методов, которые обеспечивают получение достаточно точного решения при росте погрешности в исходных данных. К таким методам относится, например, метод наименьших квадратов.
4. МЕТОД НАИМЕНЬШИХ КВАДРАТОВ ДЛЯ ИДЕНТИФИКАЦИИ РАСПРЕДЕЛЕННОГО ПАРАМЕТРА
Наряду с регуляризацией Тихонова, для идентификации распределенных параметров модели теплофизического процесса предлагается применить идею метода наименьших квадратов (МНК).
Сформулируем обратную задачу нахождения коэффициентов конвективного теплообмена для одномерного уравнения теплопроводности с граничными условиями 3-го рода.
Для простоты изложения метода предполагаем нагрев симметричным, поэтому а1 = а2 = а. В этом случае минимальный объем необходимой для решения задачи информации соответствует измерению температуры тела в какой-либо одной точке одномерной области о < х < 1. Предположим, известна температура тела на границе с внешней средой:
Задача состоит в нахождении а = а(т) как функции, зависящей от времени. Предполагая функцию а(т) непрерывной, с целью аппроксимации искомой функции воспользуемся полиномом п-й степени
степень полинома будет определяться по принципу невязки.
Решив задачу Дирихле методом конечных разностей, получим температуры Т(тг,х^-), которые будут необходимы для вычисления производной в граничном условии
в котором неизвестной является только величина конвективного теплообмена а(т).
После замены конечными разностями теплового потока получаем систему уравнений граничных условий в г моментах времени:
где ак — значение полинома в момент времени к и ак = а0 + а1 ■ к ■ Ат + ... + апкпАтп, к = 1,г, г >> п; или в матричной форме
Т (г, 0) = / (т).
а(т) = ао + аі ■ т + ... + ап ■ тn,
= а(т) [Тгр-(т) - Т(т’ 0)] ’
x=0
(Т(тг,Жі) - Т(т,хо)) = ai[Tv - Т(т,хо)], i = l,r,
/ Т(ті,xi) - Т(ті,хо) \
Ах
\ Т(тг,хі) - Т(тг,хо) /
/ ai 0 ... 0 \ / То.с. - Т(ті,хо) \
_ 0 a2 ... 0 То.с. - Т(т2,хо)
(2)
V 0 0 ... ar ) \ То.с. - Т(тг,хо) /
Научный отдел
Выделим компоненты неизвестных параметров полинома, представленного в матричной форме:
( а1 0 ... 0 \
0 а2 ... 0
V 0 0 ... аГ )
= «о
1 0 .. .0 / 1Ат 0. .. 0
0 1 .. .0 + а1 ■ 0 2Ат . .. 0
0 0 .. .1 0 0. .. гАт
/ 1Атп 0 0 \
+ ...
+ аг
0
2пАтп ... 0
\ 0 0 ... гпАтп /
Для удобства записи преобразований введем следующие обозначения:
Аг =
( 1Атг 0 0 2г Атг
\
0
0
Р =
Ах
ггАтг )
а =
а1 0 ... 0
0 а2 ... 0
00
Н =
аг )
/ Т (т1 ,Ж1) - Т (т1 ,хо) \ Л Т (т2 ,Ж1) - Т (т2 ,хо)
\ Т(тг,Ж1) - Т(тг,хо) )
( То.с. - Т(т1 ,хо) \
То.с. - Т(т2 ,хо)
V То.с. - Т(тг,хо) /
Рк и Нк — к-е элементы столбцов.
Тогда уравнение (2) в матричной форме будет иметь вид
Р = аН.
Для нахождения функции методом МНК введем меру отклонения в виде суммы квадратов разности измеренных температур от расчетных по модели
Б(а) = ^(Рк - акНк) к = 1
Используя введенные обозначения, параметр а можно записать как а(т) = ао ■ Ао + а1 ■ А1 + ... + ап ■ Ап, тогда функционал Б (а) будет иметь вид
Г
Б (а) = ^ (Рк - ак Нк )2 = (Р - (ао ■ Ао + а1 ■ А1 + ... + ап ■ Ап) ■ Н )т,
к = 1
(Р - (ао ■ Ао + «1 ■ А1 + ... + апАп) ■ Н) = РтР - 2 ^ акНтР + ^ ^ ага^-НтАг+^'Н.
к=о
г=о ^'=о
Задача сводится к нахождению коэффициентов полинома, которые бы минимизировали функционал Б (а). Пусть минимум достигается при а; = а(т'). Тогда для значений коэффициентов полинома должно быть выполнено условие
йБ(а; )
= 0, к = 0, п;
(3)
Вычислим векторы первых производных суммы квадратов отклонений по компонентам ао,а1,... ,ап: дБ _ д [(Р - аН)т(Р - аН)]
даг
даг
= 2аоНтАгН + 2а1 НтАг+1Н + ... + 2апНтАг+ПН - 2НТАгР = 0,
2
Известия Саратовского университета. 2008. Т.8. Сер. Математика. Механика. Информатика, вып. 1 Получили систему из n +1 линейных алгебраических уравнений с n + 1 неизвестными
Е ак Н1 Н = Н1 Ао Р,
к=о
п
Ё ак Нт Ак+1 Н = Нт А1Р,
к=о
п
Ё ак Нт Ак+п Н = Нт Ап Р. к=о
Из линейной независимости строк следует, что детерминант не обращается в ноль, и решение системы единственное.
Легко убедиться, что уравнения (3) определяют минимум, так как матрица вторых производных (матрица Гессе) положительно полуопределена [4].
При выборе степени аппроксимирующего полинома повышение его степени целесообразно до такого значения, при котором отклонение измеренных температур от расчетных по модели остается больше погрешности измерений. В противном случае (при более высоких степенях) погрешность измерений прямо редуцируется в погрешность определения распределенного параметра по методу МНК. И метод теряет свои фильтрующие и стабилизирующие свойства.
Библиографический список
1. Ткаченко В.Н., Иванова А.А. Анализ температурных полей криволинейной МНЛЗ на основе математического моделирования // Материалы 3-й междунар. науч.-практ. конф. «Прогресивні технології у металургії сталі: XXI сторіччя». Донецк: ДонНТУ, 2007. С. 242-249.
2. Мацевитый Ю.М. Обратные задачи теплопроводности. В 2 т. Т.2. Приложения. Киев: Наук. думка, 2003.
3. Самарский А.А., Вабищевич П.Н. Численные методы решения обратных задач математической физики. М.: Едиториал УРСС, 2004.
4. Гантмахер Ф.Р. Теория матриц. М.: Наука, 1966. 576 с.
5. Иванова А.А. Математическая модель процесса затвердевания непрерывного слитка в зоне вторичного охлаждения. // Труды ИПММ НАН Украины. Донецк: ИПММ, 2006. Вып.12. С. 76-84.
6. Ткаченко В.Н. Моделирование тепловых процессов в автоматизированных системах обработки информации // Вісник Донецького національного університету. Сер. А. Природничі науки. 2002. №2. С. 379-383.
УДК 519.713:681.3
СИНТЕЗ ПСЕВДОСЛУЧАЙНЫХ КОНТРОЛИРУЮЩИХ ТЕСТОВ ДЛЯ ДИСКРЕТНОГО УСТРОЙСТВА
И.В. Уколов
Саратовский государственный университет,
кафедра математической кибернетики и компьютерных наук
E-mail: [email protected]
В статье рассматривается псевдослучайный метод построения контролирующих тестов для дискретных устройств, применимый как к комбинационным, так и к последовательностным устройствам. Используется энтропийный подход для отыскания оптимального распределения вероятностей входных наборов, что позволяет сокращать среднюю длину генерируемых тестов. Для определения упомянутых вероятностей решается многомерная задача поиска максимума выходной энтропии с применением генетического алгоритма. Рассмотрена задача моделирования исправного устройства и неисправных модификаций с использованием системы моделирования Active-HDL. Приведены результаты построения тестов для схем из международного каталога ISCAS’89. Выполнено их сравнение с данными, полученными другими авторами с использованием иных генетических алгоритмов.
© И.В. Уколов, 2008
Synthesis of Pseudorandom Test Patterns for the Discrete Device
I.V. Ukolov
In article the pseudorandom method of test patterns generation for discrete devices applicable for both to combinational devices and to sequential devices is considered. The entropy approach for search of optimum distribution of probabilities of input vectors is used that allows reducing average length of generated tests. For definition of the mentioned probabilities the multivariate problem of search of a maximum of output entropy with application of genetic algorithm is solved. The problem of simulation of the correct device and faulty devices is considered based on simulation system Active-HDL. Results of tests generation for circuits from international catalogue ISCAS’89 are produced. Their comparison with the data collected by other authors with use of other genetic algorithms is done.
