话不多说，直接贴代码和效果图，会用之人自然会觉得有用。

效果图：

代码：

clear;

clc;

load('ALGO_GPS.neu');  % 加载数据

data = ALGO_GPS;

dyr = data(:,1);

resi_u_gps = detrend(data(:,6))*1000;  % 去除趋势项，单位：mm

e = data(:,9)*1000;  % 误差项

errorbar(dyr,resi_u_gps,e,'-o','Color',[0.2 0.8 0.8]','LineWidth',2,'MarkerEdgeColor','k',...

'MarkerFaceColor','g','MarkerSize',8);  % 带误差棒的图，('Color',[0.2 0.8 0.8]' 这个设置为颜色配比，根据自己喜好设置颜色配比，范围为(0,1))

grid on;

hl = legend('GPS','Location','NorthWest');

set(hl,'Orientation','horizon');  % 图例水平排列

set(hl,'Box','off')  % 去除图例边框

xlim([2002.0 2012.0]);

数据：

NOTE：数据文件和代码放在同一个目录下。

转载本文请联系原作者获取授权，同时请注明本文来自周锋科学网博客。

链接地址：http://blog.sciencenet.cn/blog-350278-712402.html

上一篇：大地测量部分时间转换Matlab工具箱

下一篇：看到这张图片后，感觉每天更应该认真科研了