You can not select more than 25 topics
Topics must start with a letter or number, can include dashes ('-') and can be up to 35 characters long.
25 lines
1016 B
25 lines
1016 B
%%***********************【改进FCME检测算法函数】*************************%%
|
|
%%=============================参数说明===================================%%
|
|
% data 待检测的频域信号
|
|
% Pf 虚警概率
|
|
% J 检测出来的干扰点集合
|
|
%%========================================================================%%
|
|
function J = FCME_pro(data_Freuq,Pf)
|
|
T = sqrt(-1*(4/pi)*log(Pf)); %%门限因子
|
|
th = 3.432; %%倍数关系
|
|
L = length(data_Freuq);
|
|
fi = [];
|
|
[y_sort,index] = sort(abs(data_Freuq),'ascend');
|
|
%%求每个频点的幅值
|
|
fi = abs(data_Freuq);
|
|
I = index(1:round(L/2)); %%无干扰集合
|
|
E = mean(fi(I));
|
|
Aaim = T*th*E; %%门限
|
|
J = [];
|
|
%%比较求得干扰集合
|
|
for mm = 1:L
|
|
if fi(mm) > Aaim
|
|
J = [J,mm];
|
|
end
|
|
end
|
|
end
|
|
|