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);