????clear?all
????clc
????
????S=5;?????%信噪比
????M=7;?????%陣元數
????N=3;?????%信號數
????nd=100;????%對信號的采樣點個數
????f=[1000000?1000000?1000000];????????%信號頻率矩陣
????t_max=5/1e4;???????%對觀測時間的限制
????ta=[0?40?60];%信號入射角度
????ta=ta*pi/180;
????dbc=0.5;???????%陣元間距d與信號波長之間的比值
????ta_max=pi/2;???%最大張角
????s1=10^（S/20）;????%把信噪比由分貝的單位轉化為對信號振幅的單位
????s2=1;
????j=sqrt（-1）;
????%phape=randn（Nnd）;?%生成服從正態分布的隨機矩陣（N*nd）階??
????%phap=phape/max（max（phape））;
????t=linspace（0t_maxnd）;???????????%等差元素向量
????s=exp（j*2*pi*kron（f‘t））;%+j*phap）;??%輸入信號（N*nd）階
????s0=s1*s;??????????????????????????%（N*nd）階
????n=s2*randn（Mnd）;???%噪聲模型?生成服從正態分布的隨機矩陣（M*nd）階
????a1=sin（ta）;??????????
????a2=（0:（M-1））‘;??????
????a3=kron（a1a2）;??????????%信號到達各陣列的相位差
????a=exp（-1*j*2*pi*dbc*a3）;?%導向矢量陣
????X=a*s0+n;????????????????%信號模型（M*nd）階
????R=X*X‘/nd;
%?????[USV]=svd（R）;?????%奇異值分解，計算主特征值并保留秩，并保留相應的特征向量矩陣
%?disp（S）;
%?Un=U（:N+1:M）;
%?Gn=Un*Un‘;
%?seaching_doa=-90:0.1:90?;%現在的搜尋范圍為-90到90
%?for?i=1:length（seaching_doa）
%?????a_theta=exp（-j*（0:M-1）‘*pi*sin（pi*seaching_doa（i）/180））;???%角度
%?????Pmusic（i）=1./abs（（a_theta）‘*Gn*（a_theta））;?%功率譜函數
%?end
%?plot（seaching_doa10*log（Pmusic）‘r‘）
%?xlabel（‘入射角/degree‘）;
%?ylabel（‘空間譜/dB‘）;
%?legend（‘Music?Spectrum‘）;
%?title（‘經典MUSIC估計‘）;
%?grid?on;
?????%確定初始值?
?????????th1=（-90:0.1:90）‘;
?????????th=th1/180*pi;
?????????q=zeros（1N）;
?????????Q=zeros（1length（th））;
for??????????i=1:N????????????%?其中N表示信號個數
?????????????p1=sin（q）;???
?????????????p2=（0:（M-1））‘;
?????????????p3=kron（p1p2）;
?????????????p4=exp（-1*j*2*pi*dbc*p3）;%第K（K>1）個信號的導向矢量
????????for??k=1:length（th）
?????????????pmu_a1=sin（th（k））;
?????????????pmu_a2=（0:（M-1））‘;
?????????????pmu_a3=pmu_a1*pmu_a2;?
?????????????pmu_a=exp（-1*j*2*pi*dbc*pmu_a3）;
?????????if?????????i==1
????????????????????P=pmu_a*pinv（pmu_a‘*pmu_a）*pmu_a‘;???%初始信號
????????????????????qk=trace（P*R）;????%求矩陣的跡
????????????????????Q（k）=qk;
?????????else???
????????????????????pmu_b=zeros（Mi）;
????????????for?????m=1:（i-1）?????????????%表示其前i-1個數
????????????????????pmu_b（:m）=p4（:m）;
????????????end
????????????????????pmu_b（:i）=pmu_a（:1）;???%表示a（i）
????????????????????P=pmu_b*pinv（pmu_b‘*pmu_b）*pmu_b‘;
????????????????????qk=trace（P*R）;
????????????????????Q（k）=qk;
?????????end???
????????end
????????????????????b=max（Q）;
?????????for????????k=1:length（th）??????
????????????????if?（b-Q（k））==0
?????????????????????c=k;
????????????????????q（i）=th（c）;
????????????????end
?????????end
end

?????%迭代
?????????????d=zeros（1N）;
for?p=1:N
?????for?????i=1:N
?????????????p1=sin（q）;???
?????????????p2=（0:（M-1））‘;
?????????????p3=kron（p1p2）;
?????????????p4=exp（-1*j*2*pi*dbc*p3）;%第K（K>1）個信號的導向矢量
????????for??k=1:length（th）
?????????????pmu_a1=sin（th（k））;
?????????????pmu_a2=（0:（M-1））‘;
?????????????pmu_a3=pmu_a1*pmu_a2;?
?????????????pmu_a=exp（-1*j*2*pi*dbc*pmu_a3）;
?????????????p4（:i）=pmu_a（:1）;
?????????????P=p4*p

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----
?????文件????????3668??2011-09-01?22:49??AP_ML.m