ИЗВЕСТИЯ
ТОМСКОГО ОРДЕНА ТРУДОВОГО КРАСНОГО ЗНАМЕНИ ПОЛИТЕХНИЧЕСКОГО
ИНСТИТУТА имени С. М. КИРОВА
1969
Том 168
ПРИМЕНЕНИЕ ЛИНЕЙНОГО ИНТЕРПОЛИРОВАНИЯ В ЗАДАЧАХ РЕГИСТРАЦИИ РАСПРЕДЕЛЕННЫХ ПАРАМЕТРОВ
В. Ф. ДЯДИК, И. Э. НААЦ
(Представлена научным семинаром вычислительной лаборатории ТПИ)
Дискретные методы находят самое широкое распространение в регистрации технологических параметров ввиду быстродействия, компактной формы представления текущей информации.
Системы дискретной регистрации могут включать в себя устройства, осуществляющие формирование интерполирующей функции по дискретным отсчетам.
В настоящей статье рассматриваются основные вопросы применения принципа линейного интерполирования для регистрации кривых распределения параметра в дискретном методе измерения.
Принцип линейного интерполирования состоит в том, что показания двух смежных датчиков соединяются отрезком прямой и в целом по длине поля [а, Ь] носит кусочно-линейный характер. Уравнение интерполирующей прямой для любой пары отсчетов (у£-, }'гм) записывается в виде
Р\ (*) = У'1 + (У 1+\ — уд —-' (!)
Х1+1 — х-,
где х0, х1, х2,..., х% — узлы интерполирования (¿ = 0,1,2...) и я + 1 — число датчиков.
Значение параметра в любой точке между отсчетами по длине агрегата определяется, таким образом, выражением (1). При этом абсолютная погрешность о определяется выражением
8 <!/(*)-Л (*)1- (2)
Ниже рассматривается случай линейного интерполирования с равномерным шагом к = х-1+\ --- х-ь для любой пары точек. Остаточный член интерполяционной формулы (1) для промежутка интерполирования (Х1,Х[+1) имеет вид [1]:
(3)
%_^.
Делая подстановку I — -— , получаем
Л
(л: — Х[) (х — х-1+\) = кЧ (Ь — 1).
Произведение ¿(^—1) принимает наибольшее по модулю значе-* 1 1
ние при I = —, равное — •
В (3) точка ; принадлежит интервалу (л*,-, ¿7+1). Введем величину Л/2 = тах|/"(Е)|
а ^ х ^ Ь.
Тогда
Л12Л2
(4)
8
и величина к определяется выражением
<5>
Соотношение (5) позволяет определить расстояние между двумя соседними узлами при линейном интерполировании для заданной величины ошибки интерполирования о, а необходимое число отсчетов на отрезке интерполирования [а, Ь]
п+ 1 = + (6)
п
Формула (4) определяет погрешность интерполирования при условии, что отсчеты у0, уи у2)..., уп и узлы интерполирования х{), хи х2,...,хп известны точно. Однако при решении задач регистрации распределенных параметров эти величины, как правило, известны приближенно.
Возникающая при этом погрешность называется неустранимой, так как она не может быть уменьшена за счет увеличения числа узлов, в отличие от погрешности интерполирования (погрешности метода), определяемой формулой (4).
Отсчеты у'ь известны с погрешностью измерения значений параметра в точках отсчетов — а погрешность фиксации узлов интерполирования—Необходимо определить неустранимую погрешность в зависимости от значений о{ и а2.
Эту задачу можно рассматривать как задачу определения величины погрешности при косвенных измерениях по известным погрешностям прямых измерений |2].
Согласно теореме сложения частных погрешностей для уравнения (1) имеем
=ти + ¡^и++(7)
\ ду0 I \ дух / \ дх01 ~ \ дхх
Вычислив значения частных производных, получим: =2(*) =
, (У1 — Уо)2
2 2 1 - 7- (х - х0) + — (* - х0)2 п А2
Ъ!
Первое слагаемое выражения (8) есть составляющая неустранимой погрешности, обусловленная неточностью отсчетов, второе слагаемое — неточностью фиксации узлов интерполирования. Выражение (8) определяет значение неустранимой погрешности для любого х и промежутке между двумя отсчетами. Для определения диапазона изменения з на интервале интерполирования определяются наибольшие и наименьшие значения а. С этой целью исследуются выражения, стоящие в квадратных скобках формулы (8) на экстремум,
и определяются значения з (х) в точках отсчетов. Соответственно находим для узлов (х0, х{)
(*о) = +
К2
(9)
а (X,) = у
К
1У1 —Уо)-'
,
-2 I ^ » ^ о
т -
(У1 - Уо)2
А2
Максимальная неустранимая погрешность имеет место в точках отсчета, минимальная в у"2 раз меньшая — в середине интервала
Если '32 = 0, то выражения (9) соответственно имеют вид
^ (Х0) - аА; ау(хг) = о,; , (10)
т. е. максимальная неустранимая погрешность уравнения (1), обусловленная неточностью измерения отсчетов, равна самой погрешности измерения.
Минимальная в ]/2 меньше погрешности измерения. При равенстве нулю имеет место неустранимая погрешность, обусловленная неточностью фиксации узлов интерполирования, при условии, что сами отсчеты уп заданы точно. Соответствующие значения этой погрешности определяются:
а (г\- У1—У о ~
сх \Хо) — -:-
к
М*1) =(И)
к
21 /2 к "
Анализ выражений (11) показывает, что неустранимая погрешность, обусловленная неточностью фиксации узлов интерполирования,
зависит от величины —-—— градиента интерполируемой функции
к
между двумя соседними отсчетами. Максимальное значение этой погрешности — так же в точках отсчетов, минимальное — в середине интервала.
Выводы
!. Рассмотрены основные вопросы применения принципа линейного интерполирования для регистрации кривых распределения параметра в дискретном методе измерения.
2. Дана методика расчета необходимого числа отсчетов.
3. Определена неустранимая погрешность интерполирования, обусловленная неточностью отсчетов и неточностью фиксации узлов интерполирования.
ЛИТЕРАТУРА
1. И. С. Березин, Н. П. Жидков. Методы вычислений. Т. 1, Физматгиз, М., 1962.
2. М. Ф. Маликов. Основы метрологии, ч. 1. Машгиз, М., 1949.
6. Заказ 2258. 81