единицы измерения для временных интервалов температуры

У меня есть временная серия температур воздуха, измеренная в градусах С, а частота сигнала - ежемесячно. Я вычисляю спектры мощности временных рядов как:единицы измерения для временных интервалов температуры

L = length(temp); 
NFFT = 2^nextpow2(L); % Next power of 2 from length of y 
Y = fft(temp,NFFT)/L; 
Fs = 1; % one sample per month 
f = Fs/2*linspace(0,1,NFFT/2+1); % frequency 
Pxx = 2*abs(Y(1:NFFT/2+1)); % power 

% Plot single-sided amplitude spectrum. 
plot(f,Pxx,'-k'); 
xlabel('Frequency (c/month)'); % cycles per unit time

, что приводит к следующему разделу.

Что единицы должны быть для YAxis? Должно ли это быть^{o} C?

Диапазон временных рядов составляет около 20, так что, думаю 10 будет работать как амплитуда ...

источник

2016-01-13 Emma Tebbs

Вы видели http://stackoverflow.com/questions/1523814/units-of-a-fourier-transform-fft-when-doing-spectral-analysis-of-a-signal? Кроме того, я думаю, что единицы оси x должны составлять 1 месяц, а не C/месяц (например, температура поднимается и опускается каждый день, независимо от того, сообщается ли она в F или C, это еще 1-дневный цикл). – tom10

Спасибо. С/месяц относится к циклам в месяц. В этом случае это 0.083 цикла в месяц, что составляет примерно один цикл в год. –

Извините, что вы придирчивы, но так как ваш другой блок - «C», из трех элементов по оси x, чтобы не произнести полностью, «циклы «вероятно, не лучший выбор. – tom10

Допустим, первоначальные оси были Т (по вертикали) и т (по горизонтали); после преобразования они равны Y и X.

Вы правильно поняли, что единицы X являются обратными единицами t. Это необходимо для таких формул, как exp(2*pi*i*t*X): продукт t*X должен быть безразмерным.

Parseval's identity может использоваться для вывода единиц по вертикальной оси. В соответствии с этим,

(T единиц) (т единиц) = (Y единиц) (блоки X)

, следовательно, Y единиц = (T единиц) (T единиц). В вашем случае это градус-раз-секунда, но это лучше выражено как степень на частоту единица.

Действительно, функция, созданная преобразованием Фурье, равна амплитуде плотности, выражая, сколько колебаний существует в каждом конкретном диапазоне частот.

источник

2016-01-14 01:15:11

Спасибо за ответ. Таким образом, ось y в моем случае может быть | P (f) | (degC), где f обозначает «частоту» - таким образом, представляет degC на единицу частоты, и поэтому на пике демонстрируется 10 градусов C на этой конкретной частоте? Или я должен просто написать его как «плотность амплитуды» или | P (f) | ? –

«Амплитудная плотность» легче проглотить, чем странные единицы. Кстати, то, что я написал, основано на понимании того, что переменная вашей горизонтальной оси непрерывна (как предлагает линейный график). Если бы у вас было это дискретно, со столбцом для каждой частоты, то единицы были бы просто градусами. Это то же самое различие, что и между дискретными и непрерывными случайными величинами: вероятностями и плотностью вероятности. –

единицы измерения для временных интервалов температуры

ответ

Смежные вопросы