分块压缩感知模型

function Phiw=aweight（Phiblock_sizeM）

N=block_size;
q=N;
for i=1:N
    for j=1:N
        w（ij）=1+1/（1+（i+j）/q）;
    end 
end
wn=w（:）;
Nb=N*N;

Phiw=zeros（MNb）;
for m=1:M
    for n=1:Nb
        Phiw（mn）=Phi（mn）*wn（n）;
    end
end