我想从太阳的时间序列数据中确定大小频率分离。下面绘制了数据(时间序列和幂序列)的摘录。
幂级数在 MATLAB 中计算如下:
n = length(t_obs);
dt = diff(t_obs(1:2));
y = fft(d_obs, n);
P = y .* conj(y)/n;
f = (0:n/2)/(n*dt);
f = f(1:(n/2));
P = P(1:(n/2));
plot(f, P);
我无法理解的是:
- 我怎样才能获得大的频率分离和小小的分离为了和无需从图中手动读取它(有很多数据集)。如果你想展示一个实现这个的例子(尽管解释或提示就足够了),我精通 Python 和 MATLAB。
- powerseries上y轴的单位是什么?
时间序列图:
幂级数图:
以供参考:
我希望你能帮忙。