给我一段MUSIC算法求随机信号功率谱的matlab代码

以下是一段基于Matlab的MUSIC算法求随机信号功率谱的代码示例：

% 生成随机信号
N = 1000; % 信号长度
fs = 1000; % 采样频率
t = (0:N-1)/fs; % 时间序列
x = sin(2*pi*50*t) + sin(2*pi*150*t) + randn(1,N); % 信号

% 构建阵列
N = 8; % 阵列元素个数
d = 0.5; % 阵列元素间距
theta = -90:0.5:90; % 信号方向范围
A = zeros(N,length(theta));
for n = 1:N
    A(n,:) = exp(1i*2*pi*d*(n-1)*sind(theta));
end

% 计算协方差矩阵和特征值分解
Rxx = x*x'/N; % 协方差矩阵
[V,D] = eig(Rxx); % 特征值分解
[~,I] = sort(diag(D),'descend'); % 特征值排序
V = V(:,I);

% 计算信号子空间和噪声子空间
Un = V(:,N+1:end);
U = V(:,1:N);

% 计算MUSIC谱
Pmusic = zeros(1,length(theta));
for th = 1:length(theta)
    a = exp(1i*2*pi*d*(0:N-1)'*sind(theta(th)));
    Pmusic(th) = 1/(a'*(Un*Un')*a);
end

% 绘制功率谱曲线
figure;
plot(theta,10*log10(abs(Pmusic)/max(abs(Pmusic))));
xlabel('信号方向（度）');
ylabel('功率谱（dB）');
title('随机信号功率谱估计');

以上代码将生成一个随机信号，并在一个8个元素的均匀线阵上进行信号方向估计，最终绘制出随机信号的功率谱曲线。其中，通过计算协方差矩阵和特征值分解得到信号子空间和噪声子空间，再利用MUSIC谱估计算法计算出信号源的方向，并绘制出功率谱曲线。