MATLAB实现稀疏表示程序及Gabor原子分解

f=MVP';                                %输入原始信号
k=5;                                   %设置迭代次数 
[a,N]=size(f);                  
signal_reconstruct=zeros(1,N);         %重建信号初始化
signal_r=f;                            %残余信号初始化
%设置原子库的参数
a_base=2;                             %伸缩因子定为2
j_min=0;                              %离散化指数j的最小值
j_max=log2(N);                        %离散化指数j的最大值
u_base=1/2;                           %位移因子定为1/2
p_min=0;                              %p的最小值为0，最大值在函数里定义
v_base=pi;                            %频率定为π
k_min=0;                              %k的最小值为0，最大值在函数里定义
w_base=pi/6;                          %相位定为π/6
i_min=0;                              %i的最小值为0
i_max=12;                             %i的最大值为12
m=0;                                  %原子个数初始值为0 
%将残余信号投影到最佳原子上，使残余信号最小化
tic
ab=zeros(1,k)
for n=1:k
[proj,scale,translation,freq,phase]=gabor_atom(signal_r,N,a_base,j_min,j_max,u_base,p_min,v_base,k_min,w_base,i_min,i_max);
%主函数调用，函数详见附录
t=0:N-1;
t=(t-translation)/scale;
g=(1/sqrt(scale))*exp(-pi*t.*t).*cos(freq*t+phase);            %原子
g=g/sqrt(sum(g.*g));                                    %使得原子范数保持1
signal_reconstruct=signal_reconstruct+proj*g;                %MP算法重构信号
signal_r=signal_r-proj*g;                                 %计算残余信号
m=m+1;