【求助】关于小波变换图！

lvxiaohong

都完全可以，你把你要做的问题告诉我把！问题的目的！

lsapan

我主要是想通过树立统计和小波分析来揭示本市近些年的洪涝和强降水频数的演变规律，其他的我已经弄好了，现在就还差小波这一方面，现在好象教流行小波分析的，所以我也引进来用用。
最主要是要找到洪涝频数是否具有某种周期性（所以我想看一下他的剖面图的效果，我想铁上去的）或者是找到哪个年代是洪涝多发期等规律吧，在由此做出长期预报:下个10年或者20年是否是多发年，大概有几次，强度如何等等。
不过，目前还是前一步啦，如果大家有什么好的建模方法的话，还请推荐一、二啊！！！！！

lvxiaohong

看小波变换的剖面图我认为没有多大的意义，因为，他是在不同的尺度上的，相当于积分在大的尺度上这些值就高。但是，用他的某个尺度，看他的趋势是可以的。对这些我不太懂，用一些统计的方法应该可以吧。还有好像目前小波域的Markov模型用来预报可能可以，你去看看。下午你的问题，等高线最简单你可以在画出后，逐个编辑，否则你的用set语句，逐个设置这些Patch对象。[c,h] = contour(w)；clabel(c,h), 然后逐个设置，h(1),h(2),...这些都是patch对象。

lsapan

真是很不好意思，又要麻烦你们了，因为我急着要交差，图不能用彩色的，所以要把负数曲线用虚线画，正数部分用时线画，儿标志的数据的小数位数好象也太多了点，最好把x轴的标志改为从（1961-2000）这样便于也看清。lvxiaohong 教我的方法我弄了好半天，找不到门，麻烦在说的细点，我画不出来，很着急的！
  
以下为具体的程序：
>> b=[3 5 2 4 1 2 3 5 8 6 4 5 5 9 12 6 6 2 1 2 1 4 11 3 1 0 5 6 12 6 1 11 3 6 7 11 5 14 8 9]
>> y=b/215*100
>> x=linspace(1961,2000,40);
>>w=cwt(y,[2:2:18],'mexh');
>>[c,h]=contour(w,20)
>>clabel(c,h)
………………

lvxiaohong

你看等高线数据若取为10个的话，那么就是有10个Patch对象。其中一部分对象是负值，一部分是正值。你逐个设置。至于改坐标轴的范围，很简单，画图的时候可以写明范围，也可以用独立的程序设置。我可以给你编一个这样的程序，你自己先试试。晚上我再来，如果再不懂，到时候再说

lsapan

我对于patch不熟悉，其中 PATCH(X,Y,Z,C)，我用x为年数，y为a ，z为w，可是sizs不对，另外其中得c好象时指颜色得嘛，从哪里给它变成实线或虚线？在画图时怎样改变x坐标得值，我是很不懂得那种哦，要具体得例子 才看得懂 ，呜呜 ……

lsapan

我明天就要交了，能不能告诉我具体得式子？

lvxiaohong

用C画完后，然后用set(h(1),'linestyle',':');如果设置其他patch的话，就设置其她的patch句柄。改变坐标的范围，很简单的，在用contour的时候，写上x,y坐标的范围即可。其实改变成任何范围都行，当然需要自己编程实现。主要是设置她的tickdir and ticklabel  属性。你自己看帮助一定不会有问题的。