gaussianpointdetection.m

强杂波背景下的的高斯点目标的检测Matlab源代码

%%
% White Ambient Noise and No Jamming or Pn is flat 
clear all;clc;
Ph = 1;
f = [-1:0.001:1];
NumF = length(f);
eplison = 1/8;
% calculate the Pn
for i = 1:NumF
    Pn(i) = exp(-abs(f(i)));
end
lmda = [0:0.001:3];
NumLm = length(lmda);
eplisonlada = zeros(1,NumLm);   %  能量值
for j = 1:NumLm
    for i = 1:NumF
        temp(i) = sqrt(Pn(i)/lmda(j))-Pn(i);
%         temp(i) = Pn(i)/sqrt(lmda(j))-Pn(i);
        if temp(i)>=0
            eplisonSf(i) = temp(i);
        else
            eplisonSf(i) = 0;
        end
    end
    f = -1:0.001:1;
    eplisonlada(j) = trapz(f,eplisonSf);
end
figure;
plot(lmda,eplisonlada); grid on;  % the energy of signal varies with the different lamda
axis([0 3 0 4]);
ylabel('Energy');xlabel('Lamda');

% find the optimal lamda
lamdaMax = exp(1);
for j =1:NumLm
    if lmda(j)<= lamdaMax
        errEnergy = abs(eplisonlada(j)-1/8);
        if errEnergy<=0.0001
            lamdaOptiaml = lmda(j);   % The optimal lamda
        end
    end
end
for i = 1:NumF
    temp1(i) = sqrt(Pn(i)/lamdaOptiaml)-Pn(i);
    if temp1(i)>=0
         eplisonSffinal(i) = temp1(i);
    else
         eplisonSffinal(i) = 0;
    end
end
figure 
plot(f,eplisonSffinal);grid on;    % the optimal Es(f) varies with the f
axis([-1 1 0 0.16]);
ylabel('Optimal Es(F)');xlabel('F');

%  calculate the value of d^2
for i = 1:NumF
    temp2ford(i) =  eplisonSffinal(i)/(eplisonSffinal(i) + Pn(i));
end
d2 = trapz(f,temp2ford)