| 查看: 2855 | 回复: 11 | |||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | |||
[求助]
M-T图求助
|
|||
各位大侠，小弟求助个技术问题，怎样绘制M-t图，用matalab或GMT皆可 时间较好表示，大小怎么来画呢？ |
回复此楼» 猜你喜欢
地质学复试调剂
已经有1人回复
-
地球物理学和空间物理学论文润色/翻译怎么收费?
已经有72人回复
sissiautumn
铜虫 (初入文坛)
- 应助: 0 (幼儿园)
- 金币: 547.3
- 散金: 300
- 红花: 1
- 帖子: 49
- 在线: 20.2小时
- 虫号: 1963608
- 注册: 2012-08-30
- 专业: 地震学
【答案】应助回帖
|
给你一个例子。可以自己修改 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之间的。范围可以自己改 |
5楼2012-10-20 15:23:39
2楼2012-09-16 20:52:27
zsyanlei
铜虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 721.1
- 散金: 300
- 帖子: 174
- 在线: 27.9小时
- 虫号: 1005242
- 注册: 2010-04-25
- 性别: GG
- 专业: 地震学
3楼2012-09-28 11:31:22
4楼2012-09-28 16:53:33
