function?[wused]?=?BCS_fast（PHItsigma2eta）
[NM]?=?size（PHI）;
PHIt?=?PHI‘*t;
PHI2?=?sum（PHI.^2）‘;
ratio?=?（PHIt.^2）./PHI2;
[maxrindex]?=?max（ratio）;
alpha?=?PHI2（index）/（maxr-sigma2）;
%?compute?initial?mu?Sig?S?Q
phi?=?PHI（:index）;
Hessian?=?alpha?+?phi‘*phi/sigma2;
Sig?=?1/Hessian;
mu?=?Sig*PHIt（index）/sigma2;
left?=?PHI‘*phi/sigma2;
S?=?PHI2/sigma2-Sig*left.^2;
Q?=?PHIt/sigma2-Sig*PHIt（index）/sigma2*left;
%
for?count?=?1:1000

????s?=?S;?q?=?Q;
????s（index）?=?alpha.*S（index）./（alpha-S（index））;
????q（index）?=?alpha.*Q（index）./（alpha-S（index））;
????theta?=?q.^2-s;

????%?choice?the?next?alpha?that?maximizes?marginal?likelihood
????ml?=?-inf*ones（1M）;
????ig0?=?find（theta>0）;%尋找theta>0的索引
????%?index?for?re-estimate
????[irefoowhich]?=?intersect（ig0index）;
????if?~isempty（ire）
????????Alpha?=?s（ire）.^2./theta（ire）;
????????delta?=?（alpha（which）-Alpha）./（Alpha.*alpha（which））;
????????ml（ire）?=?Q（ire）.^2.*delta./（S（ire）.*delta+1）-log（1+S（ire）.*delta）;
????end
????%?index?for?adding
????iad?=?setdiff（ig0ire）;%尋找theta>0但未被選中的索引
????if?~isempty（iad）
????????ml（iad）?=?（Q（iad）.^2-S（iad））./S（iad）+log（S（iad）./（Q（iad）.^2））;
????end
????is0?=?setdiff（[1:M]ig0）;%尋找theta<0的索引
????%?index?for?deleting
????[idefoowhich]?=?intersect（is0index）;
????if?~isempty（ide）
????????ml（ide）?=?Q（ide）.^2./（S（ide）-alpha（which））-log（1-S（ide）./alpha（which））;
????end

????[ML（count）idx]?=?max（ml）;%最大化的位置
????%?check?if?terminates?
????if?count?>?2?&?abs（ML（count）-ML（count-1））?????????break;
????end

????%?update?alphas
????which?=?find（index==idx）;
????if?theta（idx）?>?0
????????if?~isempty（which）?%?re-estimate?不是空集則重估
????????????Alpha?=?s（idx）^2/theta（idx）;
????????????Sigii?=?Sig（whichwhich）;?mui?=?mu（which）;?Sigi?=?Sig（:which）;
????????????delta?=?Alpha-alpha（which）;
????????????ki?=?delta/（1+Sigii*delta）;
????????????mu?=?mu-ki*mui*Sigi;
????????????Sig?=?Sig-ki*Sigi*Sigi‘;
????????????comm?=?PHI‘*（phi*Sigi）/sigma2;
????????????S?=?S?+?ki*comm.^2;
????????????Q?=?Q?+?ki*mui*comm;
????????????%
????????????alpha（which）?=?Alpha;
????????else?%?adding
????????????Alpha?=?s（idx）^2/theta（idx）;
????????????phii?=?PHI（:idx）;?Sigii?=?1/（Alpha+S（idx））;?mui?=?Sigii*Q（idx）;
????????????comm1?=?Sig*（phi‘*phii）/sigma2;
????????????ei?=?phii-phi*comm1;
????????????off?=?-Sigii*comm1;
????????????Sig?=?[Sig+Sigii*comm1*comm1‘?off;?off‘?Sigii];
????????????mu?=?[mu-mui*comm1;?mui];
????????????comm2?=?PHI‘*ei/sigma2;
????????????S?=?S?-?Sigii*comm2.^2;
????????????Q?=?Q?-?mui*comm2;
????????????%
????????????index?=?[index;idx];
????????????alpha?=?[alpha;Alpha];
????????????phi?=?[phiphii];
????????end
????else
????????if?~isempty（which）?%?deleting
????????????Sigii?=?Sig（whichwhich）;?mui?=?mu（which）;?Sigi?=?Sig（:which）;
????????????Sig?=?Sig-Sigi*Sigi‘/Sigii;?Sig（:which）?=?[];?Sig（which:）?=?[];
????????????mu??=?mu-mui/Sigii*Sigi;?mu（which）?=?[];
??

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

?????文件???????5479??2009-12-12?13:22??米波一維兩球高頻\Algorithms\BCS_fast_rvm.m

?????文件???????1133??2009-12-12?13:22??米波一維兩球高頻\Algorithms\CoSaMP.m

?????文件??????10040??2017-03-11?16:50??米波一維兩球高頻\Algorithms\ExCoV.m

?????文件???????8090??2009-12-12?13:22??米波一維兩球高頻\Algorithms\ExCoVapp.m

?????文件???????7322??2009-12-29?09:36??米波一維兩球高頻\Algorithms\ExCoVls.m

?????文件??????21647??2009-12-12?13:22??米波一維兩球高頻\Algorithms\GPSR_BB.m

?????文件????????666??2018-05-08?16:56??米波一維兩球高頻\Algorithms\gui_fig51.m

?????文件??????14446??2009-12-12?13:54??米波一維兩球高頻\Algorithms\hard_l0_Mterm_mod.m

?????文件???????2763??2009-12-12?13:22??米波一維兩球高頻\Algorithms\l1qc_logbarrier.m

?????文件???????4179??2009-12-12?13:22??米波一維兩球高頻\Algorithms\l1qc_newton.m

?????文件???????1597??2009-12-13?10:55??米波一維兩球高頻\Algorithms\SBL.m

?????文件??????????0??2018-05-08?16:59??米波一維兩球高頻\Algorithms\Untitled6.m

?????文件???????3392??2017-06-30?11:13??米波一維兩球高頻\BCS_fast.m

?????文件????????676??2017-11-20?11:16??米波一維兩球高頻\cx1.m

?????文件????????664??2018-01-20?19:47??米波一維兩球高頻\fft_1d.m

?????文件????????538??2018-01-20?10:33??米波一維兩球高頻\fft_isar.m

?????文件??????36057??2018-01-20?09:51??米波一維兩球高頻\mibo1d_qiu_gaopin_imag.txt

?????文件??????36057??2018-01-20?09:51??米波一維兩球高頻\mibo1d_qiu_gaopin_real.txt

?????文件???????2646??2018-01-20?16:42??米波一維兩球高頻\OneD.m

?????文件????????490??2017-12-27?09:56??米波一維兩球高頻\Solve_OMP.m

?????文件???????4436??2018-05-04?08:57??米波一維兩球高頻\TwoD.m

?????目錄??????????0??2018-05-08?16:57??米波一維兩球高頻\Algorithms

?????目錄??????????0??2018-01-20?19:47??米波一維兩球高頻

-----------?---------??----------?-----??----

???????????????162318????????????????????23