music算法的matlab程序

clc
clear
close all

% MUSIC algorithm
% standar
c = 3e8;

f0 = 24.15e9;
lambda = c/f0;
d = lambda*0.5;

N_t = 3;
theta_t = [-7 -3.5 0 3.5];
P = 8;

N = 128;
p = 0:P-1;
% amp = randn（1N_t）;
% phi = rand（1N_t）;
% s = （amp‘.*exp（1i*phi‘）） * exp（-1i*2*pi*0.1*（0:N-1））;
s = randn（N_tN）;
% n = （randn（PN） + 1i*randn（PN））;
A = exp（-1i*2*pi*f0/c*p‘*d*sind（theta_t（1:N_t）））;

X = A*s;
X = awgn（X33‘measured‘）;

Rx = （X*X‘）/N;
[EVD] = eig（Rx）;
Un = EV（:1:P-N_t）;

% G = Rx（:1:N_t）;
% H = Rx（:N_t+1 : end）;
% P1 = inv（G‘*G）*G‘*H;
% Q = （[P1‘ eye（P-N_t）]）‘;
% Un = Q;
ang = -30:0.5:30;
A = exp（-1i*2*pi*f0*p‘*d*sind（ang）/c）;
P = 1./diag（A‘*Un*Un‘*A）;
plot（angdb（abs（P）））

% smooth music
% c = 3e8;
% 
% f0 = 24.15e9;
% lambda = c/f0;
% d = lambda*0.5;
% 
% N_t = 4;
% theta_t = [-7 -3.5 0 3.5];
% P1 = 128;
% 
% N = 128;
% p1 = 0:P1-1;
% 
% P = 8;
% p = 0:P-1;
% 
% amp = randn（1N_t）;
% phi = rand（1N_t）;
% s = （amp‘.*exp（1i*phi‘）） * exp（-1i*2*pi*0.1*（0:N-1））;
% % s = randn（N_tN）;
% % n = （randn（PN） + 1i*randn（PN））;
% A = exp（-1i*2*pi*f0/c*p1‘*d*sind（theta_t（1:N_t）））;
%