在matlab中4维的数据如何进行Bootstrap resampling

维欻? · 发表于 2013-7-23 11:15:10

liuzf406 发表于 2013-7-23 10:36

登录/注册后可看大图

物理意义倒是好说，关键是算法不会啊。如图，我计算了两个不同太阳活动阶段合成的hgt距平，现在我想用boo ...

你的意思是程序不会？算法就是：对空间一个点拿出来就是一个时间序列，对这个时间序列进行有放回的重采样
利用随机数就可以完成啊

liuzf406 · 发表于 2013-7-23 11:50:54

维欻? 发表于 2013-7-23 11:15

登录/注册后可看大图

你的意思是程序不会？算法就是：对空间一个点拿出来就是一个时间序列，对这个时间序列进行有放回的重采样 ...

是啊，对一维的还好说，对这么多维的我就晕了。唉！

维欻? · 发表于 2013-7-23 14:41:50

本帖最后由维欻? 于 2013-7-23 14:43 编辑

liuzf406 发表于 2013-7-23 11:50

登录/注册后可看大图

是啊，对一维的还好说，对这么多维的我就晕了。唉！

一维的会那就好办了
假如你的数据为data(360,180,20,120) %lon*lat*lev*time
那么你就针对每一个点 x=data(i,j,k,:) 做bootstrap得到一个y
最后整合成一个newdata(i,j,k,:)=y;
循环1000次就可以了
这样就可以了
如果你的数据量比较大的话很很很费内存的，建议在大机器上做，PC是做不了的，嵌套这么多的循环matlab会很慢

liuzf406 · 发表于 2013-7-23 14:55:36

维欻? 发表于 2013-7-23 14:41

登录/注册后可看大图

一维的会那就好办了
假如你的数据为data(360,180,20,120) %lon*lat*lev*time
那么你就针对每一个点 x ...

谢谢指导啊！
是啊，数据量比较大，对于hgt还好说，2.5X2.5的半球数据,但对于气温和降水,可能就比较大了,0.5X0.5. 无论如何,现在知道思路了,我试试做.

谢谢楼上各位的指导!

飘渺居 · 发表于 2015-11-17 20:56:44

请问楼主搞定了嘛？弱弱的问下Fortran可以实现Bootstrap resampling 吗？

MlleWXY · 发表于 2016-1-8 05:53:57

		自动登录	找回密码
密码			立即注册

[程序设计] 在matlab中4维的数据如何进行Bootstrap resampling