MIMO_DOAMIMO-雷达DOA估计中克拉美罗界MATLAB代码.zip

克拉美罗下界 Cramer-Rao Lower Bound (CRLB)可以用于计算无偏估计中能够获得的最佳估计精度，因此经常用于计算理论能达到的最佳估计精度，和评估参数估计方法的性能（是否接近CRLB下界）。 本次上传内容为：MIMO雷达DOA估计中克拉美罗界MATLAB代码

%%%%%%%Cramer-Rao bound on DOA for a single target
clear all
clc
M=10;
snr=1;
i=sqrt（-1）;
syms theta;                           %定义一个符号变量theta
%atheta = zeros（M1）;  
for m=1:M
    atheta（m1）=exp（-i*（m-1）*pi*sin（theta））;
end
atheta;  
dtheta=diff（atheta）;                  %求导引矢量关于角度theta的一阶微分
theta1 = 0;
theta2 = （12*pi）/180;
theta3 = linspace（theta1theta2）;
for k=1:100
    theta=theta3（k）;
    athetanum=eval（atheta）;
    dthetanum=eval（dtheta）;                    
    b = ones（M）;
    c = eye（M）;
    CRB1（k）=1/（2*snr*（M*dthetanum‘*b*dthetanum+athetanum‘*b*athetanum*（（norm（dthetanum））.^2）...
                       -M*（（（abs（athetanum‘*b*dthetanum））.^2）/（athetanum‘*b*athetanum））））;  %CRB on DOA for coherent signals 
    CRB2（k）=1/（2*snr*（M*dthetanum‘*c*dthetanum+athetanum‘*c*athetanum*（（norm（dthetanum））.^2）...
                       -M*（（（det（athetanum‘*c*dthetanum））.^2）/（athetanum‘*c*athetanum））））;  %CRB on DOA for orthogonal signals
    CRBthe1（1k）=CRB1（k）;
    CRBthe2（1k）=abs（CRB2（k））; 
    
end
theta3 = theta3*（180/pi）;
CRBthe11=CRBthe1*（180/pi）;
CRBthe22=CRBthe2*（180/pi）;
figure
semilogy（theta3 CRBthe11‘b--‘theta3CRBthe22‘r‘）
grid on
xlabel（‘θ[deg]‘）ylabel（‘CRB[deg]‘）
axis（[0 12 0 100]）
legend（‘β=1‘‘β=0‘2）
title（‘CRB on DOA for M=10L=1SNR=0dB.for orthogonal and coherent signals‘）

     

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----
     文件        1391  2020-07-07 08:55  MIMO_DOA.m