- 积分
- 4994
- 贡献
-
- 精华
- 在线时间
- 小时
- 注册时间
- 2013-12-3
- 最后登录
- 1970-1-1
|

楼主 |
发表于 2014-5-25 11:14:41
|
显示全部楼层
clc;
clear;
close all;
data1=importdata('E:\workshop\biweight\Z-score=1.5\biweight_remove_ch3.dat');
data2=importdata('E:\workshop\biweight\Z-score=1.5\biweight_retain_ch3.dat');
[m,n]=size(data1);
[s,t]=size(data2);
lat1=data1(:,1);
bias1=data1(:,3);
lat2=data2(:,1);
bias2=data2(:,3);
%plot(lat1,bias1,'r.','MarkerSize',1);
plot(lat2,bias2,'g.','MarkerSize',1);
hold on;
%plot(lat2,bias2,'g.','MarkerSize',1);
plot(lat1,bias1,'r.','MarkerSize',1);
XLabel('Latitude','FontSize',20);
set(gca,'Xlim',[-90,90]);
YLabel('O-M','FontSize',20);
set(gca,'Ylim',[-30,30]);
set(gca,'xticklabel',{'-30','-25','-20','-15','-10','-5','0',...
'5','10','15','20','25','30'});
Title( 'ch3 '); |
|