%?MUSIC_based_on_ULA
%?Author:?PENGWEILAI
%?Date:?2016-04-20
%?理想均勻線陣模型，9個陣元，三個信號（兩個相干，一個不相干），從103070入射。半徑波長比為0.5
%信噪比為20，采樣2048

close?all;?clear;?clc
%%?Necessary?parameter?and?signal?generating??S1=10S2=30S3=70
K??????=?3;??????????????????????????%?信號個數
M??????=?9;??????????????????????????%?陣元數
N??????=?2048;???????????????????????%?快拍數
Ratio??=?0.5;??????????????????????????%?半徑波長比

f1?????=?10e4;
f2?????=?f1;
f3?????=?3e4;
fs?????=?400e6;
t??????=?（1:N）/fs;
theta1?=?10*pi/180;??????????????????%?俯仰角??
theta2?=?30*pi/180;
theta3?=?70*pi/180;
?
FAI1???=?-2*pi*Ratio*sin（theta1）*（0:M-1）;???????%相移
FAI2???=?-2*pi*Ratio*sin（theta2）*（0:M-1）;?
FAI3???=?-2*pi*Ratio*sin（theta3）*（0:M-1）;?
???????
var_n??=?1;??????????????????????????????????????????????????%?噪聲功率
SNR1???=?10;?????????????????????????????????????????????????%?信噪比
A11?????=?sqrt（var_n*10.^（SNR1/10））;??????????????????????????%?信號幅度
s1?????=?A11*exp（-i*（2*pi*f1*t+2*pi*rand））;???????????????????%?s1（n）?和?s2（n）?只要不成比例即可，否則相干
s2?????=?A11*exp（-i*（2*pi*f2*t+2*pi*rand））;??
s3?????=?A11*exp（-i*（2*pi*f3*t+2*pi*rand））;

A1??????=?exp（i*FAI1）‘;??????????