24小时热门版块排行榜    

查看: 249  |  回复: 0

xxping2007

新虫 (初入文坛)

[求助] 小波系数方差求Hurst指数的程序

我不太懂小波系数方差法,找了一篇文献上面有编程原理,参见附件。然后又找了一段程序,但是求随机一组数据H总是负的,理论上是大于0小于1.但不知道是否是程序问题还是哪有问题。另外我想画出H随时间的图,怎么画呢?谢谢啦。
程序如下
data=100*rand(1,4096);
[d,a]=wavedec(data,10,'db2');%利用小波“db3”对信号data进行10次分解,求得的近似系数存放在数组d中,细节系数存放在数组a中)
[cd1,cd2,cd3,cd4,cd5,cd6,cd7,cd8,cd9,cd10]=detcoef(d,a,[1,2,3,4,5,6,7,8,9,10]);%提取小波分解细节系数分量
yvals=[var(cd1),var(cd2),var(cd3),var(cd4),var(cd5),var(cd6),var(cd7),var(cd8),var(cd9),var(cd10),];
xvals=[1,2,3,4,5,6,7,8,9,10];
logx=xvals;
logy=log2(yvals);
p2=polyfit(logx,logy,1);%一次线性拟合,“1”代表1次
figure;
plot(logx,logy);
hurst=(p2(1)-1)/2 %p的第一个就数就是斜率即p(1).Hurst exponent(指数) is the slope(斜率) of the linear  r(r=2H+1)
回复此楼

» 本帖附件资源列表

  • 欢迎监督和反馈:小木虫仅提供交流平台,不对该内容负责。
    本内容由用户自主发布,如果其内容涉及到知识产权问题,其责任在于用户本人,如对版权有异议,请联系邮箱:xiaomuchong@tal.com
  • 附件 1 : 基于自相似检测DDoS攻击的小波分析方法.pdf
  • 2015-01-24 11:09:17, 707.24 K

» 猜你喜欢

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 xxping2007 的主题更新
信息提示
请填处理意见