clc;clear;close all
dataTable = readtable('perf erosion .xlsx','Sheet','perf4','Range','A2:B901');
data_x1=dataTable{:,1};
data_y1=dataTable{:,2};
plot(data_x1/60,data_y1,'k-','LineWidth',2,'MarkerSize',10);
hold on
dataTable = readtable('perf erosion .xlsx','Sheet','perf4','Range','C2:D901');
data_x2=dataTable{:,1};
data_y2=dataTable{:,2};
plot(data_x2/60,data_y2,'r-','LineWidth',2,'MarkerSize',10);
hold on
dataTable = readtable('perf no erosion.xlsx', 'Sheet', 'perf4', 'Range', 'A2:B908');
data_x3=dataTable{:,1};
data_y3=dataTable{:,2};
plot(data_x3/60,data_y3,'k--','LineWidth',2,'MarkerSize',10);
hold on
dataTable = readtable('perf no erosion.xlsx', 'Sheet', 'perf4', 'Range', 'C2:D908');
data_x4=dataTable{:,1};
data_y4=dataTable{:,2};
plot(data_x4/60,data_y4,'r--','LineWidth',2,'MarkerSize',10);
hold on
xlabel('\fontname{宋体}时间 \fontname{Times New Roman}/ min', 'FontSize',22),
ylabel('\fontname{Times New Roman}perf fracture width / mm', 'FontSize',22)
axis([0 100 0 10])
legend('\fontname{宋体}射孔簇1和4 ;考虑冲蚀', '\fontname{宋体}射孔簇2和3 ;考虑冲蚀', 'Location', 'northwest');
% 创建虚拟坐标轴并添加第二个图例
ah = axes('position',get(gca,'position'),'visible','off');
legend(ah, '\fontname{宋体}射孔簇1和4 ;不考虑冲蚀', '\fontname{宋体}射孔簇2和3 ;不考虑冲蚀', 'northeast');