用MUSIC算法计算随机信号的功率谱密度，matlab代码

MUSIC算法可以用于估计随机信号的功率谱密度，下面是一个简单的MATLAB代码示例：

% 生成随机信号
N = 1024;
x = randn(N,1);

% 设置MUSIC算法参数
n = 8; % 阵元数
d = 0.5; % 阵元间距
theta = linspace(-90,90,181); % 角度搜索范围

% 构造阵列坐标
pos = (0:n-1)*d;

% 构造信号空间矩阵
S = zeros(n,n);
for k = 1:n
    a = exp(-1i*2*pi*pos(k)*sind(theta)/343);
    S = S + a'*a;
end

% 计算MUSIC谱
[V,D] = eig(S);
[~,idx] = sort(diag(D),'descend');
V = V(:,idx);
En = V(:,2:end);
Pmusic = zeros(size(theta));
for k = 1:length(theta)
    a = exp(-1i*2*pi*pos(1)*sind(theta(k))/343);
    Pmusic(k) = 1/(a'*En*En'*a);
end

% 绘制MUSIC谱
plot(theta,Pmusic);
xlabel('Angle (degree)');
ylabel('Spectrum');

在上面的代码中，我们生成了一个长度为N的随机信号，并且设置了MUSIC算法的参数，包括阵元数、阵元间距和角度搜索范围。接着我们构造了信号空间矩阵，并且计算了MUSIC谱。最后我们绘制了MUSIC谱。