小波分析：WAVELETjfind——>sufer绘图

言深深

对小波分析已经有几个帖子了，群里有@iamyangboy在
http://bbs.06climate.com/forum.php?mod=viewthread&tid=1617
提出了问题，貌似自己解决了。
昨天晚上，@xiaoting在进行分析时出现一些小问题，经过讨论后将结果发出来，如右不妥，欢迎拍砖指导！
原包中（wavelet.rar）已经说明很详细，请允许再次啰嗦一下：
程序分析生成三个文件，其中：
nhr_mse中存放的是小波的模
nhr_mod中存放的是标准化http://bbs.06climate.com/forum.php?mod=viewthread&tid=935后的小波的模
nhr_re中存放的是小波的实部

大家出来之后的小波文件是dat格式的，如果是会写gs文件的话可以之际操作绘图，要是不行的话可以用surfer实现哦。
1.生成的dat文件转换成txt之后方便查看
这个是用@xiaoting提供的数据 data.txt (200 Bytes, 下载次数: 148)

2011-10-17 10:13 上传

点击文件名下载附件 进行计算的
  parameter(n=50,m=17)    !注意这边是50哦O(∩_∩)O~
   dimension  wtr(n,m),wtm(n,m),wte(n,m)
       open(21,file='nhr_re.dat',access='direct',form='binary',recl=10)
       open(61,file='nhr_mse.dat',access='direct',form='binary',recl=n*m*4)
       open(31,file='nhr_mod.dat',access='direct',form='binary',recl=n*m*4)
        read(21,rec=1) ((wtr(j,i),j=1,n),i=1,m)
        read(31,rec=1) ((wtm(j,i),j=1,n),i=1,m)
        read(61,rec=1) ((wte(j,i),j=1,n),i=1,m)

       open(22,file='nhr_re.txt')
       open(62,file='nhr_mse.txt')
       open(32,file='nhr_mod.txt')
        write(22,*) ((wtr(j,i),j=1,n),i=1,m)
        write(32,*) ((wtm(j,i),j=1,n),i=1,m)
        write(62,*) ((wte(j,i),j=1,n),i=1,m)
end


2.查看生成的txt文件中可以看到小波很好的分成m行，n列，这就是小波的实部和模方啦，下面我们用sufer绘图
parameter(m=17,n=50)
real x(m,n)
open(1,file="nhr_re.txt")
do i=1,m
read(1,*) (x(i,j),j=1,n)
enddo
close(1)
open(2,file="pic.txt")
do i=1,m
do j=1,n
write(2,"(2i8,f10.2)") j+1958,i,x(i,j)
enddo
enddo
close(2)
end

这个小程序可以帮助你将数据方程我们需要的三列格式，这样suferhttp://bbs.06climate.com/forum.php?mod=viewthread&tid=946能够很好的识别，结果如下

传说中的谁

深深的沙发卖不掉，嘿嘿不错，来学习了，虽然我不知道小波分析是干嘛用的。也许以后用得着的说，收下了。

言深深

传说中的谁 发表于 2011-10-17 10:59

                               
登录/注册后可看大图

深深的沙发卖不掉，嘿嘿不错，来学习了，虽然我不知道小波分析是干嘛用的。也许以后用得着的说 ...

what ever 3q3q

善人/jw

好东西好东西~~
论坛里有好多人在问呢~~
这样大家应该就OK了~~

阿华田

深深发力了啊

言深深

善人/jw 发表于 2011-10-17 13:22

                               
登录/注册后可看大图

好东西好东西~~
论坛里有好多人在问呢~~
这样大家应该就OK了~~

就是学习学习吧，呵呵

言深深

阿华田 发表于 2011-10-17 13:43

                               
登录/注册后可看大图

深深发力了啊

为什么你们都有幸运卡！！！汗···

xiaoting

谢谢深深版主的正解

我上午又练习了一遍
当n=50时算出来的m=18哦，n=40时算出来的m=17
还有我查看txt时数据不是18行 50列而是180行 5列

言深深

xiaoting 发表于 2011-10-17 15:06

                               
登录/注册后可看大图

谢谢深深版主的正解

我上午又练习了一遍

貌似m需要人工调一下的吧？试试看

chliu2010

呵呵，值得分享