|
|
Если шум, то его нужно фильтровать. Можно так.
Сигнал по прежнему умножаете на окно. Берете пр. Фурье от произведения. Находите квадраты амплитуд всех гармоник 50 Гц. Сколько их убралось, наиболее заметных, в Вашей полосе,от 50 Гц до Fdiskr/2.
Находите удвоенную сумму квадратов амплитуд гармоник.
Извлекаете корень из суммы. Все.
Не удовлетворит - остается по сути только увеличение размера выборки,точнее числа периодов в ней.
Что-то можно отыграть выбросив окно вообще, но немного видимо, но будет сложнее, и вряд ли будет стоит овчинка выделки.
Дальше для случая с окном.
По нахождению квадратов амплитуд.
Предварительно нужно взять пр. Фурье от окна. Просуммировать квадраты модулей отсчетов с номерами N-2,N-1,0,1,2. 5 штук. Центр, он же максимум и по паре справа и слева от него. Полученную сумму запомнить. Назову ее оконная нормировка. Вообще число членов в сумме зависит от окна, от ширины его главного лепестка в частотной области. Тем шире, тем больше членов в сумме. Например 7 - центр и по 3 с боков. Или 3 - центр и по одному с боков. 5 - некая усредненная цифра, подходящая в большинстве случаев.
Теперь с квадратами амплитуд гармоник. Находите максимум в районе очередной гармоники. Суммируете квадраты модулей 2 отсчетов справа от него, 2-х слева и квадрат модуля ссамого максимума. Тоже 5 членов в сумме. Делите на оконную нормировку. Все. Делаете то же самон для следующей гармоники. И т.д. Складываете результаты для всех гармоник. Удваиваете. Извлекаете корень.
Если бы оконная нормировка вычислялась по 7, значит и кв. амплитуды тоже по 7.
Вообще более правильно делать аппроксимацию в районе максимума. Аппроксимировать нужно спектром окна. Но это сложнее.
Описанная процедура проста, поэтому и приведена.
E-mail: 
info@telesys.ru
