liushuangq05 发表于 2011-3-18 13:03:23

日本海啸的matlab效果(投影坐标)

活跃下matlab版气氛,这里给一个这两天做的日本9级地震的海啸模拟效果给大家,完全用matlab来做,使用worldmap函数和pcolorm函数.效果如下:

mhkmars 发表于 2011-3-18 14:31:23

1# liushuangq05
很有意思，不知道lz使用的那个工具箱做的呢？觉得很有趣，想学习一下

turbinesv 发表于 2011-3-19 12:58:18

2# mhkmars

不错，不知道用的什么理论做的。
只是图好看呢还是有实际效果。

wxm27811 发表于 2011-3-20 14:45:04

很羡慕。。。。

alexqxp 发表于 2011-3-21 20:26:59

喜欢贴近时事的例子。学习下相关函数。海啸相关数据没有，可以的话楼主发程序大家看看。
>> worldmap('japan')
>> load coast
>> plotm(lat,long)

alexqxp 发表于 2011-3-21 20:45:06

本帖最后由 alexqxp 于 2011-3-21 20:48 编辑

学习下pcolorm函数。
h = pcolorm(Z)
h = pcolorm(Z,npts)
h = pcolorm(lat,lon,Z)
h = pcolorm(lat,lon,Z,PropertyName,PropertyValue,...)

h = pcolorm(Z) projects the data grid Z on a graticule grid the size of Z between the latitude and longitude limits of the current map axes. The handle h of the displayed surface can be returned.
h = pcolorm(Z,npts) results in a graticule grid defined by npts, which is a two-element vector of the form . The default npts is .
h = pcolorm(lat,lon,Z) allows three other methods of defining the graticule grid. If lat and lon are matrices, they represent the actual graticule vertices as might be returned by meshgrat. If vectors, they are the representative coordinates of the rows and columns, respectively, of such a grid. If they are two-element vectors, they are treated as latitude and longitude limits, and a graticule mesh the size of the default npts is calculated.
h = pcolorm(lat,lon,Z,PropertyName,PropertyValue,...) allows the input of property name/property value pairs to control the surface object properties. Any property supported by the standard MATLAB function surface except XData, YData, and ZData can be altered in this manner.
RemarksThis function warps a data grid to a graticule mesh, which itself is projected according to the map axes property MapProjection. The fineness, or resolution, of this grid determines the quality of the projection and the speed of plotting it. There is no hard and fast rule for sufficient graticule resolution, but in general, cylindrical projections need fewer graticule points in the longitudinal direction than do complex curve-generating projections.
Examples
load topo
axesm miller
pcolorm(topo,)
demcmap(topo)
tightmap

根据例子来看，将海面高度数据替换topo，或许就可以将经纬度对应的海面高度映射到日本近海。随着时间高度发生变化，最后形成随时间变化的海啸图。当然，仅仅是分析猜测。欢迎楼主指正。

memo1221 发表于 2011-3-22 10:43:59

很有意思，看来我还太年轻

liushuangq05 发表于 2011-3-22 12:16:29

alexqxp版主说的很对,不过没那么复杂的.worldmap('hi',,,'lineonly');就可以把我的区域显示出来了.然后利用pcolorm来覆盖一下就可以.海啸模拟用的是线性球坐标浅水波方程,模拟出来的结果还是比较好的,和华东片区的7个验潮站的数据基本吻合.
对模拟的各时刻提取各点的最大值,就可以做浪高分布图出来.

alexqxp 发表于 2011-3-22 20:08:06

看了下楼主的相关帖子，似乎做灾害模拟仿真的。把楼主的帖子收集了点放在下面，学习了下。几年的前后帖子连起来看，楼主也是经过了五六年的学习和提高才走到今天这一步的，看到一个技术人员的成长经历让人感慨。
matlab做的海啸传播模拟图
http://forum.simwe.com/viewthread.php?tid=961219&highlight=
这个是ansys做的。汶川地震静态位移模拟及命令流
http://forum.simwe.com/viewthread.php?tid=892713&highlight=

liushuangq05 发表于 2011-3-22 22:11:33

呵呵,让我申请进入simwe论坛的起初主要是自学ansys时问题多多,然后baidu到simwe的,后来觉得matlab版块更有意思.

zhengxing2000 发表于 2011-3-28 15:22:19

早就知道matlab很强大，没有想到matlab这么强大。
现在在做混沌，做不出程序，害怕搞理论找不到工作，很是纠结。

longyuff 发表于 2011-3-28 17:13:04

哇，很厉害哟。

sdojiangsdo0729 发表于 2011-3-28 19:29:40

顶一下呀，以后得好好学习下matlab了

zzs10105218 发表于 2011-4-17 11:23:20

呵呵，看来可以学习的东西还是很多的。

zhutou0517 发表于 2011-4-22 09:56:05

对楼主顶礼膜拜！

chuanhailc 发表于 2011-4-22 22:30:27

各种牛人啊 ！

geoyang 发表于 2014-7-15 09:01:35

alexqxp 发表于 2011-3-22 20:08
看了下楼主的相关帖子，似乎做灾害模拟仿真的。把楼主的帖子收集了点放在下面，学习了下。几年的前后帖子连 ...

楼主，你给的链接现在打不开了啊，能不能更新下:loveliness:
页: [1]
查看完整版本: 日本海啸的matlab效果(投影坐标)