淺談壓縮感知(十七):測量矩陣之有限等距常數RIC的計算


有限等距常數(RestrictedIsometry Constant, RIC)是與有限等距性質(Restricted IsometryProperty, RIP)緊密結合在一起的一個參數。

一、RIC定義:

在前面的一篇RIP文章中其實已經提到了,這里在貼出該定義:

CS滿足條件:S階RIP性質只要要求0<δS<1就可以了,而RIC是指滿足RIP的最小δS。

二、RIC計算:

RIC與特征值的關系:

三、MATLAB代碼

function y=RIPText(x,t)
% calculete the kth Restricted Isometry Constant of measurement matrix x.
% Modified on March 13th, 2012.
% Reference: Wei Dai, Olgica Milenkovic. Subspace Pursuit for Compressive Sensing
% Signal Reconstruction. IEEE Trans. ona Information Theory. vol.55, no.5,
% 2230-2249. May, 2009
%calculate the Restricted Isomentry Constant of matrix constructed by
%random n columns of the original matrix x.
%created by Li Bo in March 16th, 2011
%Harbin Institute of Technolgy
n=size(x,2);%the numbers of column of original matrix x
Num=1:n;%create a row of numbers from one to n with iterval 1
Com=combntns(Num,t);% all the combination of the selected t columns from total n columns
ma=size(Com,1);% the number of combinations
delta=zeros(1,ma);% a vector used to store the restricted isometry constant candidate
for i=1:ma
    b=x(:,Com(i,:));% new matrix constructed by the t selected columns 
    e=eig(b'*b-eye(t));%calculate all the eigen values of matrix b'*b-eye(t)
    delta(i)=(max(abs(e)));%select the largest absolute eigen value as restricted isomentry candidate
end
y=max(delta);% select the largest one as restricted isometry constant

用此代碼來驗證矩陣其實是不現實的,原因解釋如下,在程序運行過程中容易導致內存不足等問題。

四、參考文章

http://blog.csdn.net/jbb0523/article/details/44622819


免責聲明!

本站轉載的文章為個人學習借鑒使用,本站對版權不負任何法律責任。如果侵犯了您的隱私權益,請聯系本站郵箱yoyou2525@163.com刪除。



 
粵ICP備18138465號   © 2018-2025 CODEPRJ.COM