clear all;
clc;
close all;
%% ROC曲线
% y = -10:0.1:10;
a1 = 2.2;
a2 = 2;
% f11 = 1/sqrt(2*pi)*exp(-(y-a1).^2/2);
% f10 = 1/sqrt(2*pi)*exp(-(y).^2/2);
% figure;
% plot(y,f11);
% hold on;plot(y,f10,'r');
%% 第一个传感器的ROC曲线
Pf1 = 0:0.01:0.99;
N = length(Pf1);
for i = 1:N
P = 1-Pf1(i);
x = norminv(P,0,1);
xx = x:0.001:10;
yy = 1/sqrt(2*pi)*exp(-(xx-a1).^2/2);
Pd1(i) = trapz(xx,yy);
end
figure;
plot(Pf1,Pd1);
%% 第二个传感器的ROC曲线
Pf2 = 0:0.01:0.99;
N = length(Pf2);
for i = 1:N
P = 1-Pf2(i);
x = norminv(P,0,1);
xx = x:0.001:10;
yy = 1/sqrt(2*pi)*exp(-(xx-a2).^2/2);
Pd2(i) = trapz(xx,yy);
end
hold on;plot(Pf2,Pd2,'r');
%% 与规则
Pd3 = Pd1.*Pd2;
Pf3 = Pf1.*Pf2;
hold on;plot(Pf3,Pd3,'g');
%% 或规则
Pd4 = 1-(1-Pd1).*(1-Pd2);
Pf4 = 1-(1-Pf1).*(1-Pf2);
hold on;plot(Pf4,Pd4,'y');
legend('传感器1','传感器2','与融合规则','或融合规则');
%% 表决融合规则
a = [2 2.3 2.1 2.05 2.11 2.21 3];%2.14 2.26 2.16
for i = 1:length(a)
Pf(i,:) = 0:0.01:0.99;
N = length(Pf(i,:));
for j = 1:N
P = 1-Pf(i,j);
x = norminv(P,0,1);
xx = x:0.001:10;
yy = 1/sqrt(2*pi)*exp(-(xx-a(i)).^2/2);
Pd(i,j) = trapz(xx,yy);
end
end
NNN = length(a);
for k = 1:NNN;
pp = 1:length(a);
PD(k,:) = zeros(1,100);
PF(k,:) = zeros(1,100);
for i = k:length(a)
pallie = combntns(pp,i);
NN = size(pallie);
PDD(i-k+1,:) = zeros(1,100);
PFF(i-k+1,:) = zeros(1,100);
for j = 1:NN(1)
Pdd = ones(1,100);
Pff = ones(1,100);
for oo = 1:NNN
[a00,a000] = find(pallie(j,:)==oo);
a0000 = size(a000);
u(oo) = a0000(2);
end
for o = 1:NNN
if(u(o)==0)
Pdd = Pdd.*(1-Pd(o,:));
Pff = Pff.*(1-Pf(o,:));
else
Pdd = Pdd.*Pd(o,:);
Pff = Pff.*Pf(o,:);
end
end
PDD(i-k+1,:) = PDD(i-k+1,:) + Pdd;
PFF(i-k+1,:) = PFF(i-k+1,:) + Pff;
end
PD(k,:) = PD(k,:) + PDD(i-k+1,:);
PF(k,:) = PF(k,:) + PFF(i-k+1,:);
end
end
figure;
plot(PF(1,:),PD(1,:));
hold on;plot(PF(2,:),PD(2,:),'r');
hold on;plot(PF(3,:),PD(3,:),'g');
hold on;plot(PF(4,:),PD(4,:),'y');
hold on;plot(PF(5,:),PD(5,:),'r--');
hold on;plot(PF(6,:),PD(6,:),'g--');
hold on;plot(PF(7,:),PD(7,:),'y--');
% hold on;plot(PF(8,:),PD(8,:),'r.');
% hold on;plot(PF(9,:),PD(9,:),'g.');
% hold on;plot(PF(10,:),PD(10,:),'y.');
legend('1','2','3','4','5','6','7');%,'8','9','10'
% %% 与规则
% Pd3 = Pd(1,:).*Pd(2,:).*Pd(3,:).*Pd(4,:).*Pd(5,:).*Pd(6,:).*Pd(7,:);
% Pf3 = Pf(1,:).*Pf(2,:).*Pf(3,:).*Pf(4,:).*Pf(5,:).*Pf(6,:).*Pf(7,:);
% hold on;plot(Pf3,Pd3,'g');
%
