Отправлено SM 09 июля 2004 г. 19:12 В ответ на: И вот:-) отправлено -=ВН=- 09 июля 2004 г. 18:09

Про сигнал и СПМ. Я думал на эту двойку сначала, еще вчера, но... Честная сумма квадратов сигнала, деленная на кол-во выборок сошлась с интегрированием его СПМа, отсюда я заключил, что матлаб эту двойку учел за меня, что собственно подтвердилось графиком, в котором на оси ясно написано, что там db/герц, а не на пару герц а потом и чтением документации (пришлось заставиться :)).

Метод Уэлча - это усреднение модифицированных периодограм, в общем куча ПФов выполненных с наложением окон во времени (я их накладывал через каждую выборку) и с усреднением результатов. ф-ция pwelch матлаба, наверняка и в маткаде нечто есть. В общем сигнал обрабатывал вот так:

s=[1 1 1 zeros(1,57)); % кусок сигнала
s=repmat(s,1,25); % весь сигнал
pwelch(s,rectwin(60),1,60,50*60); % посмотреть на графике
[pxx,w]=pwelch(s,rectwin(60),1,60,50*60); % получить в переменную
sqrt(sum(pxx*50)) % весь rms 50 - это сетка частот
sqrt(sum(s.*s)/1500) % то-же, но по-человечески
rmsp=sqrt((pxx(2)+pxx(4)+pxx(6))*4) % rms в полосе

4 герца взялись не из длины окна, а из того, что я изначально предложил пропускать сигнал через фильтр (идеальный), пропускающий +-2 герца вокруг 3-х гармоник. И рмс ему считать без всяких фурьёв. Вот и 4 для каждой полосы. Вот и 4*3=12 для шума. Окно там не причем, в реальности фильтр и IIR можно сделать.

Кстати похоже что моя ошибка - это рассматривать данный сигнал как случайный процесс... Или просто как-то неправильно рассматривать.

• Ответ: — -=ВН=-   (09.07.2004 19:38, 820 байт)
  • Ок. Большое спасибо за разъяснения (+) — SM   (09.07.2004 19:47, 278 байт)
    • —    (14.11.2004 19:52, 281 байт)
• ну не [1 1 1, а [1.5 1.5 1.5 естессно в первой строке — SM   (09.07.2004 19:13, пустое)