music1.m

music算法

clear all;
close all;
clc;

N=100;
fs=1e3;
f0=100;
fbo=1000/20;
f1=30;
M = 2; % Use 4-ary modulation.
xx = randint(1,N/5,M); % Random digital message
s1=dmod(xx,f0,fbo,fs,'psk',M);

t=0:0.1:pi;
%rs=int('ss(t+tao/2)*conj(ss(t-tao/2)*exp(j*2*pi*f0*t))','tao','-100','100');
zeta1=pi/6;
a=[1,exp(-j*pi*sin(zeta1)),exp(-j*2*pi*sin(zeta1)),exp(-j*3*pi*sin(zeta1)),exp(-j*4*pi*sin(zeta1)),exp(-j*5*pi*sin(zeta1))]';
xx=a*s1;
x=awgn(xx,10);
rx=x*x';
[v,d]=eig(rx);
v5=v(:,2:6);
vv=v5*v5';
i=1;
for zeta1=0:pi/30:2*pi;
a1=[1,exp(-j*pi*sin(zeta1)),exp(-j*2*pi*sin(zeta1)),exp(-j*3*pi*sin(zeta1)),exp(-j*4*pi*sin(zeta1)),exp(-j*5*pi*sin(zeta1))]'
p(1,i)=1/(a1'*vv*a1); 
i=i+1;
end
i=1:61;
plot((i-1)*6,p(i));