f1=0.1.f2=0.13,给我一段matlab代码，使用Levinson-Durbin递推求该随机信号功率谱

好的，以下是使用Levinson-Durbin递推求随机信号功率谱的Matlab代码：

% 定义随机信号
f = [0.1, 0.13];
a = [1, -0.8, 0.64, -0.512, 0.4096];
N = length(a)-1;
r = zeros(1,N+1);
for k = 0:N
    r(k+1) = a(1:N-k+1)*a(k+1:N+1)';
end

% 使用Levinson-Durbin递推求解自回归模型系数和功率谱
P = zeros(1,N+1);
E = r(1);
a_hat = zeros(N+1,N+1);
a_hat(1,:) = 1;
for n = 1:N
    k = -a_hat(n,1:n)*flipud(a(1:n+1))'/E;
    a_hat(n+1,1:n+1) = [a_hat(n,1:n)+k*flipud(a_hat(n,1:n))', k];
    E = E*(1-k^2);
    P(n+1) = E;
end
P(1) = a_hat(N+1,:)*r';

% 绘制功率谱曲线
f_axis = linspace(0, 0.5, 1000);
H = freqz(1,a_hat(N+1,:),2*pi*f_axis);
P = P/P(1);
H = H.*conj(H);
plot(f_axis, H);
xlabel('Frequency (cycles/sample)');
ylabel('Power Spectrum');
title('Power Spectrum of Random Signal');

这段代码中，我们首先定义了一个随机信号，然后计算该信号的自相关函数。接着使用Levinson-Durbin递推算法求解自回归模型系数和功率谱，最后绘制功率谱曲线。需要注意的是，我们在计算功率谱时将其归一化，使得最大功率为1。