|
|
【答案】应助回帖
给你一个例子。可以自己修改
clc
clear
quakedate=datenum(1976,7,28); %起止时间
enddate=datenum(2010,7,27); %起止时间
data=load('%数据文件%');
loc=[118.1,39.6]; %作图范围
dis=0.5;
range=[loc(1)-dis,loc(1)+dis,loc(2)-dis,loc(2)+dis];
%% mag<4
mag=4;
minmag=3;
choose1=data(:,2)>=range(1);
choose2=data(:,2)<=range(2);
choose3=data(:,3)>=range(3);
choose4=data(:,3)<=range(4);
choose=choose1&choose2&choose3&choose4;
tempLoc=data(choose, ;
tempMag=tempLoc(tempLoc(:,1)<=mag, ;
tempMag=tempMag(tempMag(:,1)>=minmag, ;
year=tempMag(:,4);
mon=tempMag(:,5);
day=tempMag(:,6);
plotdate=datenum(year,mon,day);
Mag=tempMag(:,1);
h=stem(plotdate,Mag);
datetick('x','yyyy')
axis([quakedate enddate 2 7])
这个画的是震级3~4之间的。范围可以自己改 |
|