在网上搜到一段利用稀疏表示来求解DOA的matlab代码。其中有一处细节没看明白，按照代码逻辑，猜测是设定了一个指标，求解后相应参数若高于指标，则说明求解效果不佳，需要进一步细化。但关于该指标，没看明白，也没有找到相关文献。求指点。
相关第一处：
“
array_num = 20; %阵列数目
free_deg = 2*snap*(array_num-target_num); %自由度
PFA = 0.001;
h0 = free_deg / (gammaincinv(1-PFA,free_deg/2)*2); %这是什么指标？用到了不完全伽马函数求逆
”
用自由度推出的h0是什么含义？
相关第二处：
“
    Rone = 1/snap*y*y'; %测量值的协方差矩阵
    AFin = exp(-J*pi*(0:array_num-1)'*sin([thetaCoa(1) thetaCoa(2)]*pi/180)); %代入计算角度后的阵列流型
    PM = eye(array_num) - AFin *inv(AFin'*AFin)*AFin';%Afin不是方阵，没法求逆，大概是差值吧，什么含义？
    alpha1 = (array_num-target_num)*10^(-SNR/10)/real(trace(PM*Rone)); %这里看不懂了

    if alpha1>h0
        est_theta(i,=thetaCoa; %第i次的估计值
    else
”
PM，alpha1的含义分别是什么？和h0的关系又是什么？

谢谢了！ 返回小木虫查看更多