急求！！请小伙伴帮忙做个实验，matlab功率谱分析这组数据得到图是啥样的？

LemonTree007

请小伙伴帮忙做个实验，matlab功率谱分析这组数据得到图是啥样的？

LemonTree007

数据在这儿，不用下载。
182
189
175
181
198
204
126
168
143
185
214
205
217
218
92
76
101
152
178
193
197
220
211
222
169
158
202
224
107
102
42
11
25
115
122
133
156
122
156
76
40
49
38
38
90
111
127
65
92
124
140
64
77
26
71
85
53
54
75
98
67
72
117
109
152
102
112
142
129
75
65
53
59
73
58
80
41
29
67
50
124
53
62
26
25
61
57
76
84
104
105
121
110
71
90
81
87
41
12
67
67
25
44
58
97
18
102
82
103
84
99
99
68
85
33
155
83
117
111
38
56
64
41
18
40
38
40
48
27
51
68
76
90
101
85
74
72
74
45
31
68
89
102
102
68
35
13
24
20
17
18
37
23
50
81
126
73
62
91
125
84
11
20
97
119
130
50
17
6
40
71
117
123
68
29
53
47
104
57
138
140
107

帅帅f

% 导入数据
data = [182 189 175 181 198 204 126 168 143 185 214 205 217 218 92 76 101 152 178 193 197 220 211 222 169 158 202 224 107 102 42 11 25 115 122 133 156 122 156 76 40 49 38 38 90 111 127 65 92 124 140 64 77 26 71 85 53 54 75 98 67 72 117 109 152 102 112 142 129 75 65 53 59 73 58 80 41 29 67 50 124 53 62 26 25 61 57 76 84 104 105 121 110 71 90 81 87 41 12 67 67 25 44 58 97 18 102 82 103 84 99 99 68 85 33 155 83 117 111 38 56 64 41 18 40 38 40 48 27 51 68 76 90 101 85 74 72 74 45 31 68 89 102 102 68 35 13 24 20 17 18 37 23 50 81 126 73 62 91 125 84 11 20 97 119 130 50 17 6 40 71 117 123 68 29 53 47 104 57 138 140 107];

% 采样频率
Fs = 1000;

% 计算窗口长度
winlen = 2^nextpow2(length(data));

% 使用pwelch函数计算功率谱密度
[p,f] = pwelch(data,[],[],[],Fs);

% 绘制功率谱图
plot(f,p);
xlabel('频率 (Hz)');
ylabel('功率')

LemonTree007

帅帅f 发表于 2023-3-13 16:54
% 导入数据
data = [182 189 175 181 198 204 126 168 143 185 214 205 217 218 92 76 101 152 178 193 19 ...

谢谢，你的代码是ncl吗？我是用的Python，出来的图感觉有问题，所以想知道正确的图是怎样的

帅帅f

LemonTree007 发表于 2023-3-20 16:12
谢谢，你的代码是ncl吗？我是用的Python，出来的图感觉有问题，所以想知道正确的图是怎样的

matlab代码