forked from dbgomez94/high-throughput-satellites
-
Notifications
You must be signed in to change notification settings - Fork 0
/
HTS_Joint_Distribution.m
51 lines (45 loc) · 1.55 KB
/
HTS_Joint_Distribution.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
%% Just HTS
load('HTS_Value_Analysis_Results_Nominal_V2.mat');
[mu_NPV,sig_NPV] = normfit(NPV);
[mu_RIOC,sig_RIOC] = normfit(Roic);
Roic=R./(C_aioc+C_ioc);
Roic_max=max(Roic);
Roic_min=min(Roic);
seg_x=(Roic_max-Roic_min)/100;
NPV_max=max(NPV);
NPV_min=min(NPV);
seg_y=(NPV_max-NPV_min)/100;
x_axis=Roic_min:seg_x:Roic_max;
y_axis=NPV_min:seg_y:NPV_max;
figure(1)
histogram2(Roic,NPV,x_axis,y_axis,'Normalization','pdf','FaceColor','flat');
h = colorbar;
colormap bone
set(get(h,'title'),'string','PMF');
view(2)
xlabel('Discounted Return on Invested Capital (ROIC)','fontname','times new roman')
ylabel('Net Present Value ($ million)','fontname','times new roman')
zlabel('Probability Mass Function (PMF)','fontname','times new roman')
%% Comparison of HTS and WB
figure(2)
yyaxis left
histogram2(Roic,NPV,x_axis,y_axis,'Normalization','pdf','FaceColor','flat');
h = colorbar;
colormap jet
set(get(h,'title'),'string','PMF');
view(2)
xlabel('Discounted Return on Invested Capital','fontname','times new roman')
ylabel('Net Present Value ($ million)','fontname','times new roman')
zlabel('Probability Mass Function (PMF)','fontname','times new roman')
hold on
clear;
load('WB_Value_Analysis_Results.mat');
yyaxis right
histogram2(Roic,NPV,x_axis,y_axis,'Normalization','pdf','FaceColor','flat');
h = colorbar;
colormap jet
set(get(h,'title'),'string','PMF');
view(2)
xlabel('Discounted Return on Invested Capital','fontname','times new roman')
ylabel('Net Present Value ($ million)','fontname','times new roman')
zlabel('Probability Mass Function (PMF)','fontname','times new roman')