话不多说,直接贴代码和效果图,会用之人自然会觉得有用。
效果图:
代码:
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工具箱
下一篇:看到这张图片后,感觉每天更应该认真科研了