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
|
function [] = energy_vs_stability( possible_energy, stability, index, stability_max )
%Plots energy vs. stability
% Using data sent from mode_match, this function plots energy vs.
% stability of each solution for eyed comparisons.
n = size(stability,2);
%Color code data points
if n == 3
c = [1, 0, 0; 0, 1, 0; 0, 0, 1];
elseif n == 2
c = [1, 0, 0; 0, 1, 0];
elseif n == 1
c = [1, 0, 0];
else
c = 1:n;
end
%Reorganize stability to coincide with energies
for i=1:n
fixed_stability(i,1) = stability(i);
end
for i = 1:n
energy(i,1) = possible_energy(index(i));
end
%Graph scatter plot
figure(n+1)
textCell = arrayfun(@(x,y) sprintf('(%3.2f, %3.2f)',x,y),energy,fixed_stability,'un',0);
scatter(energy,fixed_stability,25,c,'filled')
xlabel('Energy')
ylabel('Stability')
ylim([0,stability_max])
dx =0.0000000005;
dy =0.0000000005;
%Label data points relative to solution number
for i = 1:n
n_solution = num2str(i);
text(energy(i)+dx, fixed_stability(i)+dy, n_solution, 'FontSize',10);
end
end
|
