function?Pmusic=fMUSIC（datafsflfhanglesxnumnPSDsPSD）
%fMVM?水聲線列陣的頻域MVDR波束形成
%???Pmusic=fMUSIC（datafsflfhangles[x[nPSDsPSD]]）?
%
%???data????--?陣列接收的數據，每行為一次快拍，每列為一個陣元輸出
%???fs??????--?時間采樣頻率（Hz）
%???fl??????--?系統通頻帶的低端頻率（Hz）
%???fh??????--?系統通頻帶的高端頻率（Hz）
%???angles??--?欲形成波束的主極大方向（°），如?[-90:2:90]
%???x???????--?線列陣各陣元的坐標（米）。如?[0:0.05:0.3]
%??????????????如果沒有給定，默認為均勻線列陣，間距為中心頻率對應的半波長
%???num?????--?目標數目
%?%???nPSD????--?通帶內噪聲PSD（線性）
%?%???nPSD????--?通帶內信號PSD（線性）
%
%輸出參數：
%???Pmusic??--??寬帶非相干合成的MVDR波束功率輸出
%?Example:
%???B?=?fMUSIC（data50310[-90:0.2:90][0:.05:0.3]）


????%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
????%?Copyright（c）?2006-2026?by?Zhang?Qunfei%
????%?$Version?1.0$?$Date:?2008/7/13$??????%
????%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

C=1490;??%水中的聲速
%%?對各列進行fft變換
[NM]=size（data）;??%?N=10240;??%FFT點數超過數據點數沒有意義
%?coef=std（data）;??data=data*diag（1./coef）;??%陣元靈敏度歸一化
if?nargin<5?||?M>N??help（mfilename）??end
sects=length（x）+2;??%分段數，為確保空間自相關矩陣為滿秩，應大于陣元數
overlap=0.5;??%分段重疊率
N=fix（N/（sects-（sects-1）*overlap））;
nl=fix（fl*N/fs）;
nh=fix（fh*N/fs）;
data1=[];????%每行是特定頻率的快拍，每列是一個陣元的頻譜
for?i=0:sects-1
????iStart=fix（i*N*（1-overlap））+1;
????temp=fft（data（iStart:iStart+N-1:））;
????data1=[data1;temp（nl:nh:）];??%取出‘fl‘kHz－‘fh‘kHz的譜線
end
%?psd=mean（abs（data）2）;??psd=1./psd;
%?data=data.*（psd*ones（1M））;??%頻域預白化，對低信噪比有利，對高信噪比不利
f=nl*fs/N:fs/N:nh*fs/N;
%?%%?按信噪比加權
%?if?nargin==8
%?????[pq]=size（data）;
%?????m=length（nPSD）;
%?????n=length（sPSD）;
%?????if?m<3?||?n<3
%?????????error（‘PSD每行至少3個點。‘）;
%?????end
%?????Spsd=resample（sPSDpn）;
%?????Npsd=resample（nPSDpm）;
%?????RSN=（Spsd./Npsd）‘*ones（1q）;
%?????data=data.*RSN;
%?end????
%%?確保x為一個行向量
if?nargin<6??%不給定陣元坐標，默認為中心頻率半波長配置的均勻線列陣
????f0=（fl+fh）/2;??lamda=C/f0;
????x=（0:M-1）*lamda/2;
else
????[mn]=size（x）;
????if?m~=1
????????if?n==1??x=x‘;
????????else?error（‘陣元坐標x必須是一個行向量!?‘）
????????end
????end
????x=x-x（1）;
????if?length（x）~=M?error（‘陣列結構與數據結構不符！‘）??end
end
%%?＝＝＝頻域波束形成＝＝＝?
invLamda=f‘/C;??%波長的倒數
BF=zeros（nh-nl+1length（angles））;
m=nh-nl+1;
for?i=1:m
????temp=data1（i+（0:sects-1）*m:）;
????Rf1=temp‘*temp;
????a=exp（-j*2*pi*invLamda（i）*x‘*sin（pi/180*angles））;
????[VD]=eig（Rf1）;?D=diag（D）;
????[DI]=sort（D‘descend‘）;
????V=V（:I）;??V=V（:num+1:end）;
????BF（i:）=1./sum（abs（a‘*V）.^22）;
????BF（i:）=BF（i:）/max（BF（i:））;
end
Pmusic=sum（BF）;??%非相干疊加

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

?????文件???????2811??2008-07-14?07:23??LFMdetection\fMUSIC.m

?????文件???????1933??2013-06-19?21:22??LFMdetection\LFM_Ambiguity.m

?????文件????????299??2013-03-14?17:10??LFMdetection\LFM_MatchedFilter.m

?????文件???????1338??2013-03-20?11:28??LFMdetection\LFM_S.m

?????文件????????113??2014-09-20?20:57??LFMdetection\readme.txt

?????文件???????2665??2013-10-27?19:52??LFMdetection\Signal_cumulate.m

?????目錄??????????0??2014-09-20?20:57??LFMdetection

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

?????????????????9159????????????????????7