将功率谱密度命令的结果存储在矩阵中 - matlab

Question

我需要一些有关 for 循环和功率谱密度的 matlab 代码的帮助。

我正在分析一个信号，需要能够测量每 30 秒段的平均功率。我已经用函数'buffer'分割信号，这样矩阵segment_mat的每一列中都有30秒的数据，但我希望能够使用for循环将功率谱密度的结果存储在矩阵中循环遍历数据并将每个段的平均功率存储在向量avpow 中。有什么想法我可能会出错吗？

for j=1:120
    Hs=spectrum.periodogram({'Hamming'});
    Power=psd(Hs,segment_mat(:,j)',fs',fs);
    avpow=mean(Power(j))';
end

Answer 1

您的问题是您在 for 循环的每次迭代中都覆盖了结果。

只需修改您的代码：

Power(j) = ...
avpow(j) = ...

您将存储计算结果。如果速度有问题，请考虑预先分配。