|
|
Оцените:
Суть метода в кратце такая (в предположении что шума нет, только синус неизвестной частоты, хотя можно и амплитуду добавить и оценить ее):
F=sin(wt) удовлетворяет дифф. уравнению:
d^2(F)/dt^2 + k*F = 0;
где k - это квадрат от w=2*pi*f
Далее не буду излагать основу метода, но конечный результат такой:
Минимизируем SUM[ (x'' + k*x)^2 ], по эйлеру требуем нуля производной по k.
и получаем оценку: k = - Sum(x'' * x)/Sum(x*x)
тоесть точность оценки определяется только точностью аппроксимации цисленной второй производной, что описано во всех вычметодах.
Если сигнал с шумом, то при моделировании оказалось что начиная с определенной длинны прямоугольного фильтра дисперсия оценки стабилизируется и меняется слабо, оценка проводилась на периоде частоты и синус был представлен 250 точками.
Если интересно, то могу посмотреть как называется статья на эту тему и сообщить.
E-mail: 
info@telesys.ru
