|
|
  |
Измерение точки Ликвидуса. |
|
|
|
Apr 29 2008, 04:19
|
Участник

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

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

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

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

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

|
Цитата(stoker @ Apr 29 2008, 12:37)  Полагаю в данном случае нужно проводить измерения с помощью синхронного детектора. Можете поискать в инете "Lock-in", они позволяют мерить полезный сигнал в присутсвие больших шумов. Например в даташите на микруху LTC1043 есть апликейшн где применяется метод синхронного детектирования, думаю изменить под вашу задачу не составит особого труда. Простите, а каким боком здесь вообще синхронное детектирование? Что Вы предлагаете детектировать, и синхронно с чем? Цитата(stoker @ Apr 29 2008, 12:37)  ...Если можно сменить способ измерения - ставьте платиновые датчики, благо их много всяких. Что есть платиновые датчики, и почему с их применением изменится способ измерения? И вообще, что Вы о нём знаете, и зачем его нужно непременно менять?
--------------------
Самонадеянность слепа. Сомнения - спутник разума. (с)
|
|
|
|
|
Apr 29 2008, 09:02
|
Местный
  
Группа: Свой
Сообщений: 375
Регистрация: 8-11-05
Пользователь №: 10 593

|
Цитата(VadimPlotnikov @ Apr 29 2008, 08:19)  У кого-какие мысли. Медианный фильтр изначально предназначен для удаления импульсных помех, т.е. редких по сравнению с длиной фильтра, выбросов значительной амплитуды, а не для фильтрации сигнала на фоне шумов. Если у Вас есть точка перегиба и нужно найти ее координаты, попробуйте приближение ваших экспериментальных данных полиномом третьего порядка. Коэффициенты полинома рассчитываются методом наименьших квадратов. Третий порядок полинома выбран потому, что это полином наименьшего порядка, имеющий точку перегиба.
|
|
|
|
|
Apr 29 2008, 09:57
|

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

|
Цитата(Stanislav @ Apr 29 2008, 12:47)  Простите, а каким боком здесь вообще синхронное детектирование? Что Вы предлагаете детектировать, и синхронно с чем?  Измерять сопротивление термочувствительного элемента, на заданной частоте или в узкой полосе, таким образом получая большее отношение сигнал - шум. Вообще техника довольно известная, сам использовал. Вот пример: http://www.cappels.org/dproj/dlmom/dlmom.html думаю тут все понятно. Цитата Что есть платиновые датчики, и почему с их применением изменится способ измерения? И вообще, что Вы о нём знаете, и зачем его нужно непременно менять?  Платиновые датчики - достаточно линейны, большой выбор сопротивлений, прощще проводить измерение сопротивления с помощью синхронного детектирования чем мерять напряжение на термопаре, которое может быть сравнимо с шумами наведёнными на провода, хотя и здесь можно подмешать опорную частоту.
|
|
|
|
|
Apr 29 2008, 10:30
|

Гуру
     
Группа: Свой
Сообщений: 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)  ...Платиновые датчики - достаточно линейны, большой выбор сопротивлений, прощще проводить измерение сопротивления с помощью синхронного детектирования чем мерять напряжение на термопаре, которое может быть сравнимо с шумами наведёнными на провода. Питание терморезистора переменным током, с последующим СД, конечно, может быть полезно для уменьшения фликкер-шума, и других помех. Однако, мы не знаем всех условий задачи: например, платиновое термосопротивление может быть по каким-то причинам неприменимо, нормирующий усилитель не может быть установлен в непосредственной близости от датчика, и т.д. Вопрос выбора и правильного включения датчика заслуживает отдельной темы; здесь же нужно обработать уже имеющиеся данные.
--------------------
Самонадеянность слепа. Сомнения - спутник разума. (с)
|
|
|
|
|
Apr 29 2008, 11:12
|
Профессионал
    
Группа: Свой
Сообщений: 1 526
Регистрация: 8-04-05
Пользователь №: 3 960

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

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

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