資源簡介
用于STAP的仿真 運用采樣矩陣求逆算法 STAP的其他算法可以再此基礎上改進得到
代碼片段和文件信息
%opt2d.m:?全自由度空時自適應處理
%--------------------------------------------------------------------------
%start?????????:?2004.11.04??AM?10:21???yunrisheng
%Latest?change?:?
%--------------------------------------------------------------------------
clear
tic
%雜波仿真參數
N?=?12;????????????????????????%?陣元個數
M?=?10;????????????????????????%?相干脈沖數
CNR?=?30;??????????????????????%?雜噪比
beta?=?1;??????????????????????%?雜波折疊系數(beta?=?2*v*T/d)
sita_a?=?-90:.9:90.;?????????%?雜波單元個數???????????????
sita?=?sita_a*pi/180;
[NN?N_bin]?=?size(sita);
%目標參數
sita_t?=?-25;??????????????????%?目標DOA
omiga_t?=?0.4;?????????????????%?目標Doppler
SNR?=?0;???????????????????????%?信噪比
%空間導向矢量和時間導向矢量
%空間頻率和Dopple頻率滿足?omiga_d?=?beta?*?omiga_s
omiga_s?=?pi*sin(sita);??????
omiga_d?=?beta*omiga_s;???????
aN?=?zeros(NN_bin);
bN?=?zeros(MN_bin);
aN?=?exp(-j*[0:N-1]‘*omiga_s)./sqrt(N);
bN?=?exp(-j*[0:M-1]‘*omiga_d)./sqrt(M);
%目標空時信號
aN_t?=?zeros(N1);
bN_t?=?zeros(M1);
aN_t?=?exp(-j*pi*[0:N-1]‘*sin(sita_t*pi/180))/sqrt(N);
bN_t?=?exp(-j*pi*[0:M-1]‘*omiga_t)/sqrt(M);
S_t?=?zeros(M*N1);
S_t?=?kron(aN_tbN_t);
%計算雜波協方差矩陣
R?=?zeros(M*NM*N);?????????????????????
S?=?zeros(M*NN_bin);???????????????????
ksai?=?10^(CNR/10)*(randn(1N_bin)+j*randn(1N_bin))/sqrt(2);???????????????%服從正態分布的隨機幅值,方差為1
for?ii?=?1:N_bin
????S(:ii)?=?kron(aN(:ii)bN(:ii));??
????R?=?R?+?ksai(ii).*(S(:ii)*S(:ii)‘);???????
end
%干擾協方差矩陣,雜噪比為30dB
R?=?R?+eye(M*N);??
- 上一篇:IEEE33節點潮流程序
- 下一篇:卡爾曼濾波matlab代碼
評論
共有 條評論