| 查看: 894 | 回复: 1 | ||
nobodyxu木虫 (正式写手)
|
[求助]
小波重构出现问题
|
|
各位大侠你们好，帮我看看下面小波去噪的程序，在matlab运行过程中出现 ??? Error using ==> horzcat All matrices on a row in the bracketed expression must have the same number of rows. Error in ==> xiaoboquzhouqi at 12 cs=[cs,coef_soft{j}]; 也就是说处理之后的小波分解系数在并入一个矩阵时，时列的个数不一样，就不能进行重构。帮我看看怎么处理，谢谢了！ 局部放电试验所采集的信号中往往混有白噪声、周期干扰信号去除。此处采用常用db系列小波中的db6小波进行9尺度的多分辨分解后，根据白噪声能量特性，估算各尺度的阈值大小，采用硬值进行处理，后进行重构。Matlab程序如下： function sd=liu_denoise(mix_signal) %此函数用于去除白躁信号＆周期性干扰信号 %输入参数mix_signal为采集到的信号波形 p=0.6745; w_dept=9; w_name='db6'; coef=cell(1,w_dept); thr=zeros(1,w_dept+1); [c,l]=wavedec(mix_signal,w_dept,w_name); %对混合信号S进行db6的9尺度一维分解 coef(1)={appcoef(c,l,w_name,w_dept)};%计算尺度为9的一维分解低频系数 cs=[cs,coef_soft{j}]; thr(1)=median(abs(coef{1}))/p*sqrt(2*log(length(coef{1})));%计算1尺度上的阈值 coef_soft(1)={wthresh(coef{1},'h',thr(1))};%对小波系数进行阈值为thr(1)的硬阈值处理 cs=[coef_soft{1}]; for j=2:w_dept+1 coef(j)={detcoef(c,l,w_dept-j+2)};%计算尺度为9到2的各尺度高频小波系数 coef1(j)={detcoef(c,l,w_dept-j+2)}; thr(j)=median(abs(coef{j}))/p*sqrt(2*log(length(coef{j})));%计算9到2各尺度上的阈值 coef_soft(j)={wthresh(coef{j},'h',thr(j))};%对小波系数进行阈值为thr(j)的硬阈值处理 cs=[cs,coef_soft{j}]; end sd=waverec(cs,l,w_name); %根据小波系数[cs,l]对信号进行重构 |
回复此楼» 猜你喜欢
一志愿华中农微生物，288分，三年实验经历
已经有7人回复
-
269求调剂
已经有3人回复
-
268分085602化学工程调剂
已经有16人回复
-
071000生物学调剂求助
已经有4人回复
-
296求调剂
已经有5人回复
-
求调剂
已经有3人回复
-
297求调剂
已经有25人回复
-
一志愿中国科学院上海有机所，有机化学356分找调剂
已经有10人回复
-
367求调剂
已经有9人回复
-
086003调剂求助
已经有17人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- gabor小波提取纹理特征的几个问题求助！
已经有10人回复
- 压缩感知过时了？
已经有18人回复
- 小波变换后系数问题
已经有6人回复
- 哪位大神来指导下数字信号处理怎么学，感激不尽
已经有30人回复
- 怎么用MATLAB生成一个三维的Shepp-Logan头骨模型并保存下来，谢谢了
已经有5人回复
- 【求助】菜鸟求助 matlab 近红外光谱小波变换问题
已经有10人回复
- 【求助】图像去噪是否可以作为博士研究的方向
已经有33人回复
2楼2014-05-05 10:04:23
