matlab功率谱密度函数

MATLAB中可以使用函数`pwelch`计算功率谱密度。该函数的基本用法如下：

[Pxx, f] = pwelch(x, window, noverlap, nfft, fs)

其中，`x`是输入信号，`window`是窗函数，`noverlap`是重叠窗口的长度，`nfft`是FFT点数，`fs`是采样频率。函数的返回值`Pxx`是功率谱密度，`f`是频率向量。

下面是一个简单的示例代码：

% 生成信号
fs = 1000; % 采样频率为1000Hz
t = 0:1/fs:1-1/fs; % 时间向量
x = sin(2*pi*100*t) + randn(size(t)); % 信号

% 计算功率谱密度
[Pxx, f] = pwelch(x, [], [], [], fs);

% 绘制结果
plot(f, 10*log10(Pxx)); % 将功率转换为分贝单位
xlabel('Frequency (Hz)')
ylabel('Power/Frequency (dB/Hz)')

这段代码生成一个包含100Hz正弦波和高斯白噪声的信号，并计算其功率谱密度。最后将结果绘制出来。