基于壓縮感知算法數(shù)據(jù)壓縮與復(fù)原代碼美賽思路即將發(fā)布科研交流程序注釋

1、一部分%-% initializations%-n=length(x); %m=length(measvec); % xrecon=zeros(1,n); % max likelihood % xreconstruction mrecon=zeros(1,n); % average % srecon=zeros(1,n); % std_dev % sreconpdf_N_to_M=zeros(n, r, model_order); % N to M , pdf_M_to_N=zeros(m, l, model_order); % M to N xx=1:model_order; xx=xx-(

2、model_order+1)/2; xx=xx/max(xx); % xx=xx*boundx; % values over which pdf is sampled % pdf_xx=xx;delta=xx(2)-xx(1); % 其中，xx是抽樣的序列符號，即：抽多少個樣、抽樣間隔是由xx決定的。二部分%-% compute signal node prior%-pdf_prior= (k/n)*normpdf(xx,0,sig_1); % if (sig_0 epsilon) % epsilon pdf_prior= pdf_prior + (1-k/n)*normpdf(xx,0,si

3、g_0);else in2=find(abs(xx)1) % damping pdf_res=gmean(pdf_res,pdf(i,:),gamma_pdbp,0); %gamma_pdbp=0: true BP % end mtt, stt, maxtt=meanvarmaxpdf(pdf_res), xx); % computes statistics xrecon(i)=maxtt; % max likelihood mrecon(i)=mtt; % mean srecon(i)=stt; % std_dev around mean pdf(i,:)=pdf_res; % store

4、result其實(shí)，此部分就對應(yīng)于公式（5）。值得注意的是函數(shù)meanvarmaxpdf，如何進(jìn)行MMSE估計的。Meanvarmaxpdf函數(shù)的公式如下（所用的變量表達(dá)式全部與meanvarmaxpdf.m中一致）：1、歸一化pdf：pdf=pdf/sum(pdf);2、計算pdf中最大值mm和其所在位置ind。3、求期望m=xx.*pdf. (因?yàn)閤x為橫軸上的值，pdf為其對應(yīng)的概率，所以乘積是期望)4、求xx2的期望e2=（xx.*xx.*pdf）5、求方差v2=e2-m*m（即對應(yīng)公式）6、求標(biāo)準(zhǔn)差sig=以上所求的各個變量均用來估計xx和判斷迭代終止條件。六部分這部分對應(yīng)于迭代終止判

5、斷，沒有相應(yīng)公式可以寫，只是判斷終止條件是否滿足， % END OF ESTIMATE OF SIGNAL for j=1:ln % to send next message pdf_tmp=reshape(pdf_M_to_N(neighbors(j),self_index(j),:),1, model_order); pdf_tmp=pdf_tmp+epsilon; % for stability pdf_tosend=divpdf(pdf_res,pdf_tmp); pdf_tosend=pdf_tosend/sum(pdf_tosend); tmptmp=pdf_tosend; if

6、(it1) % MDBPF prevpdf=reshape(pdf_N_to_M(i, j,:),1,model_order); tmptmp=gmean(pdf_tosend,prevpdf,gamma_mdbpf,0); end pdf_N_to_M(i, j, :)=tmptmp/sum(tmptmp); end %of for j end end if (it =2) % display and break on last iteration dispvec_anderrors(mrecon, phi, phisign, measvec, dispind, x); if (it=ite

7、r) break; end end其中，函數(shù)divpdf.m也很簡單，就是返回pdf_res除以pdf_tmp的值。七部分下面是程序的最后一部分，即信息從變量節(jié)點(diǎn)傳到測量節(jié)點(diǎn)。若為第一次迭代，則程序從初始化pdf_N_to_M(i, rr, :)=pdf_prior(:);后就進(jìn)步本循環(huán)。 %- % BACKWARD ITERATION - from signal to measurements %- for i=1:m neighbors=phi(i,:); neighbors=setdiff_shri(neighbors,0); ln=length(neighbors); phisigni

8、=phisign(i,1:ln); self_index=self_indexM(i,:); pdf_res_all=; for jj=1:ln % process neighbors tmptmp=reshape(pdf_N_to_M(neighbors(jj),self_index(jj),:),1, model_order); if (phisigni(jj) epsilon) pdf_res_all=convpdf_fft(pdf_res_all, y_noise, epsilon); end for j=1:ln %To send next message tmptmp=reshap

9、e(pdf_N_to_M(neighbors(j),self_index(j),:),1, model_order); mvi=measvec(i); if (phisigni(j) 0) & (b2) tmptmp=reverse(tmptmp,model_order); end pdf_res=unconvpdf_fft(pdf_res_all, tmptmp, epsilon); pdf_res=abs(ifft(pdf_res); pdf_res=shiftpdf_fft(pdf_res, mvi, delta, model_order); pdf_res=pdf_res/sum(pdf_res); if (phisigni(j) 1) % MDBPF prevpdf=reshape(pdf_M_