怎样在matlab中自己编写一个小波的程序

怎样在matlab中自己编写一个小波的程序

a1.a2,d2,d1是ca1,ca2,cd1,cd2这些小波系数的重构。ca1,ca2,cd1,cd2是小波系数，它们的数据点数随分解层次的增大而减少，这就难以与原始信号对比分析，通常会经过重构变为与原始信号个数相同的a1.a2,d2,d1，从物理意义上讲，只有a1.a2,d2,d1才是有实际量纲的信号，ca1,ca2,cd1,cd2是没有量纲和物理意义的。

1、基于低频融合策略的小波图像融合程序（Imfus.rar）；

2、嵌入式零树小波（EZW）算法的图像压缩编解码程序（EZW.rar）；

3、多级树集合分裂（SPIHT）算法的图像压缩编解码程序（SPIHT.rar）；

4、一种SPIHT改进算法（Mod-SPIHT.rar）。

其中都包含了基本的小波图像分解与重构程序，使用的是基本的Haar小波，根据Mallat算法编写。

哈哈！傻眼了吧，理论和实际不可兼得啊！对于你这个问题有如下建议：

小波去噪的试验，十个有九个都喜欢用正、余弦函数，但由于小波函数的几何特征，其实不易得到满意效果，你只能选择线性强的小波基，即对称性强并且光滑的小波，嘿嘿，貌似能做DWT的所有小波基中只有sym8的对称性最强（参看http://zhidao.baidu.com/question/585223273699457565.html?from=pubpage&msgtype=2），这样你将会得到其降噪效果好于其它小波基的谬论，实际应用中如果原始信号本身特征就不规律对称则有可能其它小波基的处理效果会好于sym8，这需要试错的方法试验，是你这个原始信号对称的前提假设造成的，但不可认为sym8适合你这个特殊例子就得出结论，你这试验只适合你的假设，其实并不适合所有情形，个人认为其实不具说服力和可比性。如同我家钥匙开不了你家门，是不能得出我家钥匙不如你家钥匙的结论滴。

当你的噪声特征与信号的特征的频率相近时，小波也无能为力，它不是神，也一样分不出啥是噪声，所以其一是可增加分解层数，你这个信号只有100个数，5层已经很高了，再增大也没啥用了，可能会过多显示小波基的特征，造成扭曲失真（如果用SWT会好很多，但需要自己编制函数）；其二是参数SCAL可以改为伸缩的sln，而不是固定的one，这样分解层数和SCAL都将起作用，你可以试着改改玩，效果还行。

分解5层

分解8层

你可以试试只分解一层的状况，频率最低的几处噪声会保留下来哦！

---