%%Capon算法DOA估計
%MUSIC算法
%%錯誤：相位s_phase沒有與j相乘，這樣就不是改變相位，而是改變幅度了。
%%問題：頻率相同，是不是相干信號？
%%仿真現像：同頻能估計出所有的入射方向，但是有偏差；不同頻，估計出所有的入射方向，偏差較小；
%%U（1:D）為信號子空間，U（D+1:end）為噪聲子空間。
%%如果信號子空間維數小于信號源個數，就是有信號子空間擴散到了噪聲子空間，這樣會導致某些導...
%%向矢量與噪聲子空間不完全正交，而無法正確估計信號源方向。
%%如果信號子空間維數大于信號源個數，噪聲子空間擴散到了信號子空間，擴散的部分與剩下的噪聲...
%%子空間正交，而出現假峰。
%%由上可知，準確獲知信號源的個數對正確估計信號的方向起著重要作用。此外，相干信號會引起信...
%%號子空間的擴散

clear?all;clc
c=3e8;%光速
s_angle=pi*[-30?10?50]/180;%入射角度
M=8;%陣元數
D=3;%信號源數
s_freq=[3.02e3?3.01e3?3e3];%信號頻率為30MHz
s_wavelength=c./s_freq;
s_sample=7e3;%采樣頻率為70MHzs
d=s_wavelength（2）/2;%陣元間距為半波長
s_phase=[0.1?0?-0.3]*（2*pi）;%初相
s_snapshoot=100;%快拍數
N=normrnd（01Ms_snapshoot）;%不相關噪聲，Mxs_snapshoot
SNR=[10?10?10];%信噪比
%Amp=sqrt（10.^（SNR/10））;%各復信號幅度，其中噪聲功率為1
Amp=sqrt（2*10.^（SNR/10））;%各實正弦信號幅度，其中噪聲功率為1

%入射信號S（Dxs_snapshoot）每一行，為同一信號不同時間的采樣
%%（1）
n=[1:s_snapshoot];
for?k=1:D
????S（k:）=Amp（k）*sin（2*pi*n*s_freq（k）/s_sample+s_phase（k））;
????%S（k:）=Amp（k）*exp（j*（2*pi*n*s_freq（k）/s_sample+s_phase

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----

?????文件???????2490??2010-07-12?20:11??capon.m

-----------?---------??----------?-----??----

?????????????????2490????????????????????1