- With yyaxis, I don't think there is a way to set which axes' lines show up in front of the other. However, a workaround is to create two separate axes and coordinate their positions and x-limits using linkprop/linkaxes. See below. Since the right axes is created first, its lines show up under the left axes' lines.
- Store the line handles when the lines are plotted, and specify which handles to include in the legend in the legend call. See below.
- I don't know why MATLAB freezes with exportgraphics. If you use the workaround in #1, you'll have two axes, and to get them both in one exported image, you'll have to export the figure. See below. Maybe that will avoid the freezing. Try it and see.
load('data1.mat')
load('data2.mat')
load('data3.mat')
clf
ax_right = axes();
hP = plot(D.Time, D.Quantumumolm2s,Color=[0.8 0.8 0.8], LineWidth=1);
ylabel('P')
ax_right.YAxisLocation = 'right';
ax_right.XTick = [];
ax_left = axes('Color','none');
hold on
hG = line(Gk_tt.Time(1:15),Gk_tt{1:15,1},"LineWidth",1,"Color",[0, 0.5, 0],"Marker",".","MarkerSize",14,LineStyle="-");
hR = line(Rk_tt.Time(1:15),Rk_tt{1:15,1},"LineWidth",1,"Color",[0.6350, 0.0780, 0.1840],"Marker",".","MarkerSize",14,LineStyle="-");
hN = stem(NK.Time(1:15),Nk_cal1(1:15,1),"LineWidth",1,"Color",[0.9290 0.6940 0.1250],"Marker",".","MarkerSize",14,LineStyle="-");
line(GM.Time,GM{:,1},"LineWidth",1,"Color",[0, 0.5, 0],"Marker",".","MarkerSize",16,LineStyle="-")
line(RM.Time,RM{:,1},"LineWidth",1,"Color",[0.6350, 0.0780, 0.1840],"Marker",".","MarkerSize",16,LineStyle="-")
stem(RM.Time,NM(:,1),"LineWidth",1,"Color",[0.9290 0.6940 0.1250],"Marker",".","MarkerSize",16,LineStyle="-")
xlabel('Time')
ylabel('GK,RK and NK')
ax_left.YAxis.Color = 'k';
ax_right.YAxis.Color = 'k';
legend([hP,hG,hR,hN],{'P','GK','RK','NK'},Location='northwest')
datetick
fontsize([ax_left ax_right],16,"points")
linkprop([ax_left ax_right],'Position');
linkaxes([ax_left ax_right],'x')
xlim([datetime('April 20, 2019'), datetime('June 5, 2019')])
yl = ylim(ax_left);
tMark = datetime(2019,5,5); % Time to mark by a vertical line
text(tMark + days(0.5), mean(yl)*2 , {'Start', 'layer'}, 'VerticalAlignment', 'bottom', 'HorizontalAlignment', 'left', 'FontSize', 14, 'Color', 'blue');
plot([tMark tMark], yl,"Color","b",LineWidth=1,LineStyle="-");
set(gcf, 'WindowState', 'maximized')
exportgraphics(gcf(),"m.png",'Resolution',600)
