Алгоритм поиска пиков

[Yaroslav11]

Добрый день!!! Помогите отсортировать данные. Пишу программу на Borland c++ 3.1 поиск пиков в аппаратурном спектре гамма излучения с плоховатым разрешением -15% пики широкие. Выбрал метод поиска пиков по второй производной. d2y=y[-1]-2y[0]+y[1]; Нахожу минимальные значения. Проблема в том что в области пиков имеется несколько минимумов, которые прошли по критерию 3 стандартных отклонения (из-за широких пиков). Далее нахожу вес каждого пика по формуле квадрат y[xmax] /дисперсию всего спектра. Необходимо отсортировать (ложные пики выкинуть) полученные данные centr[i], Ves[i] по следующим критериям:
1)Если centr[i] < 0.6*centr[i+1] то пики из одной области, сравниваем их вес
если Ves[i]>Ves[i+1] то Ves[i+1]=0 и centr[i+1]=0 выкидываем пик, если нет то отбрасываем centr[i].
2) Если centr[i] > а*centr[+1] то это пик.а=0,6 или около того
И так нужно сравнить соседние найденные пики.
Пример данных первый столбец канал второй вес:
55 20,1
81 1,97
189 2,61
220 5,19
247 2,23
Должны остаться 55 20,1;220 5,19
Помогите. Подскажите если кто знает более эффективный метод поиска пиков. Буду благодарен. Спасибо!