function?[qXSNewInd]=AddOneElement（OmegaQXSYopt）
%AddOneElement?adds?one?entry?to?the?co-support?updates?the?signal?
%estimate?and?the?accumulated?orthogonal?set.
%??[qXSNewInd]=AddOneElement（OmegaQXSYopt）
%??========================================================================
%??Input:
%??Omega?-?analysis?dictionary.
%??Q?????-?matrix?of?the?orthogonal?set?accumulated?so?far?（in?its?rows）.
%??X?????-?the?current?estimate?for?the?signal.
%??S?????-?the?current?estimate?for?the?co-support.
%??Y?????-?the?observed?（possibly?noisy）?signal.
%??opt???-?the?type?of?optimality?criterion?used?for?adding?elements?to?the?
%??????????co-support:?0?-?THR?1?-?BG?2?-?OBG.?
%??========================================================================
%??Output:
%??q??????-?the?new?row?vector?added?to?the?orthogonal?set.
%??X??????-?the?new?estimate?for?the?signal.
%??S??????-?the?new?estimate?for?the?co-support.
%??NewInd?-?the?new?element?added?to?the?co-support.
%??========================================================================
%??Tomer?Peleg
%??Department?of?Electrical?Engineering
%??Technion?Haifa?32000?Israel
%??tomerfa@tx.technion.ac.il
%
%??October?2012
%??========================================================================
d=size（Omega2）;
inds=find（S==0）;
q=[];
NewInd=[];
if?~isempty（inds）
????Q_temp=Omega（inds:）;
????if?opt==2
????????if?~isempty（Q）
????????????%?modified?gram-schmidt?for?all?possible?directions
????????????for?j=1:size（Q1）
????????????????Q_temp=Q_temp-（Q_temp*Q（j:）‘）*Q（j:）;
????????????end
????????????Q_temp=Q_temp./repmat（sqrt（diag（Q_temp*Q_temp‘））1d）;
????????end
????end
????%?project?on?each?of?the?possible?directions
????if?opt==0
????????Proj=（Q_temp*Y）.^2;
????else
????????Proj=（Q_temp*X）.^2;
????end
????%?choose?the?direction?with?the?smallest?energy?reduction
????[dummypos]?=?min（Proj）;
????q=Q_temp（pos:）;
????if?opt<2
????????if?~isempty（Q）????????
????????????for?j=1:size（Q1）
????????????????q=q-（q*Q（j:）‘）*Q（j:）;
????????????end
????????????q=q/sqrt（q*q‘）;
????????end
????end
????if?nargout>1
????????X=X-q‘*（q*X）;
????end
????if?nargout>2
????????S=abs（Omega*X）<1e-4;
????end
????if?nargout>3
????????NewInd=inds（pos）;
????end
end

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----

?????文件?????185727??2001-07-12?13:02??images\barbara.png

?????文件?????177762??2001-07-12?12:42??images\boats.png

?????文件??????34985??2005-09-11?14:44??images\house.png

?????文件?????151199??2009-10-04?09:08??images\lena.png

?????文件??????40181??2002-08-29?16:48??images\peppers256.png

?????文件???????8601??2011-11-01?15:21??images\PWCImage.png

?????文件???????6294??2012-10-28?16:19??AnalysisKSVD.m

?????文件????????703??2012-10-28?15:54??readme.txt

?????文件???????3115??2012-10-28?15:54??Contents.m

?????文件???????3611??2012-10-28?15:20??TestKSVDNaturalImages.m

?????文件???????3405??2012-10-28?15:15??TestKSVDPWCImage.m

?????文件???????1645??2012-10-28?15:12??TestKSVDRandomOmega.m

?????文件???????2283??2012-10-28?15:11??TestKSVDOmegaDIF.m

?????文件???????4049??2012-10-28?15:11??TestPursuitOmegaDIF.m

?????文件???????2243??2012-10-28?15:01??GenerateAnalysisSignals.m

?????文件???????1346??2012-10-28?15:00??GenerateOmegaDIF.m

?????文件????????966??2012-10-28?14:59??GenerateRandomRow.m

?????文件???????1063??2012-10-28?14:59??ComputeOrthoSet.m

?????文件???????2304??2012-10-28?14:59??AddOneElement.m

?????文件???????2064??2012-10-28?14:58??ErrorBGP.m

?????文件???????1898??2012-10-28?14:58??RankBGP.m

?????文件???????1268??2012-10-28?14:58??DisplayOmega.m

?????文件???????2117??2012-10-27?09:19??dictdist.m

?????文件???????7680??2011-05-16?09:08??col2imstep.mexw32

?????文件???????7168??2011-05-16?09:08??im2colstep.mexw32

?????文件???????3924??2011-05-16?08:48??col2imstep.c

?????文件???????1001??2011-05-16?08:48??col2imstep.m

?????文件???????1178??2011-05-16?08:48??countcover.m

?????文件???????3388??2011-05-16?08:48??im2colstep.c

?????文件???????1341??2011-05-16?08:48??im2colstep.m

............此處省略16個文件信息