
- •Оглавление
- •Глава 1. Оценивание контролируемых параметров по экспериментальным данным 12
- •Глава 2. Критерий качества задачи оценивания параметра 31
- •Глава 3. Выбор алгоритма обработки экспериментальных данных в автоматизированных системах управления и анализ их свойств 63
- •Глава 4. Формирование модели измерения задачи оценивания по экспериментальным данным 109
- •Глава 5. Формирование модели измерения в задаче оценивания параметра по экспериментальным данным на начальном этапе разработки нового электрооборудования 150
- •Глава 6. Создание устройств формирования модели измерения 176
- •Введение
- •Глава 1. Оценивание контролируемых параметров по экспериментальным данным
- •1.1. Общий анализ этапов структурирования эмпирической информации
- •1.2. Концептуальная модель процесса оценивания контролируемых параметров
- •1.3. Совершенствование алгоритмов обработки экспериментальных данных
- •1.4. Современная концепция оценивания измеряемого параметра
- •1.5. Классификация результатов измерения по характеру оценивания погрешности
- •1.6. Концептуальная модель задачи алгоритмизации оценивания результатов измерения
- •1.7. Ретроспективный анализ этапов развития теории оценивания
- •Контрольные вопросы
- •Глава 2. Критерий качества задачи оценивания параметра
- •2.1. Определение качества задачи оценивания измеряемого параметра
- •2.2. Формирование критерия качества гомоморфной математической модели измерения
- •2.3. Информационная мера степени изоморфности модели
- •2.4. Расчет информационной меры изоморфности
- •2.5. Оценка информационного объема и риска модели измерения
- •2.6. Асимптотическое оценивание пропускной способности математической модели измерения
- •2.7. Асимптотический метод выделения признаков модели измерения
- •Типовые задачи и примеры их решения
- •Контрольные вопросы
- •Глава 3. Выбор алгоритма обработки экспериментальных данных в автоматизированных системах управления и анализ их свойств
- •3.1. Особенности формирования алгоритмов оценивания в автоматизированных системах управления
- •3.2. Общий анализ алгоритмов оценивания по критерию минимума риска
- •3.3. Общий алгоритм оценки измеряемого параметра
- •3.4. Оптимальный одношаговый алгоритм
- •3.5. Модификации алгоритма обработки экспериментальных данных
- •3.6. Моделирование алгоритма обработки экспериментальных данных
- •3.7. Исследование сходимости алгоритма
- •3.8. Определение весовых коэффициентов алгоритма
- •3.9. Определение начальных условий алгоритма Язвинского при оценке результатов измерения
- •Типовые задачи и примеры их решения
- •Контрольные вопросы
- •Глава 4. Формирование модели измерения задачи оценивания по экспериментальным данным
- •4.1. Анализ задачи формирования модели измерения
- •4.2. Принципы построения модели измерения параметра контролируемого объекта
- •4.3. Этапы решения задачи формирования модели измерения
- •4.4. Общая постановка задачи формирования модели измерения
- •4.5. Выбор критерия близости
- •4.6. Способы преобразования переменных модели измерения
- •4.7. Общий анализ формирования модели контролируемого объекта
- •4.8. Решение задачи формирования модели контролируемого объекта
- •4.9. Оптимизация алфавита классов и словаря признаков
- •4.10. Взаимосвязь размерности алфавита классов и качества и эффективности модели измерения
- •4.11. Взаимосвязь размерности вектора признаков и вероятности правильности формирования модели измерения
- •4.12. Формализация задачи оптимального взаимосвязанного выбора алфавита классов и словаря признаков
- •4.13. Формирование оптимального алфавита классов и словаря признаков в условиях ограничений
- •Контрольные вопросы
- •Глава 5. Формирование модели измерения в задаче оценивания параметра по экспериментальным данным на начальном этапе разработки нового электрооборудования
- •5.1. Исследование условий формирования модели измерения задачи оценивания измеряемого параметра
- •5.2. Геометрический способ формирования модели измерения
- •5.3. Лингвистический алгоритм формирования модели измерения на начальных этапах оценивания
- •5.4. Обоснование выбора критерия расхождения для формирования модели измерения по экспериментальным данным
- •5.5. Метод формирования модели измерения
- •5.7. Синтез метода и алгоритма формирования стратегии постановки начальных экспериментов
- •Типовые задачи и примеры их решения
- •6. Найденный квант исключается из множества г, т.Е.
- •Контрольные вопросы
- •Глава 6. Создание устройств формирования модели измерения
- •6.1. Состав устройства формирования модели измерения
- •6.2. Общий анализ устройств формирования моделей измерения
- •6.3. Общий анализ моделирования устройств формирования модели измерения
- •6.4. Структура устройства формирования модели измерения
- •6.5. Моделирование контролируемого объекта
- •6.6. Моделирование средств измерения параметров контролируемого объекта
- •6.7. Моделирование каналов измерения
- •Типовые задачи и примеры их решения
- •6.8. Модель алгоритма формирования модели измерения
- •6.9. Модуль оценки качества и эффективности устройства формирования
- •6.10. Модуль управления моделью устройства формирования
- •6.11. Использование принципов опытно-теоретического метода при моделировании устройства формирования
- •6.12. Моделирование в задачах создания и оптимизации устройства формирования
- •Заключение
- •Библиографический список
- •Системы реального времени Синтез алгоритмов оценивания технологического параметра
- •346428, Г. Новочеркасск, ул. Просвещения, 132
Типовые задачи и примеры их решения
Задача 3.1. Синтезировать алгоритм обработки сигнала для определения места повреждения высоковольтных кабельных линий.
Описание предметной области. Одним из путей повышения эксплуатационной надежности электрических сетей является быстрое выявление и восстановление поврежденных участков линий электропередачи. Протяженность распределительных сетей 0,4 – 35 кВ достигает нескольких миллионов километров. В условиях промышленных предприятий и городов до 80% электроэнергии доставляется потребителям по кабельным линиям. По данным эксплуатации кабельные линии повреждаются в 2–3 раза чаще, чем другие элементы цепи электроснабжения. Ежегодно в кабельных сетях отыскивается около полумиллиона мест повреждения. Время ликвидации таких повреждений составляет от нескольких часов до нескольких суток.
В комплексе мероприятий, связанных с выявлением повреждения кабельной линии, поиск места повреждения непосредственно на трассе линии остается одной из наиболее трудоемких операций. Большое разнообразие видов повреждений (замыкание одной, двух, трех жил между собой и на оболочку, обрыв жил, заплывающие пробои) и параметров поврежденной кабельной линии (переходное сопротивление, длина линии до места повреждения и за ним, марка кабеля и глубина его залегания) вызвало создание различных методов и устройств для непосредственного отыскания места повреждения на трассе кабельной линии. Однако ряд повреждений, прежде всего металлические однофазные замыкания и заплывающие пробои при неразрешенной оболочке, а также повреждения в кабелях с отдельным металлическим покрытием жил остаются еще трудно выявляемыми и требуют дополнительных теоретических оценок и получения на их основе новых технических решений как в части разработки новых способов поиска повреждений, обладающих большей разрешающей способностью в различных ситуациях, так и в части совершенствования аппаратуры.
Для совершенствования поисковой аппаратуры следует иметь необходимый расчетный аппарат при выборе параметров поисковой аппаратуры с учетом характеристик поврежденной кабельной линии, оценить помехозащищенность, чувствительность и в конечном итоге – надежность выявления повреждения и массогабаритные показатели аппаратуры, что позволяет повысить эффективность работы оператора и в целом снизить затраты и время на поиск места повреждения кабельной линии.
Примером поискового устройства являются акустические приемники, которые выпускаются в составе универсальных серийных кабелеискателей и предназначенные как для индукционного, так и для акустического метода определения места повреждения кабельных линий [В. В. Платонов, В. Ф. Быкадоров. Определение мест повреждения на трассах кабельной линии. – М.: Энергоатомиздат, 1993. – 256 с.]. Эти кабелеискатели различаются в основном типом акустического датчика, поэтому рассмотрим особенности выполнения кабелеискателей "Поиск" П-6 с электродинамическим микрофоном и КАИ-77 (КАИ-80) с пьезоэлектрическим датчиком и с цифровым измерителем расстояния.
Кабелеискатель "Поиск" содержит акустический и индукционный датчики с усилителями, которые через переключатель вида поиска могут подключаться к выходному усилителю с телефоном. В качестве акустического датчика может быть использован виброизмерительный пьезоэлектрический преобразователь типа ДН–3–М1, анализ которого приведен ниже в задаче 3.3, а в качестве индукционного датчика – катушка с сердечником, настроенная в резонанс на поисковую частоту.
Кабелеискатель КАИ-77 содержит пьезоэлектрический датчик со встроенным усилителем и встроенную в корпус прибора индукционную катушку, которые через усилитель приемника подключаются к головным телефонам. Предусмотрен выход на стрелочный индикатор. Кабелеискатель имеет ограниченную чувствительность по акустическому каналу.
Повышенной избирательностью к заданному диапазону частот обладает кабелеискатель с электродинамическим датчиком и электромеханической положительной обратной связью (рис 3.2). Подвижный магнит NiSn, подвешенный в поле неподвижного магнита N3S3, при колебаниях грунта наводит в неподвижных катушках L1 и L2 противоположные по знаку ЭДС. Благодаря встречному включению катушек, полезный сигнал удваивается, а сигнал помехи вычитается. Суммарный полезный сигнал через полосовой фильтр Ф1 и основной усилитель У1 посту
|
Рис. 3.2. Структурная схема избирательного кабелеискателя с электродинамическим датчиком и электромеханической положительной обратной связью |
Индукционно-акустический кабелеискатель с цифровым измерителем расстояния позволяет быстро локализовать место повреждения за счет определения расстояния до источника электроакустических колебаний. Принцип работы измерителя основан на определении времени запаздывания акустического сигнала относительно индукционного сигнала, которое практически пропорционально расстоянию от оператора до места повреждения. Структурная схема кабелеискателя (рис. 3.3) – двухканальная, но сигнал индукционного канала используется в качестве опорного не только для измерителя расстояния», но и для контроля трассы и наличия разрядов.
|
Рис. 3.3. Структурная схема индукционо-акустического кабелеискателя ПА-02 |
Поиск места повреждения на трассе ведется следующим образом. Генератор акустических сигналов подключается к поврежденной кабельной линии. Оператор входит в зону поиска (измеритель расстояния выключен) и по индукционному каналу фиксирует трассу кабеля. Перемещаясь по трассе и периодически устанавливая акустический датчик на грунт, оператор находит зону приема акустического сигнала с помощью головных телефонов. Одновременно возникает синхронное мигание светодиодов индукционного и акустического каналов, что свидетельствует о возможности включения измерителя расстояния. Показания двухразрядного индикатора (от О,9 до 9,9 м) позволяют вести ускоренный целенаправленный поиск. Над местом повреждения прибор показывает глубину залегания кабеля.
Решение задачи
1. Системный анализ решаемой задачи и этап формализации постановки решаемой задачи. При традиционном подходе к задаче оценивания принимают соответствующие гипотезы о приближенных нормальности, независимости, однородности шумов, сопровождающих процесс измерения, и т.д.. Это позволяет при выборе алгоритма обработки информации использовать метод максимального правдоподобия [1], обеспечивающий минимизацию риска оценки результатов измерения.
Так для выборки Y = (Y1,..., Yn), каждый элемент Yi которой имеет плотность распределения f(yi; x), x x, плотность распределения выборки Y определяется выражением
f(y;
x)=
(yi;
x),
y
= (y1,...,
yn)Rn,
xx,
а функция правдоподобия данного вектора Y имеет вид
(Y; x) = f(Y; x) = (yi; x), xx.
При
регулярном семействе F
= {f(y;
x),
x
x}
для нахождения оценки максимального
правдоподобия
(Y)
пользуются стандартными приемами
дифференциального исчисления, решая
уравнение правдоподобия
(Y; x) = 0.
При этом всякий математический результат, предназначенный для непосредственного практического использования, должен обладать свойством устойчивости по отношению к небольшим отклонениям от исходных предположений. Однако в ряде случаев даже небольшие вариации исходных предпосылок о вероятностной структуре контролируемой информации могут привести к существенному ухудшению свойств используемых процедур. Поэтому необходимо постоянно оценивать условия измерения и осуществлять адаптацию к ним.
При измерении влияние одной из компонент может настолько перевесить суммарный вклад других компонент в ln (Y; x), что ее значение будет полностью определять величину оценки. Сам по себе этот эффект не вызывает беспокойства, если модель для плотности вероятности верна, однако он может возникать и при случайном ошибочном измерении, который не описывается принятой плотностью распределения вероятности f, то есть в общем случае f задает только приближенную исходную модель. Характеристики номинально оптимальных процедур обработки сигналов во многих случаях могут резко ухудшаться, даже при сравнительно малых отклонениях от исходных данных. В этом состоит основная причина, по которой приходится строить алгоритмы обработки сигналов, которые характеризуются свойством робастности [7, 64], то есть алгоритмы, обладающие высокой эффективностью в условиях, когда свойства сигналов и шума могут меняться в пределах заданных классов возможных характеристик.
Наличие априорной неопределенности относительно статистических характеристик информации измерения приводит к необходимости построения методов обработки, обладающих пониженной чувствительностью к всевозможной неадекватности между априорно принятой вероятностной моделью и свойствами реальной среды, то есть робастностью.
Важнейший подход к проблеме синтеза робастных методов обработки, обеспечивающих достаточно точное решение подобных задач в условиях неопределенности, дает принцип минимакса. Минимаксный подход требует заранее ориентироваться на ситуацию, наименее благоприятную для задачи оценивания, т.е. обеспечивает получение некоторого гарантированного решения. Сущность подхода состоит в том, что, сохраняя высокую устойчивость найденных решений, он позволяет избежать значительного снижения точности оценивания в случае, когда условия оценивания близки к предполагаемым. Один из наиболее эффективных путей построения минимаксных алгоритмов дает метод максимального правдоподобия (или метод М-оценивания), впервые предложенный П. Хьюбером в 1964г [7].
Для
того, чтобы как - то уменьшить нежелательную
чувствительность критерия, основанного
на (Y),
вместо функции (Y;
x),
соответствующей исходной модели,
используют ее ограниченную модификацию
(Y;
x)
вида
(Y;
x)
=
где a и b – постоянные.
При справедливости заданной модели и не слишком малых значениях a и b, эффективность критерия будет падать на краях области изменения x. Ограниченность (Y; x) влечет за собой нечувствительность (робастность) к влиянию небольшого числа ложных наблюдений. Изменяя длину интервала [– a, b], можно устанавливать компромисс между степенью робастности и снижением эффективности для принятой модели. Робастные методы оказываются полезными в тех случаях, когда при синтезе оптимальных процедур приходится учитывать возможность снижения их эффективности из - за отклонения характеристик входных данных от априорно предполагаемой модели. Робастные алгоритмы обработки сигналов обладают высокой эффективностью при номинальных условиях и приемлемой эффективностью в условиях, когда свойства сигналов и шума могут меняться в пределах заданных классов возможных характеристик.
Робастная процедура обработки сигнала для номинальной модели может уступать в эффективности процедуре, которая для этой модели оптимальна, но зато полная эффективность для данного класса моделей будет высокой или, по крайней мере, достаточной. При этом решение необходимо искать в виде процедуры, которая в наихудшем случае обладает оптимальной эффективностью, что приводит к минимаксным робастным алгоритмам.
Адаптивные алгоритмы сложнее минимаксных процедур фиксированной структуры, поэтому их целесообразно использовать в тех случаях, когда априорная неопределенность слишком велика (что не позволят обеспечить приемлемую эффективность минимаксных методов) или же имеются время или данные, необходимые для адаптации. Минимаксным АОЭИ отдают предпочтение в тех случаях, когда неопределенность классов имеет более ограниченный характер. Здесь их целесообразно использовать как робастные процедуры, предохраняющие от свойственного оптимальным методам чрезмерного падения эффективности при отклонении от исходных допущений. Минимаксные АОЭИ могут применяться и в рамках адаптивного подхода, поскольку механизм обучения в адаптивных процедурах никогда не предполагается настолько совершенным, чтобы за имеющееся ограниченное время произошла полная адаптация.
Для формализации робастного подхода при оценивании экспериментальных данных рассмотрим модель наблюдений
Y(t) = X(t) + N(t), – < t < ,
где X(t), N(t) – сигнал и шум (действительные, центрированные, ортогональные, стационарные в широком смысле случайные процессы). Пусть X(t) и N(t) имеют соответственно спектральные плотности мощности ФX и ФN. Задача оценивания экспериментальных данных (фильтрации, предсказания или сглаживания) сводится к получению оценки X(t) по заданному процессу наблюдений Y(t) в форме
=
,
где h – импульсная характеристика стационарного линейного фильтра. Качество полученной оценки сигнала характеризуется величиной квадратичного риска, который определяется выражением
E{
=
=
=e(ФX,
ФN;
H),
где H – передаточная функция, соответствующая h.
Робастный подход основан на минимизации по H величины падения эффективности фильтра, в наихудшем случае, по критерию:
.
Фильтр, удовлетворяющий данному критерию, является робастным по отношению к классам неопределенности и .
Для решения этой задачи ищется седловая точка минимаксной игры, т.е. спектральная пара (ФX,L, ФN,L) и фильтр HR, удовлетворяющие условию
=
=
.
Согласно правому равенству, HR является оптимальным фильтром для (ФS,L, ФN,L), то есть
HR() = ФX,L()/(ФX,L() + ФN,L()).
Таким образом, определение седловой точки включает в себя отыскание пары (ФX,L, ФN,L), удовлетворяющей данному выражению. Левое равенство показывает, что фильтр HR достигает наихудшей эффективности на спектральной паре (ФX,L, ФN,L), для которой он оптимален. Это наихудшее значение e(ФX,L, ФN,L, HR) определяет гарантированный уровень эффективности фильтра HR на классах и . Так, если оба класса и порождаются всеми спектрами, которые отличаются от исходного не более чем на заданную величину , то робастный фильтр имеет вид
где H0 = ФS/(ФS+ФN) – исходный фильтр, а константы k1 и k2 определяются значением и конкретным типом использованной модели.
Соответствующее минимальное значение квадратичного риска определяется по формуле
e(ФX,
ФN;
HW)
=
= eW(ФX,
ФN).
Из полученного выражения следует, что для определения структуры алгоритма обработки информации необходимо знание спектра оцениваемого сигнала ФX() и шума ФN(), сопровождающего его измерение. Однако при измерении эти составляющие сигнала получают в смеси и отделить их друг от друга затруднено.
2. Выбор алгоритма обработки экспериментальной информации. Одним из наиболее широко используемых в обработке сигналов алгоритмов является дискретный фильтр Калмана [53], основанный на линейной модели наблюдений:
Yn = HnXn + Vn, n = 0, 1, ....,
где Yn – вектор наблюдения размером r для каждого n; Hn – матрица наблюдений; Xn – вектор состояния размером n; Vn – шум наблюдений.
Оптимальные по критерию минимума среднего квадратичного риска оценки Xn и Xn+1 по наблюдениям Y0, ..., Yn определяются рекуррентно, в соответствии с известным алгоритмом калмановской фильтрации:
n+1= Fn n + Wn,
где Fn – одношаговая переходная матрица состояний; Wn – шум состояния.
Из общей теории минимаксного подхода к задачам математической статистики известно, что конечно-выборочная дисперсия минимаксными свойствами не обладает, что определяет использование робастных алгоритмов.
Для завершения формулирования требований к синтезируемому АОЭИ учтем, что процесс вычислений необходимо организовать таким образом, чтобы результат можно было не получать каждый раз заново, а только корректировать с учетом очередного поступившего результата измерения. Поэтому целесообразно использовать алгоритмы, в которых вычислитель, перед проведением любого информационного вычисления, располагает информацией о результатах всех предшествующих вычислений, достаточной памятью для хранения этих результатов и достаточным временем для их обработки, т.е. для проведения соответствующего алгоритмического вычисления необходимо использовать не пассивные (неадаптивные), а последовательные (адаптивные) алгоритмы. В данной ситуации радикальным решением является создание процедур с конечной памятью (на скользящем окне). В связи с этим целесообразно использовать рекуррентный метод вычисления в рамках критерия наименьших квадратов на окне, совпадающем с интервалом усреднения результатов измерения и перемещающемся по оси времени в виде “скользящего окна”. Сформулированным требованиям соответствует алгоритм, определенный по эволюционной теореме [53], для вычисления выборочного среднего с конечной памятью на “скользящем окне”, причем в рассматриваемом алгоритме используются усредненные весовые коэффициенты, с учетом всех значений, входящих в рассматриваемое “окно”.
В данном алгоритме в процессе адаптации определяется “окно” усреднения, а поскольку спектр оцениваемого сигнала неизвестен, то используется рекуррентный метод вычисления. Критерием качества результата решения является минимизация среднего риска вычисления на окне, совпадающем с интервалом усреднения результатов измерения и перемещающемся по оси времени при поступлении очередного значения (“скользящее окно”).
Для уменьшения смещений оценки параметра диагностирования предложено, в рамках метода, использовать усредненные весовые коэффициенты с учетом всех значений, входящих в рассматриваемое окно. Алгоритм, основанный на данном методе, состоит из следующих операций:
1. Начальной оценки по критерию минимума эмпирического риска [1], представляющую собой выборочное среднее:
.
2. Определения весового коэффициента предыдущего измерения:
.
3. Определения весового коэффициента последующего измерения:
.
4. Определения коэффициента передачи для нового измерения, на основании усреднения текущих и предыдущих весовых коэффициентов:
.
5. Нахождения оценки нового измерения, получаемой согласно эволюционной теоремы [74]:
.
6.
Принимая во внимание, что погрешности
представляют собой некоррелированные
случайные величины с нулевым средним
и постоянной дисперсией
,
определяется величина текущего окна,
при которой величина эмпирического
риска [1, 89, 78] равна
.
Разработанный алгоритм реализован в виде программы для обработки осциллограмм, полученных при диагностировании зоны повреждения силовых КЛ (“ORSK”) [110]. Робастный АОЭИ для обработки осциллограммы диагностирования зоны повреждения силовой КЛ, синтезированный в разд. 6.5, реализован в программе "ORSK"V 2.0. Так, для кабеля СБ-325мм2 на напряжение 6кВ, длиной Lкл=300м, зона поражения которого произошла на расстоянии Lп=100м, получены результаты, представленные на рис. 3.4. При этом получены несмещенные оценки измеряемой величины, что достигается за счет введения в АОЭИ [110] усредненных весовых коэффициентов и , которые определяются через выборочное среднее, но на соответствующих интервалах “скользящего окна”.
Другим достоинством модифицированного АОЭИ является отсутствие ухода оценок по интервалам измерения, что приводит к более точному отсчету данных интервалов при ОМП. Необходимо отметить увеличение числа операций в модифицированном АОЭИ на скользящем окне, которое ограничивает его применение при оценке параметров измерения в реальном масштабе времени. Вторым ограничением, накладываемым самим АОЭИ, является необходимость использования нечетного количества замеров в рассматриваемых окнах усреднения, что требуется для исключения ухода получаемых оценок по оси времени. Несмотря на указанные недостатки, для обработки экспериментальных данных, полученных в процессе диагностирования зоны повреждения силовых КЛ, предпочтительно применение модифицированного АОЭИ.
Рис. 3.4. Результаты обработки экспериментальных данных при определении расстояния до места повреждения силового кабеля модифицированным алгоритмом
Данная программа защищена свидетельством об официальной регистрации РОСПАТЕНТом [110].
Задача 3.2. Синтезировать алгоритм обработки экспериментальной информации, получаемой при исследовании акустических характеристик электротехнического оборудования.
Синтезированные в разд. 4.3, 4.4 АОЭИ, а также их робастный тип, полученный в разд. 6.5, можно использовать как при формировании высококачественной идентификационной модели (МИд) исследуемого объекта, так и при высококачественной оценке измеряемого параметра управляемого процесса, а для повышения качества экспериментальных данных на начальных этапах оценки целесообразно использовать алгоритмы, сформулированные выше.
Задача идентификации ЭСВМ включает постановку идентифицирующего эксперимента и решение задачи оценивания экспериментальных данных, полученных в эксперименте. При идентификации ЭСВМ эксперименты должны быть взаимосвязаны общей стратегией проведения. Для измерения параметров в автоматизированной системе научных исследований, контроля и испытаний при разработке и производстве элементов АСУ разработаны и внедрены алгоритмы обработки экспериментальной информации. Программная реализация данного алгоритма использована для автоматизации проектирования и наладки электроприводов корабельных вспомогательных механизмов в автоматизированной системе научных исследований на базе стенда. При этом контролировалась радиальная составляющая вибрационного ускорения ВИД в горизонтальной плоскости ax. Для этого в верхнем углу вертикальной части основания стенда, к которой непосредственно крепится двигатель, был установлен датчик вибрации. Были получены осциллограммы вибрационного ускорения для различных значений тока питания фазы А при прямом и согласованном положении ротора, которые после обработки синтезированным алгоритмом приведены на рис. 3.5. На рис. 3.5 приведены результаты обработки осциллограмм вибраций для различных значений тока питания фазы А до и после обработки предлагаемыми алгоритмами. При этом линия 1 получена для разных значений риска.
Рис. 3.5. Осциллограммы вибраций для различных значений тока питания фазы А до и после обработки предлагаемыми алгоритмами
На участке 1а риск АОЭИ выбран большим, чтобы быстро достигнуть значения оценки обрабатываемого процесса за счет снижения ее качества, а на участке 1б – равным риску кривой 2. В результате получаются оценки обрабатываемого процесса одинакового качества, но за разное время обработки, т.е. АОЭИ на участке 1 обладает большей эффективностью, что согласуется с выводом разд. 4.9. Синтезированные АОЭИ для исследования ЭСВМ обеспечили точность измерения параметра, не превышающую точность датчиков измерения –2–3%, что, следовательно, позволило более точно оценивать параметры МИд ЭСВМ.
Рис.
3.6. Пьезоэлектрический преобразователь
виброизмерительный ДН–3–М1: 1 –
основание; 2 – пьезоэлемент ЭП4К-1-3; 3 –
винт; 4 – крышка; 5 – масса инерционная;
6 – втулка; 7 - контакт; 8 – втулка; 9 –
колпачок; 10 – кольцо.
Задача 3.3. Определить начальную импульсную переходную функцию виброизмерительного пьезоэлектрического преобразователя.
Решение задачи. На рис. 3.6 изображено устройство виброизмерительного пьезоэлектрического преобразователя ДН–3–М1 (паспорт 5Ф2.781.102 ПС), который используется в испытательном стенде для оценивания параметров ВИП. При воздействии механических колебаний на основание вибропреобразователя с ускорением а на инерционную массу m действует сила Р. Пьезоэлементы вибропреобразователя испытывают деформации сжатия – растяжения, в результате чего на грани пьезоэлементов, перпендикулярной вертикальной оси, возникают электрические заряды q, пропорциональные действующей силы P, с коэффициентом пропорциональности, равному индексу пьезомодуля а11. Сигнал с кварцевых пластин снимается экранированным кабелем и поступает на измеритель шума и вибраций ВШВ – 003 – М3 (паспорт ИМ 2.745.000 ПС).
Эквивалентная схема преобразователя, соединенного кабелем с измерительной цепью, представлена на рис. 3.7. Произведем оценку параметров элементов данной схемы замещения. Так, собственное сопротивление пьезоэлемента R0 определяется удельным сопротивлением материала пластин и их поверхностным сопротивлением.
Рис. 3.7. Эквивалентная схема замещения пьезоэлектрического преобразователя: C0 – емкость между гранями пьезоэлектрика (емкость преобразователя); Ск – емкость кабеля между жилой и экраном; Свх – входная емкость измерительной цепи; R0 – сопротивление преобразователя; Rк – сопротивление изоляции кабеля; Rвх – входное сопротивление измерительной цепи
Первая составляющая для кварца (1015 – 1016 Ом) значительно больше второй, поэтому определяющим является поверхностное сопротивление, для повышения которого до R = 109 1010 Ом преобразователь приходится герметизировать, защищая его поверхность от влажности и загрязнения. Электрическое сопротивление изоляции используемого вибропреобразователя ДН–3–М1 и кабеля Rк при нормальных условиях применения не менее 10 ГОм (Под нормальными условиями применения принято: температура окружающего воздуха (205)С, относительная влажность от 30 до 80%, атмосферное давление от 84 до 196 кПа (от 630 до 795 мм рт. ст.)).
Наиболее высокое входное сопротивление Rвх (до 1013–1015 Ом) и входную емкость Свх (до 20 пФ) измерителя шума и вибраций обеспечивают МДП – транзисторы, однако они имеют высокий уровень шумов. Поэтому для этих целей используют полевые транзисторы с p – n – переходом типа КП303Г, входное сопротивление которого составляет не менее 1011Ом.
Используемый вибропреобразователь ДН–3–М1 имеет электрическую емкость при работе с кабелем длиной 2 м (C0 + Ск) 1535325 пФ, а с кабелем 5 м 1800370 пФ. Анализируемую схему замещения можно привести к схеме, в которой включены параллельно сопротивление R и емкость С, где сопротивление R равно сопротивлению параллельного соединения R0, Rк и Rвх, а емкость С = С0 + Ск + Свх. Данная схема замещения позволяет определить передаточную функцию пьезоэлектрического преобразователя с подключенной к нему измерительной цепью. Так, изменение силы P(t) создает переменный ток:
i = а11dP(t)/dt,
где а11 – индекс пьезомодуля для грани, перпендикулярной оси Х (вертикальной). Выходное напряжение преобразователя с подключенной к нему измерительной цепью составляет
вых
=
,
где
= ja11
.
После простых преобразований получаем
вых
= а11R
.
Из данного выражения можно получить передаточную функцию виброизмерительного пьезоэлектрического преобразователя в виде
Wпэп(p)
=
,
которая соответствует передаточной функции неидеального дифференцирующего звена, переходная функция которого имеет вид [111]:
hпэп(t)
=
,
а импульсная переходная функция –
kпэп(t)
= –
.
Амплитудно-частотная характеристика преобразователя, включенного в измерительную цепь, имеет вид
Hпэп()
=
,
а логарифмическая амплитудная характеристика –
20
lg
Hпэп()
= 20 lg
+ 20 lg
RС–
20 lg
.
Данная
логарифмическая амплитудная характеристика
аппроксимируется при <
отрезком низкочастотной асимптоты вида
20 lg + 20 lg RC,
имеющей наклон 20 дб/дек, а при > – отрезком высокочастотной асимптоты – горизонтальной прямой 20 lg , пересекающейся с низкочастотной асимптотой при значении 0 = . При этом фазо – частотную характеристику в области высоких частот можно представить в виде
=
arctg
=
–
RC.
Полученные характеристики позволяют определить выходной сигнал пьезоэлектрического преобразователя в виде
()
= W()
()
= a11Rexp[j(
–
RC)]
(),
то есть выходной сигнал виброизмерительного пьезоэлектрического преобразователя ДН–3–М1 по величине пропорционален входному воздействию и сдвинут на соответствующую задержку.
Анализ полученной амплитудно-частотной характеристики показал, что для расширения частотного диапазона измеряемых величин в сторону низких частот (за счет сдвига точки 0 = ) необходимо увеличивать постоянную времени T = RC. При этом расширение частотного диапазона путем увеличения емкости С осуществляется включением параллельно преобразователю конденсаторов, но из полученного выражения для амплитудно–частотной характеристики следует, что это уменьшает выходное напряжение преобразователя. Увеличение же сопротивления приводит к расширению частотного диапазона без потери чувствительности, однако это требует улучшения качества изоляции и применения усилителей с высокоомным входом.
Подводя итог проведенного исследования виброизмерительного пьезоэлектрического преобразователя ДН–3–М1, можно сделать вывод, что, несмотря на пропорциональность преобразования, данный тип датчика имеет достаточно большую временную задержку. Так, оценки параметров виброизмерительного пьезоэлектрического преобразователя ДН–3–М1 позволили определить постоянную времени RC < 1 с для датчика с усилителем напряжения. Для датчиков с усилителями заряда постоянные времени составляют 10 – 100 с. Однако уже при таких постоянных времени оказывается возможным проводить квазистатическую градуировку пьезоэлектрических датчиков.