用matlab写一段时延功率谱算法

以下是一个简单的时延功率谱算法的Matlab代码：

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

% 时延范围
tau_max = 100;

% 时延功率谱矩阵
P = zeros(tau_max, length(x));

% 计算时延功率谱
for tau = 1:tau_max
for n = 1:length(x)-tau
P(tau, n) = x(n)*x(n+tau);
end
end

% 对时延功率谱进行平均
P = mean(P, 2);

% 绘制时延功率谱图像
plot(P);
xlabel('时延');
ylabel('功率谱密度');
title('时延功率谱');