Matlab 恒虚警率检测CFAR

function [res, th_lr, th_ud]= cfar_filter(data,varargin)
% Function: 恒虚警率检测
%
% Input:
%   data: 数据
% Output:
%   res: 检测结果
%   th_lr:左右方向计算门限
%   th_ud:上下方向计算门限
%
% Example:
%     res = cfar(data);
%
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

guard_r_lr = find_Para(varargin,'guard_r_lr',2); % 左右方向保护单元半径
training_r_lr = find_Para(varargin,'training_r_lr',10); % 左右方向训练单元半径
guard_r_ud = find_Para(varargin,'guard_r_ud',1); % 上下方向保护单元半径
training_r_ud = find_Para(varargin,'training_r_ud',2); % 上下方向训练单元半径
Pfa_lr = find_Para(varargin,'Pfa_lr',1e-10); % 左右方向虚警率
Pfa_ud = find_Para(varargin,'Pfa_ud',1e-4); % 上下方向虚警率
dim = find_Para(varargin,'dim',3); % 维度 1一维(左右方向) 2一维十字 3宽十字 4窗口
type = find_Para(varargin,'type',1); % 类型 1均值 2最小值 3最大值 4排序

row = size(data,1);
col = size(data,2);
res = false(size(data));
th_lr = zeros(size(da