Аппроксимация спектра как суммы гауссовых пиков
От: dims12 http://www.relativity.ru
Дата: 09.12.03 21:15
Оценка:
Привет!

У меня есть спектр в виде таблицы значений. Если его отобразить на графике, то там будут несколько пиков разного местоположения и вышины. Я бы хотел вычислить эти пики, то есть, предположить, что каждый пик — это "гауссоида" и рассчитать параметры все гауссоид для спектра.

Есть ли какой-нибудь общепринятый алгоритм на эту тему или придется самому итерациями подбирать?

Димс.
Re: Аппроксимация спектра как суммы гауссовых пиков
От: vasiap  
Дата: 10.12.03 09:09
Оценка:
Здравствуйте, dims12, Вы писали:

D>Привет!


D>У меня есть спектр в виде таблицы значений. Если его отобразить на графике, то там будут несколько пиков разного местоположения и вышины. Я бы хотел вычислить эти пики, то есть, предположить, что каждый пик — это "гауссоида" и рассчитать параметры все гауссоид для спектра.


D>Есть ли какой-нибудь общепринятый алгоритм на эту тему или придется самому итерациями подбирать?


D>Димс.


Рекомендую посмотреть нелинейный метод наименьших квадратов.

Василий.
Re[2]: Аппроксимация спектра как суммы гауссовых пиков
От: zabivan  
Дата: 11.12.03 10:55
Оценка:
Здравствуйте, vasiap, Вы писали:

V>Рекомендую посмотреть нелинейный метод наименьших квадратов.


А ссылку не подбросите?
Re: Аппроксимация спектра как суммы гауссовых пиков
От: Ignoramus  
Дата: 11.12.03 21:27
Оценка:
Здравствуйте, dims12, Вы писали:

D>У меня есть спектр в виде таблицы значений. Если его отобразить на графике, то там будут несколько пиков разного местоположения и вышины. Я бы хотел вычислить эти пики, то есть, предположить, что каждый пик — это "гауссоида" и рассчитать параметры все гауссоид для спектра.


D>Есть ли какой-нибудь общепринятый алгоритм на эту тему или придется самому итерациями подбирать?


Я бы сначала обработал НЧ-фильтром чтобы сгладить спектр и выделить настоящие пики на фоне шумовых выбросов, затем бы нашел все локальные максимумы (мат ожидания), и двигаясь от них в стороны, вычислил бы дисперсию.

Кстати, почему именно "гауссоиды"?
Re[3]: Аппроксимация спектра как суммы гауссовых пиков
От: vasiap  
Дата: 12.12.03 14:21
Оценка:
Здравствуйте, zabivan, Вы писали:

Z>Здравствуйте, vasiap, Вы писали:


V>>Рекомендую посмотреть нелинейный метод наименьших квадратов.


Z>А ссылку не подбросите?


Увы, навскидку не скажу. Я в молодости пользовался CERNLIB (www.cern.ch-> далее поиск по сайту). Только с ее загрузкой,
насколько я помню могут быть проблемы- она свободна для научных
организаций- в смысле IP адрес научной организации должен быть зарегистрирован в ЦЕРНе. Если есть знакомые из Дубны или подобных
организаций- попросите их помочь. Эта библиотека на фортране и великолепно работает. Так же в ЦЕРНе есть т.н. ROOT на С++.(распространяется на тех же условиях). В нем нелинейный МНК тоже должен быть реализован. Но как он работает- не имею представления.

а так, простите за банальность, только гугл.

Василий.
 
Подождите ...
Wait...
Пока на собственное сообщение не было ответов, его можно удалить.