plot(x,polyval(t1,x),'--b',x,polyval(t2,x),'--r',x,polyval(t3,x),'--k','LineWidth',1);
axis([1960,2010,200,2500])
hold on
p2=plot(x,GTI,'-b','LineWidth',1.5);
p21=plot(x,pGTI,'-r','LineWidth',1.5);
p22=plot(x,sGTI,'-g','LineWidth',1.5);
hold on
[AX,H1,H2]=plotyy(x,GTI,x,sp,'plot','bar'); %画双轴,AX(1)左轴,AX(2)右轴,H为曲线本身
set(AX(2),'TickDir','out','YDir','reverse','Ylim',[0,0.5*10^6],'YTick',[0:0.1*10^6:0.5*10^6],'YTickLabel',[0 0.1 0.2 0.3 0.4 0.5],'YColor','b','FontSize',12); %设置右边轴为倒立
set(AX(1),'TickDir','in','Ylim',[200,3000],'YTick',[200:500:3000],'FontSize',12,'YColor','k');
set(get(AX(1),'Xlabel'),'String','Time');
set(get(AX(1),'Ylabel'),'String','GTI');
set(get(AX(2),'Ylabel'),'string','');
legend([p2,p21,p22,H2],'GTI','PGTI','SGTI','');
hold off