genbounddata1.m

来自「一种基于压缩感知技术的图像重建程序」· M 代码 · 共 48 行

% Generates data for controlled L^p norm example when:
% M constant, P constant, N varies

P = 0.25:0.05:1;
M = [1000 2000 4000];
Nfrac = 0.05:0.05:0.75;
numTrials = 20;
errArray = zeros(length(P), length(M), length(Nfrac), numTrials);
Pnorm = zeros(length(P), length(M), length(Nfrac), numTrials);

for pp = 1:length(P)
    for mm = 1:length(M)
        N = floor(Nfrac .* M(mm));
        for nn = 1:length(N)
            for jj = 1:numTrials
                % Generate Random Dictionary
                Phi = MatrixEnsemble(N(nn),M(mm));

                % Generate sparse signal
                alpha0 = GenRandVec(M(mm),P(pp));

                % generate the vector S
                S = Phi * alpha0;

                % Solve the BP problem
                alpha = SolveBP(Phi, S);

                errArray(pp,mm,nn,jj) = twonorm(alpha - alpha0);
                Pnorm(pp,mm,nn,jj) = pnorm(alpha0,P(pp));
            end
        end
    end
end

save BoundData1.mat M P Nfrac numTrials errArray Pnorm

%
% Copyright (c) 2006. Yaakov Tsaig
%  

%
% Part of SparseLab Version:100
% Created Tuesday March 28, 2006
% This is Copyrighted Material
% For Copying permissions see COPYING.m
% Comments? e-mail [e-mail]
%