经典谱估计周期图法matlab

以下是经典谱估计周期图法的Matlab代码示例：

% 输入信号
x = randn(1, 1000);

% 周期图法参数设置
N = 256; % FFT长度
M = 64; % 子段长度
L = N - M + 1; % 子段重叠长度
K = fix((length(x)-M)/L) + 1; % 子段数

% 初始化周期图
Pxx = zeros(1, N);

% 计算周期图
for k = 0:K-1
    % 提取子段
    xk = x(k*L+1:k*L+M);
    
    % 加窗
    xk = xk .* hamming(M)';
    
    % FFT计算
    Xk = fft(xk, N);
    
    % 周期图累加
    Pxx = Pxx + abs(Xk).^2;
end

% 平均周期图
Pxx = Pxx / K;

% 绘制频谱图
f = (0:N-1) / N;
plot(f, 10*log10(Pxx));
xlabel('Frequency');
ylabel('Power Spectral Density (dB)');
title('Periodogram Power Spectral Density');

% 显示结果
Pxx