【正文】 n a l西安文理學院本科畢業(yè)設計（論文） 第 22 頁 也能精確恢復原始信號，相對誤差小在 左右 ，而且具有很好的收斂性。 OMP 算法的 Matlab 語言實現(xiàn) 因為 OMP 算法只能對稀疏信號進行精確重構(gòu)，所以為了實現(xiàn)對 OMP 的仿真，首先需要將圖片轉(zhuǎn)換為稀疏信號，如 FFT、 DCT、小波變換等，在本文中采取小波西安文理學院本科畢業(yè)設計（論文） 第 18 頁 變換，然后將得到的稀疏信號與測量矩陣相乘得到測量值 Y， 然后通過 OMP 算法，根據(jù) Y 精確重建出原始信號最后經(jīng)過小波反變換得到原始圖像。 本文后面的仿真就是采用基于壓縮感知的正交匹配追蹤算法來重建圖像的。因而，在最小化驅(qū)動的壓縮感知理論完善工作的基礎上，希望能夠基于稀疏性自適應的貪婪迭代和基于多層超先驗建模的非凸迭代思想設計適于噪聲情形大尺度問題的快速魯棒重建算法。 (2)硬件易實現(xiàn)的確定性測量矩陣設計。壓縮感知理論的大部分文獻中，測量矩陣 ? 都是線性的且設計好的 ， 不需根據(jù)觀測信號自適應地變化。 ? 基礎理論層面： （ 1） 基于非 正交稀疏字典的壓縮感知信號重建理論。 CS 理論自從被 （美國科學學院院士） ( 創(chuàng)始人 )及 (華裔科學家， 20xx 年菲爾茨獲獎得主， 20xx 年被評為世界上最聰明的科學家 )等人提出后，在信息論、信號 /圖像處理、醫(yī)療成像、模式識別、地質(zhì)勘探、光學 /雷達成像、無線通信等領域受到 高度關注，并被美國科技評論為 20xx 年度 10 大科技進展之一 [14]。然而由于壓縮重構(gòu)算法的計算量比較大，難以達到實時性要求，因此實時高性能壓縮感知成像系統(tǒng)是未來重要的研究方向?；?ROMP算法和 SAMP算法的突出優(yōu)勢。 20xx年 Donoho等人提出了分段正交匹配追蹤 (STOMP， Stagewise OMP)算法。 Chen， Donoho和 Saunders指出，求解一個更加簡單的 1l 優(yōu)化問題會產(chǎn)生同等的解(要求 ?和 ? 不相關 )： 1min TX? . CS TA X X Y? ?? ? (式 ) 稍微的差別使得問題變成了一個凸優(yōu)化問題，于是可以方便地化簡為線性規(guī)劃問題，典型算法代表： BP算法．盡管 BP算法可行，但在實際應用中存在兩個問題：第一，即使是常見的圖像尺寸， BP算法的計算復雜度也難以忍受，在采樣點個數(shù)滿足 M cK? ， ? ?2log / 1c N K??時，重構(gòu)計算復雜度的量級在 3()ON ；第二，由于 1?范數(shù)無法區(qū)分稀疏系數(shù)尺度的位置，所以盡管整體上重構(gòu)信號在歐氏距離上逼近原信號，但存在低尺度能量搬移到了高尺度的現(xiàn)象，從而容易出現(xiàn)一些人工效應，如一維信號會在高頻出現(xiàn)振蕩。 信號重構(gòu) 如何設計快速重構(gòu)算法，從線性觀測 CSY A X? 中恢復信號，是第三步要將解決的問題，即信號的重構(gòu)問題。隨機高斯矩陣具有一個有用的性質(zhì)：對于一個 MN? 的隨機高斯矩陣 ? ，可以證明當 M≥ cKlog(N／ K)時 T CSA?? ? 在很大概率下具有 RIP性質(zhì) (其中 c是一個很小的常數(shù) )。從中可以看出，問題的關鍵是如何確定非零系數(shù)的位置來構(gòu)造出一個可解的 MK? 線性方程組。 顯然，如果觀測過程破壞了 X 中的信息，重構(gòu)是不可能的。 在構(gòu)造冗余字典方面，文獻 [7]中提出使用局部 Cosine基來刻畫聲音信號的局部頻域特性；利用 bandlet基來刻畫圖像中的幾何邊緣；還可以把其它的具有不同形狀的基函數(shù)歸入字典，如適合刻畫紋理的 Gabor基、適合刻畫輪廓的 Curvelet基等等。即在某個正交基字典里，自適應地尋找可以逼近某一種信號特征的最優(yōu)正交基，根據(jù)不同的信號尋找最適合信號特性的一個正交基，對信號進行變換以得到最稀疏的信號表示。 圖 壓縮感知理論 的信號重構(gòu)過程框圖。如果 ? 的非零個數(shù)比 N 小很多，則表明該信號是可壓縮的。 CS 理論利用到了許多自然信號在特定的基 ? 上具有緊湊的表示。此外在壓縮編碼過程中，大量變換計算得到的小系數(shù)被丟棄，造成了數(shù)據(jù)計算和內(nèi)存資源的浪費。 第三章： 正交匹配追蹤重建算法。自從 20xx 年 CS 的提出，在 IEEE 的信號處理匯刊、信號處理快報匯刊、信號處理雜志、信息論匯刊等國際知名期刊上開始涌現(xiàn)出上百篇關于 CS 理論與應用方面的文獻。 在 20xx 年 Donoho 等人針對稀疏性信號，提出了壓縮感知（ Compressive sensing,簡稱 CS）理論。 Sparse transform。該理論表明，用遠低于 Nyquist 采樣定理要求的頻率對信號進行采樣也能實現(xiàn)信號的精確重構(gòu)。s. This theory is under the condition that the signal is pressible or sparse. In this case, using far less than the required sampling frequency of the Nyquist theory to sample the signal is able to accurately reconstruct the theory breaks though the traditional Nyquist sampling theory, which overes a lot of problems such as a great number of sampling data, time wasting, data storage space wasting and so on. As a result, it reduces signal processing cost and device cost. The pressed theory has three key sides: (1) Sparse transformation, for a non sparse signal, we need to find a proper orthogonal basis on which the signal has a sparse representation。如何高效處理這些數(shù)據(jù)并且最大限度的節(jié)省存儲空間 及傳輸成本已成為目前信息領域進一步向前發(fā)展的主要瓶頸之一。信號的稀疏重建與壓縮感知理論有重大的實用價值和應用前景，已經(jīng)成為信號領域中一個新的研究方向 [1]。首先介紹了壓縮感知理論的研究背景及意義，然后介紹了國內(nèi)外研究背景和現(xiàn)狀，最后整理出全文內(nèi)容的結(jié)構(gòu)安排。 西安文理學院本科畢業(yè)設計（論文） 第 4 頁 第二章 壓縮感知理論相關知識 壓縮感知理論框架 傳統(tǒng)的信號采集、編解碼過程如圖 所示 。 西安文理學院本科畢業(yè)設計（論文） 第 5 頁 圖 縮感知理論的編解碼框圖 壓縮感知的基本理論及核心問題 壓縮感知，也被稱為壓縮傳感或壓縮采樣，是一種利用稀疏的或可壓縮的信號進行信號重構(gòu)的技術 [3]。為簡化問題，假定這些基是規(guī)范正交的。 Candes證明了只要信號在某一個正交空間具有稀疏性，就能以較低的頻率 ? ?MN?? 采樣信號，而且可以以高概率重構(gòu)該信號。 Candes和 Tao研究表明，滿足具有冪次速度衰減的信號，可利用壓縮感知理論得到恢復，并且重構(gòu)誤差滿足： 62? ( / l o g ) rrE X X C K N ?? ? ? ? (式 ) 其中 r=1/p– 1/2， 0p1. 文獻 [6]指出光滑信號的 Fourier系數(shù)、小 波系數(shù)、有界變差函數(shù)的全變差范數(shù)、振蕩信號的 Gabor系數(shù)及具有不連續(xù)邊緣的圖像信號的 Curvelet系數(shù)等都具有足夠的稀疏性，可以通過壓縮感知理論恢復信號。從冗余字典中找到具有最佳線性組合的 K 項原子來表示一個信號，稱作信號的稀疏逼近或高度非線性逼近。 西安文理學院本科畢業(yè)設計（論文） 第 8 頁 信號的觀測矩陣 如何設計一個平穩(wěn)的、與變換基 ? 不相關的 MN? 維的觀測矩陣 ? ，保證稀疏向量 ? 從 N 維降到 M 維時重要信息不遭破壞，是第二步要解決的問題，也就是信號低速采樣問題。對此，kk ff ?? ????? 112222 即有限等距性質(zhì)（ RIP）給出了存在確定解的充要條件。不相干性越強，互相表示時所需的系數(shù)越多；反之，相 關性則越強。 Donoho給出了觀測矩陣所必需具備的三 個條件 [9]，并指出大部分一致分布的隨機矩陣都具備這三個條件，均可作為觀測矩陣，如：部分 Fourier集、部分 Hadamard集、一致分布的隨機投影 (uniform Random Projection)集等，這與對 RIP性質(zhì)進行研究得出的結(jié)論相一致。因此，求解式（ )的數(shù)值計算極不穩(wěn)定而且是 NP難問題。該方法針對 BP、 MP和 OMP方法沒有考慮信號的多尺度分解時稀疏信號在各子帶位置的關系，是將稀疏系數(shù)的樹型結(jié)構(gòu)加以利用，進一步提升了重構(gòu)信號的精度和求解的速度。 ROMP算法與 OMP算法的不同之處在于，該算法首先根據(jù)相關原子挑選多個原子作為候選集，然后從候選集中按照正則化 原則挑選出部分原子，最后將其并入最終的支撐集，從而實現(xiàn)了原子的快速、有效選擇。該相機具有一種全新的相機結(jié)構(gòu)，使用數(shù)字微鏡陣列完成圖像在偽隨機二值模型上線性投影的光學計算。使其在測量舉證的作用下，用遠小于原圖像的數(shù)據(jù)量進行計算得到融合結(jié)果還原為圖像表示，可節(jié)省中間融合所需的計算量，并且能夠更好地利用原圖像中像素間的內(nèi)在聯(lián)系，是一個非常值得研究的課題?；?RIP理論，目前已研制出了一些設備，有萊斯大學研制的單像素相機和 A/I轉(zhuǎn)換器，麻省理工學院研制的編碼孔徑相機，耶魯大學研制的超譜成像儀，麻省理工學院研制的 MRI RF脈沖設備，伊利諾伊州立大學研制的 DNA微陣列傳感器。雖然基于線性投影的壓縮感知理論能夠直接應用于自 然圖像這樣的復雜高維信號，但是由于沒有考慮到自然圖像的固有特性，諸如結(jié)構(gòu)多成分性、高階統(tǒng)計性等 ,對于自然圖像壓縮采樣本身沒有特殊的指導作用。目前，基于構(gòu)造方法的自然圖像過完備字典設計具有很好的理論支撐，正則化幾何方法、幾何多尺度分析、基于信息論的 “有效編碼假設 ”為其奠定了堅實廣闊的理論基礎。 （ 3） 噪聲情形大尺度問題的快速魯棒重 建算法 [15]設計。經(jīng)過數(shù)次迭代，該信號便可以由一些原子線性表示。 OMP的重建算法是在給定迭代次數(shù)的條件下重建，這種強制迭代過程停止的方法使得 OMP需要非常多的線性測量來保證精確重建。 接下來，取 觀測向量 M 的長度為 64，即采樣率 M/N= 對其進行重建。原始信號進行了小波變換后，變得更加稀疏，高頻段的信息基本沒了，這樣圖像縮小了很多，便于保存與傳輸?；粉櫵惴ň褪腔谝陨显韽亩_到重構(gòu)和對白噪聲去噪的目的。仿真如下圖 ： O r i g i n e i m a g e N = 6 5 5 3 6 B P , s a m p = 3 9 . 8 4 3 8 % t i m e = 5 7 . 0 3 8 s e cO M P , s a m p = 3 9 . 8 4 3 8 % t i m e = 1 . 8 0 1 s e c FD R , s a m p = 3 9 . 8 4 3 8 % t i m e = 0 . 7 0 5 s e c 圖 原始圖像與 BP、 OMP、 STOMP_FDR 恢復對比圖 西安文理學院本科畢業(yè)設計（論文） 第 26 頁 圖 從圖像的重構(gòu)質(zhì)量上可以看出， OMP 的重建效果次于 BP 算法 但優(yōu)于STOMP_FDR 算法。西安文理學院本科畢業(yè)設計（論文） 第 27 頁 結(jié)束語 近年來，信號處理領域出現(xiàn)了一種新的信息采樣理論 —— 壓縮感知。因此需要研究基于最小 L1 范數(shù)算法的快速算法，做到重構(gòu)效果與時間發(fā)的統(tǒng)一。正交變換選取的合適與否直接關系到信號是否滿足稀疏度的要求，從而會影響到信號是否能被精確重建，因此選取合適的正交變換基也是重建算法需要考慮的一個重要部分。我只有在今后的學習、工作中，以鍥而不舍的精神，努力做出點成績，以 不負恩師的厚望 。 % 信號長度 M=64。 % 采樣序列 x=*sin(2*pi*f1*Ts*ts)+*sin(2*pi*f2*Ts*ts)+*sin(2*pi*f3*Ts*ts)+*sin(2*pi*f4*Ts*ts)。 % 信號頻率 1 f2=100。 % 測量矩陣 (高斯分布白噪聲 ) s=Phi*x.39。 % 迭代次數(shù) (有噪聲的情況下 ,該迭代次數(shù)為 K) for col=1:N。*s。 plot(hat_x,39。) norm(hat_x.39。) % 重建信號 plot(x,39。 % 殘差 西安文理學院本科畢業(yè)設計（論文） 第 32 頁 pos_array(times)=pos。*r_n)。 % 算法迭代次數(shù) (m=K) Psi=fft(eye(N,N))/sqrt(N)。 % 信號頻率 3 f4=400。r39。 % 信號頻率 1 f2=10