91av视频/亚洲h视频/操亚洲美女/外国一级黄色毛片 - 国产三级三级三级三级

  • 大小: 5KB
    文件類型: .rar
    金幣: 2
    下載: 0 次
    發布日期: 2024-01-31
  • 語言: 其他
  • 標簽: LFM信號??

資源簡介

LFM信號及回波的生成,時域及頻譜分析,模糊度函數分析,時頻分析,匹配濾波器能量檢測、循環相關檢測、及fmusic參數估計

資源截圖

代碼片段和文件信息

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


評論

共有 條評論