功率谱密度

最新推荐文章于 2023-11-16 10:17:19 发布

Supre_yuan

最新推荐文章于 2023-11-16 10:17:19 发布

阅读量4.8k

点赞数 4

文章标签： matlab 信号处理

本文链接：https://blog.csdn.net/weixin_43880225/article/details/120598147

版权

这篇博客介绍了如何在MATLAB中计算和绘制功率谱密度（PSD），包括使用周期图法、Welch方法以及不同窗口函数的影响。通过示例代码展示了从原始信号到不同方法计算PSD的过程，并通过 subplot 进行了可视化比较。

摘要由CSDN通过智能技术生成

功率谱密度简介

PSD（Power Spectrum Density）。PSD定义为单位频带内的信号功率。

MATLAB中PSD函数

Hpsd = dspdata.psd(Data)

最简单的一种计算psd方式。

Hpsd = dspdata.psd(Data, Frequencies)

计算功率谱密度的频率矢量。这个矢量的范围取决于谱型值。对于单侧，默认范围可以为[0, Fs/2]。对于双侧，可以为[0,Fs]。

Hpsd = dspdata.psd(...,'Fs',Fs)

使用采样率Fs。

Matlab里提供很多PSD工具，比如：周期图法、Welch法、Yule-Walker法、Burg法、协方差法。

代码

fs = 5000;    % 规定采样率为5000Hz
t = (0 : fs/10) ./ fs;    % 规定时间序列
A = [1 2];
f = [150; 140];
xn = A * sin(2 * pi * f * t) + 0.1*randn(size(t));

subplot(321)    % 原信号
plot(t, xn)
title('原始信号');
xlabel('时间（Sec）');
ylabel('幅值');

subplot(322)    % 周期图估计
Hrect = spectrum.periodogram;
psd(Hrect, xn, 'Fs', fs, 'NFFT', 1024);

subplot(323)    % hamming窗估计周期图
Hhamm = spectrum.periodogram('Hamming');
psd(Hhamm, xn, 'Fs', fs, 'NFFT', 1024);

subplot(324)    % 矩形窗周期图估计
Hs = spectrum.periodogram('rectangular');
psd(Hs, xn, 'Fs', fs, 'NFFT', 1024);

subplot(325)    % welch法周期图估计，使用150点，重合50点
Hs = spectrum.welch('rectangular', 150, 50);
psd(Hs, xn, 'Fs', fs, 'NFFT', 512);

subplot(326)    % welch法周期图估计，使用100点，重合75点
Hs = spectrum.welch('rectangular', 100, 75);
psd(Hs, xn, 'Fs', fs, 'NFFT', 512);