chen89
2011-03-22, 10:12
求一个时间序列的功率谱密度,采样频率1000,采样1024点,求功率谱密度后总是多4个数据点,程序如下:
Fs=1000; % sample frequency
nFFT=1024; %sample points
y=load('20.txt');
[Pyy,F]=psd(y,nFFT,Fs,hamming(2^10),100);
psdplot(Pyy,F,'','','Sample PSD Plot')
80.00
121.00
121.00
1930754103.34
582224618.65
162011695.99
119220659.20
107660982.24
91904073.50
49826977.29
32915149.02
17284186.79
9433419.55
5865603.87
4083406.50
5164048.98
3666514.83
2704488.81
2035118.86
1969625.62
1543491.65
930662.21
1204306.98
983375.88
451247.86
392591.08
342907.77
290368.79
173074.36
168133.75
325292.19
323148.91
155282.21
186869.43
147028.48
91356.60
57691.76
35775.04
27474.79
53656.75
58722.28
63198.95
108538.09
104919.77
64712.75
52918.32
65289.42
34623.31
35731.66
29798.90
29475.03
13006.81
11826.30
10807.90
14903.87
12853.44
17387.41
14139.86
21216.18
21095.11
24637.21
21022.16
19897.05
16796.11
13002.55
11518.88
8212.26
6924.11
9412.02
7002.95
8022.40
6841.29
6065.95
....
上面是截取的一部分结果,Pyy有516点,应该只有512点,而且前面几个数据明显不对,那位大佬指点一下,十分感激!
Fs=1000; % sample frequency
nFFT=1024; %sample points
y=load('20.txt');
[Pyy,F]=psd(y,nFFT,Fs,hamming(2^10),100);
psdplot(Pyy,F,'','','Sample PSD Plot')
80.00
121.00
121.00
1930754103.34
582224618.65
162011695.99
119220659.20
107660982.24
91904073.50
49826977.29
32915149.02
17284186.79
9433419.55
5865603.87
4083406.50
5164048.98
3666514.83
2704488.81
2035118.86
1969625.62
1543491.65
930662.21
1204306.98
983375.88
451247.86
392591.08
342907.77
290368.79
173074.36
168133.75
325292.19
323148.91
155282.21
186869.43
147028.48
91356.60
57691.76
35775.04
27474.79
53656.75
58722.28
63198.95
108538.09
104919.77
64712.75
52918.32
65289.42
34623.31
35731.66
29798.90
29475.03
13006.81
11826.30
10807.90
14903.87
12853.44
17387.41
14139.86
21216.18
21095.11
24637.21
21022.16
19897.05
16796.11
13002.55
11518.88
8212.26
6924.11
9412.02
7002.95
8022.40
6841.29
6065.95
....
上面是截取的一部分结果,Pyy有516点,应该只有512点,而且前面几个数据明显不对,那位大佬指点一下,十分感激!