Welch方法matlab代码

以下是使用MATLAB实现Welch方法的示例代码：

% 生成一个长度为N的随机信号
N = 1024;
x = randn(1, N);

% 设置分段长度和重叠长度
L = 128;
overlap = L/2;

% 计算分段数和窗函数
M = fix((N-overlap)/(L-overlap));
window = hamming(L);

% 对每个分段进行傅里叶变换并平方
Pxx = zeros(1, L/2+1);
for i = 1:M
    segment = x((i-1)*(L-overlap)+1:i*(L-overlap)+L);
    segment = segment .* window';
    X = fft(segment);
    Pxx = Pxx + abs(X(1:L/2+1)).^2;
end

% 平均各个分段的傅里叶变换结果，得到谱密度估计值
Pxx = Pxx / M / L;
f = (0:L/2) / L * 2;

% 绘制谱密度估计曲线
plot(f, 10*log10(Pxx))
xlabel('Normalized frequency')
ylabel('Power/frequency (dB/Hz)')

上述代码中，我们首先生成一个长度为N的随机信号x。然后设置分段长度L和重叠长度overlap，并计算分段数M和窗函数window。接下来，我们对每个分段进行窗函数加权，然后进行傅里叶变换，并将变换结果的平方加到Pxx中。最后，我们平均各个分段的傅里叶变换结果，得到谱密度估计值Pxx，并绘制谱密度估计曲线。