-
Notifications
You must be signed in to change notification settings - Fork 42
/
plot_resultsV.m
143 lines (128 loc) · 5.62 KB
/
plot_resultsV.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
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
function [ ] = plot_resultsV( selectionParams, targetSa, IMs, simulatedSpectra, SaKnown, knownPer )
%% Produce a number of plots of target and selected response spectra
% Variables used here
% SaKnown : As before, it contains the response spectra of all the
% available ground motions (N*P matrix) - N ground motions,
% P periods
% SaKnownV : Same as SaKnown except for V component
% simulatedSpectra : a matrix of simulated response spectra defined
% at selectionParams.TgtPer and selectionParams.TgtPerV
% sampleBig : Same as SaKnown, but is only defined at PerTgt, the
% periods at which the target response spectrum properties
% are computed
% sampleBigV : Same as sampleBig except for V component
% sampleSmall : The response spectra of the selected ground motions,
% defined at TgtPer
% sampleSmallV : Same as sampleSmall except for V component
% meanReq : Target mean for the (log) response spectrum
% covReq : Target covariance for the (log) response spectrum
% Extract data for H and V components of GM
nPer = length(selectionParams.TgtPer);
tgtSpecMed = exp(targetSa.meanReq(1,1:nPer));
tgtSpecMedV = exp(targetSa.meanReq(1,(nPer+1):end));
tgtSpecStd = targetSa.stdevs(1,1:nPer);
tgtSpecStdV = targetSa.stdevs(1,(nPer+1):end);
simSpec = simulatedSpectra(:,1:nPer);
simSpecV = simulatedSpectra(:,(nPer+1):end);
SaKnownV = selectionParams.SaKnownV;
% Plot simulated response spectra
nSTD = 1.96;
figure
subplot(1,2,1);
loglog(selectionParams.TgtPer, tgtSpecMed, '-r', 'linewidth', 3)
hold on
loglog(selectionParams.TgtPer, tgtSpecMed .* exp(+nSTD*tgtSpecStd), '--r', 'linewidth', 3)
for i=1:selectionParams.nGM
loglog(selectionParams.TgtPer, simSpec(i,:), 'k');
end
loglog(selectionParams.TgtPer, tgtSpecMed .* exp(-nSTD*tgtSpecStd), '--r', 'linewidth', 3)
axis([min(selectionParams.TgtPer) max(selectionParams.TgtPer) 1e-2 5])
axis square
xlabel('T (s)')
ylabel('S_a (g)')
hLeg = legend('Median','2.5 and 97.5 percentile', 'Simulated spectra');
set(hLeg,'location','southwest','box','off');
title('Horizontal')
subplot(1,2,2);
loglog(selectionParams.TgtPerV, tgtSpecMedV, '-r', 'linewidth', 3)
hold on
loglog(selectionParams.TgtPerV, tgtSpecMedV .* exp(+nSTD*tgtSpecStdV), '--r', 'linewidth', 3)
for i=1:selectionParams.nGM
loglog(selectionParams.TgtPerV, simSpecV(i,:), 'k');
end
loglog(selectionParams.TgtPerV, tgtSpecMedV .* exp(-nSTD*tgtSpecStdV), '--r', 'linewidth', 3)
axis([min(selectionParams.TgtPerV) max(selectionParams.TgtPerV) 1e-2 5])
axis square
xlabel('T (s)')
title('Vertical')
% Plot selected response spectra
figure
subplot(1,2,1);
loglog(selectionParams.TgtPer, tgtSpecMed, 'b', 'linewidth', 3)
hold on
loglog(selectionParams.TgtPer, tgtSpecMed .* exp(+nSTD*tgtSpecStd), '--b', 'linewidth', 3)
loglog(knownPer,SaKnown(IMs.recID,:).*repmat(IMs.scaleFac,1,size(SaKnown,2)),'k');
loglog(selectionParams.TgtPer, tgtSpecMed .* exp(-nSTD*tgtSpecStd), '--b', 'linewidth', 3)
axis([min(selectionParams.TgtPer) max(selectionParams.TgtPer) 1e-2 5])
axis square
xlabel('T (s)');
ylabel('S_a (g)');
hLeg = legend('Median','2.5 and 97.5 percentile','Selected ground motions');
set(hLeg,'location','southwest','box','off');
title('Horizontal')
subplot(1,2,2);
loglog(selectionParams.TgtPerV, tgtSpecMedV, 'b', 'linewidth', 3)
hold on
loglog(selectionParams.TgtPerV, tgtSpecMedV .* exp(+nSTD*tgtSpecStdV), '--b', 'linewidth', 3)
loglog(knownPer,SaKnownV(IMs.recID,:).*repmat(IMs.scaleFacV,1,size(SaKnownV,2)),'k');
loglog(selectionParams.TgtPerV, tgtSpecMedV .* exp(-nSTD*tgtSpecStdV), '--b', 'linewidth', 3)
axis([min(selectionParams.TgtPerV) max(selectionParams.TgtPerV) 1e-2 5])
axis square
xlabel('T (s)');
title('Vertical')
% Target, initial, and finally selected medians
figure
subplot(1,2,1);
loglog(selectionParams.TgtPer, tgtSpecMed,'k','linewidth',1)
hold on
loglog(knownPer, exp(IMs.stageOneMeans),'r*--', 'linewidth',1)
loglog(knownPer,exp(mean(log(SaKnown(IMs.recID,:).*repmat(IMs.scaleFac,1,size(SaKnown,2))))),'b--','linewidth',1)
axis([min(selectionParams.TgtPer) max(selectionParams.TgtPer) 1e-2 5])
axis square
xlabel('T (s)');
ylabel('Median S_a (g)');
hLeg = legend('Target', 'Stage 1 selection', 'Final selection');
set(hLeg,'location','southwest','box','off');
title('Horizontal')
subplot(1,2,2);
loglog(selectionParams.TgtPerV, tgtSpecMedV,'k','linewidth',1)
hold on
loglog(knownPer, exp(IMs.stageOneMeansV),'r*--', 'linewidth',1)
loglog(knownPer,exp(mean(log(SaKnownV(IMs.recID,:).*repmat(IMs.scaleFacV,1,size(SaKnownV,2))))),'b--','linewidth',1)
axis([min(selectionParams.TgtPerV) max(selectionParams.TgtPerV) 1e-2 5])
axis square
xlabel('T (s)');
title('Vertical')
% Target, initial, and finally selected standard deviations
figure
subplot(1,2,1);
semilogx(selectionParams.TgtPer,tgtSpecStd,'k','linewidth',1)
hold on
semilogx(knownPer, IMs.stageOneStdevs,'r*--','linewidth',1)
semilogx(knownPer,std(log(SaKnown(IMs.recID,:).*repmat(IMs.scaleFac,1,size(SaKnown,2)))),'b--','linewidth',1)
axis([min(selectionParams.TgtPer) max(selectionParams.TgtPer) 0 1])
axis square
xlabel('T (s)');
ylabel('Standard deviation of lnS_a');
hLeg = legend('Target', 'Stage 1 selection','Final selection');
set(hLeg,'location','northeast','box','off');
title('Horizontal')
subplot(1,2,2);
semilogx(selectionParams.TgtPerV,tgtSpecStdV,'k','linewidth',1)
hold on
semilogx(knownPer, IMs.stageOneStdevsV,'r*--','linewidth',1)
semilogx(knownPer,std(log(SaKnownV(IMs.recID,:).*repmat(IMs.scaleFacV,1,size(SaKnownV,2)))),'b--','linewidth',1)
axis([min(selectionParams.TgtPerV) max(selectionParams.TgtPerV) 0 1])
axis square
xlabel('T (s)');
title('Vertical')