Отправлено misyachniy 16 марта 2005 г. 18:35 В ответ на: Можно ли средствами матлаб определить частоту дискретизации по следующим данным. отправлено misyachniy 16 марта 2005 г. 11:43

Полезный сигнал в полосе от 0.05 до 200 Герц
Снято 2048 отсчетов все числа положительные.
Максиальное значение 12439.
Минимальное 50
Если посмотреть на график сигнала

vovka=fopen('c:\close.bin','r')
m=fread(vovka,inf,'*uint16');
j=1:2048;
plot(j, m(j))

То визуально видно что период сигнала 20 отсчетов.
20 отсчетов помноженное на 50 дает предположение, что частота
дискретизации 1000 отсчетов/секунду.

Разложил в ряд Фурье

f=fft(m);
plot(j, f(j))

На сколько я понимаю совет SM - значения выборок(амплитуд) сигналов от 1025 до 2048 не учитываются.
График получился симметричным относительно средины отсчетов 1024/1025.

При рассмотрении графика выборок отмечаю следующие особенности:
В точке 1 значение 3E6 - 3 миллиона может и многовато - не знаю, это по моему мнению постоянная составляющая.

Затем самый первый видимый минимум (интересно почему не максимум?)
-4.8274e+005 выборка 105
затем
-2.9687e+005 выборка 211
-2.7608e+005 выборка 316

Вторая и третья гармоники?

По рассчету предложенному SM получаю 50 * 1024/105 = 487 выборок в секунду.
Что не сходится с видимым.

Если 50 * 2048/105 = 975 выборок в секунду, что похоже на
результат видимый на первом графике.

Прошу считать мой опус - написанный инженером-электронщиком, а не математиком DSP ;-)
И не бить больно, а дать толковый совет.

• Поработал над сигналом еще, похоже на правду(+) — misyachniy   (17.03.2005 10:55, 849 байт, картинка)
• Частотный интервал между соседними спектр. отсчетами=Fdiskr/N. N - длина временнОй реализации, над которой берется Фурье. N=2048. — -=ВН=-   (16.03.2005 20:08, пустое)
• ну я имел в виду модуль фурья - то есть abs(fft(m)) и постоянка откуда? — SM   (16.03.2005 19:49, пустое)
  • Вот - сигнала у меня нет, сам сгенерил нечто (+) — SM   (16.03.2005 19:57, 281 байт)
    • И в догонку - я имел в виду (+) — SM   (16.03.2005 20:24, 471 байт)