black = [17 18 12]/256 ;

green = [101 183 145] /256;

red = [215,48,39]/256;

blue = [69,117,180]/256;

fid = fopen('E:\code\matlab\论\ramp structures from turbulence to climate\HawaiiCO2DataSet.txt');

textscan(fid,'%s',149169,'delimiter','\n');%start from 1990

Data = textscan(fid,'%s%f%f%f%f%f%f%f%f%f%f%f%f%f%f%s%s');

CO2Data = Data{1,8};

CO2Data(find(CO2Data == -999.99)) = nan;

CO2Data = [CO2Data; nan];

%-----------------------------画出二氧化碳的值，包括nan的缺测值---------------

figure

width =18; %图片大小

height = 16;

set(gcf, 'unit','centimeters','Position', [0 0 width height]); %参数含义为：xmin,ymin,width,height

plot(CO2Data,'color',red)

set(gca,'xtick',1:43824:length(CO2Data),'xticklabel',{'1990','1995','2000','2005','2010','2015'})

set(gca, 'Fontname', 'Arial', 'Fontsize', 9,'LineWidth', 1,'Fontweight','normal');

xlabel('Year','Fontname', 'Arial', 'Fontsize', 10,'Fontweight','normal')

ylabel('CO_2(\mumol mol^-^1) ','Fontname', 'Arial', 'Fontsize', 10,'Fontweight','normal')

title('CO_2 Concentration','Fontname', 'Arial', 'Fontsize', 10,'Fontweight','normal')

print('CO2_Climate','-deps','-r600');%eps格式 精度为600 dpi 颜色为cmyk系统