Matlab: Pwelch оконное и усреднение сигнала
У меня есть вопрос о функции pwelch в Matlab. Я записал сейсмический сигнал в течение 20 минут, используя частоту дискретизации 100 выборок в секунду. Теперь я хочу определить оценку спектральной плотности мощности с помощью функции pwelch. Я хочу разделить мой сигнал на 1-минутные окна, т.е. 6000 выборок на окно. Перекрытие должно составлять 50% (по умолчанию).
1-й вопрос: дает ли следующий код то, что я хочу?
Fs = 100; Nyq = Fs/2; N = length(signal); f = 0:(Nyq/(N-1)):Nyq; Nf = length(f); win = 6000; pwelch(signal,win,[],Nf,Fs);2-й вопрос: почему в результирующем графике отображается только один сигнал, а не 20 (по одному на каждое окно)?
- Третий вопрос: получу ли я такой же результат, если разделю данные на 20 частей с 6000 точками для каждой части, нанесу эти части вместе на одной периодограмме и возьмем среднее значение?
Заранее спасибо.