aboutsummaryrefslogtreecommitdiff
path: root/Measurements/At_fhp.m
diff options
context:
space:
mode:
authorEugeniy E. Mikhailov <evgmik@gmail.com>2021-01-29 16:23:05 -0500
committerEugeniy E. Mikhailov <evgmik@gmail.com>2021-01-29 16:23:05 -0500
commit3983eb46023c1edd00617729ba929057fda8d0bd (patch)
tree816ad084f355000656c43da9160f1c257bbb1ddc /Measurements/At_fhp.m
downloadl1magic-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.m28
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);