NCL specx_anal计算功率谱的横坐标？

一顾留白 · 发表于 2022-3-6 21:46:14

登录后查看更多精彩内容~

您需要登录才可以下载或查看，没有帐号？立即注册

x

我想看一下海温变化的次季节周期，在ncl官网找到了计算的函数和步骤（specx_anal (ucar.edu)），但是对横坐标还是不太理解，我参考了spec3的脚本和图像（NCL: Spectral Analysis and Complex Demodulation (ucar.edu)），图中纵坐标是方差，横坐标是频率，但是单位是cycles/month，那我如果想看0-90天的情况是不是横坐标需要设成1/90-1?

function specx_anal ( x

: numeric, iopt [1] : integer, jave [1] : integer, pct [1] : numeric )
另外参数的不同似乎对结果影响很大，iopt=0求距平，=1去趋势+求距平有何差异呢？ jave选取是否和我想得到的周期有关呢？数值应该如何选取？

你是最美七月天 · 发表于 2022-3-7 00:26:02

感谢分享！学习一下！

stepdance · 发表于 2022-3-7 08:09:03

线性趋势对应于长周期低频分量，差异就在于你的研究目的是否需要包含这部分分量

一顾留白 · 发表于 2022-3-7 09:59:48

啊理解了，谢谢~~

一顾留白 · 发表于 2022-3-8 09:52:20

来自问自答一波
1.横坐标单位是频率，取决于输入数据的单位，输入数据是每天的，那横坐标就是1/day
2.jave表示是否平滑，按需选择，理论上来说与我想要的周期无关

雨子哥 · 发表于 2022-4-2 15:39:49

楼主是用NCL做的谱分析吗？学习一下！

一顾留白 · 发表于 2022-4-3 15:19:29

是用ncl做的，官网有现成的脚本，很方便，就是要理解一下原理

雨子哥 · 发表于 2022-4-4 20:26:55

请问楼主做过横坐标的刻度缩放吗？我现在要画的图需要横坐标不等间距显示，不知您是否做过呢？

一顾留白 · 发表于 2022-4-5 14:08:33

我之前是把横坐标设置为1/fre，这样应该单位就是day（看你的输入数据是按天还是按月），画图的时候调整坐标轴设置就好了

ghs · 发表于 2022-4-10 09:26:07

我记得我当时用的是tmXBValues，tmXBLabels这两个来设置的先选择刻度线的位置，在确定刻度线的下标

		自动登录	找回密码
密码			立即注册

NCL specx_anal计算功率谱的横坐标？

登录后查看更多精彩内容~

回帖奖励 +5 金钱

回帖奖励 +5 金钱

回帖奖励 +5 金钱

回帖奖励 +5 金钱