声振论坛

 找回密码
 我要加入

QQ登录

只需一步,快速开始

查看: 3309|回复: 5

[综合讨论] 求教:傅立叶变换求时间序列平均周期P的matlab程序?

[复制链接]
发表于 2006-10-23 19:27 | 显示全部楼层 |阅读模式

马上注册,结交更多好友,享用更多功能,让你轻松玩转社区。

您需要 登录 才可以下载或查看,没有账号?我要加入

x
求教:傅立叶变换求时间序列平均周期P的matlab程序?
哪位大哥有相关程序啊,最好带注释,先谢谢了!!
回复
分享到:

使用道具 举报

 楼主| 发表于 2006-10-23 19:53 | 显示全部楼层
Y = fft(wolfer);

Y的幅度是功率谱,画出功率谱和频率的对应关系就得出了周期图,去掉第一点,因为他只是所有数据的和,画图有:

N = length(Y);
Y(1) = [];
power = abs(Y(1:N/2)).^2;
nyquist = 1/2;
freq = (1:N/2)/(N/2)*nyquist;
plot(freq,power),
grid on
xlabel(’cycles/year’)
title(’Periodogram’)
这段代码是什么意思?
period = 1./freq;
plot(period,power),
axis([0 40 0 2e7]),
grid on
ylabel(’Power’)
xlabel(’Period(Years/Cycle)’)
发表于 2006-10-24 19:21 | 显示全部楼层
不知道平均周期是一个什么样的概念?
发表于 2007-1-5 21:21 | 显示全部楼层
请问这个程序是计算什么的,可以算地震波的卓越周期吗
发表于 2007-1-8 13:12 | 显示全部楼层
平均周期P是个什么参数?你是要问平均次数嘛?
发表于 2007-7-24 06:51 | 显示全部楼层
原帖由 liutuo 于 2006-10-23 19:53 发表
Y = fft(wolfer);

Y的幅度是功率谱,画出功率谱和频率的对应关系就得出了周期图,去掉第一点,因为他只是所有数据的和,画图有:

N = length(Y);
Y(1) = [];
power = abs(Y(1:N/2)).^2;
nyquist = 1/2;  ...

由DFT的变换对的定义第一点为所有数据的和,但为什么要去掉呢?
您需要登录后才可以回帖 登录 | 我要加入

本版积分规则

QQ|小黑屋|Archiver|手机版|联系我们|声振论坛

GMT+8, 2024-11-18 20:29 , Processed in 0.066974 second(s), 18 queries , Gzip On.

Powered by Discuz! X3.4

Copyright © 2001-2021, Tencent Cloud.

快速回复 返回顶部 返回列表