关于采样频率定义的问题

叶惜

最近在做傅里叶变换和功率谱分析，对于里面的频率Fs产生了疑惑。采样频率的定义查到是一秒钟采样的个数，但是我现在有一段资料，原始资料是日平均资料，进行分析的时候的资料是我做了季节平均的了。做功率谱分析的时候我就产生了疑惑，首先第一个问题是，采样频率应该是算季节平均之后的吧，第二问题是，怎么算呢？是24（小时）*3600（秒）*90（天）HZ，是这样吗?我被这个问题困扰了很久，感觉绕不出去了，也查不到答案，请大家帮帮忙，感谢了。

叶惜

不好意思啊，第一次发帖紧张选错类别了，应该是matlab的，自己好像改不了{:eb303:}{:eb303:}

二爷名声在外

麻烦楼主稍微详细一点点…第一个问题没太看明白，第二个问题是对的

言深深

变换分析是针对数据而言的，变换方式并不会因为你是月资料，或者日资料有所变化，换而言之，变换方法的本身是不识别你的资料属于什么类型的。因此···这频率就是你输入的频率【此处频率等价于周期，因为我们知道二者其实是一回事】

叶惜

二爷名声在外 发表于 2014-7-20 21:48
麻烦楼主稍微详细一点点…第一个问题没太看明白，第二个问题是对的

就是说采样频率是按照一天来算呢，还是按照季节来算呢？不过你说第二个问题是对的话，意思就是应该按照季节来算的吧。

叶惜

言深深 发表于 2014-7-20 22:06
变换分析是针对数据而言的，变换方式并不会因为你是月资料，或者日资料有所变化，换而言之，变换方法的本身 ...

那我做功率谱分析之后得到的频率也可以用这样换算成年吗？比如说，我得到的某个峰值对应的频率是100000hz，然后我把他就能转换成100000s，是这样吗？

言深深

叶惜 发表于 2014-7-21 09:33
那我做功率谱分析之后得到的频率也可以用这样换算成年吗？比如说，我得到的某个峰值对应的频率是100000hz ...

不是吧?你的序列，假如你的序列长度是100，那么你的周期就是1-100之间，计算的频率就是2pi/100-2pi之间，这和秒什么的没有关系的。

叶惜

言深深 发表于 2014-7-21 11:48
不是吧?你的序列，假如你的序列长度是100，那么你的周期就是1-100之间，计算的频率就是2pi/100-2pi之间， ...

要的就是这个，为什么计算的频率就是2pi/100-2pi之间呢？应该就是这样这里不明白啊
你看一下我做的图。就像知道这个图的横坐标要怎么转成周期呢。序列是一个季度一次的，共176个数据

用的是这样的程序。
Fs = 7776000;   
h = spectrum.welch;                  % Create a Welch spectral estimator.
psd(h,x,'Fs',Fs);                    % Calculate and plot the one-sided psd.
hpsd = psd(h,x,'ConfLevel',.98);     % psd with confidence level
figure,plot(hpsd)

言深深

叶惜 发表于 2014-7-21 14:13
要的就是这个，为什么计算的频率就是2pi/100-2pi之间呢？应该就是这样这里不明白啊
你看一下我做的图。 ...

频率和周期的乘机等于2pi，这个你知道的吧？就是这么定义的啊