поиск пиков на спектре гамма изл. (вейвлет преобразование)

[Volosat1y]

наткнулся на интересный файлик в интернете:
A package for gamma-ray spectrum analysis and routine neutron activation analysis pdf 327 KB

хочу воспроизвести (в смысле закодить) алгоритм поиска пиков (GWTM), помогите разобраться )

[/math]

[незваный гость]

Так в чем вопрос? Какого рода помощь Вы хотите получить?

[Volosat1y]

если я правильно понял, то коэффициенты $a_n$. эт коэффициенты зависомости $Y = a_1x^1+a_2x^2+...+a_Nx^N (полиномиальная регрессия), где $x_i , y_i$ - эт значения энергии и центройда/ширины соответственно и N - степень полинома.
если я ошибаюсь, пожалуйста поправьте =)

следующий пункт 2.2 поиск пиков.

в формуле (3)
поскольку A(i)- число каунтов в i-ом канале. (или, выражаясь по-русски, значение i-ой составляющей вектора спектра), то N в данном случае это количество каналов в спектре.

тогда что такое a, b и С?

дальше сказанно что а зависит от FWHM, и еще чуть-ниже что $ a = \frac {FWHM} {2.353} а это вообще есть формула $ \sigma =\frac {FWHM} {2\sqrt{ln(2)}} $, что вводит меня в больший ступор поскольку сигма зависит от A(i), и эту зависимость можно получить пересчитав "Energy and width calibration" в терминах каналов (а не энергии).

вопрос: как считается этот вейлет - формула (3) или (5).

Добавлено спустя 56 минут 52 секунды:

что-то мне подсказывает что а берется с определеным шагом из какого-то интервала. и вейвлет считается для всех этих значений.
C задаваемая константа.
а вот что делать с b так и не могу понять =)

[Volosat1y]

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