小波分析中遇到的问题，求指导

tszhaoweiwen

本人想画小波分析的图，然后在论坛找到大家分享的教程（搜索下应该是比较热门和经典教程）。http://bbs.06climate.com/forum.php?mod=viewthread&tid=1617&extra=&page=1（感谢分享）
1、按照楼主给的示例数据，1959至2010，一共是52年的资料，smjfind.f中应该是M=LOG10(52.0*DT/S0)/(LOG10(2.0)*DJ)吧，算出的m=18。sm.txt的内容如下：
   2.00   2.38    2.83    3.36    4.00    4.76    5.66    6.73    8.00    9.51
  11.31  13.45  16.00  19.03  22.63  26.91  32.00  38.05  45.25
   2.07   2.46    2.92    3.47    4.13    4.91    5.84    6.95    8.26    9.83
  11.69  13.90  16.53  19.66  23.38  27.80  33.06  39.31  46.75

2、WAVELETjfind.f 中应该是    parameter(n=52,dt=1.0,s0=2*dt,dj=0.25,m=18)吧。
    怎么没有读取源数据gz.txt的语句呢？data.grd是在哪儿生成的呢？
出错啊，万望指点，谢谢！



希望得到画过小波分析的朋友帮助 ，谢谢

topmad

end of file encountered 一般指的是超出了文件的尾部

假如文件有100行 你按行读取的时候  do i=1,102 过了100的时候 就会出现这个问题

tszhaoweiwen

topmad 发表于 2013-8-19 18:52

                               
登录/注册后可看大图

end of file encountered 一般指的是超出了文件的尾部

假如文件有100行 你按行读取的时候  do i=1,102 过 ...

现在是怎么没有读取源数据gz.txt的语句呢？data.grd是在哪儿生成的呢？

兰州小高压

这有个word里面有很详细的说明

tszhaoweiwen

兰州小高压 发表于 2013-9-21 19:05

                               
登录/注册后可看大图

这有个word里面有很详细的说明

谢谢 分享

likailing2008

个人认为用这段程序比较好http://bbs.06climate.com/forum.p ... 34&fromuid=8223

二爷的屁股

顶一顶，赚点金子

xiaopingzjpq

兰州小高压 发表于 2013-9-21 19:05

                               
登录/注册后可看大图

这有个word里面有很详细的说明

谢谢分享，我来研究下下。

xiaopingzjpq

兰州小高压 发表于 2013-9-21 19:05

                               
登录/注册后可看大图

这有个word里面有很详细的说明

您好，兰州小高压，请问您的小波程序出图后，画出来的图像很纵坐标是经纬度的形势，请问怎样设置，可以画出横坐标是年份的图啊，谢谢。

qfmy

谢谢分享！