![](/user_photo/_userpic.png)
книги из ГПНТБ / Кавалеров, Г. И. Введение в информационную теорию измерений
.pdfрактеристик. Некоторые частные варианты оптимизации были рассмотрены ранее три изложении материала гл. 2—5. Ниже излагается более общая 'постановка зада чи параметрического синтеза ИИС с известной структу рой.
Применительно к решению задачи измерения опре деленной совокупности величин в известных условиях выбирается один или несколько вариантов структурной схемы и 'ключевая характеристика, обеспечение которой особенно важно в данных условиях. Примером ключевой характеристики может быть объем аппаратуры в литрах для некоторых задач специальной техники, себестои мость массовых приборов и т. д. Основные метрологиче ские характеристики отдельных блоков, из которых воз можна комплектация различных вариантов системы в пределах выбранной структурной схемы, также можно считать заданными, так как они выбираются из нормаль ных рядов агрегатированных средств, о которых будет сказано ниже. Кроме того, требуется знание зависимости ключевой характеристики от принятых показателей точ ности. Задача синтеза при этом формулируется либо как минимизация (или максимизация) ключевой характери стики при ограничении на результирующую погрешность, либо как задача минимизации погрешности при задан ных ограничениях на значение ключевой характеристики.
Из сказанного вытекает возможность решения сле дующей частной задачи оптимального синтеза ИИС. Пусть структура системы определена с точностью до ти пов отдельных звеньев и их взаимного расположения. Такая определенность, как правило, имеется у опытного проектировщика после получения технического задания. Требуется оптимально распределить требования к ме трологическим характеристикам отдельных блоков
водном из двух вариантов:
1)минимизируя затраты на изготовление ИИС (клю чевая характеристика) при заданном ограничении на точность системы в целом в смысле принятого метроло гического показателя качества измерений;
2)минимизируя суммарную погрешность (ключевая характеристика — метрологический показатель качества всей РІИС) при заданном ограничении на затраты.
Рассмотрим решение этой задачи применительно к из мерительной информационной системе с одним входом и одним выходом,
?33
Итак, рассматривается ИИС, состоящая из іі звеньев, каждое из которых описывается реальным оператором Л(г)р; і = 1, .. ., п. В связи с отличием реальных операто ров от идеальных всеми звеньями вносятся погрешности. Обозначим погрешность, вносимую t-м звеном в услови ях нормировки погрешностей, через уь г=1, .. ., п, истин ное значение измеряемой величины (на входе системы) через X и идеальный оператор всей системы через Ап- Тогда желательный сигнал на выходе системы
v(,t)=Aü[x(t)], |
(6-1) |
а результирующая погрешность |
|
y { t)= z {t)—v{4), |
(6-2) |
где z(t) — реальный сигнал на выходе системы.
Для решения задачи необходимо знание зависимости стоимости конкретных видов звеньев от принятого метро логического показателя качества. Зададим эти зависи мости в виде Сі= С,ДТі (У,)], где МЛ—принятый норми руемый метрологический показатель качества t-го звена системы или совокупность метрологических характери
стик этого звена; С,- — функция, зависящая |
от типа зве |
на. Суммарная стоимость такой системы |
|
с , = £ сд'УДУ/)], |
(6-3) |
;=і |
|
т. е. при прочих равных условиях зависит от распреде ления погрешностей между звеньями системы.
В то же время «суммарная» точность системы |
|
V s= WQPl,Wt,...,Wn,X). |
(6-4) |
Последняя запись отражает тот факт, что погреш ность уі, которой аттестуется звено, не совпадает с по грешностью, вносимой этим звеном в реальных условиях измерений. Последняя в первую очередь зависит от x(t), но связана также с величиной Уі(і) определенными со отношениями.
В рамках первого варианта поставленной выше зада чи можно минимизировать (6-3) при заданном ограни чении сверху на функцию от (6-2), т. е. искать
С * = min S Сі Г ЗД )]. |
(6-5) |
{ч'(П )}€'-г=І |
|
224
Здесь Li — область, в которой выполняется неравенство
^ (T b |
Ѵя, X * ) < ¥дод, |
(6-6) |
где 'Гдоп — допустимое значение нормируемого метроло гического показателя; х* — истинное значение измеряе мой величины, соответствующее условиям нормирования и поверки всех звеньев и системы в целом.
Второй вариант задачи предполагает создание систе мы с максимальной точностью при заданном ограниче нии на стоимость, т. е. отыскание
W* |
= min ЧГ(У„ У., .... Уя, X*). |
(6-7) |
В данном случае L2— область, в которой выполняет |
||
ся неравенство |
|
|
|
ЕС<[ЧГ*(У*)]<СД0П, |
(6-8) |
где Сдои— максимально допусти,мое значение стоимости системы.
Исследование конкретных вариантов схем ИИС по зволяет определить вид зависимостей (6-3) и (6-4), а также область допустимых решений Li и Ь2. Для вы пуклых областей Li и L2 при условии сепарабельности ключевой характеристики может быть использован ме тод динамического программирования. При условии, что функции, характеризующие ключевые характеристики, являются выпуклыми, возможно применение итерацион ных методов поиска минимума. Однако эти вопросы вы ходят за рамки настоящей работы. Некоторые примеры приведены в [Л. 6-2].
Кроме того, в качестве примера применения изложен ной методики рассмотрим следующий простейший слу чай. Пусть система состоит из п блоков, каждый из ко торых характеризуется коэффициентом усиления а, и аддитивной погрешностью ij-u приведенной к выходу дан ного преобразователя. Тогда дисперсия погрешности на выходе всей системы равна:
а2. |
п |
Г |
п |
\ |
п— 1 |
/ п |
\ |
|
£< |
П < ) |
+ * 2 |
п НХ |
|||||
■^вых |
||||||||
|
/=і |
\/=; +1 |
) |
А=1 |
\/=ft+i |
) |
||
|
|
X t |
\ |
ві |
( П aJ |
|
15-301 |
225 |
Пусть далее стоимость і-го блока в зависимости от
овыразится формулой
Сі — Л' ехр (— Bp*).
Тогда в случае некоррелированности у,- н уи при ■іфіг имеем оптимальное распределение погрешностей между блоками как решение следующей системы урав нений:
2 2
где Я — неопределенный множитель Лагранжа; а2ДОп — допустимое значение дисперсии погрешности на выходе системы.
Таким образом, изложенный подход позволяет гово рить об оптимизации измерительной системы с учетом экономико-метрологических факторов. Одним из путей дальнейшего развития описанной задачи является сведе ние требований по точности также к экономическому фактору путем учета функции потерь потребления. По добная постановка возможна для электроизмерительных систем, в отношении которых заранее точно определен объект .применения. Однако эти вопросы также выходят за рамки настоящего исследования.
Другим важным аспектом создания и использования агрегатироваиных средств электрических измерений является обоснованный выбор рядов типоиомипалов раз личных групп приборов и устройств. В связи с зависи мостью себестоимости от серийности излишне мелкое дробление невыгодно, так как при этом, хотя и имеет место удовлетворение потребности заказчиков по но менклатуре в соответствии с запросами, каждая серия очень мала, а число серий велико.
С другой стороны, малое число номиналов укрупняет серии, т. е. выгодно заводам, но заставляет потреби теля выбирать приборы с «запасом» по параметрам. Другими словами, если потребителю требуется, напри мер, некоторое быстродействие, не соответствующее вы пускаемым номиналам, то он должен выбрать прибор с большим быстродействием (ближайший номинал), т. е. купить более дорогой прибор. Метод расчета оптимума с учетом специфических особенностей измерительной т?х-
226
Никй |
отличается от принятого в смежных областях |
[Л. 6-3, |
6-4]. |
Формализуем задачу выбора гипономиналов следую щим образом. Охарактеризуем каждый прибор совокуп ностью п номиналов и. Компоненты вектора и есть ос новные характеристики (параметры) прибора, из числа которых нужно в первую очередь выделить точность, бы стродействие и пределы измерения. Совокупность всех векторов и порождает некоторое /і-мерное пространство Еп, в котором для каждого конкретного вида приборов можно оговорить из 'физических соображений область функционирования WnczEn. В этой области каждому типономиналу us соответствует подобласть WVsczWn.
Если возможна замена приборов с типономиналом u.s приборами с типономиналом и,., то имеется возможность сформулировать задачу унификации. Под рядом типономииалов номенклатурной группы будем понимать лю
бой набор типономиналов |
{«.,}, |
5= |
1, 2, . . . |
такой, что |
U Wnu zdWn. Предполагается, |
что |
на. множестве пара- |
||
s . |
распределение |
потребности |
||
метров задано «-мерное |
в приборах wn(ui, Uz, . .., ип) по параметрам и. Удовлет ворение потребности возможно при различных разбиени ях области Wn, однако эффективность различных раз биений неодинакова. Таким образохм необходимо задать ся критерием эффективности и выбрать такое разбиение, которое соответствует экстремальному значению крите рия.
Пусть выбраны значения параметров { и t = l, 2, . ..,
.. ., п. Тогда критерий оценки качества разбиения Е бу дет некоторым функционалохм вида
Е[т(иь . . ип), {Ui}\, |
(6-9) |
причем плотность вероятностей w(ui, ..., ип) не зависит от выбранного разбиения. Наибольшие трудности пред ставляет выбор вида функционала Е. Наибольшей общ ностью, по-видимому, обладает запись в форме, которую мы -поясним на примёре одномерной задачи
ті
Е== 2 Yz'jC (Н-г j +1, Тг'з)~Ь
/=1
ті "і, i+i |
— Ui)w(ui)dui, |
(6-10) |
+S j |
15* |
227 |
где щ — г'-й |
параметр; |
ііц — нижняя |
граница |
/-й |
зоны |
|||
(подобласти |
разбиения) |
jio |
параметру |
w(Uj) — рас- |
||||
пределение |
потребности |
по |
щ; |
4 a = |
L |
"i,j+1 |
|
|
| |
w{ui)dui — |
|||||||
число приборов, входящих в |
/-іо |
|
“іІ |
по |
пара |
|||
зону |
||||||||
метру г; Ші — число подобластей |
(зон) |
разбиения по |
||||||
параметру и,-; С(щ-, ,-+ь |
гг,-,-)— стоимость |
в рублях |
всей |
годовой .партии приборов, входящих в у'-ю зону параме тра і; <p(Uitj+i—Ui) —-функция потерь от несоответствия
номинального параметра |
|
3-+і и текущих значений по |
|||
требности по |
параметру |
щ; L — общее |
количество |
||
средств измерения, выпускаемых подотраслью в год. |
|||||
Нетрудно видеть, что |
первое слагаемое |
-в |
формуле |
||
(6-10) отражает |
затраты |
на |
производство, |
а |
второе — |
ущерб от «запасов» потребителей приборов по параме тру щ. Напомним, что имеется в виду под запасами. До пустим, потребителю требуется некоторое значение нор мируемого показателя точности, которого нет в ряду типономиналов. Тогда он вынужден покупать более точ ное средство измерения, т. е. иметь ненужный ему запас по точности, за который приходится переплачивать. По добная же ситуация возможна по быстродействию и дру гим параметрам.
В многомерном случае формула (6-10) принимает вид:
п т
|
|
|
Xi |
Yz'jC1(ui,j -ц, Yi'j) + |
|
|
|
|
«=i ;=1 |
|
|
|
|
|
n |
m t |
|
|
|
|
+ |
S |
E |
l («X. j +, — Щ ) W |
(U i) d u i , |
(6-11) |
|
где |
/=I |
/=> |
|
|
|
|
|
|
|
|
|
|
|
“ l . i + l " 2 . І + 1 |
|
" n . J + I |
|
|
||
Тгj L j |
I |
... |
j |
wn (ul, и |
un)du1du2...dun. |
|
" l . J |
" 2 . І |
|
|
|
|
|
(6- 12)
К сожалению, в данных формулах имеется по край ней мере одна величина, в общем виде не поддающаяся оценке. Действительно, вид функции ущерба ф(«,-, 3-+і— Ui) не может быть определен как в большинстве кон кретных применений измерительных приборов, так и для средств измерения универсального назначения.
228
Поэтому возможен следующий подход. В качестве критерия качества разбиения используется первое сла гаемое в формуле (6-11). Требуется минимизировать этот функционал. Так как согласно пятилетнему плану заданы контрольные цифры (в стоимостном выражении), то, получив в результате минимизации по всем разбие ниям некоторую цифру суммарной стоимости всех при боров, выпущенных подотраслью, в виде
пт
C j^rn in £ |
S |
С6' 13) |
{“t*} i=i |
i=i |
|
следует изменением величины L в результате повторных |
||
пересчетов получить СХ= |
С0, |
где С0— контрольная го |
довая цифра на конец пятилетки в рублях.
Следует иметь в виду, что существуют так называе мые ряды предпочтительных чисел. Поэтому практиче ски задача разбиения сводится к отысканию разрешен ных и запрещенных значений в этих рядах. Кроме того, все исходные характеристики могут бъіть заданы только в численной форме. Отсюда естественно вытекает целе сообразность численного решения поставленной задачи с помощью ЭЦВМ. Именно по этому пути пошли при на значении рядов средств АСЭТ (см. гл. 1). В результате этой работы была разработана методика назначения типономиналов средств АСЭТ, отдельные аспекты кото рой изложены в [Л. 6-5, 6-6].
Изложенные задачи не исчерпывают возможные ва рианты оптимизации. Имеется целый ряд других мето дов, позволяющих найти оптимальные или квазиопти мальные решения в ряде частных аспектов построения измерительной аппаратуры, к рассмотрению которых мы и перейдем в дальнейших параграфах этой главы.
6-2. ОПТИМАЛЬНОЕ КВАНТОВАНИЕ
Рассмотрим прежде всего частную задачу оптими зации ■—выбор шкалы 'Квантования по уровню, обеспе чивающей наилучшее в определенном смысле согласова ние квантователя с измеряемой величиной.
Задача оптимального квантования по уровню может быть сформулирована следующим образом: пусть из-
229
ізестны вероятностные свойства сигнала, подлежащего квантованию; требуется осуществить его квантование по уровню при заданном числе N уровнен 'квантования та ким образом, чтобы 'построенный квантователь был наплучшим в смысле некоторого точностного критерия сре ди всех квантователей с N уровнями.
Впервые эту задачу начали изучать специалисты по теории связи более десяти лет тому назад '(см., напри мер, [Л. 6-7—6-12)].
Побудительной причиной такого исследования яви лось развитие импульсных методов передачи сообщений. Оптимальное квантование по уровню может рассматри ваться как одни из способов сжатия данных, т. е. сокраще ния полосы частот, необходимой для их передачи. Действи тельно, если сравнить некоторый оптимальный квантова тель е 'квантователем с постоянным шагом, обеспечиваю щим ту же точность, то очевидно, что число уровней последнего превосходит число уровней оптимального квантователя. Поэтому (во всяком случае, если не при меняется статистическое кодирование) можно получить некоторую «экономию», тратя на передачу одного отсче та, взятого на оптимальном квантователе, меньшее число двоичных единиц по сравнению со случаем равномерно го квантования по уровню.
Указанный выше факт — возможность более компакт ного представления сигнала при заданной точности вос произведения— лежит в основе применения методов оптимального квантования в измерительной практике.
Теория, излагаемая в этом параграфе, не может ис пользоваться при создании измерительных приборов уни версального назначения: ее предпосылкой является предположение, что «вид» измеряемого сигнала априо ри известен. В то же время построение специализирован ных АЦП, входящих в измерительную систему, постоян но используемую в одинаковых условиях с одинаковой целью, зачастую может быть осуществлено с помощью
этой теории (см., |
например, {Л. 6-13]). Важным «благо |
приятствующим» |
фактором для использования опти |
мальных способов |
квантования является наличие ЭВМ |
в измерительной системе. Это позволяет повысить точ ность восстановления квантуемого сигнала оптимальным выбором «точек отнесения» внутри кванта, не усложняя конструкции АЦП: в вычислительную машину вводится только код номера кванта, в который попал измеряемый
230
сигнал, вы'бор же точки отнесения производится программой воостановлейия.
Приведем точную формулировку одной из задач тео рии оптимального 'Квантования. Пусть фиксирован не который критерий точности, основанный на усреднении функции штрафов ср(х, z), т. е. качество устройства, пре образующего случайный сигнал X в случайный же сиг
нал Z, оценивается величиной ср(А', Z). Будем предпола гать, что функции ф(х, z) удовлетворяет условиям, накладываемым на функции штрафов (см. гл. 1) .Требует ся по заданному числу уровней квантования N и одно мерной плотности вероятностей квантуемого сигнала w і ( х ) определить параметры квантователя с N уровня
ми, минимизирующие <р(Х, Z), т. е. |
найти вектор (лщ .. .. |
||
..., xN- ü z 1........Zn), |
доставляющий |
минимум |
функции |
Ч; (хі, |
..., A'jv-i; Zi, ..., |
zN) = |
|
N |
A'h |
|
(6-14) |
= 2 |
[ T (X, zk) w, (X) dx. |
Здесь хи-—'пороги квантования; zk — точки отнесения;
.Vo — минимальное, а хк — максимальное значение изме ряемой величины, т. е. если значение квантуемого сигна ла лежит между порогами хк-і и Xh, то значение кванто ванного сигнала есть zk-
Как показано в работах А. И. Величкина [Л. 6-14, 6-15], эта задача, кажущаяся весьма частной, является основной в данной теории. Так, рассмотрение более об щих ситуаций (учет запаздывания в квантователе, нали чия помех, отличных от шума квантования по уровню) сводится к минимизации функций вида (6-14). Мы вкратце рассмотрим эти результаты несколько ниже.
Прежде чем перейти к изложению способов решения задачи оптимизации, обсудим один вопрос, связанный с возможностью практического применения методов оп тимального квантования. Априорные данные об измеряе мом сигнале получаются, как правило, эксперименталь но. Естественно, они не вполне точны. Поэтому необходи
мо |
исследовать |
корректность рассматриваемой задачи, |
т. |
е. выяснить, |
непрерывно ли зависит ее решение от |
априорных данных (ниже мы встретимся еще с одной ситуацией, в которой важно знать ответ на этот вопрос).
231.
Введем следующие обозначения:
|
|
N |
|
|
|
Ф '1'<**>(«,) = |
I] |
9 (х, z{) да (х) dx |
|
||
= |
¥ (* „ |
І — 1 X 1-1 |
Z^, |
_2д,), |
(6-15) |
. JC. |
|||||
|
|
’N-l' |
|
|
|
R (да) = |
min r £ ‘>- {Zi} (да). |
|
|||
|
{*<}.{*<} |
|
|
|
|
Функционал r j ^ ’ |
каждой |
плотности |
вероят |
||
ностей да (x), |
сконцентрированной |
на отрезке |
[х0, Хд,], |
ставит в соответствие величину среднего штрафа, соот ветствующего способу квантования с ворогами {х,} и точками отнесения {г,}, а функционал R x ( w ) — величину
среднего |
штрафа при оптимальном квантовании. |
|||||
В |[Л. 6-13] показано, что функционал |
R n (w ) |
является |
||||
непрерывным в том смысле, что для любого |
еі>0 най |
|||||
дется б2>0 такое, |
что если |
|
|
|
||
|
|
XN |
I да, (x) — Ш, (x) |
I dx < |
|
|
|
|
j |
e2, |
(6-16) |
||
|
|
XI. |
|
|
|
|
то |
|і?іѵ(ші)—Rx(wi) ] <ei, |
|
|
(6-17) |
||
где |
Ei и 62— любые произвольно |
малые положительные |
||||
числа. |
|
пусть |
|
|
|
|
|
Действительно, |
|
|
|
||
|
|
|
Ь0— max 9 (x,z). |
|
(6-18) |
ге[*„. xn\.
Положим, E2=ei/Lo. Тогда для любого способа кван тования с N уровнями
Nх\
^J ?(Х, Zn) [да, (х) — да, (л:)] dx
6=1 XН-1 |
|
|
N |
х ь |
|
< А> J J |
Г I |
(Х) — Щ (х) I dx - |
XN |
|
да2(х) I dx < L0 = в,. (6-19) |
—Z-o j* I О», (л:) - |
||
x0 |
|
|
232