声振论坛

 找回密码
 我要加入

QQ登录

只需一步,快速开始

查看: 1556|回复: 1

[稳定性与分岔] 平均周期的疑问?

[复制链接]
发表于 2009-5-25 11:31 | 显示全部楼层 |阅读模式

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

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

x
由以下程序求出的平均周期怎么和时间序列的长度相等?按道理是不相等的,不知道问题出在哪里??
大家帮忙解释解释  ,谢谢了


function T_mean=period_mean_fft(data)
%该函数使用快速傅里叶变换FFT计算序列平均周期
%data:时间序列
%T_mean:返回快速傅里叶变换FFT计算出的序列平均周期
data=load('6.txt');
Y = fft(data);     %快速FFT变换
N = length(Y);    %FFT变换后数据长度
Y(1) = [ ];           %去掉Y的第一个数据,它是data所有数据的和
power = abs(Y(1:N/2)).^2;  %求功率谱
nyquist = 1/2;
freq = (1:N/2)/(N/2)*nyquist; %求频率
subplot(121)
plot(freq,power); grid on     %绘制功率谱图
xlabel('频率')
ylabel('功率')
title('功率谱图')
period = 1./freq;                %计算周期
subplot(122)
plot(period,power); grid on  %绘制周期-功率谱曲线
xlabel('周期')
ylabel('功率')
title('周期—功率谱图')
[mp,index] = max(power);       %求最高谱线所对应的下标
T_mean=period(index)        %由下标求出平均周期
回复
分享到:

使用道具 举报

发表于 2009-8-21 09:23 | 显示全部楼层

平均周期的疑问?'6.txt'是什么那?

data=load('6.txt');
'6.txt'是什么那?
发我邮箱看下
yuantongyi@tom.com
您需要登录后才可以回帖 登录 | 我要加入

本版积分规则

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

GMT+8, 2024-12-19 17:07 , Processed in 0.055780 second(s), 18 queries , Gzip On.

Powered by Discuz! X3.4

Copyright © 2001-2021, Tencent Cloud.

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