小波去噪matlab程序

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

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

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

分解5层

分解8层

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

按照你上面的代码，a3 dd3 dd2 dd1是应该都为行向量，它们的长度a3和dd3长度相等，dd2约等于其一倍，dd1又是dd2的约一倍左右，c1=[a3 dd3 dd2 dd1]是4个行向量按顺序依次组成总小波系数数组，这代码没啥问题，waverec和wavedec函数都是依靠CL组构的，它对C处理没啥难以理解的，L数组存放的就是a3 dd3 dd2 dd1长度值，靠他在C中取出最高阶的a和其他各阶的d，a和d的小波系数本来就不等长，这样才能实现DWT理论中各层分解中小波函数伸长的效果。就是因为它们不等长，用元胞数组 *** 作太麻烦，所以才把系数混排成一行，再加一个长度数组L共同 *** 作来访问数据。所以不知道你对这有啥疑问？

---