㈠ 用matlab軟體怎麼進行小波分析
它自己有小波工具箱,wav開頭的,裡面有很多相關函數,一般直接用就是了。初學建議看幫助裡面小波工具箱的demo,可能更好理解。
㈡ 用matlab中工具箱進行小波去噪步驟
matlab讀取excel文件比較方便,建議你把數據放到xls文件中保存,然後在matlab中用xlsread這個函數讀取出來。版
讀取出的數據應該是一權個一維數組了,用plot畫出圖的話,就是常見的曲線。
然後做小波分解:選用你覺得合適的小波基,例如haar,然後用這個小波基做小波分解,再把高頻部分去掉,然後用低頻部分還原,就得到了去噪後的信號。
其實你這個問題估計也可以用神經網路或者其它曲線擬合一類手段來解決。具體的情況要根據數據特徵來判斷。
以上。
專業路過的老狼
㈢ 小波分析在層序地層劃分中的應用
1.小波分析簡介
20世紀80年代後期至今,一種著名的、在各行各業有重要應用價值的數學理論和方法技術在科學技術界得到了廣泛的重視和採用,它就是被譽為「數學顯微鏡」的小波分析(李世雄,1994)。小波分析的主要功能和特點是,它具有多分辨分析或多尺度分析功能,可以把信號分解成各種不同的尺度成分;它具有很強的局部分析功能,同時具有時間(或空間)域和頻率域的局部分析性質,它可自動地通過伸縮、平移聚焦到信號的任一細節對其加以分析(侯遵澤,1998)
(1)小波分析基本原理。小波(wavelet),即小區域的波,是一種特殊的長度有限、平均值為0的波形。它有兩個特點:一是「小」,即在時域都具有緊支集或近似緊支集;二是正負交替的波動性。如果用小波和構成傅里葉分析基礎的正弦波做比較的話,傅里葉分析所用的正弦波在時間上沒有限制,從負無窮到正無窮,但小波則傾向於不規則與不對稱。
傅里葉分析是把信號分解到一組相互正交的正弦波上的,也就是基函數,我們可以把基函數看成是度量信號某些特徵的一把「尺子」,傅里葉分析度量的就是信號的頻譜特徵,但是如果這把「尺子」過於規則,有時候就不能十分精確地表達信號蘊含的信息,而在小波分析中,「尺子」換成了規則程度更低的小波函數,從而可以更加有效地表達信號中信息的成分。
小波變換對不同頻率在時域上的取樣步長是調節性的,即在低頻時小波變換的時間解析度較差,而頻率解析度較高;在高頻時小波變換的時間解析度較高,而頻率解析度較低(圖2-13),這正符合低頻信號變化緩慢而高頻信號變化迅速的特點(胡昌華,1999)。這就構成了利用小波變換進行信號分析的基礎。
圖2-13 數字信號的小波變換
(2)一維連續小波變換。小波變換實際上是求取信號在各小波函數上的投影值。每個小波函數均由一個母小波函數經過尺度伸縮與時間平移得來的。信號分析的一般思路就是分解與組合,尋找一組最能代表信號特徵的函數形式,將信號用這些量來逼近,或者寫成這些量的線性組合的形式。
小波分析的思想來源於伸縮和平移方法:對波形的尺度伸縮就是在時間軸上對信號進行壓縮與伸展,而時間平移就是指小波函數在時間軸上的波形平行移動。
(3)離散小波變換。由於連續小波變換的伸縮和平移系數是相互獨立的,所以通過伸縮和平移得到的各個小波函數之間有一定的相似性,但由於這兩個系數之間的獨立,就引入了信息的冗餘。在解析度一定的情況下,描述了多餘的信息,使得反映信號特徵的一些參數相互重疊,給我們的分析帶來不便,但這些特點可以用在本身就有自相似性的信號上,可以讓我們更清楚地看到信號自身的自相似性。
此外,由於冗餘信息的存在,也使得小波逆變換的重構過程不唯一,也就是說,由同一母小波生成的不同的小波變換函數可能重構成同一個信號。為了減少冗餘信息,就引入了離散小波變換的概念,其中的伸縮和平移系數是可數的,重構過程用求和的形式給出。如果伸縮和平移系數滿足一定的對應關系,則稱為二進小波變換(尺度以2的冪的形式給出)。離散小波變換主要是建立在二進制小波變換的基礎上的。
測井曲線數據也恰好是離散數據,符合離散變換的要求。在利用小波分析進行層序地層劃分時,主要是對測井曲線進行多尺度分解,得到不同尺度下的小波變換圖,利用其表現出來的特徵來劃分不同級次的層序。
2.利用小波分析進行層序地層劃分
利用小波分析方法是層序地層劃分方法上的一種新的嘗試,其目的是盡量減少層序劃分過程中的主觀因素,依靠地層自身表現出來的客觀特徵來識別層序、准層序組以及准層序。在我們研究的沉積岩地層中,沉積物的特徵可以反映沉積時水體的特徵。隨著沉積水深的變化,沉積物中多種特徵都會相應的發生變化,如放射性物質含量、有機質含量或其他微量元素的含量等,這種變化就會在相應的測井曲線上反映出來。而沉積水深變化受到了多種因素的影響,有長期和短期的旋迴,是多個不同周期旋迴的疊加,因此測井曲線應該是沉積地層中某種隨水深變化的特徵的多種頻率變化的響應的疊加。也就是說,測井曲線中包含著沉積水深不同周期變化的信息,是多個沉積水深變化周期相互疊加的響應。而小波分析能夠將信號分解為不同頻率不同周期的旋迴,因此,小波分析的特點恰好可以和測井曲線的特點相對應,利用小波分析的方法可以比較准確地將測井曲線中相互疊加的反映水深變化的不同周期的信息分別識別出來,識別出的這些信息就可以用來進行沉積旋迴的劃分。
同時,小波分析方法還可以幫助解決傳統研究方法所不能解決的一些難題,如大段單一岩性地層中的沉積旋迴識別。大段單一岩性尤其是大段泥岩、頁岩,並不是一個小的沉積旋迴里沉積的產物,相反,應是一個相當長時期沉積下來的,但是通過傳統的岩性劃分方法卻很難將其劃分開,這就給准層序甚至准層序組的劃分造成了困難。小波分析方法可以較好地解決這一問題,利用這種方法可以從測井曲線的細微變化中識別沉積間斷和沉積旋迴。
(1)測井曲線的選擇。不同的曲線具有不同的地質含義,進行相同的變換可能會得到不同的結果。但在研究中通過對GR、AC、COND、電阻率等多條曲線進行小波變換後對比發現,不同測井曲線所得出的變換結果盡管形態上不完全一樣,但在旋迴的劃分上卻比較一致(圖2-14)。圖中曲線a是COND測井曲線經過db5小波變換後的結果,曲線b是同一井段AC曲線變換後的結果。出現這個結果是由於雖然不同的曲線代表著不同的地層響應,會呈現出不同的特徵,但地層中各種參數的變化主要受沉積環境的影響,會隨著沉積環境的旋迴變化呈現出基本一致的旋迴特徵。這也從一個方面反映了小波變換在沉積旋迴劃分中的客觀性。因此,可以選擇目標井的測量精度較高、質量較好的曲線來進行小波變換,進而進行沉積旋迴的劃分。
圖2-14 對COND和AC曲線進行小波變換結果對比
(2)小波的選擇。同傅里葉分析不同,小波分析的基(小波函數)不是唯一存在的,所有滿足小波條件的函數都可以作為小波函數,那麼小波函數的選取就成了十分重要的問題,實際選取小波的標准主要有以下三種。
1)自相似性原則:對二進小波變換(因為在正交小波變換中,取樣的方式就是按照小波函數取樣的,所以不存在這個問題),如果選擇的小波對信號有一定相似性,則變換後的能量就比較集中,可以有效減少計算量。
2)判別函數:針對某類問題,找出一些關鍵性的技術指標,得到一個判別函數,將各種小波函數代入其中,得到一個最優原則。
3)支集長度:大部分應用選支集長度在5~9之間的小波。因為支集太長會產生邊界問題,支集太短不利於信號能量的集中。
但在實際應用中,因為大部分信號的信息量太大,很難找到相應的模式,因此主要是依靠經驗來選取。根據前人研究經驗及作者對不同函數所做結果的對比,本書採用的是Daubechies小波,階數為5。
Daubechies小波是由著名小波學者Ingrid Daubechies所創造,她發明的緊支集正交小波是小波領域的里程碑,使得小波的研究由理論轉到可行。這一系列的小波簡寫成dbN,其中N表示階數。
(3)工作流程。測井曲線能比較准確地反映井旁地層的電性、物性等特徵,但往往會受到測井儀器、鑽井泥漿等其他非地層因素的干擾,且不同頻率的旋迴相互疊加,對正確識別和劃分沉積旋迴造成一定的影響。小波分析能真正消除干擾信號,放大真實信息,按不同頻率反映出測井曲線中包含的真正旋迴特徵,以此劃分不同級別層序單元,同時還可以在劃分高精度沉積旋迴的基礎上,與Fischer圖解相結合劃分出體系域。
MATLAB軟體的小波分析工具箱是一種比較常用的工具。MATLAB是Math works公司於1982年推出的一套高性能的數值計算和可視化軟體。MATLAB的推出得到了廣大專家學者的廣泛關注,其強大的擴展功能為各個領域的應用提供了基礎。各個領域的專家學者相繼推出了MATLAB工具箱,包括信號處理、神經網路、圖像處理、小波分析等。其中小波分析工具箱可以滿足對測井曲線進行小波變換的需要。
圖2-15 小波分析流程圖
在對測井曲線進行小波變換時,首先需要對所研究層段的頂底界面進行准確的標定,然後將需要變換的該深度段的測井曲線數值建立單獨的文本文件作為原始文件。將原始文件導入後保存成.m格式的信號文件。選擇MATLAB軟體的wavelet(小波分析)工具箱進行離散一維小波變換,小波類型選擇db,階數為5,最大級數定為12,將上述參數選好後進行分析,即可得到一組12條不同級次的db5小波變換曲線(圖2-15)。此外對其進行連續一維小波變換,可以得到小波的頻譜分析圖,選擇合適的最大顯示值,根據頻譜圖上圖形的閉合方向可以區分出層序的界面和層序單元(圖2-6,圖2-7)。
(4)單井分析實例。牛100井位於牛庄窪陷西部,地層以砂泥岩互層為主,岩性變化較快(圖2-16)。利用小波分析方法對AC、R25兩條測井曲線進行了一維連續小波變換,分別得到其小波變換譜系圖,對AC曲線進行了一維離散變換,得到不同階數的小波,根據與地震、測井及錄井岩性資料的對比,選用d11,d9,d7三個層的小波變換曲線分別進行層序、准層序組和准層序的劃分。
將傳統劃分方法所得的結果與小波分析方法所得結果進行對比可以比較明顯的看出,在層序和准層序組的劃分上兩種方法劃分的層序單元基本一致,可以相互驗證。在准層序級別上的劃分,小波分析方法的優勢就比較明顯地體現了出來,劃分的旋迴數較多,精度也有提高。這也正是小波分析作為「數學顯微鏡」的特點所決定的。
從圖2-16中小波分析得到的d11曲線可以看出,這一段地層可以劃分為兩個大的旋迴,對應兩個層序,譜系圖上的特徵也比較明顯。其中每個大的旋迴又可以分為三個次一級的旋迴,在d9及譜系圖上可以找到相關界面,相當於每個層序劃分出三個准層序組,每個准層序組在測井曲線及錄井資料上也有較明顯的反旋迴特徵。在進行准層組的劃分時,小波分析方法可以劃分出肉眼不易識別的旋迴,從而提高了劃分精度。整段地層一共可以劃分為21個准層序,代表不同的沉積旋迴。其旋迴特徵在d7曲線上有較好體現,從譜系圖上也可以找到各界面的標志。從測井曲線和岩性上看,基本上每一個准層序都是一個反旋迴,代表著一期的水體變換,這也完全符合層序地層學的基本原理。
圖2-16 牛100井小波分析資料的層序地層劃分
王62井位於牛庄窪陷東部,與牛100井相比,劃分出的各層序單元的厚度發生了明顯的變化,但數目基本一致,這也證明了小波分析劃分層序地層的結果是比較可靠的。通過對AC曲線的小波變換得到AC曲線的小波變換譜系圖和小波變換曲線,如圖2-17所示。從譜系圖和d11曲線上可以將整段地層劃分為兩個大的旋迴,分別對應層序Ⅲ和層序Ⅳ。其中每個層序又可以劃分為3個准層序組,在d9曲線上可以看到相應的旋迴出現,譜系圖上可以找到界面的標志(圖2-17)。王62井這一段地層一共可以劃分成20個准層序,缺失第一個准層序。各准層序在岩石類型、顏色和測井曲線上基本上可以看出反旋迴特徵,符合層序地層劃分方法。
通過牛100井、王62井的劃分可以看出,小波分析方法在砂泥岩互層的地層中有較好的應用效果,可以提高層序劃分的精度和准確性。在層序劃分中有比較好的可重復性,使得全區的劃分結果比較客觀和統一,減少了人為因素造成的干擾。
㈣ 求助matlab中的小波工具包
打開matlab軟體,進入軟體主界面在軟體的左下方找到start按鈕,點擊選擇toolbox,然後選擇wavelet進入wavemenu界面,選擇一維小波中的wavelet1-D並進入 7.選擇要處理的信號,界面出現loaded信號,這就是沒有去噪前的原9.分析後在左邊欄目中出現s,a*,d*,其中s為原信號,a*為近似信號,d*為細節信號11.閾值方法常用的有4種fixed(固定閾值),rigorsure,heusure,minmax根據需要選擇,一般情況下rigorsure方式去噪效果較好12.oft(軟閾值),hard(硬閾值)一般選擇軟閾值去噪後的信號較為平滑13.在雜訊結構中選擇unscaled white noise,因為在工程應用中的雜訊一般不僅僅含有白雜訊14.在雜訊結構下面的數值不要隨意改,這是系統默認的去噪幅度16.在此窗口下點擊file-save denoised singal,保存輸出去噪後的信號18.去噪結束後,把去噪後信號(.mat格式)拖至matlab主界面的workspace中,與原信號一起打包,以便以後計算統計量不會用就查幫助文檔啊!waverec函數是不需要你自己加零延拓的,上面的代碼完全不知所謂,waverec函數的使用是要依賴wavedec函數得到的CL組構的,CL組構中存放小波系數的數組C本身就已經延拓了,而且你不知道它對數據延拓了多少,延拓的方式有多種根本不是你這樣直接加零就行的。我發現你很有才,經常提問和編出一些匪夷所思的問題和代碼,不恥下問的精神是好的,但我個人是很不提倡這種做法的,有時間在這打字提問,不如找幾本基礎參考書看看,不了解就查嗎,不明白就往明白搞嗎,但看你這些「新奇」的問題和代碼真很抓狂,自己對於這些基本問題都懶得琢磨,打著勤奮好學,不恥下問的幌子,太沒勁了!哦,看錯了,waverec函數是可以用上面的代碼的,我看成wrcoef函數了,wrcoef函數可以實現waverec、upwlev和upcoef三個函數的功能之和,所以比waverec函數應用簡單,不需要你將其他分量置零,用它實現小波工具箱功能最方便。
㈤ matlab小波分析工具箱時信號載入出現了問題,請高手看下,謝謝~
這是多變數,也就是多個一維信號同時小波分析,如果是分析一個一維信號,直接使用工具箱的第一個按鈕即可
㈥ matlab中小波分析工具箱中wrcoef和waverec的區別是什麼
這兩個函數是與wavedec這個函數緊密相連的。
wavedec就是小波分解,將一個信號分解成指定回層數n,並返回答各層的小波系數。
waverec——它的作用與wavedec相反,即將給定的小波系數一次性完全重建出信號。
wrcoef——這個也是輸入小波系數,重建信號。但是它與上面有些區別,區別在於它重建的是原信號在指定層次的,高頻或者低頻分量。也就是說,這個信號不是原本的信號,而且某個層次上的逼近。
㈦ 如何使用matlab中的工具箱
如果是系統自帶的,你可以直接用,如果是外部的或者是自編的你需要先把文件夾拷貝到tools文件夾下,再設置路徑。x0dx0aMatlab常用工具箱介紹(英漢對照)x0dx0aMatlab Main Toolbox——matlab主工具箱x0dx0aControl System Toolbox——控制系統工具箱x0dx0aCommunication Toolbox——通訊工具箱x0dx0aFinancial Toolbox——財政金融工具箱x0dx0aSystem Identification Toolbox——系統辨識工具箱x0dx0aFuzzy Logic Toolbox——模糊邏輯工具箱x0dx0aHigher-Order Spectral Analysis Toolbox——高階譜分析工具箱x0dx0aImage Processing Toolbox——圖象處理工具箱x0dx0aLMI Control Toolbox——線性矩陣不等式工具箱x0dx0aModel predictive Control Toolbox——模型預測控制工具箱x0dx0aμ-Analysis and Synthesis Toolbox——μ分析工具箱x0dx0aNeural Network Toolbox——神經網路工具箱x0dx0aOptimization Toolbox——優化工具箱x0dx0aPartial Differential Toolbox——偏微分方程工具箱x0dx0aRobust Control Toolbox——魯棒控制工具箱x0dx0aSignal Processing Toolbox——信號處理工具箱x0dx0aSpline Toolbox——樣條工具箱x0dx0aStatistics Toolbox——統計工具箱x0dx0aSymbolic Math Toolbox——符號數學工具箱x0dx0aSimulink Toolbox——動態模擬工具箱x0dx0aSystem Identification Toolbox——系統辨識工具箱x0dx0aWavele Toolbox——小波工具箱x0dx0ax0dx0a例如:控制系統工具箱包含如下功能:x0dx0a連續系統設計和離散系統設計x0dx0a狀態空間和傳遞函數以及模型轉換x0dx0a時域響應(脈沖響應、階躍響應、斜坡響應)x0dx0a頻域響應(Bode圖、Nyquist圖)x0dx0a根軌跡、極點配置x0dx0ax0dx0a較為常見的matlab控制箱有:x0dx0ax0dx0a控制類:x0dx0ax0dx0a控制系統工具箱(control systems toolbox)x0dx0a系統識別工具箱(system identification toolbox)x0dx0a魯棒控制工具箱(robust control toolbox)x0dx0a神經網路工具箱(neural network toolbox)x0dx0a頻域系統識別工具箱(frequency domain system identification toolbox)x0dx0a模型預測控制工具箱(model predictive control toolbox)x0dx0a多變數頻率設計工具箱(multivariable frequency design toolbox)x0dx0ax0dx0a信號處理類:x0dx0a信號處理工具箱(signal processing toolbox)x0dx0a濾波器設計工具箱(filter design toolbox)x0dx0a通信工具箱(communication toolbox)x0dx0a小波分析工具箱(wavelet toolbox)x0dx0a高階譜分析工具箱(higher order spectral analysis toolbox)x0dx0ax0dx0a其它工具箱:x0dx0a統計工具箱(statistics toolbox)x0dx0a數學符號工具箱(symbolic math toolbox)x0dx0a定點工具箱(fixed-point toolbox)x0dx0a射頻工具箱(RF toolbox)x0dx0ax0dx0a1990年,MathWorks軟體公司為Matlab提供了新的控制系統模型化圖形輸入與模擬工具,並命名為Simulab,使得模擬軟體進入了模型化圖形組態階段,1992年正式命名為Simulink,即simu(模擬)和link(連接)。matlab7.0里的simulink為6.0版本,matlab6.5里的simulink為5.0版本。x0dx0ax0dx0aMATLAB的SIMULINK子庫是一個建模、分析各種物理和數學系統的軟體,它用框圖表示系統的各個環節,用帶方向的連線表示各環節的輸入輸出關系。x0dx0a啟動SIMULINK十分容易,只需在MATLAB的命令窗口鍵入「SIMULINK」命令,此時出現一個SIMULINK窗口,包含七個模型庫,分別是信號源庫、輸出庫、離散系統庫、線性系統庫、非線性系統庫及擴展系統庫。x0dx0a1.信號源庫x0dx0a包括階躍信號、正弦波、白雜訊、時鍾、常值、文件、信號發生器等各種信號源,其中信號發生器可產生正弦波、方波、鋸齒波、隨機信號等波形。x0dx0a2.輸出庫x0dx0a包括示波器模擬窗口、MATLAB工作區、文件等形式的輸出。x0dx0a3.離散系統庫x0dx0a包括五種標准模式:延遲,零-極點,濾波器,離散傳遞函數,離散狀態空間。x0dx0a4.線性系統庫x0dx0a提供七種標准模式:加法器、比例環節、積分環節、微分環節、傳遞函數、零-極點、狀態空間。x0dx0a5.非線性系統庫x0dx0a提供十三種常用標准模式:絕對值、乘法、函數、回環特性、死區特性、斜率、繼電器特性、飽和特性、開關特性等。x0dx0a6.系統連接庫包括輸入、輸出、多路轉換等模塊,用於連接其他模塊。x0dx0a7.系統擴展庫x0dx0a考慮到系統的復雜性,SIMULINK另提供十二種類型的擴展系統庫,每一種又有多種模型供選擇。x0dx0a使用時只要從各子庫中取出模型,定義好模型參數,將各模型連接起來,然後設置系統參數,如模擬時間、模擬步長、計算方法等。SIMULINK提供了Euler、RungeKutta、Gear、Adams及專用於線性系統的LinSim演算法,用戶根據模擬要求選擇適當的演算法。x0dx0ax0dx0a當然,不同版本的Matlab/Simulink內容有所不同。x0dx0ax0dx0a另外,Simulink還提供了諸如航空航天、CDMA、DSP、機械、電力系統等專業模塊庫,給快速建模提供了很大的便利。
㈧ matlab中的小波工具箱怎麼用,希望能詳細介紹
將原始數據文件夾到裝有matlab的電腦
打開matlab軟體,進入軟體主界面
在軟體的左下方找到start按鈕,點擊選擇toolbox,然後選擇wavelet
進入wavemenu界面,選擇一維小波中的wavelet1-D並進入
5.將數據文件(.Mat格式)托到matlab軟體主界面的workspace
6.在wavemenu主界面中選擇file-load signal或者import from workspace—import signal
7.選擇要處理的信號,界面出現loaded信號,這就是沒有去噪前的原
始信號
8.右上角選擇用於小波分析的小波基以及分解層數並點擊analyse開始分析
9.分析後在左邊欄目中出現s,a*,d*,其中s為原信號,a*為近似信號,d*為細節信號
10.然後點擊denoise去噪
11.閾值方法常用的有4種fixed(固定閾值),rigorsure,heusure,minmax根據需要選擇,一般情況下rigorsure方式去噪效果較好
12.oft(軟閾值),hard(硬閾值)一般選擇軟閾值去噪後的信號較為平滑
13.在雜訊結構中選擇unscaled white noise,因為在工程應用中的雜訊一般不僅僅含有白雜訊
14.在雜訊結構下面的數值不要隨意改,這是系統默認的去噪幅度
15.點擊denoise開始正式去噪
16.在此窗口下點擊file-save denoised singal,保存輸出去噪後的信號
17.去噪結束
18.去噪結束後,把去噪後信號(.mat格式)拖至matlab主界面的workspace中,與原信號一起打包,以便以後計算統計量
19.Matlab編程計算相關統計量以及特徵量
20.得出統計量和特徵量後結束
㈨ MATLAB 小波分析
[c l]=wavedec(x,5,'db6'); %%%x是要分析的數據%%%
d1=wrcoef('d',c,l,'db6',1);
d2=wrcoef('d',c,l,'db6',2);
d3=wrcoef('d',c,l,'db6',3);
d4=wrcoef('d',c,l,'db6',4);
d5=wrcoef('d',c,l,'db6',5);
a5=wrcoef('a',c,l,'db6',5);
plot(d1); %%%輸出d1的圖像
合並的話,直接用x1=a5+d5+d4+d3就可以吧