Погрешности математического моделирования. Информационная коррекция погрешностей измерений Математическая модель погрешности мостовой измерительной схемы

В общем виде модель погрешности 0,95(t) может быть представлена в виде 0,95(t) = 0 + F(t), где D0 -- начальная погрешность СИ; F(t) -- случайная для совокупности СИ данного типа функция времени, обусловленная физико-химическими процессами постепенного износа и старения элементов и блоков. Получить точное выражение для функции F(t) исходя из физических моделей процессов старения практически не представляется возможным. Поэтому, основываясь на данных экспериментальных исследований изменения погрешностей во времени, функцию F(t) аппроксимируют той или иной математической зависимостью.

Простейшей моделью изменения погрешности является линейная:

где v скорость изменения погрешности. Как показали проведенные исследования, данная модель удовлетворительно описывает старение СИ в возрасте от одного до пяти лет. Использование ее в других диапазонах времени невозможно ввиду явного противоречия между определенными по этой формуле и экспериментальными значениями частоты отказов.

Метрологические отказы возникают периодически. Механизм их периодичности иллюстрирует рис.1, а, где прямой линией 1 показано изменение 95%-ного квантиля при линейном законе.

Рис. 2.

При метрологическом отказе погрешность D0,95(t) превышает значение Dпр=D0+nD3, где D3 -- значение запаса нормируемого предела погрешности, необходимого для обеспечения долговременной работоспособности СИ. При каждом таком отказе производится ремонт прибора и его погрешность возвращается к исходному значению D0. По прошествии времени Тр= ti - ti-1 опять происходит отказ (моменты tt, t2, t3 и т.д.), после которого вновь производится ремонт. Следовательно, процесс изменения погрешности СИ описывается ломаной линией 2 на рис.1, а, которая может быть представлена уравнением

где n -- число отказов (или ремонтов) СИ. Если число отказов считать целым, то это уравнение описывает дискретные точки на прямой 1 (рис.2, а). Если же условно принять, что п может принимать и дробные значения, то формула (2) будет описывать всю прямую 1 изменения погрешности D0,95(t) при отсутствии отказов.

Частота метрологических отказов увеличивается с ростом скорости v. Она столь же сильно зависит от запаса нормируемого значения погрешности D3 по отношению к фактическому значению погрешности средства измерений D0 на момент изготовления или окончания ремонта прибора. Практические возможности воздействия на скорость изменения v и запас погрешности D3 совершенно различны. Скорость старения определяется существующей технологией производства. Запас погрешности для первого межремонтного интервала определяется решениями, принятыми производителем СИ, а для всех последующих межремонтных интервалов -- уровнем культуры ремонтной службы пользователя.

Если метрологическая служба предприятия обеспечивает при ремонте погрешность СИ, равную погрешности D0 на момент изготовления, то частота метрологических отказов будет малой. Если же при ремонте лишь обеспечивается выполнение условия D0 (0,9... 0,95) Dпр, то погрешность может выйти за пределы допустимых значений уже в ближайшие месяцы эксплуатации СИ и большую часть межповерочного интервала оно будет эксплуатироваться с погрешностью, превышающей его класс точности. Поэтому основным практическим средством достижения долговременной метрологической исправности средства измерений является обеспечение достаточно большого запаса D3, нормируемого по отношению к пределу Dпр.

Постепенное непрерывное расходование этого запаса обеспечивает на некоторый определенный период времени метрологически исправное состояние СИ. Ведущие приборостроительные заводы обеспечивают D3 = (0,4...0,5) Dпр, что при средней скорости старения v = = 0,05АП /год позволяет получать межремонтный интервал Тр= D3 = 1/Т/v = 8... 10 лет и частоту отказов р= 0,1... 0,125 год-1.

При изменении погрешности СИ в соответствии с формулой (1) все межремонтные интервалы Тр = 1/Т будут равны между собой, а частота метрологических отказов р будет постоянной в течение всего срока эксплуатации. Однако проведенные экспериментальные исследования показали, что на практике это не выполняется.

Математические модели

Построенные выше физические модели крайне важно описать с помощью символов в виде математических формул и уравнений. Эти символы – параметры объектов (они же обозначают физические величины) – связаны между собой в виде выше сформулированных физических законов.

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

Следовательно, о физических величинах можно говорить как о параметрах, характеризующих и качественно, и количественно построенные физические модели.

Процесс создания математической модели можно также разделить на 3 этапа:

Этап 1. Составление формул и уравнений, описывающих состояние, движение и взаимодействия объектов в рамках выбранных физических моделœей.

Этап 2. Решение и исследование сугубо математических задач сформулированных на первом этапе. Основным вопросом здесь является решение так называемой прямой задачи, ᴛ.ᴇ. получение теоретических следствий и численных данных. На этом этапе важную роль играет математический аппарат и вычислительная техника (компьютер).

Этап 3. Выяснение того, согласуются ли результаты анализа и вычислений с результатами измерений в пределах точности последних. Отклонение результатов расчётов от результатов измерений свидетельствует:

Либо о неправильности применённых математических методов;

Либо о неверности принятой физической модели;

Либо о неверности процедуры измерений.

Выяснение источников ошибок требует большого искусства и высокой квалификации исследователя.

Бывает, что при построении математической модели некоторые её характеристики или связи между параметрами остаются неопределёнными вследствие ограниченности наших знаний о физических свойствах объекта. К примеру: иногда оказывается, что число уравнений, описывающих свойства объекта и связи между объектами, меньше числа параметров (физических величин), характеризующих объект. В этих случаях приходится вводить дополнительные уравнения, характеризующие объект и его свойства, иногда даже пытаются угадать эти свойства, для того, чтобы задача была решена, а результаты соответствовали результатам опытов в пределах заданной погрешности. Подобного образа задачи называются обратными.

Проблема достоверности наших представлений об окружающем мире, ᴛ.ᴇ. проблема соответствия модели объекта и реального объекта͵ является ключевой проблемой в теории познания. Сегодня общепринято, что критерием истинности наших знаний является опыт. Модель адекватна объекту, в случае если результаты теоретических исследований (расчёт) совпадают с результатами опыта (измерений) в пределах погрешности последнего.

Погрешности имеют место не только при измерениях, но и при теоретическом моделировании. Для теоретических моделœей, в соответствии с природой возникновения, будем различать:

Погрешности, возникающие при разработке физической модели;

Погрешности, возникающие при составлении математической модели;

Погрешности, возникающие при анализе математической модели;

Погрешности, связанные с конечным числом разрядов чисел при вычислениях.

В последнем случае, к примеру, число π в рамках символической записи как отношение длины окружности к диаметру представляет собой точное число, но попытка записать его в численном виде (π=3,14159265…) вызывает погрешность, связанную с конечным числом разрядов.

Перечисленные погрешности возникают всœегда. Избежать их невозможно, и их называются методическими . При измерениях методические погрешности проявляют себя как систематические.

Пример : погрешности физической и математической модели маятника, возникающие при измерении периода колебаний маятника в виде тела, подвешенного на нити.

Физическая модель маятника :

Нить – невесома и нерастяжима;

Тело – материальная точка;

Трение отсутствует;

Тело совершает плоское движение;

Гравитационное поле – однородное (ᴛ.ᴇ. g =const во всœех точках пространства, в которых находится тело);

Влияние других тел и полей на движение тела отсутствует.

Очевидно, что реальное тело не должна быть материальной точкой, оно имеет объём и форму, в процессе движения или со временем тело деформируется. Вместе с тем, нить имеет массу, она обладает упругостью и также деформируется. На движение маятника влияет движение точки подвеса, обусловленное действием вибраций, всœегда имеющих место. Также на движение маятника влияет сопротивление воздуха, трение в нити и способ ее крепления, внешние магнитное и электрическое поля, неоднородность гравитационного поля Земли и даже влияние гравитационного поля Луны, Солнца и окружающих тел.

Перечисленные факторы, в принципе, бывают учтены, однако сделать это достаточно трудно. Для этого потребуется привлечь почти всœе разделы физики. В конечном счете, учет этих факторов значительно усложнит физическую модель маятника и ее анализ. Не учет перечисленных, а также множества других, не упомянутых здесь факторов, существенно упрощает анализ, но приводит к погрешностям исследования.

Математическая модель маятника :

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

, (1), где L – длина нити; φ – отклонение тела от положения равновесия.

При φ<<1 обычно считают, что sin φʼʼφ, и тогда уравнение движения записывается:.(2)

Это – линœейное дифференциальное уравнение, ĸᴏᴛᴏᴩᴏᴇ должна быть решено точно. Данноерешение имеет вид , где . Отсюда следует, что период колебаний маятника Т 0 =2p/w 0 не зависит от амплитуды φ 0 . При этом, это решение нельзя считать точным решением задачи о колебаниях маятника, представленного простейшей физической моделью, поскольку исходное уравнение (1) было другим.

Можно уточнить решение. В случае если разложить sin φ в ряд и учесть хотя бы первые два члена разложения, ᴛ.ᴇ. считать, что sinφʼʼφ+φ 3 /6, то решение дифференциального уравнения существенно усложнится. Приближенно его можно записать в виде , где . Отсюда следует, что в данном приближении период колебаний маятника Т =2p/w зависит от амплитуды колебаний по параболическому закону.

Τᴀᴋᴎᴍ ᴏϬᴩᴀᴈᴏᴍ, погрешность математической модели (уравнение (2)), связанная с заменой sin φ на φ, приводит к погрешности результата расчета периода колебаний маятника. Оценка этой погрешности должна быть получена из решения задачи во втором приближении.

Проблема построения и анализа математической модели объекта исследования с заданной точностью, а также оценка погрешности расчётов в ряде случаев очень сложна. Требуется высокая математическая культура исследователя, необходим тщательный математический анализ и самой модели, и применяемых методов решения.

К примеру, не имеет смысла требование решения уравнения (1) с точностью, существенно превышающей точность построения физической модели. В частности, в предыдущем примере нет смысла делать замену sinφʼʼφ+φ 3 /6 вместо sinφʼʼφ, в случае если нить заметно деформируется или сопротивление воздуха велико.

Применение ЭВМ значительно увеличило возможности построения и исследования математических моделœей в технике, однако не следует думать, что совершенное знание математики, численных методов и языков программирования позволит решить любую физическую и прикладную задачу. Дело в том, что даже самые изящные и точные методы расчетов не могут исправить ошибки, допущенные при построении физической модели. Действительно, в случае если длина L не постоянна, или если размеры тела сопоставимы с длиной нити, или трение велико и колебания маятника быстро затухают, то даже абсолютно точное решение уравнения (1) не позволит получить точное решение задачи о колебаниях маятника.

Общая характеристика понятия “измерение” (сведения из метрологии)

В метрологии определœение понятия “измерение” даёт ГОСТ 16.263-70.

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

Измерение включает в себя следующие понятия:

Объект измерения;

Цель измерения;

Условия измерения (совокупность влияющих величин, описывающих состояние окружающей среды и объектов);

Метод измерения, ᴛ.ᴇ. совокупность приёмов использования принципов и средств измерений (принцип измерения – совокупность физических явлений, положенных в основу измерения);

Методика измерения, ᴛ.ᴇ. установленная совокупность операций и правил, выполнение которых обеспечивает получение необходимых результатов в соответствии с данным методом.

Средства измерения:

▪ измерительные преобразователи,

▪ измерительные приборы,

▪ измерительные установки,

▪ измерительные системы,

▪ измерительно-информационные системы;

Результаты измерений;

Погрешность измерений;

Понятия, характеризующие качество измерений:

▪ достоверность (характеризуется доверительной вероятностью, ᴛ.ᴇ. вероятностью того, что истинное значение измеряемой величины находится в указанных пределах);

▪ правильность (характеризуется значением систематической погрешности);

▪ сходимость (близость друг к другу результатов измерений одной и той же величины, выполняемых повторно одними и теми же методами и средствами и в одних и тех же условиях; отражает влияние случайных погрешностей на результат);

▪ воспроизводимость (близость друг к другу результатов измерений одной и той же величины, выполняемых в разных местах, разными методами и средствами, но приведенных к одним и тем же условиям).

Погрешности теоретических моделей - понятие и виды. Классификация и особенности категории "Погрешности теоретических моделей" 2017, 2018.

Структурная схема измерительного устройства обычно представляется соединением звеньев, каждое из которых отражает математическую модель отдельного элемента или части ИУ.

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

На рис. 3.1 показаны обозначения наиболее распространенных звеньев стационарных ИУ: безынерционного звена с нелинейной функцией преобразования у = f(x) (рис. 3.1, о), безынерционного звена с линейной пропорциональной функцией преобразования у = kx (рис. 3.1, б), инерционного линейного звена с передаточной функцией W(p) (рис. 3.1, в).

Рис. 3.1.

а - нелинейное безынерционное; 6 - линейное безынерционное; в - линейное инерционное; г - сумматор; д - сравнивающее устройство: е - с т входами и п выходами; ж - с двумя входами и одним выходом; з - с одним входом и двумя выходами

Здесь же показаны структурные схемы сумматора (рис. 3.1, г) и сравнивающего устройства (рис. 3.1, Э), для которых, соответственно, можно записать z-хл-у и z = x-y. Эти устройства также можно рассматривать в качестве звеньев структурной схемы ИУ, так как в общем случае звено может иметь несколько (т) входов и несколько (??) выходов (см. рис. 3.1, е). Например, мостовые схемы с двумя активными плечами имеют два входа и один выход (рис. 3.1, ж), а дифференциальные И У - один вход и два выхода (рис. 3.1, з ).

Звенья могут быть пронумерованы, а их передаточные функции, уравнения, или характеристики указаны вне структурной схемы. На рис. 3.2 показан пример такой структурной схемы. В этом случае ИУ состоит из двух безынерционных звеньев (линейного звена 1 (с характеристикой у х = 5х{) и нелинейного звена 2 (с характеристикой у 2 = У / (10 + у {))), одного линейного инерционного звена 3 с передаточной функцией W 3 (p) = 2/(р + 3), сравнивающего устройства и сумматора, определенным образом связанных между собой.

Рис. 3.2.

Справа от рисунка показана соответствующая система уравнений. Эта система уравнений и структурная схема И У эквивалентны. Однако, поскольку математическое описание одинаковых измерительных преобразований может быть различным, возможны разные варианты структурной схемы одного и того же ИУ.

При анализе статического режима измерений считается, что измеряемая величина и информативные параметры измерительного сигнала во всех точках структурной схемы И У не изменяются во времени. В этом случае в уравнении звена 3 следует принять dy^/dt = 0. Тогда система уравнений, описывающих структурную схему ИУ, примет вид

Исключая промежуточные переменные, найдем статическую характери стику ИУ

Видно, что вследствие нелинейности статической характеристики звена 2 общая статическая характеристика ИУ оказалась нелинейной. При малых изменениях измеряемой величины х можно вместо нелинейной характеристики звена 2 у 2 = */i/(10 + у х) использовать линеаризованную характеристику этого звена у 2 = yj 10. В этом случае вместо прежнего

результата получим у - - х.

В возмущенном режиме измерений структурная схема ИУ дополняется элементами, отражающими воздействие влияющих величин и помех (см. рис. 3.9). Их также можно считать постоянными или переменными, рассматривая соответственно возмущенный статический и возмущенный динамический режимы измерений.

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

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

где т и σ – математическое ожидание и среднеквадратическое отклонение.

Распределение Гаусса неоднократно подтверждалось экспериментальными данными в диапазоне значений, соответствующих размаху ±3σ. В соответствии с этим распределением, погрешность совмещения в конкретной точке εх в направлении Х воспринимается как случайная величина, распределенная по нормальному закону, со следующими характеристиками:

(3.16)

где rx – коэффициент корреляции между величинами смещений соседних единичных участков в направлении X ; С2 x – число сочетаний из Х по 2, рассчитываемое из выражения

Из соотношений (3.15) и (3.16) выводится аналитическая запись плотности вероятности распределения величин:

Графики зависимости погрешностей совмещения от координат точек по одной оси, вытекающие из соотношения (3.18), показаны на рис. 3.59.

Рис. 3.59. Диаграмма погрешностей совмещения слоев в направлении Х

При наличии статистических данных могут быть найдены числовые характеристики распределения (3.18) для участка длиной L с шагом сетки h . Они находятся из соотношений:

(3.19)

где ML , σ L – соответственно математическое ожидание и дисперсия деформации участка длиной L ; – число сочетаний из L / h по 2.

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

Случайным процессом X(t) называется процесс (функция), значение которого при любом фиксированном значении t = tQ является случайной величиной X(t ). Конкретный вид процесса (функции), полученный в результате опыта, называется реализацией .

Рис. 4. Вид случайных функций

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

Математическим ожиданием случайной функции X(t) называется неслучайная функция которая при каждом значении аргумента t равна математическому ожиданию соответствующего сечения:

где p(x, t) - одномерная плотность распределения случайной величины х в соответствующем сечении случайного процесса X(t).

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

Корреляционная функция - неслучайная функция R(t, t") двух аргументов t и t", которая при каждой паре значений аргументов равна ковариации соответствующих сечений случайного процесса:

Корреляционная функция, называемая иногда автокорреляционной, описывает статистическую связь между мгновенными значениями случайной функции, разделенными заданным значением времени т = t"-t. При равенстве аргументов корреляционная функция равна дисперсии случайного процесса. Она всегда неотрицательна.

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

Математическое ожидание постоянно;

Дисперсия по сечениям является постоянной величиной;

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

Важной характеристикой стационарного случайного процесса является его спектральная плотность S(w), которая описывает частотный состав случайного процесса при w>О и выражает среднюю мощность случайного процесса, приходящуюся на единицу полосы частот:

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

При построении математической модели погрешности измерений следует учитывать всю информацию о проводимом измерении и его элементах.

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

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

Применение методов математической статистики к обработке результатов измерений правомочно лишь в предположении о независимости между собой отдельных получаемых отсчетов;

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

Таким образом, условия непрерывности и независимости для результатов измерений и их погрешностей соблюдаются приближенно, а иногда и не соблюдаются. В математике под термином "непрерывная случайная величина" понимается существенно более узкое, ограниченное рядом условий понятие, чем "случайная погрешность" в метрологии.

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

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

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

2) исследование математических задач, к которым приходят на первом этапе. Основным вопросом здесь является решение прямой задачи, т.е. получение численных данных и теоретических следствий. На этом этапе важную роль приобретают математический аппарат и вычислительная техника (компьютер).

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

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

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