% [D, x] = random_bp_problem(d, nAtoms)
%
% Generate a random signal representation problem.
% The signal x is d-dimensional.
% The dictionary D is d-dimensional and contains nAtoms atoms.
%
function [D, x] = random_bp_problem(d, nAtoms)

D = zeros(nAtoms, d);
for i = 1 : nAtoms,
    D(i,:) = randus(d);
end
x = randus(d);

%%%%%%%%%%%%%%%%%%%%%%%
function v = randus(d) 
   v = randn(1,d);
   r = sqrt(sum(v.^2));
   v = v/r;