реклама на сайте
подробности

 
 
 
Reply to this topicStart new topic
> Измерение точки Ликвидуса.
VadimPlotnikov
сообщение Apr 29 2008, 04:19
Сообщение #1


Участник
*

Группа: Участник
Сообщений: 34
Регистрация: 20-03-07
Из: Красноярск
Пользователь №: 26 321



Стоит задача измерения точки Ликвидуса (переход из жидкое в твёрдое состояние), характеризуется замедлением остывания вещества . Измерение происходит термопарой с занесением данных в массив. На графике ищется точка перегиба по второй производной. В теории всё просто, а на практике проблемы. Сигнал зашумлён, и в результате вторая производная тоже зашумлена что резко снижает достоверность результата или невозможность его получения, хотя визуально на графике исходной функции видно где происходит перелом. Применение фильтров (медиана – убивает данные полностью, усреднение раз на раз) на входные данные не даёт нужного результата. Реально помогает только увеличение шага производной (т.е берём не подряд точки 1-2-3 а 10-20-30), но ухудшается результат. Читал, что хорошие результаты даёт вейвлет-преобразования но пока не разобрался как использовать.
У кого-какие мысли.
Go to the top of the page
 
+Quote Post
Stanislav
сообщение Apr 29 2008, 04:54
Сообщение #2


Гуру
******

Группа: Свой
Сообщений: 4 363
Регистрация: 13-05-05
Из: Москва
Пользователь №: 4 987



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


--------------------
Самонадеянность слепа. Сомнения - спутник разума. (с)
Go to the top of the page
 
+Quote Post
stoker
сообщение Apr 29 2008, 08:37
Сообщение #3


Местный
***

Группа: Свой
Сообщений: 340
Регистрация: 28-11-05
Из: Москва
Пользователь №: 11 469



Цитата(VadimPlotnikov @ Apr 29 2008, 08:19) *
Стоит задача измерения точки Ликвидуса (переход из жидкое в твёрдое состояние), характеризуется замедлением остывания вещества . Измерение происходит термопарой с занесением данных в массив. На графике ищется точка перегиба по второй производной. В теории всё просто, а на практике проблемы. Сигнал зашумлён, и в результате вторая производная тоже зашумлена что резко снижает достоверность результата или невозможность его получения, хотя визуально на графике исходной функции видно где происходит перелом. Применение фильтров (медиана – убивает данные полностью, усреднение раз на раз) на входные данные не даёт нужного результата. Реально помогает только увеличение шага производной (т.е берём не подряд точки 1-2-3 а 10-20-30), но ухудшается результат. Читал, что хорошие результаты даёт вейвлет-преобразования но пока не разобрался как использовать.
У кого-какие мысли.

Полагаю в данном случае нужно проводить измерения с помощью синхронного детектора. Можете поискать в инете "Lock-in", они позволяют мерить полезный сигнал в присутсвие больших шумов. Например в даташите на микруху LTC1043 есть апликейшн где применяется метод синхронного детектирования, думаю изменить под вашу задачу не составит особого труда. Если можно сменить способ измерения - ставьте платиновые датчики, благо их много всяких.
Go to the top of the page
 
+Quote Post
Stanislav
сообщение Apr 29 2008, 08:47
Сообщение #4


Гуру
******

Группа: Свой
Сообщений: 4 363
Регистрация: 13-05-05
Из: Москва
Пользователь №: 4 987



Цитата(stoker @ Apr 29 2008, 12:37) *
Полагаю в данном случае нужно проводить измерения с помощью синхронного детектора. Можете поискать в инете "Lock-in", они позволяют мерить полезный сигнал в присутсвие больших шумов. Например в даташите на микруху LTC1043 есть апликейшн где применяется метод синхронного детектирования, думаю изменить под вашу задачу не составит особого труда.
Простите, а каким боком здесь вообще синхронное детектирование? Что Вы предлагаете детектировать, и синхронно с чем? 07.gif

Цитата(stoker @ Apr 29 2008, 12:37) *
...Если можно сменить способ измерения - ставьте платиновые датчики, благо их много всяких.
Что есть платиновые датчики, и почему с их применением изменится способ измерения? И вообще, что Вы о нём знаете, и зачем его нужно непременно менять? 07.gif


--------------------
Самонадеянность слепа. Сомнения - спутник разума. (с)
Go to the top of the page
 
+Quote Post
NickNich
сообщение Apr 29 2008, 09:02
Сообщение #5


Местный
***

Группа: Свой
Сообщений: 375
Регистрация: 8-11-05
Пользователь №: 10 593



Цитата(VadimPlotnikov @ Apr 29 2008, 08:19) *
У кого-какие мысли.


Медианный фильтр изначально предназначен для удаления импульсных помех, т.е. редких по сравнению с длиной фильтра, выбросов значительной амплитуды, а не для фильтрации сигнала на фоне шумов. Если у Вас есть точка перегиба и нужно найти ее координаты, попробуйте приближение ваших экспериментальных данных полиномом третьего порядка. Коэффициенты полинома рассчитываются методом наименьших квадратов. Третий порядок полинома выбран потому, что это полином наименьшего порядка, имеющий точку перегиба.
Go to the top of the page
 
+Quote Post
stoker
сообщение Apr 29 2008, 09:57
Сообщение #6


Местный
***

Группа: Свой
Сообщений: 340
Регистрация: 28-11-05
Из: Москва
Пользователь №: 11 469



Цитата(Stanislav @ Apr 29 2008, 12:47) *
Простите, а каким боком здесь вообще синхронное детектирование? Что Вы предлагаете детектировать, и синхронно с чем? 07.gif

Измерять сопротивление термочувствительного элемента, на заданной частоте или в узкой полосе, таким образом получая большее отношение сигнал - шум. Вообще техника довольно известная, сам использовал.
Вот пример: http://www.cappels.org/dproj/dlmom/dlmom.html думаю тут все понятно.
Цитата
Что есть платиновые датчики, и почему с их применением изменится способ измерения? И вообще, что Вы о нём знаете, и зачем его нужно непременно менять? 07.gif

Платиновые датчики - достаточно линейны, большой выбор сопротивлений, прощще проводить измерение сопротивления с помощью синхронного детектирования чем мерять напряжение на термопаре, которое может быть сравнимо с шумами наведёнными на провода, хотя и здесь можно подмешать опорную частоту.
Go to the top of the page
 
+Quote Post
Stanislav
сообщение Apr 29 2008, 10:30
Сообщение #7


Гуру
******

Группа: Свой
Сообщений: 4 363
Регистрация: 13-05-05
Из: Москва
Пользователь №: 4 987



Цитата(stoker @ Apr 29 2008, 13:57) *
Измерять сопротивление термочувствительного элемента, на заданной частоте или в узкой полосе, таким образом получая большее отношение сигнал - шум. Вообще техника довольно известная, сам использовал.
Вот пример: http://www.cappels.org/dproj/dlmom/dlmom.html думаю тут все понятно.
Понятно.
Ваше предложение касается методики измерений. Вопрос темы не совсем об этом: есть уже данные измерений, нужно получить конкретный результат.

Цитата(stoker @ Apr 29 2008, 13:57) *
...Платиновые датчики - достаточно линейны, большой выбор сопротивлений, прощще проводить измерение сопротивления с помощью синхронного детектирования чем мерять напряжение на термопаре, которое может быть сравнимо с шумами наведёнными на провода.
Питание терморезистора переменным током, с последующим СД, конечно, может быть полезно для уменьшения фликкер-шума, и других помех.
Однако, мы не знаем всех условий задачи: например, платиновое термосопротивление может быть по каким-то причинам неприменимо, нормирующий усилитель не может быть установлен в непосредственной близости от датчика, и т.д.
Вопрос выбора и правильного включения датчика заслуживает отдельной темы; здесь же нужно обработать уже имеющиеся данные.


--------------------
Самонадеянность слепа. Сомнения - спутник разума. (с)
Go to the top of the page
 
+Quote Post
stoker
сообщение Apr 29 2008, 10:44
Сообщение #8


Местный
***

Группа: Свой
Сообщений: 340
Регистрация: 28-11-05
Из: Москва
Пользователь №: 11 469



Ну это на тот случай если все таки изменения можно вносить в железку. А уж если нет и данные получены, можно попробовать их загнать в матлаб, да и поиграться в нём с различными способами фильтраций, дабы выбрать оптимальный вариант обработки, в частности там есть вейвлет-анализ, правда я не совсем понимаю как его можно здесь использовать.
Go to the top of the page
 
+Quote Post
Stanislav
сообщение Apr 29 2008, 10:48
Сообщение #9


Гуру
******

Группа: Свой
Сообщений: 4 363
Регистрация: 13-05-05
Из: Москва
Пользователь №: 4 987



Цитата(stoker @ Apr 29 2008, 14:44) *
А уж если нет и данные получены, можно попробовать их загнать в матлаб, да и поиграться в нём с различными способами фильтраций, дабы выбрать оптимальный вариант обработки, в частности там есть вейвлет-анализ, правда я не совсем понимаю как его можно здесь использовать.
Вот именно. smile.gif

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


--------------------
Самонадеянность слепа. Сомнения - спутник разума. (с)
Go to the top of the page
 
+Quote Post
vvs157
сообщение Apr 29 2008, 11:12
Сообщение #10


Профессионал
*****

Группа: Свой
Сообщений: 1 526
Регистрация: 8-04-05
Пользователь №: 3 960



Цитата(VadimPlotnikov @ Apr 29 2008, 08:19) *
Стоит задача измерения точки Ликвидуса (переход из жидкое в твёрдое состояние), характеризуется замедлением остывания вещества . Измерение происходит термопарой с занесением данных в массив. На графике ищется точка перегиба по второй производной. В теории всё просто, а на практике проблемы. Сигнал зашумлён, и в результате вторая производная тоже зашумлена что резко снижает достоверность результата или невозможность его получения, хотя визуально на графике исходной функции видно где происходит перелом. Применение фильтров (медиана – убивает данные полностью, усреднение раз на раз) на входные данные не даёт нужного результата. Реально помогает только увеличение шага производной (т.е берём не подряд точки 1-2-3 а 10-20-30), но ухудшается результат. Читал, что хорошие результаты даёт вейвлет-преобразования но пока не разобрался как использовать.
У кого-какие мысли.
Попробуйте дифференцирование со сглаживанием, где k - глубина сглаживания, dx - шаг оцифровки. Обычно удается вытянуть то, что видно на глаз
Прикрепленные изображения
Прикрепленное изображение
 
Go to the top of the page
 
+Quote Post
VadimPlotnikov
сообщение May 8 2008, 05:18
Сообщение #11


Участник
*

Группа: Участник
Сообщений: 34
Регистрация: 20-03-07
Из: Красноярск
Пользователь №: 26 321



Решили следующим образом. Сначала по всем отсчета используем кривую Безье. Далее меняя шаг дефиренцирования по трём точкам набираем мат. статистику. Потом по максимуму мат. статистики выбираем перелом. Получилось довольно неплохо.
Go to the top of the page
 
+Quote Post

Reply to this topicStart new topic
1 чел. читают эту тему (гостей: 1, скрытых пользователей: 0)
Пользователей: 0

 


RSS Текстовая версия Сейчас: 6th July 2025 - 11:02
Рейтинг@Mail.ru


Страница сгенерированна за 0.01481 секунд с 7
ELECTRONIX ©2004-2016