%%?MUSIC算法
%%非相干輸入信號
%%四信號輸入八陣元的線陣
clc;
close?all;
clear?all;
i=sqrt（-1）;
j=i;
M=8;??????????????????????????????????????%陣元數
K=4;??????????????????????????????????????%信源數
lma=2;????????????????????????????????????%信號中心波長
d=lma/2;??????????????????????????????????%天線陣元間距
angle1=-40;???????????????????????????????%四個輸入信號的入射角
angle2=0;?
angle3=20;?
angle4=60;?
th=[angle1;angle2;angle3;angle4];
degrad=pi/180;
SNR=[10;10;10;10];
%%?構造信號和高斯白噪聲
n=1024;
nn=1:n;
v1=0.015;????????????????????????????????????????????????%四個輸入信號的頻率
v2=0.05;
v3=0.02;
v4=0.035;
A=exp（-2*pi*d*j*[0:M-1]‘*sin（th‘*degrad）/lma）;????????????%陣列因子
S=[1*cos（v1*nn）;1*sin（v2*nn）;1*sin（v3*nn）;1*sin（v4*nn）];??%輸入信號
U=randn（Mn）+j*randn（Mn）;????????????????????????????????%構造噪聲源
Ps=diag（S*S‘/n）;??????