Code-for-arXiv-1507.01639/Raman_plot3D.m at master · qftphys/Code-for-arXiv-1507.01639 · GitHub

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
%Here I make 3D sets of 2D plots to help visualize data

%sets defaults for the plots to look good
width = 5.1;     % Width in inches
height = 3;    % Height in inches
alw = 1;       % AxesLineWidth
fsz = 14;      % Fontsize
fna = 'Helvetica'; %Fontname
lw = 1.5;      % LineWidth
msz = 8;       % MarkerSize
interp = 'tex';


%Get Data from fig files.
fig10 = load('3D_Raman_10^7_Jx_100_h_0_2.fig','-mat');
fig14 = load('3D_Raman_10^7_Jx_143_h_0_2.fig','-mat');
fig03 = load('3D_Raman_10^7_Jx_30_h_0_2.fig','-mat');

d1=fig10.hgS_070000.children.children;
[d1,d2] = d1.properties;
Ev10  = d1.XData;
R10 = d1.YData;
R102 = d2.YData;

d1=fig14.hgS_070000.children.children;
[d1,d2] = d1.properties;
Ev14  = d1.XData;
R14 = d1.YData;
R142 = d2.YData;

d1=fig03.hgS_070000.children.children;
[d1,d2] = d1.properties;
Ev0  = d1.XData;
R03 = d1.YData;
R032 = d2.YData;

%hgS_070000
%hgM_070000

%set up a vector of the same length with constant value
xv = Ev10*0+1;

%Plot the 2D plots in a 3D environment
hh = figure; hold on

plot3(1.43*xv,Ev14,R14,xv,Ev10,R10,xv*.3,Ev03,R03)
plot3(1.43*xv,Ev14,R142,xv,Ev10,R102,xv*.3,Ev03,R032,'--')

hold off

xs = .3;
x2= 1;
xb = 1.43;
yb = 17;
zb = max(max([R10,R14,R03,R102,R142,R032]));

axis([xs xb 0 yb 0 zb])
daspect([1.3 4 .8*zb])
az=80;el=34;
view(az,el)


c = .2;
w =1.2;

line([xs,xs],[0,0],[0,zb],'LineWidth',w,'Color',[c c c],'LineStyle','-')
line([xs,xs],[0,yb],[0,0],'LineWidth',w,'Color',[c c c],'LineStyle','-')

line([x2,x2],[0,0],[0,zb],'LineWidth',w,'Color',[c c c],'LineStyle','-')
line([x2,x2],[0,yb],[0,0],'LineWidth',w,'Color',[c c c],'LineStyle','-')

% line([x3,x3],[0,0],[0,zb],'LineWidth',w,'Color',[c c c],'LineStyle','-')
% line([x3,x3],[0,yb],[0,0],'LineWidth',w,'Color',[c c c],'LineStyle','-')

line([xb,xb],[0,0],[0,zb],'LineWidth',w,'Color',[c c c],'LineStyle','-')
line([xb,xb],[0,yb],[0,0],'LineWidth',w,'Color',[c c c],'LineStyle','-')

grid on
set(gca,'XGrid','off');
set(gca,'YGrid','off');
%set(gca,'ZGrid','off');

set(gca,'XTick',[.3,1,1.43])
%set(gca,'YTick',[.1,.2])
%set(gca,'ZTick',[.1,.2])


title('DOS')
ylabel('\omega/J_z');
zlabel('Raman Spectra');
legend({'J_x = 1.43 J_z','J_x = 1.0 J_z','J_x=0.3 J_z'}, 'Location', 'NorthWest');
hold off;
filename = ['3D_Ramans_h_0'];
saveas(hh,filename)
print(hh, '-dpng', filename);
print(hh, '-depsc', filename);