分享一个今天刚编写的简单功率谱分析函数

风叶ele · 发表于 2019-3-1 22:06:32

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

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

x

有关功率谱分析的帖子家园里已经有很多了，但是我感觉有的帖子说得还是让人摸不着头脑。
于是今天自己按照魏凤英课本的P72所介绍的“直接傅里叶变换”的方法，编写了一个简单脚本。

由于我自己是要对年序列分析（只看年际、年代际周期），所以该函数只需要输入一个简单的年序列即可，返回谱值和周期，后期可进一步绘图。

说明如下：

POWER_SPECTRUM power spectrum analysis for year series
[ ] = powerspectrum_for_yearseries(x)
计算一个简单时间序列（年序列）的功率谱，返回周期以及相应的谱值
附有一段用于绘制谱图的代码（但被注释掉，仅供参考）
按照“直接使用傅里叶变换”的方法进行计算，详见魏凤英《现代气候统计诊断与预测》（第二版）P72
输入 x：年序列，例如1951-2018年冬季Nino3.4指数序列
输出 s：年序列x的谱值；T：对应的周期，例如10年周期、2年周期等

powerspectrum_for_yearseries.m (1.45 KB, 下载次数: 202)

sealsu · 发表于 2020-3-30 22:00:48

我测试看看

shen159876 · 发表于 2020-4-29 20:10:48

周期是不是就是拐角？

风叶ele · 发表于 2020-5-21 11:38:21

shen159876 发表于 2020-4-29 20:10
周期是不是就是拐角？

嗯，明显的峰值区，不过严格来讲还需要配合显著性检验的曲线配合来看，我这里没有画检验的线

si18 · 发表于 2020-5-24 10:00:34

楼主，显著性检验的线应该怎样画

大物小子 · 发表于 2020-5-24 11:43:59

赞受益匪浅

2daidai · 发表于 2020-11-24 15:51:05

楼主，显著性检验的线应该怎样做？

人可山 · 发表于 2020-12-10 09:20:56

谢谢，楼主的分享，刚好需要用到功率谱

Monicaaaaaaa · 发表于 2021-3-21 12:30:37

谢谢楼主的分享棒

livteng · 发表于 2021-8-10 21:53:19

您好请问一下逐日的数据这个可以用吗

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

[程序设计] 分享一个今天刚编写的简单 功率谱分析 函数

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

[程序设计] 分享一个今天刚编写的简单功率谱分析函数