-
Notifications
You must be signed in to change notification settings - Fork 0
/
run_AD.m
80 lines (71 loc) · 1.84 KB
/
run_AD.m
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
clear all
close all
I = 1;
eps_o = 0.9;
mu = 10;
theta_o = 0.1;
lambda = 10;
gamma = linspace(0, 5, 100);
gamma = gamma(2:end)';
C = linspace(0, I, 100);
C = C(2:end)';
W = zeros(length(gamma), length(C));
Y = zeros(size(W));
R = zeros(size(W));
Ud = zeros(size(W));
Ua = zeros(size(W));
Pay = zeros(size(W));
for i = 1:length(gamma)
i
for j = 1:length(C)
[W(i, j), Y(i, j), R(i, j), Ud(i, j), Ua(i, j), Pay(i ,j)] = AD_exp(gamma(i), C(j), I, eps_o, mu, theta_o, lambda);
end
end
fontsize = 20;
figure;
imagesc(C, gamma, R);
colormap(magma(50));
colorbar;
caxis([0,1]);
set(gca,'YDir','normal');
caxis([0 1]);
xlabel('$C$', 'interpreter', 'latex');
ylabel('$\gamma$', 'interpreter', 'latex');
title('Ransom ($R^*$)', 'interpreter', 'latex');
set(gca, 'FontSize', fontsize);
set(gcf,'Position', [100 100 600 500]);
figure;
imagesc(C, gamma, W);
colormap(magma());
colorbar;
%caxis([0,1])
set(gca,'YDir','normal');
caxis([0 0.5]);
xlabel('$C$', 'interpreter', 'latex');
ylabel('$\gamma$', 'interpreter', 'latex');
title('Deterrence Effort ($W^*$)', 'interpreter', 'latex');
set(gca, 'FontSize', fontsize);
set(gcf,'Position', [100 100 600 500]);
figure;
imagesc(C, gamma, Y);
colormap(magma());
colorbar;
%caxis([0,1])
set(gca,'YDir','normal');
caxis([0 0.5]);
xlabel('$C$', 'interpreter', 'latex');
ylabel('$\gamma$', 'interpreter', 'latex');
title('Backup Effort ($Y^*$)', 'interpreter', 'latex');
set(gca, 'FontSize', fontsize);
set(gcf,'Position', [100 100 600 500]);
figure;
imagesc(C, gamma, Pay);
colormap(magma());
colorbar;
set(gca,'YDir','normal');
caxis([0 1]);
xlabel('$C$', 'interpreter', 'latex');
ylabel('$\gamma$', 'interpreter', 'latex');
title('Ransom ($R^*$)', 'interpreter', 'latex');
set(gca, 'FontSize', fontsize);
set(gcf,'Position', [100 100 600 500]);