Добрый день!

Вопрос следующий:
есть смесь гармонического сигнала с шумом, параметры сигнала не известны, как найти ОСШ?

Сначало думал через спектр взять мощность гармоники и поделить на мощность шума в спектре, но если у нас возникает растекание спектра то этот метод не подходит

У кого какие предложения?

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

в том то и дело что частота семплирования не кратна и ее менять я не могу

а если к основной гармонике запетлеваться по фазе и вычесть в противофазе ? Естественно с учетом амплитудного выравнивания.

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

никаких фурье и спектров не надо = все во временной области делается.
ПС - если частота сигнала не известна, то можно ее поискать, если не спешите тем же фильтром с перестройкой полосы пропускания

Если параметры сигнала не известны, то тогда нужно сказать диапазоны их изменений. И на каком временном интервале их можно считать постоянными:
Амплитуда
Частота
Отношение сигнал-шум
Также, насколько сильны предположения об аддитивности, стационарности и белости шума.

А почему ТС боится растекания спектра?
Используйте оконную функцию и усреднение спектра, в чем проблемы?
При с/ш = 12 дБ сигнальная палка превышает средний уровень белого шума примерно на 20 дБ, так что идентифицировать сигнал особой сложности не представляет.

И таки да, современные измерители с/ш делают на основе спектрального анализа.

об этом много написано, много сложных и не очень способов. Многое зависит от желаемого результата оценки и отношения С/Ш.
я предложил максимально простое решение.

Есть возможность поиспользовать авторегрессионные модели, как это описано в книжке Марпла "Цифровой спектральный анализ" гл 9.

Зная параметры модели, можно получить оценку мощности синусоиды.

Неоднократно делал. Оценивал, правда, фактически С/(С + Ш), но при хорошем точности вычислений и C/Ш до 30-40 дБ точность получаемого С/Ш вполне нормальная. Сигнал находим как мощность на выходе узкополосного фильтра (удобно реализовать в виде: перенос сигнала на 0 частоту IN*e^(-j*Wсигн*t) ->комплексный ФНЧ или просто интегратор -> OUT*OUT"), сигнал+шум - как просто мощность (IN^2 при действительном входе или IN*IN" при комплексном). Обе мощности интегрируем по достаточно длинному окну, после чего делаем вычисляем оценку. Sorry, если какой-нить масштабный множитель забыл.

гармоника + шум (белый) оценивается просто:
- мощность гармоники - квадрат средней амплитуды
- мощность шума - дисперсия
чем больше усредняем, тем точнее оценка

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

Уже неоднократно обсуждали тему как найти период и фазу синусоиды: "Измерение частоты основной гармоники.."

Предлагаете начать по новой?

Вероятно, имелось в виду усреднение спектров