clear?all;
N?=?1000;
noise?=?（?randn（1??N）?+?1i?*?randn（1??N））?/?sqrt（2）;
%產生帶噪聲的信號樣本u（n）
signal1?=??exp（?1i?*?0.5?*?pi?*?（0?:?N-1）?+?1i?*?2?*?pi?*?rand）;
signal2?=??exp（?-1i?*?0.3?*?pi?*?（0?:?N-1）?+?1i?*?2?*?pi?*?rand）;?
un?=?signal1?+?signal2?+?noise;
%計算自相關矩陣
M?=?8;
?for?k=1:N-M
?????xs（:k）=un（k+M-1:-1:k）.‘;
?end?
R=xs*xs‘/（N-M）;
%自相關矩陣的特征值分解
[UE]?=?svd（R）;
ev?=?diag（E）;

%根據AIC準則進行信號源個數的估計
for?k?=?1:M
????dec?=?prod（ev（k:M）.^（1/（M-k+1）））;?%計算第一項中對數的自變量的分子
????nec?=?mean（ev（k:M））;%計算第一項中對數的自變量的分母
????lnv?=?（dec?/nec）^（（M-k+1）*N）;%計算第一項中對數的自變量
????AIC（k）?=?-2?*?log（lnv）?+?2?*?（k-1）?*?（2?*?M?-?k?+?1）;
end
[AminK]?=?min（AIC）;
N1?=?K?-?1;
%根據MDL準則進行信號源個數的估計
for?k?=?1:M
????