Добавил:
Upload Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

Аппроксимация характеристик преобразователей 1

.pdf
Скачиваний:
25
Добавлен:
26.03.2015
Размер:
244.65 Кб
Скачать

Лекция №

Аппроксимация характеристик преобразователей

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

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

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

Простейший вид аппроксимации – линейная аппроксимация, т.е.

аппроксимация функции преобразования линейной зависимостью.

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

В этом случае прибегают к аппроксимации более сложными функциями:

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

Наиболее часто используют следующие виды аппроксимации.

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

F(x, ai) = a0 + a1x + a2x2 + … + anxn

- представляет разложение аппроксимирующей функции в ряд Маклорена. Принимая некоторые коэффициенты ai равными 0, получим ряд новых аппроксимирующих функций.

а) если все коэффициенты ai =0, кроме a0 , поучим функцию нулевого порядка, которая применяется при линейно ступенчатой аппроксимации.

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

F(x)

F(x)

У5

 

У4

 

У3

 

У2

F(x,a)i

 

У1

 

0 X1 X2 X3 X4 X5

 

F(x)

 

F(x)

 

 

 

 

 

F(x)

У3

 

 

 

 

 

 

 

 

 

 

 

 

У2

 

 

F3(x,ai)

 

 

 

У1

 

 

F2(x,ai)

 

 

 

 

 

 

 

 

 

 

F1(x,ai)

 

x

0

x

0

X1

X2

X3 x

б) Ограничиваясь линейным членом полиномиального выражения

F(x, ai) = a0 + a1x получим линейную аппроксимацию.

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

г) При коэффициентах a2 … an не равных нулю получаем нелинейную аппроксимирующую функцию.

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

Однако линейная аппроксимация обычно сложнее в реализации.

Практическое удобство достигается при аппроксимации функции преобразования стандартной математической функцией.

Пример : При аппроксимации дробно-линейной функцией во многих случаях обеспечивается достаточная для практики точность.

F(x, ai) = (a1x + a2) / (a3x + a4)

где a1 - a4 - постоянные коэффициенты.

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

Во многих случаях на практике применяется аппроксимация полиномом. Это обусловлено тем, что:

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

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

В общем виде экспоненциальный полином как аппроксимирующая функция имеет вид:

F(x, ai) = a0 + a1еα1х + a2 еα2х + … + an еαnх

где ai – αi - постоянные коэффициенты.

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

F(x, ai) = a0 + a1еα1х

При различных значениях коэффициентов a1 – α1 характер зависимости

F(x, ai) могут иметь самый разнообразный характер.

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

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

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

Эти точки должны выбираться из условия минимизации погрешности аппроксимации во всем диапазоне интервала преобразования.

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

сумма квадратов отклонений экспериментальных точек от кривой по вертикальному

Метод наименьших квадратов

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

аппроксимируют.

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

сумма квадратов отклонений экспериментальных точек от кривой

по вертикальному отклонению должна быть наименьшей.

2 min

Пусть имеется ряд n экспериментальных точек x1,y1; x2,y2; x3,y3;…; xn,yn, а эмпирическая формула, которую мы ищем имеет вид многочлена:

y f (x) B0 B1 x B2 x2 B3 x3 ... Bm xm

Подставляя последовательно все значения x и y мы получаем n

начальных равнений

B B x

B x 2

B x 3

... B x m y

 

1

 

 

0

1 1

2 1

3 1

 

m 1

1

 

 

 

 

B B x

B x 2

B x 3

... B x m y

 

 

2

0

1 2

2 2

3 2

 

m 2

2

 

 

 

B B x

B x 2

B x 3

 

... B x m y

 

 

3

0

1 3

2 3

3 3

 

m 3

3

 

 

 

………………………………………………

B B x

B x

2 B x 3

... B x m y

n

 

n

 

 

 

0

1 n

 

 

2 n

3 n

 

m n

 

 

 

 

 

 

n должно быть значительном меньше, чем m+1/

Условие 2

 

0 приводит к следующему выражению

F B B x B x 2 B x 3 ... B x m

y

2

 

 

 

 

0

 

1

1

 

2

1

3

1

m 1

 

 

1

 

 

 

 

 

B B x

2

B x 2

B x 3

... B

x m y

2

 

 

 

0

1

 

 

 

2 2

3

2

m

2

 

2

 

 

 

 

 

 

B

B x B x 2

B x 3

... B

x m y

2

 

 

 

 

 

0

1

3

 

 

2 3

3

3

m

3

 

3

 

 

 

 

 

... B B x

n

B x 2 B x 3 ... B x m

y

n

2

min

 

 

0

 

 

1

 

2 n

 

3 n

m

 

n

 

 

 

 

Искомыми величинами в этом выражении являются коэффициенты В1,

В2, В3, …, Вm, они должны иметь такие значения при которых это выражение имеет наименьшую величину.

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

Обозначив данное выражение через А вычисляем частные производные:

F

2 B B x B x 2

B x 3

... B x m y

 

B0

0

 

1

1

2 1

 

3 1

 

m

1

1

 

 

 

 

 

 

 

2 B B x

2

B x

2 B x 3

... B x m y

 

 

0

1

 

 

2 2

3 2

 

m

2

2

 

 

2 B

B x

 

B x

2 B x 3

... B

x m y

 

 

0

1 3

 

 

2 3

3

3

 

m

3

3

 

 

... 2 B B x

n

B x 2

B x 3

... B x m

y

 

 

 

0

 

1

 

2 n

 

3 n

 

m

n

n

 

 

F

 

 

2 B B x B x 2

B x 3

... B x m y x

 

 

 

B1

0

 

1 1

2 1

 

3 1

m 1

 

1

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

2 B B x

2

B x 2

B x 3

... B x m

y

x

2

 

 

 

0

1

 

 

2 2

3 2

 

m 2

2

 

 

 

 

 

2 B

B x

 

B x 2

B x 3

... B x m y

x

 

 

0

1 3

 

 

2 3

3 3

 

m 3

3

3

 

x

 

 

... 2 B B x

n

B x 2

B x 3

... B x m

y

 

n

 

 

 

 

 

 

0

 

1

 

2 n

 

3 n

m n

 

n

 

 

 

F

 

2 B B x B x 2

B x 3

... B x m y x 2

 

 

 

 

B2

0

 

1 1

2 1

 

3 1

m 1

 

1

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

2 B B x

2

B x 2

B x 3

... B x m

y

x 2

 

 

 

0

1

 

 

2 2

3 2

 

m 2

2

 

2

 

 

 

 

2 B

B x

 

B x 2

B x 3

... B x m y

x

2

 

 

 

0

1 3

 

 

2 3

3 3

 

m 3

3

3

 

x

 

... 2 B B x

n

B x 2

B x 3

... B x m

y

 

2

 

 

 

 

 

0

 

1

 

2 n

 

3 n

m n

 

n

 

n

.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

F

2 B B x B x 2

B x 3

... B x m y x m

 

 

 

Bm

0

 

1 1

2 1

 

3 1

m 1

 

1

1

 

 

 

 

 

 

 

 

 

 

2 B B x

2

B x 2

B x 3

... B x m

y

x m

 

 

 

0

1

 

 

2 2

3 2

 

m 2

2

 

2

 

 

 

 

2 B

B x

 

B x 2

B x 3

... B x m y

x m

 

 

 

 

0

1 3

 

 

2 3

3 3

 

m 3

3

3

 

x m

... 2 B B x

n

B x 2

B x 3

... B x m

y

 

 

 

 

 

 

0

 

1

 

2 n

 

3 n

m n

 

n

 

n

Приравнивая эти выражения к нулю, и сокращая их на 2, получаем

следующие уравнения (после небольших преобразований):

 

 

nB (x

x

 

x ... x )B (x2 x2 x2

... x2 )B

... ( y

y

2

y

... y

)

 

 

 

0

 

1

2

 

3

 

n

1

1

2

3

 

n

2

 

1

 

3

n

 

(x

x

2

x ... x

 

)B (x2

x2

x2

... x2 )B

(x3 x3 x3 ... x3 )B

...

 

 

1

 

 

3

 

n

 

0

1

2

3

 

n

1

 

1

2

3

 

n 2

 

 

 

( y1 x1

y2 x2

y3 x3

... yn xn )

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(x2

x2

x2 ... x2 )B

(x3

x3

x3 ... x3 )B (x4

x4

x4 ... x4 )B ...

 

 

 

 

1

2

 

3

 

n

0

1

2

 

3

n

1

1

2

3

 

n

 

2

 

 

 

( y x2

y x2

y x2 ... y

x2 )

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1 1

 

2 2

 

3 3

 

n

n

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(xm

xm

xm

 

... xm )B (xm1

xm1

xm1

... xm1 )B

 

 

 

 

 

 

 

 

 

1

 

2

3

 

 

n

0

1

 

2

 

3

 

n

1

 

 

 

 

 

 

 

 

(xm2

xm2

 

xm2

... xm2 )B

... ( y xm y

xm y xm

... y xm )

 

 

 

 

 

 

1

 

2

 

3

 

 

n

2

 

1 1

2

2

3 3

 

 

 

n n

 

 

Эти уравнения называются нормальными уравнениями.

Число уравнений равно(m+1) и они содержат (m+1) неизвестных коэффициентов Bi по отношению к которым эти уравнения являются линейными. Решая систему уравнений (линейных) определяем значения Bi.

Правила составления нормальных уравнений:

1. В первом нормальном уравнении:

а) коэффициентом про B0 служит число n, т.е. число начальных уравнений;

б) коэффициентами при В1, В2, В3, …, Вn служит сумма всех х, причѐм каждое из слагаемых этой суммы входит в степени равной индексу при В;

2. Во втором и следующих нормальных уравнениях:

а) коэффициентами при В служит сумма всех х в степени:

во втором уравнении на единицу больше соответствующих индексов при В;

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

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