1樓:匿名使用者
t=wpdec(y,5,'db40');
%對訊號
copy
baiy進行小波包分解,層數為du5,得到zhi的t為小波樹dao,plot一下就可看到
a10=wprcoef(t,[1,0]);
%a10是對節點[1,0]進行重構後得到的訊號。貌似沒有對那一層重構這一說法吧,只能是對某層的某個節點進行重構。節點的編號你可以從小波樹中看出來
這是我的做法,不過用的是小波包分解。不知對你有沒有用
matlab中,離散小波變換函式dwt獲得的中的低頻係數與高頻係數的具體含義是什麼,獲得這兩個引數有什麼用
2樓:邵碩小童鞋
橫座標就是取樣點,縱座標是訊號的幅值,橫座標除以取樣頻率就可以得到音訊訊號的時長。
matlab中小波變換cwt函式中的尺度與小波函式具體尺度的關係一樣嗎??
3樓:
這問題看似簡單實際好有意思,有些不知如何說起,尺度就是尺度,不是時間也不是空間,所以從理論上沒有單位,你沒見過哪個文獻中尺度是有單位的吧?它是cwt中的概念,dwt通常是不提尺度的,要簡要定義,尺度表明的是小波函式伸縮的狀態,例如尺度1.2,表明將原始小波函式伸長1.
2倍後進行cwt,得到尺度1.2下訊號的cwt小波係數。「那麼這個2到底對應haar小波中尺度a多少米呢」這話本身從理論上就沒有定義過,也不知道如何定義,或根本沒有必要定義。
為啥不知如何說起呢?是因為尺度在時域是可以通過小波函式的支撐長度和訊號本身的取樣頻率換算為頻率的,例如尺度1到底對應多少hz的頻率(到底代表多少是根據你訊號的取樣頻率變化的),這之間的轉換可以在網上搜尋「小波時頻圖」有詳細**和解釋。但空域中沒有定義頻率hz這些玩意兒,一般定義的是波數,你得進行時空的轉換,可以參考上面的**理一理。
所以尺度不代表多少米或多少秒,硬要轉換可以轉為頻率和波數,但概念上它並不代表或描述頻率和波數,只是小波函式伸縮狀態的描述,可以對應頻率和波數。
一維連續小波變換cwt呼叫方式
4樓:中地數媒
(1)coefs=cwt(s,scales,『wname』)
(2)coefs=cwt(s,scales,『wname』,『plot』)
說明:該函式為一維連續小波分析函式,其中coefs為連續小波變換後的返回係數wf(a,b)矩陣,係數以行方向儲存在矩陣中。
[例6-14]對matlab中所帶有的noissin訊號進行連續小波變換,尺度a分別為0.12,0.24,0.
48,1.2,2,4,6,8,10,小波函式用db3,請求出連續小波變換後的係數。程式:
load noissin;%裝載訊號
s=noissin(1:100)
ls=length(s);%計算訊號點的個數ls
%對s進行一維連續小波變換,把返回係數存到矩陣w中
w=cwt(s,[12.12,10.24,15.48,1.2,2:2:10],『db3』,『plot』)
xlabel(『時間』);
ylabel(『變換尺度』);
title(『對應於尺度a=0.12,0.24…小波變換系數的絕對值』);
執行程式後,返回矩陣為一個9×1000矩陣。在此為節省篇幅,我們不將結果列印出來,讀者可自己上機執行該程式觀察結果。
圖6-41 一維離散小波變換dwt
matlab 離散小波變換的階數與尺度的關係???
5樓:
呵呵!問題有意思。第一段的理解按照matlab的設定是完全正確的,j在dwt中就是階次(或叫層次),dwt1,2,3.。。。
階(層)對應cwt的尺度2,4,8.。。。。這樣隨著尺度的增加,可以認為低頻分的越來越低頻,越接近真實訊號的總體趨勢。
第二段,中的尺度真的就是通常生活中尺度的意思,就是大小,例如汽車就是比自行車大,就這麼個意思。這種理解更貼近於cwt對大小的理解,訊號的波長長,頻率低,就是尺度大的意思。在dwt中都為了避免這問題產生歧義不使用尺度而直接用階次一詞,這本書的作者可能不夠專業,這類人通常將cwt的尺度用於dwt的解釋中,是很讓人抓狂的,只能說明不太嚴謹,考慮不夠周全。
「較大尺度的訊號 進行小波變換,再選取其中的低頻部分在原尺度的1/2尺度上再進行小波變換。」這話相當外行,從表達到用詞都不嚴謹,很多專有的小波概念根本不懂,易產生歧義。按他的意思用matlab中dwt的設定可以改為,較低頻(長波長,大尺度)的訊號進行dwt,得到a1和d1,再選取a1d1中的低頻部分a1在原小波函式縮短1/2後再進行dwt,得到頻率更高的a2和d2(與matlab的設定相反,後面會有說明).
這樣就避免了尺度一詞對理解造成的困擾。
你可以看到這裡使用的是小波函式的縮短,而不是伸長,這恰好與matlab的設定相反。所以造成了dwt階次與cwt的尺度對應關係與matlab的設定不同。在dwt階次與cwt尺度或頻率的設定中一般是有兩種的,一種是以matlab使用的定義,它用小波函式的伸長做dwt,dwt1,2,3.。。。
階對應cwt的尺度2,4,8.。。。,先得到的a1d1頻率與其它ad相比最高,越向下分,頻率越低。另一種是出於數學公式推導的方便,設定與之恰好相反,用小波函式的縮短做dwt,a=1/2^j,例如dwt1,2,3。。。.
階對應cwt的尺度1/2,1/4,1/8。。。先得到的a1d1頻率最低,越向下分,頻率越高。
在非matlab的小波書籍中,尤其是數學推導較多的書中,為方便泛函空間的公式推導通常使用後一種的設定,這就造成對dwt階次和頻率變化理解的困擾,可笑的是很多文章連抄都不會抄,公式的推導用後一種,可描述解釋卻用前一種matlab的設定,簡直就是糟蹋行當,悲哉!
關於連續小波變換的幾個問題,求教 80
6樓:
首先應明白連續或不連續多指數學的概念而已,應用中的訊號都是離散的,只是你的取樣足夠高就可認為是連續的,所以小波變換中關心的是點數問題,而不關心訊號是否連續。對於cwt或dwt其連續與否不是指分析訊號,而是你說的a或b的問題,但你仍可以借鑑上面對於訊號連續的理解。cwt中a是連續的,b其實就是點數,也可認為是連續的。
最早的dwt是沒有mallat演算法的,那時a是以2的冪次方變化離散,b卻是連續變化的,即二進小波變換。這種變換很雞肋,還不如直接做cwt。dwt的應用之所以遠遠多於cwt就是引入了mallat演算法,好處是終於可以分解和重構訊號了,這種方式對訊號特徵的研究非常有利。
dwt的核心思想其實就是cwt引出的伸縮和平移的概念,a以2的冪次方變化實現了小波的伸縮,b通過下抽樣實現了小波的平移。從實際應用中進行小波變換的目的和效果來看,cwt中2/4/8/16/32的小波係數結果應該對應dwt中的階次(層數)1/2/3/4/5的小波細節係數(或更準確的是重構後的小波細節,因為cwt的係數個數是不變的等於原訊號長度,但dwt細節係數是每層近似減半的,重構後才會等長,b也是姑且認為是減半的不連續吧)。
再追問吧,第二問題可能更多,我儘量精簡。哎,幹嘛要把問題寫在一起,這就是麻煩啊,你必須追問我才能再寫!
怎麼用matlab實現小波變換???急!!!
matlab小波變換跪求用Matlab將影象經過小波變換
detcoef提取 高頻係數,到相關的論壇上看看吧 找找帖子 matlab怎麼用小波包進行影象去噪 小波影象去噪的方法大概分為3類 1 基於小波變換摸極大值原理 2 基於小波變換系數的相關性回 3 基於小波閾值的去噪答。基於小波閾值的去噪方法3個步驟 1 計算含噪聲影象的小波變換。選擇合適的小波基和...
怎樣從零基礎學懂小波變換,Matlab小波變換的係數和恢復問題
嗯,我也是 bai剛剛在小波du變換方面入了點門。其實如zhi果僅僅dao是入門的話,需要了解的並版 不是太多,到圖書館權的訊號分析分類下找一本小波變換的相關書籍就可以了。需要注意的是很多所謂的經典書籍可能並不適合入門,反正我是找了幾本看起來比較實用的書學習的,互相參照著看,加上書裡面的例項很快就能...
麻煩小波變換高手分別解釋一下matlab小波基中db1,db
你可以開啟小波基來看看它們的數值,畫出圖來看的話更加直觀,附圖是我畫的,是cdf7 9小波基,和matlab裡面的bior4.4類似。db sym都是小波的名稱,或者說是 族 後面的數字可以看作是這個族裡面的長幼次序了。你用圖把它們畫出來再比較,就一目瞭然了。不妨從簡單的haar小波基入手,看看它的...