x1 = FB.*FR; FBR = BlockMM(nr,nc,Nb,m,x1); invW = BlockMM(nr,nc,Nb,m,F2B); invWBR = FBR./(invW + tau*Nb); fun = @(block_struct) block_struct.data.*invWBR; FCBinvWBR = blockproc(FBC,[nr,nc],fun); FX = (FR-FCBinvWBR)/tau; Xest = real(ifft2(FX));
(FB, FBC, F2B, FR, tau, Nb, nr, nc, m)