我在matlab中使用pwelch方法來計算一些風速測量的功率譜。所以,到目前爲止,我寫了下面的代碼爲例: t = 10800; % number of seconds in 3 hours
t = 1:t; % generate time vector
fs = 1; % sampling frequency (seconds)
A = 2; % amplitude
P = 1000;
我嘗試使用whelch方法,我發現零頻率異常 import numpy as np
import scipy.signal as signal
import matplotlib.pyplot as plt
n = 100000
s = np.ones(n)
f, psd = signal.welch(s, return_onesided=False)
plt.plot(f, psd
我想知道我的數據中是否有6個月的季節性。我想用光譜分析來做到這一點,我會在特定的頻率下觀察頻譜的Ftest。 我尋找有助於瞭解什麼x軸在譜分析裝置(即什麼頻率的意思),如果我計算它基於與頻率12. 這裏一些數據的時間序列的對象上產生。每個月我的數據被收集超過10年: data = ts(1:120, frequency = 12, start = c(1959, 1)) # data that i