基于压缩感知算法数据压缩与复原代码美赛思路即将发布科研交流程序注释

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其实，此部分就对应于公式（5）。值得注意的是函数meanvarmaxpdf，如何进行MMSE估计的。Meanvarmaxpdf函数的公式如下（所用的变量表达式全部与meanvarmaxpdf.m中一致）：1、归一化pdf：pdf=pdf/sum(pdf);2、计算pdf中最大值mm和其所在位置ind。3、求期望m=xx.*pdf. (因为xx为横轴上的值，pdf为其对应的概率，所以乘积是期望)4、求xx2的期望e2=（xx.*xx.*pdf）5、求方差v2=e2-m*m（即对应公式）6、求标准差sig=以上所求的各个变量均用来估计xx和判断迭代终止条件。六部分这部分对应于迭代终止判

5、断，没有相应公式可以写，只是判断终止条件是否满足， % 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其中，函数divpdf.m也很简单，就是返回pdf_res除以pdf_tmp的值。七部分下面是程序的最后一部分，即信息从变量节点传到测量节点。若为第一次迭代，则程序从初始化pdf_N_to_M(i, rr, :)=pdf_prior(:);后就进步本循环。 %- % 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_