matlab代码，一段含有噪声的语音信号，利用相关图法的三角窗与矩形窗进行信号的功率图估计

以下是基于相关图法的三角窗与矩形窗对含有噪声的语音信号进行功率图估计的Matlab代码：

% 读取含有噪声的语音信号
[s, fs] = audioread('noisy_voice.wav');

% 设置分析参数
N = 8192;
Nfft = 8192;
n0 = 1000;

% 提取样本信号
x = s(n0 : n0+N-1);

% 计算自相关函数
r_rect = zeros(2*N/2-1, 1);
r_tri = zeros(2*N/2-1, 1);
for k = 1 : N/2
    % 矩形窗
    x1 = x(k : N);
    x2 = x(1 : N+1-k);
    r_rect(N/2+k-1) = x1' * x2 / N;
    r_rect(N/2-k+1) = r_rect(N/2+k-1);

    % 三角窗
    w = triang(N-k+1);
    x1 = x(k : N) .* w;
    x2 = x(1 : N+1-k) .* w;
    r_tri(N/2+k-1) = x1' * x2 / sum(w.^2);
    r_tri(N/2-k+1) = r_tri(N/2+k-1);
end

% 计算功率谱
f = (0 : Nfft/2-1)*fs / Nfft / 1000;

% 矩形窗
rx_rect = r_rect;
Sxz_rect = fft(rx_rect, Nfft);
Sxdbz_rect = 10*log10(abs(Sxz_rect(1 : Nfft/2)));

% 三角窗
rx_tri = r_tri;
Sxz_tri = fft(rx_tri, Nfft);
Sxdbz_tri = 10*log10(abs(Sxz_tri(1 : Nfft/2)));

% 绘制功率谱图
figure;
subplot(2,1,1); plot(f, Sxdbz_rect); ylabel('强度 (dB)');
xlabel('频率 (kHz)'); title('矩形窗功率谱');
subplot(2,1,2); plot(f, Sxdbz_tri); ylabel('强度 (dB)');
xlabel('频率 (kHz)'); title('三角窗功率谱');

其中，语音信号文件名为'noisy_voice.wav'，可以根据实际情况进行修改。运行该代码，即可得到基于相关图法的三角窗与矩形窗对含有噪声的语音信号进行功率图估计的结果。