close

Вход

Забыли?

вход по аккаунту

?

Патент BY10526

код для вставкиСкачать
ОПИСАНИЕ
ИЗОБРЕТЕНИЯ
К ПАТЕНТУ
РЕСПУБЛИКА БЕЛАРУСЬ
(46) 2008.04.30
(12)
(51) МПК (2006)
НАЦИОНАЛЬНЫЙ ЦЕНТР
ИНТЕЛЛЕКТУАЛЬНОЙ
СОБСТВЕННОСТИ
(54)
G 06G 7/00
СПОСОБ ИДЕНТИФИКАЦИИ НЕЛИНЕЙНОГО
ДИНАМИЧЕСКОГО ОБЪЕКТА
(21) Номер заявки: a 20060756
(22) 2006.07.19
(43) 2008.02.28
(71) Заявитель: Государственное научное учреждение "Объединенный
институт машиностроения Национальной академии наук Беларуси"
(BY)
(72) Авторы: Петько Валерий Иванович;
Куконин Владимир Егорович; Выставная Елена Ивановна (BY)
BY 10526 C1 2008.04.30
BY (11) 10526
(13) C1
(19)
(73) Патентообладатель: Государственное
научное учреждение "Объединенный
институт машиностроения Национальной академии наук Беларуси" (BY)
(56) SU 290282, 1971.
RU 2004135154 A, 2006.
RU 2256950 C2, 2005.
SU 1385122 A1, 1988.
US 2005/0044053 A1.
US 4368510, 1983.
(57)
Способ идентификации нелинейного динамического объекта, заключающийся в том,
что на объект идентификации подают случайный сигнал x(t) и регистрируют его отклик
y(t), определяют дискретные значения сигнала x(n) и отклика y(n), отличающийся тем,
что из y(n) определяют максимально возможные отрицательные Ymin и положительные
L
Ymax значения отклика, разбивают их на , где L - четное число, равномерных диапазо2
нов, определяют границы и середины этих диапазонов, формируют границы и середины
диапазонов разбиения входного сигнала x(n) путем подачи на объект идентификации сигналов с постоянными мгновенными значениями и определяют среди последних значения,
соответствующие границам диапазонов отклика y(n) и их серединам, вычисляют ширину
Фиг. 1
BY 10526 C1 2008.04.30
диапазонов разбиения входного сигнала x(n), определяют размерность быстрого преобразования Фурье для идентификации объекта и коэффициенты передачи для корректировки частотных характеристик на нулевой частоте, подают на объект идентификации
суммарные сигналы с постоянными мгновенными значениями, равными серединам диапазонов входного сигнала x(n), и псевдослучайными мгновенными значениями xr(n), где
r = 1, 2, …, L, размах которых выбирают таким, чтобы размах отклика объекта идентификации не выходил за пределы соответствующего диапазона, и по ним и по соответствующим им откликам yr(n) определяют L частотных характеристик с последующей
корректировкой последних на нулевой частоте, по которым путем обратного дискретного
преобразования Фурье определяют L импульсных характеристик объекта идентификации
в серединах диапазонов разбиения входного сигнала, а после определения дискретных
значений входного сигнала x(n) формируют из последнего частичные входные сигналы
xl(n) путем представления дискретных отсчетов x(n) в виде числовых отрезков, разбиения
их на диапазоны входного сигнала и вычисления отсчетов xl(n) как части числовых отрезков, попадающие в l-й диапазон входного сигнала, и создают модель объекта идентификации в виде уравнения:
L
y' (n ) = ∑ y l (n ) ,
l=1
где y'(n) - отклик модели;
L - общее число диапазонов;
y l (n ) =
ml
∑ a lj x l (n − j) ,
j= 0
где ml + 1 - количество дискретных отсчетов в l-й импульсной характеристике;
a lj - j-й отсчет l-й импульсной характеристики;
xl(n-j) - последовательность задержанных на j дискретных отсчетов частичного входного сигнала xl(n);
n - последовательность дискретных отсчетов входного сигнала с интервалом дискретизации не менее 16 отсчетов на верхнюю частоту в спектре входного сигнала.
Изобретение относится к технической кибернетике и может быть использовано при
идентификации нелинейных динамических объектов (НДО).
Известен способ скользящей кусочно-линейной аппроксимации функции. Этот способ
пригоден для кусочно-линейной аппроксимации вход-выходной характеристики безынерционного нелинейного объекта (БНО) и не применим к НДО (SU 1742831 А1, 1992).
Наиболее близким к предлагаемому способу является способ идентификации нелинейного динамического объекта путем определения дискретных значений входного сигнала х(n) и его отклика у(n), подачи их на регрессиометр с последующим вычислением
линии регрессии отклика НДО y(t) относительно входного сигнала x(t) по формуле:
M[ y / x ] =
1
Ni
Ni
∑ y ji ,
(1)
j =1
где i = 1, 2, ... - номер уровня анализа;
yji - ордината процесса y(t), отсчитанная в момент достижения процессом x(t) значения
xi;
Ni - количество ординат процесса y(t), подсчитанное при анализе i-го уровня;
j = 1, 2, 3, ..., Ni - текущий номер ординаты процесса y(t) при анализе i-го уровня
(SU 290282, 1971).
2
BY 10526 C1 2008.04.30
Недостатком этого способа является невысокая точность, т.к. в нем определяется не
детерминированная, а статистическая связь между входным сигналом на НДО и его откликом.
Задача изобретения - повышение точности идентификации НДО путем определения
детерминированной связи между входным сигналом на НДО и его откликом.
Поставленная задача решается тем, что в способе идентификации нелинейного динамического объекта, заключающемся в том, что на объект идентификации подают случайный сигнал x(t) и регистрируют его отклик y(t), определяют дискретные значения сигнала
x(n) и отклика y(n), причем, согласно изобретению, из y(n) определяют максимально возL
можные отрицательные Ymin и положительные Ymax значения отклика, разбивают их на ,
2
где L - четное число, равномерных диапазонов, определяют границы и середины этих
диапазонов, формируют границы и середины диапазонов разбиения входного сигнала x(n)
путем подачи на объект идентификации сигналов с постоянными мгновенными значениями и определяют среди последних значения, соответствующие границам диапазонов отклика y(n) и их серединам, вычисляют ширину диапазонов разбиения входного сигнала
x(n), определяют размерность быстрого преобразования Фурье для идентификации объекта и коэффициенты передачи для корректировки частотных характеристик на нулевой
частоте, подают на объект идентификации суммарные сигналы с постоянными мгновенными значениями, равными серединам диапазонов входного сигнала x(n), и псевдослучайными мгновенными значениями xr(n), где r = 1, 2, ..., L, размах которых выбирают
таким, чтобы размах отклика объекта идентификации не выходил за пределы соответствующего диапазона, и по ним и по соответствующим им откликам yr(n) определяют L частотных характеристик с последующей корректировкой последних на нулевой частоте, по
которым путем обратного дискретного преобразования Фурье определяют L импульсных
характеристик объекта идентификации в серединах диапазонов разбиения входного сигнала, а после определения дискретных значений входного сигнала x(n) формируют из последнего частичные входные сигналы xl(n) путем представления дискретных отсчетов x(n)
в виде числовых отрезков, разбиения их на диапазоны входного сигнала и вычисления отсчетов xl(n) как части числовых отрезков, попадающие в l-й диапазон входного сигнала, и
создают модель объекта идентификации в виде уравнения:
L
y' ( n ) = ∑ y l ( n ) ,
(2)
l=1
где у'(n) - отклик модели;
L - общее число диапазонов;
ml
y l ( n ) = ∑ a lj x l ( n − j),
(3)
j= 0
где ml+1 - количество дискретных отсчетов в l-й импульсной характеристике;
аlj - j-й отсчет l-й импульсной характеристики;
xl(n-j) - последовательность задержанных на j дискретных отсчетов частичного входного сигнала xl(n);
n - последовательность дискретных отсчетов входного сигнала с интервалом дискретизации не менее 16 отсчетов на верхнюю частоту в спектре входного сигнала.
На фиг. 1 представлена схема для реализации способа.
На фиг. 2 показан процесс разбиения на диапазоны минимально возможного Ymin и
максимально возможного Ymax откликов НДО.
На фиг. 3 иллюстрируется процесс формирования частичных входных сигналов xl(n)
из входного сигнала x(n) на примере двух его временных отсчетов x(3) и x(8).
3
BY 10526 C1 2008.04.30
Устройство содержит НДО 1 с датчиком 2 входного сигнала и датчиком 3 его отклика,
безынерционные нелинейные пружины 4, 5 и твердое тело 6 массой М, стенд 7 для задания входных сигналов, блок 8 АЦП, ЭВМ 9.
Выход ЭВМ 9 соединен с управляющим входом стенда 7, выход которого механически жестко соединен с датчиком 2 входного сигнала, выход которого электрически соединен с первым входом блока 8 АЦП, другой вход которого соединен с выходом датчика 3
отклика НДО 1, выход блока 8 АЦП через интерфейс соединен с ЭВМ 9. Первый конец
пружины 4 механически жестко связан с датчиком 2 входного сигнала, а другой - с датчиком 3 отклика НДО 1, первый конец пружины 5 соединен с верхней точкой твердого тела
6, а другой - с корпусом стенда 7.
Работа устройства происходит следующим образом. В ЭВМ 9 заносятся исходные
данные для идентификации НДО 1, т.е. количество диапазонов L разбиения его отклика, а
также высшая частота fв, до которой будет производиться идентификация.
На объект идентификации 1 с помощью стенда 7 подают случайный сигнал x(t) с достаточно большой дисперсией, такой, чтобы отклик объекта y(t) достигал как можно больших значений. Входной случайный сигнал x(t) регистрируется с помощью датчика 2, а
отклик y(t) - с помощью датчика 3. Из сигналов от датчиков 2 и 3 с помощью блока 8 АЦП
получают дискретные значения сигнала x(n) и отклика y(n), которые заносятся в ЭВМ 9. В
ЭВМ 9 определяются максимально возможные отрицательные Ymin и положительные Ymax
значения отклика НДО 1.
L
По программе, заложенной в ЭВМ, Ymin и Ymax разбивают на , где L - четное число,
2
равномерных диапазонов. Для этого сначала находят величины шагов для отрицательных
и положительных значений отклика НДО 1:
2Y
h ( −) = − min ,
(4)
L
2Ymax
(5)
h=
.
L
Затем вычисляют границы диапазонов Dr (r = 1, 2, ..., L,) разбиения отклика НДО по
следующим выражениям:
Ymin ≤ Dl < Ymin + h ( − ) ,
Ymin + h ( − ) ≤ D 2 < Ymin + 2h ( −) ,
−−−−−−−−−−−−−−−−−−
L
L
( )Ymin + ( − 1)h ( − ) ≤ D L < 0,
2
2
2
L
( + 1)h ≤ D L < h ,
+1
2
2
L
( + 2) h ≤ D L < 2 h ,
+2
2
2
−−−−−−−−−−−−−−−−−−
(L)(L − 1)h ≤ D L < Lh.
(6)
По программе, заложенной в ЭВМ 9, также находят середины диапазонов разбиения
отклика Yr (r = 1, 2, ..., L) НДО 1:
4
BY 10526 C1 2008.04.30
Y1 = Ymin + 0,5h ( −) ,
Y2 = Ymin + 1,5h ( −) ,
−−−−−−−−−−−−
L
YL = Ymin + ( − 0,5)h ( −) ,
2
(7)
2
−−−−−−−−−−−−
YL = 0,5h,
2
YL
2
+1
+2
= 1,5h
−−−−−−−−−−−−
L
YL = ( − 0,5)h.
2
Формирование границ диапазонов разбиения входного сигнала x(n) и рабочих точек
осуществляется путем подачи на объект 1 сигналов с постоянными мгновенными значениями, и находят среди них значения Xг, соответствующие границам диапазонов разбиения отклика y(n) (выражения 6), и значения Xr (r = 1, 2, ..., L), соответствующие серединам
этих диапазонов (выражения 7):
X г(1) → Ymin ,
X г( 2) → Ymin + h ( −) ,
−−−−−−−−−−−−
L
( )
Xг 2
L
→ Ymin + ( − 1)h (− ) ,
2
(8)
L
( +1)
X г 2 → 0,
L
( + 2)
→ h,
Xг 2
−−−−−−−−−−−−
X (гL ) → Ymax − h,
X (гL +1) → Ymax .
X1 → Y1 ,
X 2 → Y2 ,
(9)
−−−−−−−
X L → YL .
В выражениях (8), (9) стрелка означает знак соответствия, т.е. показывает, какому отклику НДО 1 с постоянными мгновенными значениями соответствует входной сигнал с
постоянными мгновенными значениями.
Вычисляют ширину диапазонов входного сигнала dr (r = 1, 2, ..., L), соответствующих
диапазонам разбиения Dr (r = 1, 2, ..., L) отклика НДО 1:
5
BY 10526 C1 2008.04.30
d1 = X г( 2) − X г(1) ,
d 2 = X г(3) − X г( 2) ,
−−−−−−−−−−
(10)
d L = X (гL +1) − X (гL ) .
Для определения размерности дискретного преобразования Фурье (ДПФ), используемого при идентификации объекта, сначала определяют длину импульсной характеристики
(ИХ) НДО 1. Для этого на вход НДО 1 с помощью стенда 7 подают небольшой испытательный сигнал Хи с постоянными мгновенными значениями и после установления отклика Yи резко убирают входной сигнал. Время от момента снятия сигнала до момента
полного затухания отклика принимают за длину m ИХ.
Затем по заданной высшей частоте, до которой будет производиться идентификация,
вычисляют значение интервала дискретизации ∆Т для дискретной модели НДО 1 по формуле:
1
∆T =
.
(11)
4f в
После этого вычисляют размерность N быстрого преобразования Фурье (БПФ) для
идентификации НДО 1:
2m
N≥
,
(12)
∆T
и округляют полученное значение N до большего числа из ряда чисел 2s (s = 2, 3...).
С целью получения коэффициентов передачи для корректировки частотных характеристик на нулевой частоте подают на вход НДО 1 сигнал с постоянными мгновенными
d
значениями Х1 + ∆X1, причем ∆X1 < 1 , и после окончания переходного процесса регист2
рируют отклик Yl + ∆Y1 НДО 1, и по нему и значению ∆Х1 определяют коэффициент передачи K (01) НДО 1 на нулевой частоте в рабочей точке, равной X1:
(Y + ∆Y1 ) − Y1 ∆Y1
K (01) = 1
=
.
(13)
∆X1
∆X1
Для получения суммарного сигнала сначала генерируют псевдослучайный процесс
εl(n):
4 N/2
2π
(14)
ε1 (n ) = A1
cos( nk + ϕk ),
∑
N k=1
N
где ϕk - случайная фаза, равномерно распределенная в диапазоне от 0 до 2π,
А1 - множитель, используемый для регулировки уровня случайного сигнала, который
должен быть таким, чтобы отклик НДО 1 от него не выходил за пределы диапазона D1,
при А1 = 1 дисперсия процесса ε1 (n ) равна единице,
n = 0 ÷ (2N-1).
Здесь ε1(n) - периодический процесс с периодом, равным N. В формуле (14) генерируются два периода этого процесса. Первый период, равный удвоенной длине ИХ, не рабочий (переходной процесс). Второй период используется для идентификации.
Затем подают на объект с помощью стенда 7 суммарный процесс Х1 + εl(n), регистрируют отклик НДО 1 на участке n = N ÷ (2N-1) и входной сигнал на том же участке.
Вычисляют частотную характеристику W1(к) НДО 1:
N − точечное БПФ отклика НДО 1
W1 (к ) =
.
(15)
N − точечное БПФ входного сигнала
6
BY 10526 C1 2008.04.30
Корректируют W1(к) путем замены значения нулевого отсчета a0 на значение K (01) и
тем самым получают W1c (к).
Вычисляют ИХ1 НДО 1 в рабочей точке, равной X1:
ИХ1 = ОБПФ от W1c (к),
где ОБПФ - обратное быстрое преобразование Фурье.
(16)
Подают на вход НДО 1 постоянный сигнал Х2 + ∆X2, причем ∆X 2 <
d2
, и получают
2
K (02 ) .
Вновь генерируют псевдослучайный процесс ε2(n):
4 N/2
2π
ε 2 (n ) = A 2
cos( nk + ϕ k ) ,
(17)
∑
N k =1
N
и так далее вплоть до получения ИХ2, а затем эту процедуру повторяют многократно и в
результате получают ИХ1, ИX2, ..., ИXL.
Формируют частичные входные сигналы xl(n) из входного сигнала x(n) путем представления его дискретных отсчетов в виде числовых отрезков, разбиения их на диапазоны
входного сигала и вычисления отсчетов xl(n) как части числовых отрезков, представляющих собой отсчеты x(n), попадающие в l-й диапазон входного сигнала. Для этого частичные входные сигналы xl(n) формируются из входного сигнала x(n) по следующему
правилу.
Для x(n) < 0.
Если отсчет x(n) попадает в l-й диапазон разбиения (выражения 10), то отсчеты частичных входных сигналов хl+1(n), хl+2(n), ..., x L (n ) принимают значение ширины своего
2
диапазона с обратным знаком, а отсчеты частичного входного сигнала в l-м диапазоне вычисляются из соотношения:
x l (n) = x(n) + (d l +1 + d l + 2 + ... + d L ),
(18)
2
где dl+1, dl+2, …, d L значения ширины соответствующего диапазона.
2
Для x(n) > 0.
Если отсчет x(n) попадает в l-й диапазон, то отсчеты частичных входных сигналов
xl-1(n), xl-2(n), ..., x L (n ) принимают значение ширины своего диапазона, а отсчеты час2
+1
тичного входного сигнала в l-м диапазоне вычисляются из соотношения:
x l (n) = x(n) - (d l -1 + d l - 2 + ... + d L ),
(19)
+1
2
где dl-1, dl-2, …, d L
2
+1
- значения ширины соответствующего диапазона.
Строят модель НДО 1 в соответствии с выражением:
L
y' (n ) = ∑ y l (n ),
(20)
l=1
где у'(n) - отклик модели,
уl(n) - свертка частичного входного сигнала xl(n) с l-й импульсной характеристикой
НДО 1,
y l (n ) =
ml
∑ a (jl) x l (n − j),
(21)
j= 0
7
BY 10526 C1 2008.04.30
a (jl) - j-й отсчет l-й импульсной характеристики,
xl(n-j) - последовательность задержанных на j дискретных отсчетов частичного входного сигнала xl(n),
n - последовательность дискретных отсчетов входного воздействия,
ml + 1 - количество дискретных отсчетов в l-й импульсной характеристике.
На этом процесс идентификации НДО 1 заканчивается.
Причем особо следует отметить, что входной сигнал на модель x(n), в отличие от теоремы Котельникова, должен быть продискретизирован с интервалом дискретизации не
менее 16 отсчетов на верхнюю частоту в спектре входного сигнала. Это обусловлено тем,
что данный способ требует того, чтобы огибающая дискретных отсчетов x(n) по форме
мало отличалась от непрерывного воздействия x(t).
Пример реализации способа
Пусть имеем компьютерный НДО 1, состоящий из последовательно соединенных первого линейного динамического объекта (ЛДО 1), безынерционной нелинейности (БН) и
второго линейного динамического объекта (ЛДО 2). При этом ЛДО 1 описывается следующим разностным уравнением:
z(n ) =
8
∑ a j x (n − j),
j= 0
(22)
где a0 = 1; a1 = 2; a2 = 3; a3 = 4; а4 = 5; а5 = 4; а6 = 3; a7 = 2; a8 = 1.
Безынерционная нелинейность выражается следующей функцией:
q(n) = arctgz(n),
а ЛДО 2 описывается следующим разностным уравнением:
y( n ) =
4
∑ b jq(n − j),
j= 0
где b0 = 1; b1 = 2; b2 = 3; b3 = 2,5; а4 = 1.
Пусть Ymin = -14,176 и Ymax = 14,176.
Зададим количество диапазонов L = 32.
В компьютерном объекте интервал дискретизации ∆T (формула 11) обычно берут равным единице, и поэтому для компьютерного объекта задавать fв для определения ∆Т не
надо.
Вычисляем величины шагов h(-), h для отрицательных и положительных значений отклика НДО по формулам (4), (5). Они оказались одинаковыми, т.к. модули Ymin и Ymax были взяты одинаковыми. При этом величины шагов равнялись 0,886.
Вычисляем диапазоны Dr (r = 1, 2, ..., L) разбиения отклика НДО по формуле (6). Ниже
приведены три первых из них:
1) -14,176 ≤ D1 < -13,29,
2) -13,29 ≤ D2 < -12,404,
3) -12,404 ≤ D2 < -11,518.
Вычисляем середины диапазонов разбиения отклика Yr (r = 1, 2, ..., L) НДО по формуле (7). Ниже приведены три первых из них:
Yl = -13,73; Y2 = -12,85; Y3 = -11,96.
Находим постоянные сигналы Хг, соответствующие границам диапазонов разбиения
отклика НДО (выражения 6), и постоянные сигналы Xr (r = 1, 2, ..., L), соответствующие
серединам этих диапазонов (выражения 7), некоторые из которых приведены ниже:
X (r1) = -0,5079; X (r 2 ) = -0,23; X (r3) = -0,1473; X (r 4 ) = -0,1067,
X1 = -0,318; Х2 = -0,18; Х3 = -0,124.
8
BY 10526 C1 2008.04.30
Вычисляем ширину диапазонов входного сигнала dr (r = 1, 2, ..., L), соответствующих
диапазонам разбиения Dr (r = 1, 2, ..., L) отклика НДО по формуле (10), некоторые из которых приведены ниже:
dl = 0,278; d2 = 0,083; d3 = 0,0405.
Как указывалось выше, в компьютерном объекте интервал дискретизации ∆T обычно
берут равным единице, и тогда длина импульсной характеристики будет измеряться количеством ее отсчетов. При этом вместо резкого перепада входного сигнала проще использовать одиночный импульс, отклик на который даст длину ИХ. Для нашего объекта длина
ИХ получилась равной 13. Тогда для идентификации НДО в соответствии с формулой (12)
минимальную размерность БПФ для идентификации НДО можно взять равной 32. Нами
для этой цели использовалось 64 точечное БПФ.
Определяем по формуле (13) коэффициенты передачи K 0(1) , K (02 ) , ..., K (0L ) НДО по постоянному входному воздействию в рабочих точках НДО, равных X1, X2, ..., XL, первые
три из которых приведены ниже:
K 0(1) = 3,015; K (02 ) = 9,98; K 0(3) = 20,717.
При этом значения ∆Y1, ∆Y2, ∆Y3 были равными -0,221, а значения ∆X1, ∆X2, ∆X3 соответственно следующими:
∆X1 = -0,073; ∆Х2 = -0,022; ∆X3 = -0,011.
Генерируем в соответствии с формулой (14) два периода псевдослучайных процессов.
При этом значение N было равно 64, а масштабные коэффициенты для первых трех процессов были следующими:
A1 = -0,037; А2 = -0,011; A3 = -0,005.
Проведя действия в соответствии с формулами (15), (16), получаем 32 импульсные характеристики НДО. При этом они оказались симметричными между собой, т.е.
ИХ1 = ИХ32, ИХ2 = ИХ31, ИХ3 = ИХ30, ..., ИХ16 = ИХ17. Ниже приведены некоторые из них.
ИХ2 = 0,0334; 0,169; 0,445; 0,8382; 1,2785; 1,6264; 1,789; 1,6717; 1,322; 0,8832; 0,4979;
0,2084; 0,0536.
ИХ4 = 0,139; 0,6025; 1,5215; 2,82; 4,2719; 5,4214; 5,962; 5,5826; 4,4311; 2,9704; 1,6618;
0,6594; 0,1171.
ИХ10 = 0,6455; 2,6528; 6,6836; 12,405; 18,81; 23,872; 26,237; 24,552; 19,492; 13,078;
7,3399; 2,9579; 0,6085.
По правилам (18), (19) формируем из входного сигнала x(n) частичные входные сигналы xl(n). Ниже приведен фрагмент входного сигнала x(n) и некоторые частичные входные сигналы xl(n).
x(n) = 0,0146; 0,0208; 0,0264; 0,0313; 0,0357; 0,0395; 0,0427; 0,0455; 0,048; 0,05; 0,0521;
0,0538; 0,0555; 0,057; 0,0584; 0,0598; 0,061; 0,062; 0,0627; 0,0631; 0,063; 0,0623; 0,0611;
0,0591; 0,0564; 0,0529; 0,0486; 0,0437; 0,038; 0,0318.
x25(n) = 0; 0; 0; 0; 0; 0,00254; 0,0058; 0,00759; 0,00759; 0,00759; 0,00759; 0,00759;
0,00759; 0,00759; 0,00759; 0,00759; 0,00759; 0,00759; 0,00759; 0,00759; 0,00759; 0,00759;
0,00759; 0,00759; 0,00759; 0,00759; 0,00759; 0,00759; 0,00672; 0,00107; 0.
x21(n) = 0; 0,00447; 0,00447; 0,00447; 0,00447; 0,00447; 0,00447; 0,00447; 0,00447;
0,00447; 0,00447; 0,00447; 0,00447; 0,00447; 0,00447; 0,00447; 0,00447; 0,00447; 0,00447;
0,00447; 0,00447; 0,00447; 0,00447; 0,00447; 0,00447; 0,00447; 0,00447; 0,00447; 0,00447;
0,00447.
x17(n) = 0,00374; 0,00374; 0,00374; 0,00374; 0,00374; 0,00374; 0,00374; 0,00374; 0,00374;
0,00374; 0,00374; 0,00374; 0,00374; 0,00374; 0,00374; 0,00374; 0,00374; 0,00374; 0,00374;
0,00374; 0,00374; 0,00374; 0,00374; 0,00374; 0,00374; 0,00374; 0,00374; 0,00374; 0,00374;
0,00374.
По импульсным характеристикам ИX1 ÷ ИХ32, и частичным входным сигналам xl(n),
используя формулы (20), (21), строим модель НДО. Подаем на вход НДО и его модели
9
BY 10526 C1 2008.04.30
один и тот же случайный процесс большой амплитуды, но в пределах 32 диапазонов x(n).
Сравниваем отклики НДО и модели. Среднеквадратическая погрешность идентификации
при этом не превысила инструментальной погрешности ЭВМ.
Таким образом, предложенный способ идентификации НДО оказывается существенно
точнее известного за счет получения L импульсных характеристик НДО в различных рабочих точках, определения из входного сигнала x(n) L частичных входных сигналов,
свертки их с соответствующими ИХ и последующим суммированием по l (l = 1 ÷ L).
Источники информации:
1. SU1742831 А1, 1992.
2. SU 290282, 1971.
Фиг. 2
Фиг. 3
Национальный центр интеллектуальной собственности.
220034, г. Минск, ул. Козлова, 20.
10
Документ
Категория
Без категории
Просмотров
0
Размер файла
150 Кб
Теги
by10526, патент
1/--страниц
Пожаловаться на содержимое документа