版塊導航: 正在加載中...

登錄注冊

應《網(wǎng)絡(luò)安全法》要求，自2017年10月1日起，未進行實名認證將不得使用互聯(lián)網(wǎng)跟帖服務。為保障您的帳號能夠正常使用，請盡快對帳號進行手機號驗證，感謝您的理解與支持！

24小時熱門版塊排行榜

返回列表

nobodyxu

木蟲 (正式寫手)

應助: 2 (幼兒園)
金幣: 3931.5
散金: 112
紅花: 2
沙發(fā): 1
帖子: 883
在線: 199.3小時
蟲號: 1571677
注冊: 2012-01-10
性別: GG
專業(yè): 電氣科學與工程

[求助] 小波重構(gòu)出現(xiàn)問題

各位大俠你們好，幫我看看下面小波去噪的程序，在matlab運行過程中出現(xiàn)
??? 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}];
也就是說處理之后的小波分解系數(shù)在并入一個矩陣時，時列的個數(shù)不一樣，就不能進行重構(gòu)。幫我看看怎么處理，謝謝了！

局部放電試驗所采集的信號中往往混有白噪聲、周期干擾信號去除。此處采用常用db系列小波中的db6小波進行9尺度的多分辨分解后，根據(jù)白噪聲能量特性，估算各尺度的閾值大小，采用硬值進行處理，后進行重構(gòu)。Matlab程序如下：
function sd=liu_denoise(mix_signal)
%此函數(shù)用于去除白躁信號＆周期性干擾信號
%輸入?yún)?shù)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的一維分解低頻系數(shù) 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))};%對小波系數(shù)進行閾值為thr(1)的硬閾值處理
cs=[coef_soft{1}];
for j=2:w_dept+1
coef(j)={detcoef(c,l,w_dept-j+2)};%計算尺度為9到2的各尺度高頻小波系數(shù)
               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))};%對小波系數(shù)進行閾值為thr(j)的硬閾值處理
  cs=[cs,coef_soft{j}];
end
sd=waverec(cs,l,w_name); %根據(jù)小波系數(shù)[cs,l]對信號進行重構(gòu)

回復此樓

» 猜你喜歡

291 求調(diào)劑已經(jīng)有5人回復
298-一志愿中國農(nóng)業(yè)大學-求調(diào)劑已經(jīng)有11人回復
一志愿南京理工大學085701資源與環(huán)境302分求調(diào)劑已經(jīng)有7人回復
316求調(diào)劑已經(jīng)有7人回復
材料與化工考研調(diào)劑已經(jīng)有4人回復
一志愿重慶大學085700資源與環(huán)境，總分308求調(diào)劑已經(jīng)有7人回復
一志愿211 初試270分求調(diào)劑已經(jīng)有4人回復
08工學調(diào)劑已經(jīng)有11人回復
350求調(diào)劑已經(jīng)有6人回復
求調(diào)劑一志愿武漢理工大學材料工程（085601）已經(jīng)有4人回復

» 本主題相關(guān)價值貼推薦，對您同樣有幫助:

gabor小波提取紋理特征的幾個問題求助！已經(jīng)有10人回復
壓縮感知過時了？已經(jīng)有18人回復
小波變換后系數(shù)問題已經(jīng)有6人回復
哪位大神來指導下數(shù)字信號處理怎么學，感激不盡已經(jīng)有30人回復
怎么用MATLAB生成一個三維的Shepp-Logan頭骨模型并保存下來，謝謝了已經(jīng)有5人回復
【求助】菜鳥求助 matlab 近紅外光譜小波變換問題已經(jīng)有10人回復
【求助】圖像去噪是否可以作為博士研究的方向已經(jīng)有33人回復

我就是我

1樓 2012-04-22 21:18:13

已閱回復此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

竹林里的星星

新蟲 (小有名氣)

應助: 1 (幼兒園)
金幣: 65.5
帖子: 70
在線: 18.2小時
蟲號: 2842825
注冊: 2013-12-02
專業(yè): 電力系統(tǒng)

這個程序有問題，反正我搞不懂

贊一下

回復此樓

You know you love me!

2樓2014-05-05 10:04:23