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

кандидата технических наук
Пальчик, Олег Викторович
город
Рязань
год
2008
специальность ВАК РФ
05.12.04
цена
450 рублей
Диссертация по радиотехнике и связи на тему «Радиотехнические средства цифровой обработки видеосигналов триангуляционных приборов оперативной дефектоскопии на железнодорожном транспорте»

Автореферат диссертации по теме "Радиотехнические средства цифровой обработки видеосигналов триангуляционных приборов оперативной дефектоскопии на железнодорожном транспорте"

ПАЛЬЧИК Олег Викторович

□03440700

РАДИОТЕХНИЧЕСКИЕ СРЕДСТВА ЦИФРОВОЙ ОБРАБОТКИ ВИДЕОСИГНАЛОВ ТРИАНГУЛЯЦИОННЫХ ПРИБОРОВ ОПЕРАТИВНОЙ ДЕФЕКТОСКОПИИ НА ЖЕЛЕЗНОДОРОЖНОМ ТРАНСПОРТЕ

05 12 04 - Радиотехника, в том числе системы и устройства телевидения

АВТОРЕФЕРАТ

диссертации на соискание ученой степени кандидата технических наук

2 2 СЕН 2003

Владимир 2008

003446700

Работа выполнена на кафедре Радиотехнических систем Рязанского государственного радиотехнического университета

Научный руководитель канд-т технических наук Андреев В Г

Официальные оппоненты

- доктор технических наук, профессор Брюханов Ю А

- доктор технических наук, доцент Полушин П А

Ведущая организация Всероссийский научно-исследовательский конструкторско-технологический институт (ВНИКТИ), г Коломна

Защита состоится « 30 » сентября 2008 года в 14 часов в аудитории 211-1 на заседании диссертационного Совета Д212 025 04 Владимирского государственного университета

Отзывы в двух экземплярах (заверенные печатью) просим направлять по адресу 600000, г Владимир, ул Горького, 87, Совет ВлГУ

С диссертацией можно ознакомиться в библиотеке Владимирского государственного университета

Автореферат разослан « 10 » июля 2008 г

Ученый секретарь

диссертационного совета

доктор технических наук, профессор

А Г Самойлов

ОБЩАЯ ХАРАКТЕРИСТИКА РАБОТЫ

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

Сенсоры, основанные на использовании интерференции (например, коноскопические приборы), обладают высокой точностью Однако для достижения потенциальных возможностей интерферометров требуются точное позиционирование и предсказуемая форма поверхности объекта (кривизна поверхности, ее шероховатость) Кроме того, характерна проблема позиционирования сенсора, когда незначительное его смещение приводит к отражению луча в сторону от приемника Сходная проблема и у теневых сенсоров, использующих эффект дифракции Фраунгофера

Поэтому часто единственно возможным средством является применение лазерной триангуляции (лазерной триангуляционной дально-метрии), что подтверждается многочисленными публикациями в нашей стране (например, фирм «АГРОЭЛ», КТИ СО РАН), и за рубежом (Optical Metrology Center, Великобритания, MEL Mikroelektronik, Германия, LMI Technologies, США, MTI Instruments, США), посвященными промышленному применению триангуляционных систем Однако научная сторона этого вопроса, как правило, остается открытой скромный математический аппарат, недостаточное число обобщающих выводов и рекомендаций не дают возможности решать ряд специфичных задач возникающих, например, на предприятиях железнодорожной отрасли К таким задачам в первую очередь относится оценка геометрических параметров движущихся объектов (как, например, колес вагонов), а также объектов с различной отражающей способностью поверхности (колеса и обрабатываемые механически и наплавкой литьевые элементы вагонной тележки, тормозного и рессорного узлов вагонов)

Сложности связаны с отсутствием анализа возникающих погрешностей на этапе первичной обработки и низкой точностью применяемых вычислительных алгоритмов (под первичной обработкой понимаются задачи оценки положения и ширины лазерного пятна на фотоприемнике -одиночного импульса, и калибровка дальномеров) Следует отметить, что в известных исследованиях их авторы оперируют с недостаточно обоснованными моделями сигналов (как правило, гауссианой, не учитывающей аберрации изображения, формируемого на фотоприемнике), оптимальные и квазиоптимальные методы обработки, рассматриваемые, например, в фундаментальных трудах Тихонова В И , Вайнштейна Л А , Зубакова В Д, Тартаковского Г П , Репина В Г, Прэтта У , Уайлда Д Дж ,

Крамера Г, используются только для получения эталонных оценок и не применяются на практике ввиду значительных вычислительных затрат, а нашедшие применение методы чрезвычайно чувствительны к искажениям формы сигнала, не обладая адаптивными свойствами к его параметрам Отметим также большой вклад в развитие и популяризацию триангуляционного подхода авторов Плотникова С В (КТИ СО РАН, Россия) и Smith К В (The Ohio State University, USA), посвятивших данному вопросу крупные научные исследования.

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

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

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

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

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

точности оценки геометрии сканируемых объектов, упрощение разработки и настройки триангуляционных измерительных систем

Для достижения данной цели необходимо решение следующих задач

- разработка модели сигнала с фотоприемника оптического триангуляционного сенсора,

- сравнительный анализ известных (оптимальных и практических) методов первичной обработки на основе модельных сигналов,

- синтез и анализ новых алгоритмов определения положения одиночного импульса в составе аддитивной смеси,

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

- структурно-параметрическая оптимизация методов первичной обработки,

- анализ погрешностей, возникающих на этапе первичной обработки,

- разработка методики калибровки измерителя, адаптивной к параметрам сигнала,

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

- анализ эффективности внедрения результатов исследования в системы технической диагностики

Методы исследований, использованные в диссертационной работе, основаны на статистической теории радиотехнических систем, параметрическом моделировании случайных процессов, численных алгоритмах поиска экстремума, математическом моделировании Основные числовые результаты получены на основе аналитических и вычислительных математических методов

Основные положения, выносимые на защиту

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

- Косвенный метод нахождения положения одиночного импульса, основанный на процедуре циклической свертки, оптимизированной по вычислительным затратам с помощью рекурсивного преобразования, и дающий возможность в 3 5 раз улучшить оценку положения пятна по критерию среднеквадратичного отклонения (СКО) на фоне белого шума

- Двухпараметрический метод калибровки, основанный на оценке пары

параметров одиночного импульса с использованием косвенного метода, дающий возможность расширить в 1,5 2 раза динамический диапазон измерителя по сравнению с однопараметрической калибровкой путем учета формы и ширины импульса при сохранении вычислительных затрат

Научная новизна диссертации заключается в следующих результатах

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

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

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

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

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

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

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

систем

Внедрение научных результатов диссертационной работы произведено в опытно-конструкторские работы и промышленные разработки ООО «АГРОЭЛ» (Рязань), в разработки вагонной службы «Московской железной дороги» филиала ОАО «Российские железные дороги» и ЗАО «Отраслевой центр внедрения новой техники и технологий»

Апробация работы произведена в форме научных докладов, дискуссий по основным результатам диссертационной работы, которые проходили на следующих научных конференциях 38-й научно-технической конференции РГРТА, II и VI международных научно-технических конференциях «Цифровая обработка сигналов и ее применения»

Публикации. По теме диссертации опубликовано 49 печатных и 5 рукописных научных работ, в том числе 33 публикаций в центральной печати (статьи, патенты, зарегистрированные программы и алгоритмы в Отраслевом фонде алгоритмов и программ Государственного координационного центра информационных технологий Минобразования России), 21 тезисов докладов на конференциях

Структура и объем диссертации

Диссертационная работа состоит из введения, пяти глав, заключения, библиографического списка из 143 наименований и трех приложений Диссертационная работа содержит 160 страниц, в том числе 114 страниц основного текста, 4 таблицы, 70 рисунков

СОДЕРЖАНИЕ ДИССЕРТАЦИИ

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

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

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

Обозначения на рис 1 Р - лазерные пятна, Э - их ширины, Г - фокус объектива, Н - дальность, (3 - угол триангуляции, В - база измерителя, ФП - фотоприемник Точки Т и О, величины 5 и а - вспомогательные параметры, входящие в найденные выражения

Рис 1 - Оптическая схема триангуляционного измерителя

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

Выходной сигнал в = {.$[/]} триангуляционного сенсора представляет собой три основных аддитивных компоненты фоновый шум I = {/[/]} (в общем случае с ненулевой постоянной составляющей), полезная составляющая сигнала в виде одиночного видеоимпульса х = {*[/]}, где у = 0 //-1, N -число отсчетов сигнала, и дополнительный шум в пределах импульса, преимущественно оптического происхождения Отметим, что сигнал с выхода триангуляционного сенсора характеризуется низким отношением сигнал/шум (А/а обычно 5 10, где А - амплитуда импульса, ст - среднеквадра-тическое значение шума) Экспериментальное и теоретическое исследование шумов показало, что наиболее значимая их компонента с высокой степенью точности может быть описана в виде белого шума

Моделью полезного импульса является его аналитическое описание При моделировании импульса предъявляются требования к гибкости изменения параметров модели, определяющих форму модельного импульса Известно, что интенсивность в сечении луча лазера имеет колоколообразную форму Описание такого импульса наиболее целесообразно рядом Эджворта,

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

чл-г. {

ф'

У-М,

_-1±ф(4) 3'

+11ф(5)

У-Ц,

площадь реального

£ ) 3' у I, ) V где Ф("'(») - производная интеграла вероятности, Рх импульса х длительностью Д/

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

£ = У. На рис 2 тонкой линией показан увеличенный фрагмент реального сигнала в, соответствующий импульсу, жирной линией - аппроксимация импульса рядом Эджворта Как видно из рисунка, реальный импульс имеет несимметричный характер, и аппроксимация хе с достаточной степенью передала эту особенность

Далее в первой главе проводится анализ популярного метода центра тяжести с предварительной пороговой обработкой с привлечением разработанной модели Анализ показал наглядность и удобство использования модели сигнала, полученной с использованием предложенной методики Разработанная модель импульса позволила имитировать различные практические ситуации изменение ширины отраженного от объекта лазерного пучка (обусловлено флюктуациями мощности лазера или изменением отражающей способности поверхности объекта), изменение асимметрии полезного импульса, происхождение которой обусловлено особенностями лазера и конструкции измерителя

Получены характеристики оценки положения импульса - смещение и СКО в зависимости от отношения сигнал/шум для различных коэффициентов асимметрии и эксцесса

60 40 20 0 -20

Ас

/ч V /" 'аИ!

X

г

уу

460

Рис 2

] 480 500 520 } - Модель Эджворта импульса

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

лами первой главы для различных мешающих воздействий

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

=£*[*] *[0-*)тосШ] (1)

*=о

Разработанный на основе исследования свойств (1) метод, характеризуется меньшей дисперсией за счет косвенного выделения информации из набора значений циклической свертки Достоинства разработанного метода вытекают из технического приема, заключающегося в совмещении низкочастотной фильтрации и процедуры близкой к дифференцированию

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

На рис 3 приведены зависимости СКО оценки положения импульса от отношения сигнал/шум, полученные различными методами длинный пунктир

- оптимальный (оценка соответствует границе Рао-Крамера), короткий пунктир -косвенный, сплошной линией

- центр тяжести без предварительной пороговой обработки

Как видно из рис 3, косвенный метод незначительно уступает оптимальному методу и значительно выигрывает у метода центра тяжести в отсутствие пороговой обработки Аналогичный анализ, но после пороговой обработки показал стремление качества косвенной оценки к оптимальной при расширении импульса Выигрыш у метода центра тяжести составляет 2 . 4 раза

Также были получены зависимости СКО и смещения от относительного положения р = М-]е импульсной помехи, где - абсолютное положение помехи в массиве отсчетов в, и вероятностные характеристики события сбоя, при котором косвенный метод оказывается практически

Метод ЦТ

косвенный ч \ тод Оптимал ьныи метол

\ \

Рис 3 - СКО различных методов

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

В рамках задачи параметрической оптимизации косвенного метода синтезирована оптимальная дискриминационная характеристика 8 = {&[£]}> проведено исследование влияния различных нелинейных форм нечетной дискриминационной характеристики на СКО оценки, на ее смещение и на полную вероятность сбоя косвенного метода Показано, что при использовании оптимальной (прямоугольной) формы дискриминационной характеристики может быть получен выигрыш в величине СКО до 3 % и уменьшена вероятность сбоя в 40 раз по сравнению с исходной (линейной или же, в силу ограниченной длительности N - треугольной) g

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

На рис 4 показаны обратные однопараметрические калибровочные зависимости М(Н) в пикселях фотоприемника для различных значений Б ширины лазерного пятна на поверхности объекта, а на рис 5 - разница калибровочных зависимостей Дш,2(Н) = М\о{Н) - М2{Щ при минимальном (2 мм) и максимальном (10 мм) значениях Ц а также Д5,г(#) = М5{Н) - Мг(Я) при минимальном и среднем (5 мм) значениях О

ЧН) 750 500 250 0

JJ=2 ч

D--5 м» 1 Ь D -10 мм \

120

140

160

180

Л(Я)

-100

Я

-200

7Г-

\

1С-

Л,. £

120

Н

Рис. 4 - Одномерные калибровочные зависимости

140 160 180 Рис 5 - Разность одномерных калибровочных зависимостей

Рис 4-5 иллюстрируют сокращение диапазона дальностей почти в два раза (границы 140 и 180 мм) Погрешность при этом достигает 6 пикселей (см рис 5, кривая A5 2) при средней ширине пятна, что подтверждает ограниченную применимость традиционной калибровки при возможной вариации ширины L.

Расширение динамического диапазона достигается применением калибровки, адаптивной к ширине импульса Оценка ¿к ширины находится в едином алгоритмическом цикле с косвенной оценкой положения импульса как диапазон индексов между максимумом и минимумом свертки (1), что следует из формы свертки v, близкой к форме производной (рис 6)

В случае оценивания ширины импульса по свертке погрешность а£ оценки существенно меньше и наблюдается устойчивая связь истинной ширины импульса и ее оценки (см рис 7) по сравнению с характеристиками оценки ¿1 по первой производной Характеристики оценки показаны не кривыми, а заштрихованной областью значений т к оценка практически не имеет зависимости от значения £ в силу недопустимо низкого отношения сигнал/шум после дифференцирования и обладает ошибкой превышающей истинное значение

0,05

ОД

0 15 а,/Л

но производной по свертке

Рис 7 - Статистические характеристики оценки ширины импульса

Смещение оценки ¿к (далее - £) и увеличение ее рассеяния при росте уровня шума являются приемлемыми, т к при двумерной калибровке основное значением имеет оценка положения импульса, а оценка ширины импульса играет вспомогательную роль, уточняющую результат оценки дальности

В общем случае адаптивная калибровка, так же как и традиционная, сводится к построению эмпирической табличной зависимости, но, в отличие от известных подходов, для двух параметров - д/ и ¿, реализованной на неравномерной сетке Очевидно, что затруднительно построить калибровочную зависимость //(Л/,£) в виде таблицы для всех возможных дальностей и ширин пятна Поэтому возникает необходимость интерполяции и экстраполяции таблицы В рамках разрешения этой проблемы разработана процедура триангуляции, т е разбиения калибровочной зависимости на треугольники с последующей линейной интерполяцией в пределах полученных примитивов, которая выгодно отличается от известных подходов возможностью обеспечения автоматического редактированием исходных данных, исключающим нереальные результаты измерений

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

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

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

- предложена структурная схема логического блока, входящего в состав устройства, реализующего косвенный алгоритм, которая позволила сократить временные затраты на принятие решения с Л^ до 3 тактов

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

Наибольшее значение имеет процедура рекурсивного преобразования косвенного метода, основанная на представлении циклической свертки V суммой линейных у[/] = упр[/] + Уобр[/], где уо6р[/] = упр[/ + Щ, условно названных прямой и обратной соответственно

Линейная свертка упр[/] может быть выполнена нерекурсивным фильтром с импульсной характеристикой (ИХ) g Для снижения вычислительных затрат произведено рекурсивное преобразование этой свертки на основе условия соответствия ИХ заданной Са = е, где в - верхняя треугольная матрица, составленная из отсчетов g[k], а - искомый вектор коэффициентов цифрового рекурсивного фильтра (ЦРФ),

Модифицированное условие Сма, = е, служащее снижению порядка фильтра {р<Щ, не может быть решено точно, т к. См - прямоугольная (р+1)хА^-мерная матрица, полученная из в путем вычеркивания последних Ы-р столбцов, Я1 - (/7+1)-мерный вектор коэффициентов

На основе формального критерия качества - минимума длины вектора невязок \у = е-Сма1 искомые коэффициенты ЦРФ находятся как решение системы нормальных уравнений

ИХ полученного ЦРФ уже при порядке р- 2 точно совпала с исходной линейной ^[А;] = {0, 1, , ЛЧ} Обратная свертка вычисляется ЦРФ (названным «обратным фильтром» ОФ) с импульсной характеристикой £[к\ = {И-\,Н-2, ,0} Входным сигналом этого фильтра является последовательность = {¿[ЛЧ], ,?[А'-2], , з[0]} Затухающая до 0-го значения ИХ £ не может быть точно описана выходом ЦРФ, поэтому следует ограничиться порядком наиболее соответствующим форме g Так, для линейной ИХ оказывается достаточным р = 1 (ЦРФ обладает при этом экспоненциальной ИХ) Погрешность аппроксимации полностью компенсируется адаптивной калибровкой Более того, наибольшие искажения имеют место при широком импульсе и всегда можно выбрать такой отсчет индексов сигнала в, которому соответствует т е когда вклад g' незначителен

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

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

еТ={1,0, ,0}

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

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

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

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

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

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

- разработана процедура построения аналитической модели с фотоприемника триангуляционного измерителя, учитывающая особенности полезного импульса и статистические свойства типичных шумов, произведено сопоставление параметров модели с возможными условиями эксплуатации измерителя (поляризация излучения, изменения отражающей способности сканируемой поверхности и настроек измерителя), исследованы качественные характеристики популярного метода центра тяжести с использованием построенной модели и показаны оптимальная по критерию минимума СКО величина порога и минимальное отношение сигнал/шум для метода центра тяжести,

- проанализированы известные оптимальные и практические методы оценки положения импульса и выявлены их недостатки путем сравнительного анализа численно определенных характеристик оценок -СКО и смещения,

синтезирован косвенный алгоритм оценки положения одиночного видеоимпульса в составе массива данных, дающий возможность в 2 4 раз уменьшить СКО оценки по сравнению с методом центра тяжести с предварительной пороговой обработкой, исследованы статистические характеристики косвенного алгоритма при наличии белого шума и импульсных помех, получен метод оценивания ширины импульса, выполняющийся в едином алгоритмическом цикле с косвенным алгоритмом, произведена оптимизация логического блока, входящего в состав устройства, реализующего косвенный алгоритм, основанная на параллельном выполнении однотипных операций,

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

разработана процедура двухпараметрической калибровки, адаптивной к ширине пятна, позволяющая в 1,5 2 раза расширить динамический диапазон измерителя, снизив влияние ширины пятна, проанализирован вопрос интерполяции и экстраполяции трехмерной калибровочной зависимости,

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

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

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

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

Все разработанные алгоритмы подходят для реализации на ПЛИС, что при сравнительно невысокой стоимости устройств обработки позволяет достичь максимальной вычислительной эффективности данных алгоритмов

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

В приложениях приведены список аббревиатур и условных обозначений, копии актов о внедрении, фотографии некоторых из устройств первичной обработки

Основные результаты диссертации опубликованы в работах

1 Кошелев В И, Андреев В Г, ПальчикОВ Оптимизация авторегрессионных моделей с квантованными коэффициентами // Цифровая обработка сигналов и ее применения Материалы докладов II международной конференции Москва, 2000 Т III С 217-218

2 Кошелев В И, Андреев В Г, ПальчикОВ Компенсация ошибок определения коэффициентов моделирующих АРСС-фильтров // Известия вузов Радиоэлектроника 2001 Т 44 №7 С 50-55

3 Пальчик О В Эффективная процедура вычисления циклической свертки // Методы и устройства формирования и обработки сигналов в информационных системах Межвузовский сборник научных трудов. Рязань, 2004 С 77-81

4 Пальчик О В, Андреев В Г Обработка данных триангуляционного измерителя // Цифровая обработка сигналов и ее применения Материалы докладов VI международной конференции Москва, 2004 TIIC 197-200

5 Андреев В Г, Пальчик О В Повышение точности измерения геометрических параметров движущихся объектов методом лазерной триангуляции // Датчики и системы 2004 № 9 С 6-9

6 ПальчикОВ Оценка параметров самосветящихся объектов // Проблемы передачи и обработки информации в сетях и системах телекоммуникаций Материалы 13-й международной научно-технической конференции Рязань, 2004 С 78-80

7 Венедиктов А 3, Пальчик О В Проблемы калибровки лазерных триангуляционных измерителей // В мире неразрушающего контроля.

2004. № 4 (26) С 62-63

8 Многокритериальная процедура минимизации собственной засветки триангуляционных измерителей / Венедиктов А 3, Пальчик О В, Ти-решкин ВН, ДоковДС //Приборы 2004 № 10 С 39-41

9 Андреев В Г, Пальчик О В Моделирование сигнала с фотоприемника лазерного триангуляционного измерителя // Приборы и системы Управление, контроль, диагностика 2005 № 1 С 34-37

10 Андреев В Г, Пальчик О В Анализ эффективности математического моделирования лазерных триангуляционных измерителей // Приборы и системы Управление, контроль, диагностика 2005 № 2 С 32-35

11 Пальчик О В, Доков Д С Анализ методов статистической обработки и оценивания параметров одиночного видеоимпульса // Информационные технологии моделирования и управления 2005 № 1 (19) С 58-64

12 Венедиктов А 3, ПальчикОВ Измерение геометрических параметров сложных цилиндрических объектов // Датчики и системы 2005 №1 С 24-28

13 Андреев В Г, Пальчик О В Метод нахождения энергетического центра одиночного импульса // Известия вузов Радиоэлектроника 2005 Т 48 № 7 С. 47-55

14 Пальчик О В, Горкин В Н Обобщение интегральных методов оценки положения импульса методом циклической свертки // Информационные технологии моделирования и управления 2005 №3(21) С. 375-383

15 Венедиктов А 3, Пальчик О В Параметрическая оптимизация алгоритма оценивания положения одиночного видеоимпульса // Цифровая обработка сигналов и ее применения Материалы докладов VII Международной конференции Москва, 2005 Т II С 283-285

16 Венедиктов АЗ, Пальчик О В Метод триангуляции неравномерной гладкой зависимости // Методы и средства измерений физических величин Материалы XI всероссийской научно-технической конференции Н Новгород, 2005 С 20

17 Венедиктов А 3, Пальчик О В Принцип обработки сигналов виброперемещения в оптических измерителях виброперемещения // Приборы и системы Управление, контроль, диагностика 2005 № 6 С 54-56

18 ПальчикОВ Фильтрация одиночного импульса совместно с оценкой положения его экстремума // Вестник Рязанского государственного радиотехнического университета Вып 22 Рязань, 2007 С 99-102

19 Венедиктов А 3, ПальчикОВ, ГорбылевМС Исследование и оптимизация метода обработки контурных изображений с многоуровневым квантованием // Известия вузов России Радиоэлектроника 2007. С 22-32

Пальчик Олег Викторович

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

Автореферат диссертации на соискание ученой степени кандидата технических наук

Подписано в печать 2008 Формат бумаги 60x84 1/16 Бумага офсетная Печать трафаретная Уел печ л 1,0 Уч -изд л 1,0 Тираж 100 экз Заказ

Рязанский государственный радиотехнический университет 390005, г Рязань, ул Гагарина, 59/1 Редакционно-издательский центр РГРТУ

Оглавление автор диссертации — кандидата технических наук Пальчик, Олег Викторович

ВВЕДЕНИЕ.

1. МОДЕЛИРОВАНИЕ СИГНАЛОВ ЛАЗЕРНОЙ ТРИАНГУЛЯЦИОННОЙ СИСТЕМЫ.

1.1. Вводные замечания.

1.2. Построение модели сигнала с фотоприёмника.

1.3. Анализ эффективности моделирования сигнала.

1.4. Выводы.

2. СИНТЕЗ И АНАЛИЗ АЛГОРИТМА ОЦЕНИВАНИЯ ПОЛОЖЕНИЯ ОДИНОЧНОГО ВИДЕОИМПУЛЬСА.

2.1. Вводные замечания.

2.2. Обзор методов получения эталонных оценок.

2.3. Синтез косвенного метода оценки положения импульса.

2.4. Анализ косвенного метода оценки положения импульса.

2.5. Параметрическая оптимизация косвенного алгоритма.

2.6. Выводы.

3. АДАПТИВНАЯ КАЛИБРОВКА ЛАЗЕРНЫХ ТРИАНГУЛЯЦИОННЫХ СИСТЕМ.

3.1. Вводные замечания.

3.2. Двухпараметрическая калибровка триангуляционного сенсора.

3.3. Интерполяция двумерной калибровочной зависимости на неравномерной сетке.

3.4. Экстраполяция калибровочной зависимости.

3.5. Выводы.

4. ВЫЧИСЛИТЕЛЬНАЯ ОПТИМИЗАЦИЯ АЛГОРИТМОВ ПЕРВИЧНОЙ ОБРАБОТКИ В ЛАЗЕРНЫХ ТРИАНГУЛЯЦИОННЫХ СИСТЕМАХ.

4.1. Вводные замечания.

4.2. Комбинированный метод нахождения центра тяжести.

4.3. Рекурсивное преобразование косвенного алгоритма нахождения положения импульса.

4.4. Оптимизация логического блока косвенного метода.

4.5. Параметрическая оптимизация цифровых фильтров с квантованными коэффициентами.

4.6. Выводы.

5. ПРИМЕНЕНИЕ АЛГОРИТМОВ ПЕРВИЧНОЙ ОБРАБОТКИ В ЛАЗЕРНЫХ ТРИАНГУЛЯЦИОННЫХ СИСТЕМАХ ДЛЯ ЗАДАЧ ТЕХНИЧЕСКОЙ ДИАГНОСТИКИ.

5.1. Вводные замечания.

5.2. Повышение точности оценки параметров подвижных объектов.

5.3. Оптический виброметр с высокой разрешающей способностью.

5.4. Визуальный конструктор триангуляционного сенсора.

5.5. Выводы.

Введение 2008 год, диссертация по радиотехнике и связи, Пальчик, Олег Викторович

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

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

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

К таким задачам [1-4] на железнодорожном транспорте в первую очередь относится оценка геометрических параметров движущихся объектов (как, например, колёс вагонов), а также объектов с различной отражающей способностью поверхности (колёса и обрабатываемые механически и наплавкой литьевые элементы вагонной тележки, тормозного и рессорного узлов вагонов, автосцепное устройство и его детали).

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

В России наибольшее внимание триангуляционным системам и триангуляционным сенсорам уделено в работах КТИ НП СО РАН (Новосибирск) и АГРОЭЛ (Рязань). Можно отметить работы, посвящённые промышленному применению триангуляторов [5-7] и их исследованию [8, 9], экспресс диагностике движущихся объектов [10, 11], задачам стационарного контроля [12, 13] и оптимизации сенсоров [14-16] для конкретных задач. Среди публикаций зарубежных авторов можно отметить обзорные работы по триангуляционному принципу [17-19], работы посвящённые промышленному применению [20,21], ориентированные на оптимизацию конструкции измерительных сенсоров [22] или на оптимизацию методов обработки [23, 24], а также общетеоретические [25]. Большой вклад в развитие и внедрение в промышленность триангуляционного принципа внесли работы ряда зарубежных фирм, занимающихся разработкой и производством триангуляционных сенсоров: Optical Metrology Center [26] (Великобритания), MEL Mikroelektronik [27] (Германия), LMI Technologies [28] (США), MTI Instruments [29] (США). Результаты исследований по оптимизации методов цифровой обработки, которые которых могут быть использованы для обработки сигналов триангуляционных сенсоров, изложены, например, в работах [30-38].

Триангуляционный прибор, как основа измерительной системы, может быть представлен как совокупность двух функциональных узлов — датчика (в данном случае - линейного фотоприёмника и приёмной оптики) и преобразователя [39]. Датчик формирует сигнал, имеющий в своём составе искомую информацию, а преобразователь осуществляет её извлечение и приведение к требуемому виду. В преобразователе осуществляется первичный этап обработки сигнала с триангуляционного датчика.

При проведении исследований была сформулирована модифицированная классификация методов обработки, относящихся к первичному и вторичному этапам обработки. Первичная обработка заключается в оценке положения полезного импульса, формируемого как изображение лазерного пятна с поверхности объекта, и может включать в себя фильтрацию сигнала, пороговую обработку и процедуру оценки положения импульса, а также преобразование полученной величины в заданную единицу измерения. Вторичная обработка заключается в построении сечения объекта плоскостями, триангуляции нескольких датчиков и оценке по сечению геометрических параметров объекта, например, как это предлагается в работах [40-43].

Особенность предложенной классификации состоит в том, что на этапе первичной обработки не предусмотрено обнаружения полезной составляющей сигнала, оно производится косвенно - при редактировании данных на этапе вторичной обработки. Методы первичной обработки в отличие от вторичной универсальны и могут быть обобщены для большинства конструкций дальномеров. Отметим, что универсальность этих методов позволяет использовать многие из них не только для оценки геометрических параметров, но и для анализа других типов диагностических сигналов на железнодорожном транспорте [44,45]. Вторичная обработка носит узкоспециализированный характер, зависит от решаемой задачи и, как правило, не имеет жёстких временных и аппаратных ограничений. Поэтому наибольший интерес представляет исследовательская деятельность в направлении оптимизации алгоритмов именно первичной обработки сигналов измерителей, основанных на лазерной триангуляции.

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

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

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

Задачам и методам исследований посвящены многочисленные фундаментальные труды по обработке сигналов отечественных и зарубежных авторов Тихонова В.И. [46], Вайнштейна JI.A., Зубакова В.Д. [47], Тартаковского Г.П., Репина В.Г. [48], ПрэттаУ. [49], УайлдаД.Дж. [50], Крамера Г. [51], Миддлтона [52], Ван Триса [53] в которых приводятся, в том числе, и оптимальные методы.

Однако вопросу моделирования сигнала с учётом его фазовых соотношений (что необходимо для сохранения формы и положения импульса в составе входной аддитивной смеси) практически не уделено внимания. В прикладных исследованиях обычно используется описание полезного импульса в виде гауссианы [54] или гауссианы, искажённой при формировании изображения [55]. Гауссова модель берётся за основу как наиболее типичное распределение интенсивности в сечении луча лазера, однако гауссиана - это идеализация, поскольку различные типы лазеров имеют отличающиеся, характеристики. Например, главным недостатком полупроводниковых лазеров является неодинаковое распределение интенсивности по сечению лазерного луча.

В то же время гауссова модель изначально является некорректной, что обусловлено принципом формирования изображения на фотоприёмнике триангуляционного дальномера (см рис. 1 — оптическая схема дальномера). Для упрощения моделирования оптическая система считается идеальной [56] (на практике такая ситуация также может иметь место при работе в параксиальных лучах (при малом угле триангуляции) и/или при малой диафрагме объектива).

На рис. 1 показаны плоскость линзы, фокусы F и F', границы диапазона по дальности и ширине. На поверхности условного объекта (отмечена штриховкой на границах диапазона по дальности) показаны лазерные пятна Рбл,

Рд и Рср ширина D которых определяется шириной сечения луча лазера. Через центры пятен проведена прямая, совпадающая с направлением лазерного излучения (точка Т - место расположения лазерного излучателя). Угол между этой прямой и нормалью к линзе является углом триангуляции [3. Расстояние от оптического центра О до точки Т называется базой В.

Оптическая схема триангуляционного дальномера

В плоскости сопряжённой с плоскостью луча лазера построены изображения Р'бл и Р'д лазерных пятен согласно правилу построения изображения линзы [56], ширина которых определяется D и индикатрисой рассеяния. В этой же плоскости располагается фотоприёмник ФП, что сделано в соответствие с принципом Шеймпфлуга (Scheimpflug rule) [57]. Угол между плоскостями фотоприёмника и линзы обозначен как а.

Как наглядно показано ниже на рисунке, ширина изображения пятна при изменении дальности может значительно варьироваться. Изображение лазерного пятна на фотоприёмнике представляет собой одиночный импульс с формой определяемой распределением интенсивности в сечении пучка лазера (т.е. априори неизвестной). Отметим, что в отсутствие аберраций характер распределения интенсивности пятна и формы импульса одинаков, но в сопряжённую плоскость попадают только изображения центров пятен, и, таким образом, относительно моды импульса происходит искажения импульса. Как следствие, часть импульса справа от моды может оказаться шире, чем слева или наоборот, причём, чем шире пятно на поверхности объекта и больше угол проецирования изображения на фотоприёмник, тем больше указанные искажения формы импульса. В результате сложнее точно определять положение моды импульса, которая является параметром наиболее близко связанным с дальностью до объекта.

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

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

- н w sin D

Pl(H,D)= ~ —^—.4 ^—(1)

L 2 0(H) ( г г г Л v ^ >

KS(H)±D/2; у sin arctg Н а v KS(H)±D/2j j где знаки «+» и «-» служат для получения ширин левой и правой частей импульса соответственно. Прочие входящие в выражение (1) функционалы и параметры приведены ниже: н . „ . . . (А (1 tg(P)V

Q(H) = — ~ 1, S(H) = H-tg®)-B, а = arctg F

F 5 J '

Для построения идеализированной модели импульса на фотоприёмнике требуется кроме ширин PL частей импульса найти его местоположение. В соответствии с теоремой Пифагора положение моды М импульса на фотоприёмнике определяется как:

М(Н) = J{X{H) + ВУ + Y{H)2 - С,, (2) где Ci - расстояние от точки Т до начала ФП, Х{Н) и 7(77) - координаты точки на ФП относительно оптического центра О:

Y(H) =--—, Х(Н) = .

Q{H) Q{H)

Пересчёт величин (1) и (2) из миллиметров в пиксели фотоприёмника осуществляется в виде х = х ■ N / С2, где N - число элементов ФП, С2 — длина ФП в миллиметрах.

Идеализированная модель, получаемая из импульса любой формы (определяемой, например, характеристиками лазерного излучателя) с помощью преобразования своих параметров по формулам (1)-(2) показывает, что предпочтительная конструкция измерителя, согласованная с принципом Шеймпфлуга, характерна существенным изменением ширины и асимметрии импульса на различных дальностях до объекта (см. рис. 2), которой часто пренебречь нельзя.

Искажения формы импульса за счёт свойств схемы Шеймпфлуга

Рис. 2

Однако на форму полезного импульса также влияют и технические условия его формирования - флюктуации мощности лазера, изменение свойств поверхности объекта и характеристик фотоприёмника и т.п. В частности, существенное влияние оказывают изменение отражающей способности различных участков объекта, которое приводит к изменению амплитуды и эксцесса импульса: боковые стороны покрытые ржавчиной имеют более интенсивное диффузное отражение по сравнению с блестящей поверхностью катания (рис. 3) и особенности настройки коллимирующего устройства: ширина лазерного пятна изменяется в зависимости от априори неизвестной дальности до объекта (рис. 4).

Влияние отражающей способности Влияние настройки

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

Совокупное рассмотрение вопросов первичной обработки и оптимизации оптических параметров конструкции приведено в работе [58]. Однако более поздние исследования аналогичных вопросов, например, приведённые в работе [55], показали неудовлетворительность полученных в [58] результатов и, в первую очередь, неэффективность предложенной модели сигнала с фотоприёмника. Главным недостатком всех известных моделей, применяемым в задачах лазерной триангуляции, является отсутствие или недостаточность связи их параметров с физическими предпосылками поверхности объекта коллимирующего устройства

Рис. 3

Рис. 4 формирования импульса, что не даёт возможности исследовать методы первичной обработки в условиях имитирующих реальную обстановку.

К используемым на практике алгоритмам часто предъявляются требования высокой скорости обработки и вопросу разработки высокоскоростных измерителей посвящено множество работ, особенно при оценке параметров трёхмерной поверхности [59-62]. Поэтому наибольшее распространение получили сравнительно быстрые методы оценки центра тяжести и медианы импульса [37] с предварительной пороговой обработкой, которая служит для исключения фоновых компонент сигнала. Применение пороговой обработки накладывает ограничение на отношение сигнал/шум, что приводит к дополнительным временным затратам на накопление сигнала и обязательную адаптацию порога к сигналу.

К первичной обработке также относится процедура преобразования полученных данных в единицы измерения дальности. Эта процедура обычно основана на использовании градуировочной таблицы триангуляционного измерителя [39] и служит для лианеризации выходного сигнала триангуляционного измерителя и учёта погрешностей конструкции. В силу того, что данная процедура использует результаты уже проведённой ранее градуировки и служит получения действительных характеристик прибора, а не произвольной (хотя и однозначной) градуировочной зависимости выхода прибора от входа, то в работах вышедших в периодической печати она, как правило, имеет название калибровка [39].

Традиционная однопараметрическая калибровка точечных (single point) триангуляционных сенсоров (см. работы фирм Optical Metrology Center [26] и MTI Instruments [29]) обычно оказывается неприменима при вариации параметров сигнала. Однако во многих работах теоретической направленности как, например, работа [55], калибровке вообще не уделено внимания, поскольку результаты моделирования вследствие неудачно выбранных параметров сенсора показали частный результат — линейную связь выходного сигнала триангуляционного датчика и реальной дальности.

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

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

Для достижения данной цели необходимо решение следующих задач:

- разработка модели сигнала с фотоприёмника оптического триангуляционного сенсора;

- сравнительный анализ известных (оптимальных и практических) методов первичной обработки на основе модельных сигналов;

- синтез и- анализ новых алгоритмов определения положения одиночного импульса в составе аддитивной смеси;

- анализ вычислительной эффективности различных методов первичной обработки и вычислительная оптимизация методов первичной обработки;

- структурно-параметрическая оптимизация методов первичной обработки;

- анализ погрешностей, возникающих на этапе первичной обработки;

- разработка методики калибровки лазерного триангуляционного измерителя, адаптивной к параметрам сигнала;

- разработка программного обеспечения, служащего для моделирования и проектирования триангуляционных сенсоров;

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

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

Основные положения, выносимые на защиту

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

- Косвенный метод нахождения положения одиночного импульса, -основанный на процедуре циклической свёртки, оптимизированной по вычислительным затратам с помощью рекурсивного преобразования, и дающий возможность в 3.5 раза улучшить оценку положения пятна по критерию СКО на фоне белого шума.

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

Научная новизна диссертации заключается в следующих результатах:

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

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

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

- разработана процедура двухпараметрической калибровки триангуляционных измерителей, адаптивная к ширине лазерного пятна с поверхности объекта;

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

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

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

Внедрение научных результатов диссертационной работы произведено в опытно-конструкторские работы и промышленные разработки ООО «АГРОЭЛ» (Рязань), в разработки вагонной службы «Московской железной дороги» филиала ОАО «Российские железные дороги» и ЗАО «Отраслевой центр внедрения новой техники и технологий».

Апробация работы произведена в форме научных докладов, дискуссий по основным результатам диссертационной работы, которые проходили на следующих научных конференциях: 38-й научно-технической конференции РГРТА, II и VI международных научно-технических конференциях «Цифровая обработка сигналов и её применения».

Публикации. По теме диссертации опубликовано 49 печатных и 5 рукописных научных работ, в том числе 33 публикаций в центральной печати (статьи, патенты, зарегистрированные программы и алгоритмы в Отраслевом фонде алгоритмов и программ Государственного координационного центра информационных технологий Минобразования России), 21 тезис докладов на конференциях.

Структура и объём диссертации

Диссертационная работа состоит из введения, пяти разделов, заключения, библиографического списка из 143 наименований и 3-х приложений. Диссертационная работа содержит 160 страниц, в том числе 114 страниц основного текста, 4 таблицы, 70 рисунков.

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

5.5. Выводы

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

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

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

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

ЗАКЛЮЧЕНИЕ

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

- разработана процедура построения аналитической модели с фотоприёмника триангуляционного измерителя, учитывающая особенности полезного импульса и статистические свойства типичных шумов; произведено сопоставление параметров модели с возможными условиями эксплуатации измерителя (поляризация излучения, изменения отражающей способности сканируемой поверхности и настроек измерителя); исследованы качественные характеристики метода центра тяжести с использованием построенной модели и показана оптимальная по критерию минимума СКО величина порога и минимальное отношение сигнал/шум для метода центра тяжести;

- проанализированы известные оптимальные и практические методы оценки положения импульса и выявлены их недостатки путём сравнительного анализа численно определённых характеристик оценок — СКО и смещения;

- получен косвенный алгоритм оценки положения одиночного видеоимпульса в составе массива данных, дающий возможность в 3.5 раз уменьшить СКО оценки по сравнению с методом центра тяжести с предварительной пороговой обработкой;

- получен метод оценивания ширины импульса, выполняющийся в едином алгоритмическом цикле с косвенным методом оценки положения импульса;

- исследованы статистические характеристики косвенного алгоритма при наличии шума и импульсных помех;

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

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

- разработана процедура двухпараметрической калибровки адаптивной к ширине пятна, позволяющая в 1,5.2 раза расширить динамический диапазон измерителя, исключив влияние ширины пятна; также проанализирован вопрос интерполяции и экстраполяции калибровочной зависимости;

- разработана процедура триангуляции поверхности калибровки с последующей линейной интерполяцией калибровочной зависимости;

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

- разработана рекурсивная процедура вычисления циклической свёртки;

- разработана процедура параметрической оптимизации цифровых

• рекурсивных линейных фильтров произвольной структуры с квантованными коэффициентами, основанная на использовании двух шкал квантования и пересчёте ошибок грубой шкалы в многоразрядные (т.н. «точные») коэффициенты;

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

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

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

Все разработанные алгоритмы подходят для реализации на ПЛИС, что при сравнительно невысокой стоимости устройств обработки позволяет достичь максимальной вычислительной эффективности данных алгоритмов.

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

Библиография Пальчик, Олег Викторович, диссертация по теме Радиотехника, в том числе системы и устройства телевидения

1. Венедиктов А. 3., Дёмкин В. И., Доков Д. С. Лазерные методы и средства контроля геометрии деталей // В мире неразрушающего контроля — 2004.-№1 (23).-С. 67-68.

2. Венедиктов А.З., Дёмкин В.Н., Доков B.C. Измерение параметров колесных пар подвижного состава во время движения // Железные дороги мира.- 2003.- № 9.- С. 33-36.

3. Венедиктов А.З., Дёмкин В.Н., Доков Д.С., Комаров А.В. Применение лазерных методов для контроля параметров автосцепки и пружин // Новые технологии железнодорожному транспорту: Сборник научных статей с международным участием - Омск — 2000.- С. 232-234.

4. Латышев Ю.В., Кудоба Т.С., Плотников С.В., Подчернин В.М. Оценка конкурентоспособности лазерных триангуляционных измерителей расстояний // Датчики и системы 2001 - № 6 - С. 46-49.

5. Плотников С.В., Подчернин В.М., Быковская И.В. Триангуляционные измерители и их промышленное применение // Техника машиностроения.- 2003,- № 4. С. 107-108.

6. Венедиктов А.З., Пальчик О.В., Власов Д.А. Опыт внедрения оптикоэлектронных измерительных средств на ремонтных предприятиях // Современные промышленные технологии: Материалы II Всероссийской научно-технической конференции Н. Новгород - 2005 - С. 21-22.

7. Вертопрахов В.В. Влияние формы объекта и ориентации его поверхности на точность лазерных триангуляционных измерений // Автометрия.- 1995-№ в.- С. 64-68.

8. Плотников С.В., Подчернин В.М., Быковская И.В. Исследование и разработка триангуляционных измерителей и их промышленное применение // Наука производству.- 2003- № 2 (58).- С. 43-44.

9. W.N. Demkin, D.S. Dokov, V.N. Tereshkin, A.Z. Venediktov, Laser control device of spring parameter // Nondestructive Testing and Computer Simulations in Materials Science and Engineering, Alexander I. Melker, Editor, Proceedings of SPAS VoL3/A15.

10. V.N. Demkin, D.S. Dokov, A.Z. Venediktov, Measurement of wheel parameters of a rolling stock during movement // Proceedings of SPIE.- Vol. 5066.

11. Венедиктов A.3., Дёмкин B.H., Доков Д.С. Особенности трехмерного измерения геометрических размеров деталей лазерным триангуляционным способом // Лазеры. Измерения. Информация. Тезисы докладов конференции Санкт-Петербург, 2004- С. 84.

12. A.Z. Venediktov, V.N. Demkin, D.S. Dokov, Choice of optimum modes of laser triangulation meter at control of surface form // Proceedings of SPIE-Vol. 5381.-PP. 103-109.

13. Многокритериальная процедура минимизации собственной засветки триангуляционных измерителей / Венедиктов А.З., Пальчик О.В., Тирёшкин В.Н., Доков Д.С. // Приборы.- 2004. № 10,- С. 39-41.

14. Dorsch G. Hausler, Herrman J.M. Laser triangulation: fundamental uncertainty in distance measurement// Appl. Opt. 33 — 1994-PP. 1306-1314.

15. F. Blais, M. Lecavalier, J. Bisson, "Real-time Processing and Validation of Optical Ranging in a Cluttered Environment", ICSPAT 1996 — PP.1066-1070.

16. T.A. Clarke, The development of an optical triangulation pipe profiling instrument // Optical 3-D Measurement Techniques, 1995 Vol. Ill - PP. 331-340.

17. T.A. Clarke, K.T.V. Grattan, N.E. Lindsey, Laser-based triangulation techniques in optical inspection of industrial structure // Proc. SPIE, 1990 — Vol. 1332.-PP. 474-486.

18. M.C. Leu, Z. Ji, Non-linear displacement sensor based on optical triangulation. U.S. Patent 5113080, 1992.

19. H. Wang, D. Malacara, Optical triangulation: A dual channel configuration//Rev. Sci. Instrum., 1996-Vol. 67-PP. 2606-2611.

20. H. Rothe, M. Tuershmann, P.P. Mager, R. Endter Improve accuracy in laser triangulation by variance-stabilizing transformation // Opt. Eng., 1992.-Vol. 31.-PP. 1538-1545.

21. G.A.W. West, T.A. Clark, A survey and examination of subpixel measurement techniques//Proc. SPIE, 1990.-Vol. 1395.-PP. 456-463.

22. Куликов Е.И. Методы измерения случайных процессов.-М.: Радио и связь, 1986 — 272 с.

23. Русинов Л.А. Оценки положения аналитического пика при различении помех с использованием обобщенного преобразования Фурье // Известия вузов. Приборостроение.— 1978 —Т. 21 — № 10 —С. 15-20.

24. G. Jacoviti, G. Scarano, Discrete time techniques for time delay estimation // IEEE Trans, on Signal Processing, 1993- Vol. 41- PP. 525-533.

25. Королев Н.И., Кренев А.Н., Анджан С.Э., Лашков Н.И.

26. Алгоритм оценивания частоты и фазы по двум главным компонентам спектрального разложения // Радиоэлектроника — 1995.- № 3.- С. 31-38.

27. Кошелев В.И., Горкин В.Н. Алгоритмические методы повышения точности оценки доплеровской фазы сигнала в процессоре БПФ // Тезисы докладов 37-й научно-технической конференции.-Рязань: РГРТА, 2001 —С. 17.

28. Зиатдинов С.И., Аграновский А.В., Осипов Л.А. Оценка параметров импульсного сигнала в дискретных системах // Известия вузов. Приборостроение 2004.- Т. 47.- № 5 - С. 10-16.

29. Плотников С.В. Сравнение методов обработки сигналов в триангуляционных измерительных системах // Автометрия— 1995—№ 6 — С. 58-63.

30. Двинских В.А. Оценка спектра составляющей цифрового сигнала с интенсивной гармонической помехой // Журнал технической физики —2004 — Т. 74.-Вып. 3.- С. 85-86.

31. МИ 2247-93. «Рекомендация. ГСИ. Метрология. Термины и определения».-СПб.: ВНИИМ, 1993.

32. Венедиктов А.З., Пальчик О.В. Измерение геометрических параметров сложных цилиндрических объектов // Датчики и системы —2005.— №1- С. 24-28.

33. Венедикитов А.З., Пальчик О.В., Горбылев М.С. Библиотека математических функций Cyclometry. -М.: ВНТИЦ, 2005.-№ ГР50200500156.

34. Венедиктов А.З., Пальчик О.В., Горбылёв М.С. Анализ физических характеристик тел вращения // Мир измерений —2005 .— № 7 — С. 15-18.

35. Пальчик О.В., Власов Д.А. Спектральный анализ как инструмент прогресса в современном приборостроении // Составляющие научно-технического прогресса: Материалы Международной научно-практической конференции Тамбов - 2005 - С. 137-138.

36. Пальчик О.В., Андреев В.Г. Статистический анализ диагностических сигналов при контроле тормозных систем на железнодорожном транспорте // Вестник Рязанской государственной радиотехнической академии Вып. 13 — Рязань, 2003- С. 114-117.

37. Автоматизированная система опробования тормозов грузовых составов: Отчёт о НИР (закл.) / «МЖД» филиал ОАО «РЖД»; Научн. рук. Венедиктов А.З. -№ 08.00.12/03.03.03 ДУ-17. -Рязань, 2003. -53 с. -Соисполн.: Пальчик О.В., Фурцев А.И., Козлов А.И. и др.

38. Тихонов В.И. Статистическая радиотехника М.: Советское радио, 1966.- 680 с.

39. Вайнштейн JI.A., Зубаков В.Д. Выделение сигналов на фоне случайных помех М.: Советское радио, I960 — 448 с.

40. Репин В. Г., Тартаковский Г. П. Статистический синтез при априорной неопределенности и адаптация информационных систем Ы'.: Советское радио, 1977 - 432 с.

41. Прэтт У. Цифровая обработка изображений, в 2-х томах — М.: Мир,1982.

42. УайлдД.Дж. Методы поиска экстремума / Пер. с англ. А.Н. Кабалевского, Е.П. Маслова, В.Д. Спиридонова; Под ред. А.А. Фельдбаума М.: Наука, 1967.-268 с.

43. Крамер Г. Математические методы статистики: Пер. с англ. / Под ред. А.Н. Колмогорова-М.: Мир, 1975 648 с.

44. Миддлтон Д. Введение в статистическую теорию связи: Пер. с англ. в 2-х т. / Под ред. Б.Р. Левина М.: Советское Радио, 1961, 1963.

45. Ван Трис Г. Теория обнаружения, оценок и модуляции. Обработка сигналов в радио- и гидролокации и прием случайных гауссовых сигналов нафоне помех Т. 3: Пер. с англ. Под ред. В.Т. Горяинова.- М.: Советское Радио, 1977.- 664 с.

46. Ахманов С.А., Никитин С.Ю. Физическая оптика: Учеб. для вузов М.: Изд-во Моск. ун-та, 1998 — 655 с.

47. Kyung-Chan Kim, Jong-Ahn Kim, SeBaek Oh, Soo Hyun Kim, Yoon Keun Kwak, Accuracy enhancement of point triangulation probes for linear displacement measurement // Proc. SPIE, 2000.- Vol. 3945.- PP. 88-95.

48. Матвеев A.H. Оптика — M.: Высшая школа, 1985 — 351 с.

49. Donald В. Kilgus, Donald J. Svetkoff Imaging geometry and error sensitivity in triangulation-based optical receivers // Proc. SPIE Vol.2599.- 1996-PP. 106-119.

50. K.B. Smith, Modeling, performance evaluation, calibration, and path planning of point laser triangulation probes in coordinate metrology // Ph.D. Dissertation, The Ohio State Univ., Columbus, 1996.

51. J.P. Lavelle, S.R. Schuet, D.J. Schuet, High-speed 3D scanner with real-time 3D processing // Proc. SPIE, 2004.- Vol. 5393.- PP. 19-28.

52. D. Svetjkoff, Towards a high-resolution, video rate, 3d sensor for machine vision // Proc. SPIE, 1986.-Vol. 728.

53. G. Hauster, J. Herrman, Physical limits of 3-D sensing // Proc. SPIE, 1992.-Vol. 1822.-P. 150.

54. Пальчик O.B. Особенности контурной обработки изображений в лазерных триангуляционных системах // Информационные технологии в науке, проектировании и производстве: Материалы XIV Всероссийской научно-технической конференции Н. Новгород — 2005.- С. 19-20.

55. Заварыкин В.М., Житомирский В.Г., Лаичик М.П. Численные методы М., Просвещение, 1990 - 176 с.

56. D.B.T. Kilgus, D.J. Svetkoff, Distortion characteristics and mapping in triangulation imaging systems//Proc. SPIE, 1994.-Vol. 2348.-PP. 106-129.

57. Секен К., ТомисетМ. Приборы с переносом заряда / Пер. с англ. Под ред. В.В. Поспелова, Р.А. Суриса М.: Мир, 1978 - 327 с.

58. Андреев В.Г., Пальчик О.В. Повышение точности измерения геометрических параметров движущихся объектов методом лазерной триангуляции // Датчики и системы — 2004 № 9- С. 6-9.

59. Плескунин В.И., Воронина Е.Д. Теоретические основы организации и анализа выборочных данных в эксперименте. Под ред. А.В. Башарина- JL: Изд-во Ленингр. ун-та, 1979.-232 с.

60. Левин Б.Р. Теория случайных процессов и её применение в радиотехнике. -М.: Советское радио, I960 664 с.

61. R. Baribeau, М. Rioux, Influence of speckle on laser range finders // App. Opt., 1991.-Vol. 30 (20).

62. R. Baribeau, M. Rioux, Centroid fluctuations of speckled targets // App. Opt., 1991.-Vol. 30 (26).

63. БакутП.А., Мандросов В.И., Матвеев И.Н., Устинов Н.Д. Теория когерентных изображений М.: Радио и связь, 1987 - 264 с.

64. Андреев В.Г., Пальчик О.В. Моделирование сигнала с фотоприёмника лазерного триангуляционного измерителя // Приборы и системы. Управление, контроль, диагностика 2005-№ 1- С. 34-37.

65. Андреев В.Г., Пальчик О.В. Анализ эффективности математического моделирования лазерных триангуляционных измерителей // Приборы и системы. Управление, контроль, диагностика —2005—№ 2 — С. 32-35.

66. Дёмкин В.Н., Доков Д.С., Привалов В.Е. Особенности применения лазерных диодов в линейных измерениях // Письма в Журнал технической физики.- 2004,- Т.30 № 13 - С. 40-44.

67. Пальчик О.В., Горкин В.Н. Обобщение интегральных методов оценки положения импульса методом циклической свёртки // Информационные технологии моделирования и управления. 2005 - № 3(21).— С. 375-383.

68. Краснов В.Н., Сахно С.П., Тымчик Г.С. Алгоритмы поиска экстремальных значений видеосигнала ПЗС-приемников // Приборостроение—1986.-№4.-С. 77-81.

69. Чураков Е.П. Оптимальные и адаптивные системы: Учеб. пособие для вузов М.: Энергоатомиздат, 1987 - 256 с.

70. Блейхут Р. Быстрые алгоритмы цифровой обработки сигналов: Пер. с англ.- М.: Радио и связь, 1989 448 с.

71. Пальчик О.В., Горкин В.Н. Обобщение интегральных методов оценки положения импульса методом циклической свёртки // Информационные технологии моделирования и управления. — 2005 — № 3(21).- С. 375-383.

72. Андреев В.Г., Пальчик О.В. Метод нахождения энергетического центра одиночного импульса // Известия вузов. Радиоэлектроника — 2005 — Т. 48.-№7.- С. 47-55.

73. Пальчик О.В., Андреев В.Г. Обработка данных триангуляционного измерителя // Цифровая обработка сигналов и её применения: Материалы докладов VI Международной конференции Москва, 2004.-T.IL-С. 197-200.

74. Гольденберг JI.M., Матюшкин Б.Д., Поляк М.Н. Цифровая обработка сигналов: учеб. пособие для вузов — М.: Радио и связь, 1990 — 256 с.

75. Венедиктов А.З., Пальчик О.В. Параметрическая оптимизация алгоритма оценивания положения одиночного видеоимпульса // Цифровая обработка сигналов и её применения: Материалы докладов VII Международной конференции — Москва, 2005 Т. II — С. 283-285.

76. Пальчик О.В. Оценка параметров самосветящихся объектов // Проблемы передачи и обработки информации в сетях и системах телекоммуникаций: Материалы 13-й Международной научно-технической конференции — Рязань 2004 - С. 78-80.

77. Венедиктов А.З., Пальчик О.В. Проблемы калибровки лазерных триангуляционных измерителей // В мире неразрушающего контроля — 2004 — №4 (26).-С. 62-63.

78. Венедиктов А.З., Андреев В.Г., Пальчик О.В. Адаптивная калибровка систем промышленного зрения // Материалы докладов III междисциплинарной конференции с международным участием НБИТТ-21 — Петрозаводск, 2004 С. 59.

79. Шенягин В.П., Битюков В.К. Измерение длительности колоколообразного импульса // Цифровая обработка сигналов и её применения: Материалы докладов VI Международной конференции Москва, 2004 — Т. II — С. 231-233.

80. Андреев В.Г., Венедиктов А.З., Пальчик О.В., Тирёшкин В.Н.

81. Двухпараметрическая калибровка лазерных триангуляционных измерителей (депонированная рукопись) / Рязан. гос. радиотехн. академия —Рязань, 2004 — 14 е.- 8 ил.- Библиогр.: 8 назв.- Рус.- Деп. в ВИМИ, 12.10.2004, №Д08975.

82. Корнейчук Н.П. Сплайны в теории приближения —М.: Наука, 1984.-352 с.

83. Скворцов А.В. Обзор алгоритмов построения триангуляции Делоне // Вычислительные методы и программирование —2002 —№ 1— С. 14-39.

84. Майкл JI. Вычислительная геометрия и компьютерная графика на С++-М.: БИНОМ, 1997.- 304 с.

85. Jonathan R. Delaunay Refinement Algorithms for Triangular Mesh Generation, Computational Geometry: Theory and Applications 22(l-3):21-74, May 2002.

86. Шишкин A.B., Боресков А.В. Компьютерная графика. Полигональные модели.- М.: Диалог-МИФИ, 2000 — 464 с.

87. Венедиктов А.З., Пальчик О.В. Метод триангуляции неравномерной гладкой зависимости // Методы и средства измеренийфизических величин: Материалы XI Всероссийской научно-технической конференции — Н. Новгород — 2005 — С. 20.

88. Отнес Р., Эноксон Л. Прикладной анализ временных рядов / Пер. с англ. А.И. Хохлова; под ред. И.Г. Журбенко- М.: Мир, 1982 428 с.

89. Венедиктов А.З., Пальчик О.В. Способ диагностирования тормозной магистрали железнодорожного состава: Патент Российской Федерации, МКИ В60Т 17/22.-Заявл. 18.08.2003, №2250167 // Опубл. 20.04.2005 в Бюл. № 11.

90. Сверхбольшие интегральные схемы и современная обработка сигналов / Под ред. Гунна С., УайтхаусаХ., КайлайтаТ. Пер. с англ.:—М.: Радио и связь, 1989 472 с.

91. Венедикитов А.З., Тирёшкин В.Н., Пальчик О.В. Алгоритм вычисления центра тяжести оптимизированный под архитектуру микроконтроллеров IP2022.-M.: ВНТИЦ, 2004.-№ ГР50200401286.

92. Андреев В.Г., Пальчик О.В. Комбинированный метод поиска центра тяжести в лазерных триангуляционных системах // Физика и технические приложения волновых процессов: Труды III Международной научно-технической конференции —Волгоград, 2004 — С. 186.

93. Нуссбаумер Г. Быстрое преобразование Фурье и алгоритмы вычисления сверток —М: Радио и связь, 1985 248 с.

94. Брейсуэл Р. Преобразование Хартли: Пер. с англ.-М.: Мир, 1990—175 с.

95. Пальчик О.В. Эффективная процедура вычисления циклической свёртки // Методы и устройства формирования и обработки сигналов винформационных системах: Межвузовский сборник научных трудов —Рязань, 2004.-С. 77-81.

96. РабинерА., Гоулд Б. Теория и применение цифровой обработки сигналов-М.: Мир, 1978 848 с.

97. РайсД.Р. Матричные вычисления и математическое обеспечение: Пер. с англ.- М.: Мир, 1984 264 с.

98. Кошелев В.И., Андреев В.Г. Оптимизация АР—моделей процессов с полимодальным спектром // Изв. вузов. Радиоэлектроника— 1996 — Т. 39 — № 5.- С. 43-48.

99. Бакулев П.А., Кошелев В.И., Андреев В.Г. Оптимизация АРСС-моделирования эхо-сигналов // Известия вузов. Радиоэлектроника —1994-Т. 37 —№ 9-С. 3-8.

100. Введение в цифровую фильтрацию / Под ред. Р. Богнера, А. Константинидиса; Пер с англ.- М: Мир, 1976.-218 с.

101. Еремеев В.П., Мизиненко О.В. Оптимизация амплитудно частотных характеристик цифровых фильтров // Сборник научных трудов факультета радиоэлектроники и вычислительных систем —Рига: РАУ, 1998.

102. Кошелев В.И., Андреев В.Г., Пальчик О.В. Оптимизация авторегрессионных моделей с квантованными коэффициентами // Цифровая обработка сигналов и её применения: Материалы докладов II Международной конференции.- Москва, 2000.- Т. III С. 217-218.

103. Кошелев В.И., Андреев В.Г., Пальчик О.В. Компенсация ошибок определения коэффициентов моделирующих АРСС-фильтров // Известия вузов. Радиоэлектроника.- 2001.- Т. 44.- № 7.- С. 50-55.

104. Марпл-мл. C.JI. Цифровой спектральный анализ и его приложения: Пер. с англ.- М.: Мир, 1990 584 с.

105. Кошелев В.И., Андреев В.Г. Модифицированный алгоритм АР— моделирования узкополосных процессов // Цифровая обработка сигналов и её применения: Материалы докладов II Международной конференции.-М., 1999Т. III. С. 703-705.

106. Котов О.И., Лиокумович Л.Б., Марков С.И., Медведев А.В., Николаев В.М. Белосветный волоконно-оптический межмодовый интерферометр // Письма в ЖТФ.- 1999.- Т. 25.- Вып. 12.- С. 44-50.119. http://www.altera.com

107. Лазер И.М., Шубарев В.А. Устойчивость цифровых микроэлектронных устройств —М.: Радио и связь, 1988.-215 с.

108. Телевизионная астрономия / Под ред. Никонова В.Б.— М.: Наука, 2-е изд., 1984.-420 с.

109. Зайдель А.П. Ошибки измерений физических величин —М.: Наука, 1974.- 108 с.

110. Рабинович С.Г. Погрешности измерений-Л.: Энергия, 1978 —262 с.

111. Трифонов А.П., Захаров А.В., Проняев Е.В. Обработка импульсов с гауссовской случайной субструктурой при наличии шума // Цифроваяобработка сигналов и её применения: Материалы докладов III Международной конференции Москва, 2000 - Т. I - С. 80-84.

112. Пальчик О.В., Доков Д.С. Анализ методов статистической обработки и оценивания параметров одиночного видеоимпульса // Информационные технологии моделирования и управления —2005 —№ 1 (19).— С. 58-64.

113. Болотин В.В. Прогнозирование ресурсов машин и конструкций — М.: Машиностроение, 1984 312 с.

114. ДатнерБ. Анализ вибраций роликовых и шариковых подшипников: Пер. с англ.— Конструирование и технология машиностроения.— М.: Мир, 1979.- Т. 101.- № 1.- С. 65-82.

115. Виброакустическая диагностика зарождающихся дефектов / Ф.Я. Балийкий, М.А. Иванова, А.Г. Соколова, Е.И. Хомяков-М.: Наука, 1984120 с.

116. Метод неразрушающего контроля состояния подшипников / Венедиктов А.З., Дёмкин В.Н., Пальчик О.В., Белокрылов А.Г. //

117. Прикладная радиоэлектроника. Состояние и перспективы развития: Сборник научных трудов по материалам Международного Радиоэлектронного Форума — Харьков, 2002,- Часть 2,- С. 475-476.

118. ГОСТ 24347-80. «Вибрация. Обозначения и единицы».

119. Вильнер Л.Д. Виброскорость как критерий вибрационной напряженности упругих систем // Проблемы прочности.— 1970.—№ 9 — С. 42-45.

120. Руссов В.А. Спектральная вибродиагностика.- Пермь, 1996.

121. Н. Kakishima et al. Quarterly Report of RTRI, 2000,-№3.-PP. 127-130.

122. Венедиктов А.З., Пальчик О.В. Принцип обработки сигналов виброперемещения в оптических измерителях виброперемещения // Приборы и системы. Управление, контроль, диагностика 2005.- № 6 — С. 54-56.

123. Пальчик О.В., Тирёшкин В.Н. Программа визуального конструирования лазерных триангуляционных измерителей v. 1.0 — М.: ВНТИЦ, 2004.- № ГР50200400896.

124. Пальчик О.В. Программное средство проектирования лазерных триангуляционных измерителей // Современные проблемы информатизации в технике и технологиях: Сборник трудов X Международной научной конференции-Воронеж, 2005-Вып. 10 — С. 195-196.

125. Архангельский А.Я. Программирование в C++Builder-М.: ЗАО «Издательство БИНОМ», 2000 1152 с.