задачи 348 + 337

[STyulb]

На самом деле - это одна и та же задача, но разные подходы. 337 будет работать быстро, так как связан с фурье, но нужны, скорее всего, большие изменения в коде программы. Такой способ применяют по всему миру (это стандарт) и поэтому результат гарантирован.
348 будет работать медленно и результат не гарантирован (по логике должно быть все правильно). Способ 348 распишу, так как в самой задаче пустота, а не постановка задачи. В задачах 348 можно будет закрыть.

1. Записываем точный период, модуль, луч, какие дисперсии перебирать (перебор дисперсий пункт 5).
2. Перебираем все дни и выкидываем плохие (это у нас реализовано)
3. Для "хороших" дней делаем нормирование по ступеньке, а потом усреднение с заданным периодом (один период, а не двойной) и перебираем меры дисперсии на данном периоде. Запоминаем зависимость SNR/DM за данный день. Берем следующий день и проделываем те же процедуры. В конце суммируем (или не в конце, а по ходу дела) все SNR/DM за "хорошие" дни. Ожидается, что появится пик на правильном DM.
4. Определение SNR сходное с тем, что в задаче по поиску транзиентов с большими DM. После усреднения по заданному периоду на данном DM перебираем все точки среднего профиля и ищем минимальное значение, медианное значение и максимальное значение среди точек среднего профиля данного дня. Называем SNR_max=2*(max-mediana)/(median-min). Правильный коэффициент между 2 и 3. Его нужно подбирать отдельно. но для данной задачи это не имеет никакого значения. В каталог для каждого "хорошего" дня записывается значение SNR_max. Минимальные SNR_max будут равны двойке.
5. У нас в усредненных спектрах мощности очень много одиночных пиков. Если для каждого пика сложение периода делать в пределах DM=1-1000, то время счета будет бесконечным. Из уже найденных пульсаров ясно, что большая часть найденных пульсаров имеет меры дисперсии меньше, чем 50. Поэтому при поиске пульсаров по одиночным пикам, найденным в спектрах мощности, лучше всего перебирать меры дисперсии сотнями (задается в пункте 1). Например на мерах дисперсии 1-100 перебрать все одиночные пики. Затем выкинуть из списка перебираемых пиков все найденные новые пульсары. Оставшиеся пики пересчитать на DM=101-200, затем вновь выкинуть из списка все найденные пульсары и т.д. Такой способ резко сократит счет.