122 lines
3.3 KiB
Matlab
122 lines
3.3 KiB
Matlab
clear;
|
|
close all;
|
|
clc;
|
|
|
|
load test_Pd_Pfa_vamp_cal_stat.mat;
|
|
|
|
Fontsize = 18;
|
|
plot_width = 800;
|
|
plot_height = 600;
|
|
Linewidth = 2;
|
|
Markersize = 8;
|
|
|
|
|
|
%% plot
|
|
figure(1);
|
|
loglog(P_fa, P_fa_CROD(1,:), '-o', ...
|
|
'Linewidth', Linewidth, ...
|
|
'MarkerSize', Markersize);
|
|
hold on;
|
|
grid on;
|
|
loglog(P_fa, P_fa_CROD(2,:), '-+', ...
|
|
'Linewidth', Linewidth, ...
|
|
'MarkerSize', Markersize);
|
|
loglog(P_fa, P_fa_CROD(3,:), '-d', ...
|
|
'Linewidth', Linewidth, ...
|
|
'MarkerSize', Markersize);
|
|
loglog(P_fa, P_fa_CROD(4,:), '-s', ...
|
|
'Linewidth', Linewidth, ...
|
|
'MarkerSize', Markersize);
|
|
loglog(P_fa, P_fa_CROD(5,:), '-^', ...
|
|
'Linewidth', Linewidth, ...
|
|
'MarkerSize', Markersize);
|
|
legend('SNR = 0', 'SNR = 2', 'SNR = 4', 'SNR = 6', 'SNR = 8');
|
|
xlabel('P_{fa} Set');
|
|
ylabel('Actual P_{fa}');
|
|
ylim([8e-5,1]);
|
|
set(gca, 'FontSize', Fontsize);
|
|
title('cVAMPro P_{fa}')
|
|
set(gcf, 'position', [200, 300, plot_width, plot_height]);
|
|
set(gca,'fontsize',20,'fontname','Times');
|
|
|
|
|
|
figure(3);
|
|
semilogx(P_fa_CROD(1,:), P_d_CROD(1,:), '-o', ...
|
|
'Linewidth', Linewidth, ...
|
|
'MarkerSize', Markersize);
|
|
hold on;
|
|
grid on;
|
|
semilogx(P_fa_CROD(2,:), P_d_CROD(2,:), '-+', ...
|
|
'Linewidth', Linewidth, ...
|
|
'MarkerSize', Markersize);
|
|
semilogx(P_fa_CROD(3,:), P_d_CROD(3,:), '-d', ...
|
|
'Linewidth', Linewidth, ...
|
|
'MarkerSize', Markersize);
|
|
semilogx(P_fa_(3,:), P_d_CROD(4,:), '-s', ...
|
|
'Linewidth', Linewidth, ...
|
|
'MarkerSize', Markersize);
|
|
semilogx(P_fa_CROD(3,:), P_d_CROD(5,:), '-^', ...
|
|
'Linewidth', Linewidth, ...
|
|
'MarkerSize', Markersize);
|
|
legend('SNR = 0', 'SNR = 2', 'SNR = 4', 'SNR = 6', 'SNR = 8');
|
|
xlabel('P_{fa}');
|
|
ylabel('P_{d}');
|
|
xlim([8e-5,1]);
|
|
set(gca, 'FontSize', Fontsize);
|
|
title('cVAMPro ROC')
|
|
set(gcf, 'position', [200, 300, plot_width, plot_height]);
|
|
set(gca,'fontsize',20,'fontname','Times');
|
|
|
|
|
|
figure(5);
|
|
semilogy(SNR, P_fa_CROD(:,1), '-o', ...
|
|
'Linewidth', Linewidth, ...
|
|
'MarkerSize', Markersize);
|
|
hold on;
|
|
grid on;
|
|
semilogy(SNR, P_fa_CROD(:,3), '-+', ...
|
|
'Linewidth', Linewidth, ...
|
|
'MarkerSize', Markersize);
|
|
semilogy(SNR, P_fa_CROD(:,5), '-d', ...
|
|
'Linewidth', Linewidth, ...
|
|
'MarkerSize', Markersize);
|
|
semilogy(SNR, P_fa_CROD(:,7), '-s', ...
|
|
'Linewidth', Linewidth, ...
|
|
'MarkerSize', Markersize);
|
|
semilogy(SNR, P_fa_CROD(:,9), '-^', ...
|
|
'Linewidth', Linewidth, ...
|
|
'MarkerSize', Markersize);
|
|
legend('Pfa = 1e-4', 'Pfa = 1e-3', 'Pfa = 1e-2', 'Pfa = 1e-1', 'Pfa = 1');
|
|
xlabel('SNR');
|
|
ylabel('P_{fa}');
|
|
set(gca, 'FontSize', Fontsize);
|
|
title('cVAMPro')
|
|
set(gcf, 'position', [200, 300, plot_width, plot_height]);
|
|
set(gca,'fontsize',20,'fontname','Times');
|
|
|
|
|
|
figure(6);
|
|
plot(SNR, P_d_CROD(:,1), '-o', ...
|
|
'Linewidth', Linewidth, ...
|
|
'MarkerSize', Markersize);
|
|
hold on;
|
|
grid on;
|
|
plot(SNR, P_d_CROD(:,3), '-+', ...
|
|
'Linewidth', Linewidth, ...
|
|
'MarkerSize', Markersize);
|
|
plot(SNR, P_d_CROD(:,5), '-d', ...
|
|
'Linewidth', Linewidth, ...
|
|
'MarkerSize', Markersize);
|
|
plot(SNR, P_d_CROD(:,7), '-s', ...
|
|
'Linewidth', Linewidth, ...
|
|
'MarkerSize', Markersize);
|
|
plot(SNR, P_d_CROD(:,9), '-^', ...
|
|
'Linewidth', Linewidth, ...
|
|
'MarkerSize', Markersize);
|
|
legend('Pfa = 1e-4', 'Pfa = 1e-3', 'Pfa = 1e-2', 'Pfa = 1e-1', 'Pfa = 1');
|
|
xlabel('SNR');
|
|
ylabel('P_{d}');
|
|
set(gca, 'FontSize', Fontsize);
|
|
title('cVAMPro')
|
|
set(gcf, 'position', [200, 300, plot_width, plot_height]);
|
|
set(gca,'fontsize',20,'fontname','Times'); |