-
Notifications
You must be signed in to change notification settings - Fork 0
/
Make_figure6.m
56 lines (45 loc) · 1.4 KB
/
Make_figure6.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
%Make figure 6 using the data file figure6data, which is generated using
%the code make_data_figure6
load('Data/figure6data.mat')
figure(1); clf; hold on
%make colormaps
colmap_trapped = zeros(7,3);
colmap_trapped(:,1) = linspace(0.5,1,7);
colmap_escape = zeros(4,3);
colmap_escape(:,end) = flip(linspace(0.3,1,4));
colmap = [colmap_trapped; colmap_escape];
colmap_full = parula(18);
colmap = [colmap_full(1:7,:); colmap_full(end-3:end, :)];
for i = 1:11
t = cell2mat(data(i,2));
xu = cell2mat(data(i,3));
lambda = cell2mat(data(i,4));
ispinned =cell2mat(data(i,5));
subplot(1,2,1); hold on
if ispinned(end,1) == 1
plot(t, xu, 'color', colmap(i,:), 'linewidth', 3)
else
plot(t, xu, 'color', colmap(i,:), 'linewidth', 3)
end
subplot(1,2,2); hold on
if ispinned(end,1) == 1
plot(t, lambda,'color', colmap(i,:), 'linewidth', 3)
else
plot(t, lambda,'color', colmap(i,:), 'linewidth', 3)
end
end
subplot(1,2,1);
box on
set(gca, 'XScale', 'log')
xlim([1e-4, 10])
xlabel('$t$', 'interpreter', 'latex');
ylabel('$x_2(t)$', 'interpreter', 'latex');
subplot(1,2,2);
box on
set(gca, 'XScale', 'log')
xlim([1e-5, 10])
ylim([0,hysparmax*1.2])%this cuts of the small error resulting from one sided approx and finite N
xlabel('$t$', 'interpreter', 'latex');
ylabel('$\lambda(t)$', 'interpreter', 'latex');
fig = gcf;
fig.Position(3:4) = [1200, 420];