diff options
author | Eugeniy E. Mikhailov <evgmik@gmail.com> | 2021-01-29 16:23:05 -0500 |
---|---|---|
committer | Eugeniy E. Mikhailov <evgmik@gmail.com> | 2021-01-29 16:23:05 -0500 |
commit | 3983eb46023c1edd00617729ba929057fda8d0bd (patch) | |
tree | 816ad084f355000656c43da9160f1c257bbb1ddc /Measurements/At_fhp.m | |
download | l1magic-3983eb46023c1edd00617729ba929057fda8d0bd.tar.gz l1magic-3983eb46023c1edd00617729ba929057fda8d0bd.zip |
Initial import from https://statweb.stanford.edu/~candes/software/l1magic/v1.11
Additional Clean up of Mac dirs and tex generated files
Diffstat (limited to 'Measurements/At_fhp.m')
-rw-r--r-- | Measurements/At_fhp.m | 28 |
1 files changed, 28 insertions, 0 deletions
diff --git a/Measurements/At_fhp.m b/Measurements/At_fhp.m new file mode 100644 index 0000000..a7879af --- /dev/null +++ b/Measurements/At_fhp.m @@ -0,0 +1,28 @@ +% At_fhp.m +% +% Adjoint of At_fhp (2D Fourier half plane measurements). +% +% Usage: x = At_fhp(b, OMEGA, n) +% +% b - K vector = [mean; real part(OMEGA); imag part(OMEGA)] +% +% OMEGA - K/2-1 vector denoting which Fourier coefficients to use +% (the real and imag parts of each freq are kept). +% +% n - Image is nxn pixels +% +% x - N vector +% +% Written by: Justin Romberg, Caltech +% Created: October 2005 +% Email: jrom@acm.caltech.edu +% + +function x = At_fhp(y, OMEGA, n) + +K = length(y); + +fx = zeros(n,n); +fx(1,1) = y(1); +fx(OMEGA) = sqrt(2)*(y(2:(K+1)/2) + i*y((K+3)/2:K)); +x = reshape(real(n*ifft2(fx)), n*n, 1); |