Идентификация модели измерительной системы

Математическая модель ИС может быть сформулирована аналитически (т. е. выведенная из законов физики и химии так называемая физическая или физико-химическая модель) либо экспериментально; в последнем случае ИС рассматривается как черный ящик — это так называемая статистическая модель.

Точность физической модели могут определять:

а) независимые переменные, то есть аргументы функции
F(x1,…xm-1,xm,T,a1,…,ak,z);

б) коэффициенты Идентификация модели измерительной системы - student2.ru функции F, вид которой задан;

в) значения коэффициентов.

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

Идентификация математических моделей — это обширное поле деятельности теории управления. Знание физических явлений и свойств элементов перед их применением для построения ИС играет в метрологии значительно более важную роль; поэтому ряд методов идентификации, присущих пункту (а), здесь не применяется. И цель идентификации, и требуемая высокая точность ее обусловили то, что для метрологии в целом характерны специфичные подходы, и из общих методов для нее приемлемы лишь некоторые. Так, использование статистического метода планирования эксперимента весьма незначительно. Во многих случаях требования точности при идентификации модели для нужд метрологии соответствуют специальным исследованиям в области физики. Примером могут служить исследования эффекта Джозефсона с целью использования этого явления для определения эталона напряжения.

Определение независимых переменных модели

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

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

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

Корреляционный анализ

Для исследования взаимосвязей между случайными переменными используется функция взаимной корреляции. Нормированная функция взаимной корреляции переменных U и V, реализация которых во времени непрерывна, u(t), v(t), t Идентификация модели измерительной системы - student2.ru T, описывается выражением

Идентификация модели измерительной системы - student2.ru (3.103)

где

Идентификация модели измерительной системы - student2.ru (3.104)

— функция взаимной корреляции, а

Идентификация модели измерительной системы - student2.ru (3.105)

и соответственно Идентификация модели измерительной системы - student2.ru — функции корреляции переменных U, V. Для случайных переменных, стационарных в узком смысле, функция корреляции зависит только от интервала между аргументами Идентификация модели измерительной системы - student2.ru , т. е. от Идентификация модели измерительной системы - student2.ru , следовательно,

Идентификация модели измерительной системы - student2.ru (3.106)

Нормированная функция взаимной корреляции принимает значения от —1 до +1. Если переменные независимы, то для них Идентификация модели измерительной системы - student2.ru . Однако в общем случае свойство Идентификация модели измерительной системы - student2.ru означает, что между переменными U, V нет корреляции (это не равнозначно их независимости). Лишь для нормальных распределений отсутствие корреляции означает независимость случайных переменных. Значение Идентификация модели измерительной системы - student2.ru = ± 1 указывает на то, что между переменными существует линейная связь, а именно:

если Идентификация модели измерительной системы - student2.ru то U = kV + a, (3.107a)

если Идентификация модели измерительной системы - student2.ru то U(t + Идентификация модели измерительной системы - student2.ru )=kV(t)+a, (3.107б)

где k, a = const. Если ρ = 1, то k > 0; ρ = - 1 означает, что k < 0.

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

V = U + Z, (3.108)

где Z — независимая случайная переменная со свойствами белого шума, E(Z)=0, E(Z—Z)2 =σz2; тогда после подстановки (3.108) в формулу (3.104) получаем выражение

Идентификация модели измерительной системы - student2.ru (3.109)

и соответственно выражение для коэффициента корреляции

Идентификация модели измерительной системы - student2.ru (3.110)

Чем больше доля случайной составляющей в зависимости V=f(U), характеризуемой случайной переменной Z, тем коэффициент корреляции меньше 1. Мерой пропорции здесь является отношение дисперсий. Если случайная составляющая значительно превышает изменения переменной U, то Идентификация модели измерительной системы - student2.ru , т. е. коэффициент корреляции уменьшается до нуля.

Зная распределение, можно вычислить доверительный интервал коэффициента R, вычисленного по результатам пробных измерений. Однако удобнее использовать распределение Стьюдента с п—2 степенями свободы, которое получается после замены переменной R на t:

Идентификация модели измерительной системы - student2.ru (3.111)

По таблицам этого распределения можно оценить, как число элементов пробного измерения (при ρ = 0) влияет на доверительный уровень значения R, вычисленного из пробы. Если результат определения R(n) приходится на область р<0,5, то предположение, что ρ = 0, более вероятно, чем противоположное. Чтобы результат определения R из эксперимента был достаточно надежным в статистическом смысле, количество измерений должно быть достаточно большим.

Регрессионный анализ

Метод регрессии позволяет построить количественную модель на основе экспериментальных исследований. Данные для построения модели в случае двух переменных имеют форму соответствующих им множеств переменных X и Y, а именно {X, Y}, элементами которых являются Идентификация модели измерительной системы - student2.ru при j = 1, ..., k. Обычно из теоретического анализа известно, какая из переменных независима (входная величина), а какая зависима (выходная). Применительно к модели ИС об этом говорят сами обозначения X, У.

Регрессионной моделью называется связь между условным значением математического ожидания зависимой переменной и независимой переменной, а именно:

Идентификация модели измерительной системы - student2.ru . (3.112)

Это регрессия переменной Y относительно переменной X, т. е. регрессия первого рода.

Определение условного математического ожидания у = m(х) может производиться последовательно для каждой точки, если значения Идентификация модели измерительной системы - student2.ru (j=1,..,k)многократно реализованы. На основании оценок Идентификация модели измерительной системы - student2.ru для каждой точки получается множество точек Идентификация модели измерительной системы - student2.ru , а не функциональная зависимость. В определенных случаях за точку Идентификация модели измерительной системы - student2.ru можно принимать небольшой интервал и вычислять в этом интервале среднее значение (это так называемая задача точечной оценки). Для отдельных точек зависимости Идентификация модели измерительной системы - student2.ru можно также определить доверительные интервалы и погрешности оценки. Функциональную зависимость можно установить в том случае, когда существуют данные теоретического характера, например двумерное распределение вероятностей р(х,у),но тогда экспериментальные данные излишни.

Задача определения модели регрессионным методом, так называемой регрессии второго рода, формулируется несколько иначе. Исходными данными при этом служат результаты экспериментов Идентификация модели измерительной системы - student2.ru , j=1, ..., k,и отыскивается модель линейного вида — так называемая линейная регрессия, например

Идентификация модели измерительной системы - student2.ru , (3.113)

либо нелинейного вида — нелинейная регрессия, например

Идентификация модели измерительной системы - student2.ru , (3.114)

f(x) = W Идентификация модели измерительной системы - student2.ru (x) = b0 + bl x+ ... +b Идентификация модели измерительной системы - student2.ru x Идентификация модели измерительной системы - student2.ru . (3.115)

Считается, что модель известна с точностью до коэффициентов, а также до случайной составляющей ε, охватывающей все погрешности модели. Принимается, что Е(ε) = 0.

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

Идентификация модели измерительной системы - student2.ru (3.116)

которое и дало название метод наименьших квадратов (МНК). Требованию минимального удаления точек от кривой регрессии соответствуют оценки Идентификация модели измерительной системы - student2.ru коэффициентов модели α, β либо общей модели (3.115). Способ построения модели будет проиллюстрирован на примерах; систематизированное изложение метода можно найти в справочниках по математической статистике.

Пример регрессионного метода. Идентифицируется модель ИС для измерения содержания углерода в ванне кислородного конвертера для выплавки стали.

Содержание углерода X как измеряемая координата состояния объекта определяется путем измерения скорости обезуглероживания Y (потоков углерода, отводимых в виде газообразных продуктов сгорания) и нахождения по ней величины X. Полученные экспериментально зависимости представлены на рис. 3.18. Каждая из кривых представляет собой характеристику одной плавки, а данные, соответствующие одной кривой, сильно скоррелированы. В качестве независимых результатов можно принимать лишь точки, относящиеся к разным плавкам. Результаты измерений модели ИС могут быть представлены так, как это показано на рис. 3.19,а,б. В первом случае значения независимой переменной одинаковы, во втором - разные. Используя метод регрессии первого рода, получают точечные оценки регрессии, как это показано на рис. 3.19, в. Если значения независимой переменной различны (как на рис. 3.19, б), то область X делится на интервалы и отыскиваются точечные оценки, как это показано на рис. 3.19, г, для центрального значения интервала либо для средней точки в интервале. При идентификации модели ИС для измерения содержания углерода в ванне конвертера получены точечные оценки регрессии первого рода, показанные на рис. 3.20.

Идентификация модели измерительной системы - student2.ru

Рис. 3.18. Экспериментальные характеристики конвертерной выплавки стали.

С — содержание углерода в ванне; Vω — скорость обезуглероживания

Идентификация модели измерительной системы - student2.ru (3.117)

Результат не определяет функцию. В тех же самых точках независимой переменной можно определить доверительный интервал оценки на основе пробного эксперимента. Оценка дисперсии составляет

Идентификация модели измерительной системы - student2.ru (3.118)

а доверительный интервал при заданном доверительном уровне равен

Идентификация модели измерительной системы - student2.ru (3.119)

Значение t зависит от распределения Идентификация модели измерительной системы - student2.ru . Методом регрессии первого рода нельзя получить модель ИС в аналитической форме. Кроме того, доверительный интервал определяется по точкам, и для получения достаточной точности каждого значения Идентификация модели измерительной системы - student2.ru число измерений п Идентификация модели измерительной системы - student2.ru должно быть значительным. В многомерных моделях число измерений растет в пропорции Идентификация модели измерительной системы - student2.ru , где п Идентификация модели измерительной системы - student2.ru — число повторений в данной точке х Идентификация модели измерительной системы - student2.ru , n Идентификация модели измерительной системы - student2.ru — число точек в области каждой из независимых переменных, N — число независимых переменных (N+1—размерность модели).

Идентификация модели измерительной системы - student2.ru

Рис. 3.19. Подбор точек измерения в методе регрессии первого рода

Идентификация модели измерительной системы - student2.ru

Рис. 3.20. Результаты регрессии первого рода характеристики системы для измерения содержания углерода в ванне конвертера. Идентификация модели измерительной системы - student2.ru - скорость обезуглероживания; Идентификация модели измерительной системы - student2.ru - оценка содержания углерода в ванне

При малом числе точек получается низкая точность и вследствие этого — малая пригодность модели. Аналогичную форму модели удается получить по результатам анализа регрессии первого рода при использовании аппроксимации точек произвольным методом к произвольно выбранной функциональной зависимости. Так, для данных, представленных на рис. 3.19, в локальных пределах можно аппроксимировать результаты линейными зависимостями. Результат такой аппроксимации показан на рис. 3.19, г.

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

Идентификация модели измерительной системы - student2.ru , (3.120)

где Идентификация модели измерительной системы - student2.ru - коэффициенты модели; Идентификация модели измерительной системы - student2.ru - произвольные однозначные функции аргумента x априори известного вида (здесь предполагается, что модель известна с точностью до коэффициентов); ε — случайный процесс, который мо­делируется стохастической (а не функциональной) зависимостью Y от X, описываемой моделью (3.120); реализации Идентификация модели измерительной системы - student2.ru неизвестны, но имеют одинаковую дисперсию Идентификация модели измерительной системы - student2.ru = const[14] и не коррелированы.

Линейная относительно коэффициентов форма модели соответствует линейной регрессии. Большое значение имеет предположение об известной до коэффициентов форме модели; при этом все отклонения рассматриваются как случайные погрешности. Существенно, что число степеней свободы, используемое для оценки случайных погрешностей, больше числа случайных переменных. В МНК далее постулируется, что модель (3.120) лучше всего описывает результаты экспериментов, когда погрешности в смысле квадратичной метрики имеют наименьшее значение, а именно:

Идентификация модели измерительной системы - student2.ru (3.121)

Этому условию должны удовлетворять оценки коэффициентов модели, т. е. Идентификация модели измерительной системы - student2.ru .

Возвращаясь к примеру ИС для измерения содержания углерода, отметим, что полученная там функциональная зависимость, которая проведена через точки, определенные в результате анализа регрессии первого рода, является нелинейной. На основе теоретического анализа постулируется модель

Y = 1 – exp[-k(X-p)], (3.122)

в которой k,p = const — коэффициенты. Данные рис. 3.18 подтверждают эту зависимость. Для получения зависимости, линейной относительно коэффициентов, преобразуем модель (3.122) путем подстановки

U = - ln(1 - Y).

Тогда получается U = b0 + b1X + n

либо X*=a0+a1U+e , (3.123)

если характеристикаφ (3.5) определяется как обратная функция модели f (3.10). Выше обозначено

Идентификация модели измерительной системы - student2.ru

B результате получаются значения

Идентификация модели измерительной системы - student2.ru для Идентификация модели измерительной системы - student2.ru (3.124a)

Идентификация модели измерительной системы - student2.ru для Идентификация модели измерительной системы - student2.ru (3.124б)

С целью проверки допущений исследована дисперсия погрешности ε для результатов анализа регрессии первого рода. Результат, показанный на рис. 3.21, свидетельствует о том, что полученные оценки коэффициентов модели не самые эффективные. Чтобы получить более эффективные оценки, необходимо преобразовать модель в форму, в которой дисперсия случайной составляющей e будет постоянна. Если в диапазоне измерений < Идентификация модели измерительной системы - student2.ru > дисперсия e изменяется согласно зависимости

Идентификация модели измерительной системы - student2.ru (3.125)

причем

Идентификация модели измерительной системы - student2.ru (3.126)

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

υ = y/g(x), Идентификация модели измерительной системы - student2.ru (3.127а, б)

Идентификация модели измерительной системы - student2.ru (3.128а, б)

вместо уравнения (3.120) получаем

Идентификация модели измерительной системы - student2.ru (3.129)

Идентификация модели измерительной системы - student2.ru

Рис. 3.21. Зависимость стандартного отклонения погрешности содержания углерода от измеряемой величины

Тогда условие Идентификация модели измерительной системы - student2.ru в точках оценки регрессии первого родa записывается так:

Идентификация модели измерительной системы - student2.ru (3.130)

где

Идентификация модели измерительной системы - student2.ru (3.131)

называется весом (весовым коэффициентом) данной точки.

Выражение (3.130) является условием метода взвешенных наименьших квадратов — сокращенно МВНК. Оценки коэффициентов модели ИС для измерения содержания углерода в ванне конвертера, полученные по МВНК, имеют другие значения в сравнении с результатами метода регрессии первого рода (3.124), а именно:

Идентификация модели измерительной системы - student2.ru для Идентификация модели измерительной системы - student2.ru (3.132a) Идентификация модели измерительной системы - student2.ru для Идентификация модели измерительной системы - student2.ru (3.132б)

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

Наши рекомендации