The Scaled Gradient Projection Method: An Application to Nonconvex Optimization
|
|
- Julie Young
- 5 years ago
- Views:
Transcription
1 2332 PIERS Proceedings, Prague, Czech Republic, July 6 9, 2015 The Scaled Gradient Projection Method: An Application to Nonconvex Optimization M. Prato 1, A. La Camera 2, S. Bonettini 3, and M. Bertero 2 1 Dipartimento di Scienze Fisiche, Inormatiche e Matematiche Università di Modena e Reggio Emilia, Via Campi 213/b, Modena 41125, Italy 2 Dipartimento di Inormatica, Bioingegneria, Robotica e Ingegneria dei Sistemi Università di Genova, Via Dodecaneso 35, Genova 16145, Italy 3 Dipartimento di Matematica e Inormatica Università di Ferrara, Via Saragat 1, Ferrara 44122, Italy Abstract The scaled gradient projection (SGP) method is a variable metric orward-bacward algorithm designed or constrained dierentiable optimization problems, as those obtained by reormulating several signal and image processing problems according to standard statistical approaches. The main SGP eatures are a variable scaling matrix multiplying the gradient direction at each iteration and an adaptive steplength parameter chosen by generalizing the well-nown Barzilai-Borwein rules. An interesting result is that SGP can be exploited within an alternating minimization approach in order to address optimization problems in which the unnown can be splitted in several blocs, each with a given convex and closed easible set. Classical examples o applications belonging to this class are the non-negative matrix actorization and the blind deconvolution problems. In this wor we applied this method to the blind deconvolution o multiple images o the same target obtained with dierent PSFs. In particular, or our experiments we considered the NASA unded Fizeau intererometer LBTI o the Large Binocular Telescope, which is already operating on Mount Graham and has provided the irst Fizeau images, demonstrating the possibility o reaching the resolution o a 22.8 m telescope. Due to the Poisson nature o the noise aecting the measured images, the resulting optimization problem consists in the minimization o the sum o several Kullbac-Leibler divergences, constrained in suitable easible sets accounting or the dierent eatures to be preserved in the object and the PSFs. 1. INTRODUCTION Intererometry is a standard way exploited in astronomical imaging to obtain high angular resolution starting rom two or more telescopes with smaller diameters. Famous examples o astronomical intererometers are the Very Large Telescope Intererometer (VLTI), the Navy Prototype Optical Intererometer (NPOI) and the Center or High Angular Resolution Astronomy (CHARA) array. A urther Fizeau intererometer which recently provided its irst images is that located at the Large Binocular Telescope (LBT, Mount Graham, Arizona) and called Large Binocular Telescope Intererometer (LBTI [1]). With its two 8.4-m primary mirrors, separated by a 14.4-m center-tocenter distance, and mounted on a common alt-az mount, LBT can be considered the very irst Extremely Large Telescope (ELT). The adaptive secondary mirrors provide high-order adaptive optics corrections which produce images with Strehl Ratio (SR), i.e., the ratio o pea diraction intensity o an aberrated versus perect waveorm (see, e.g., [9]), greater than 0.9. LBTI collects the light rom both primary mirrors and combines the two beams on the image plane (Fizeau imaging). In order to obtain a uniorm coverage o the requency plane, several LBTI images at dierent rotation angles are needed. Each o these measured images will be aected by several sources o noise, whose statistical nature includes both Poisson (due to, e.g., radiation rom the object, bacground, dar current) and Gaussian (read-out-noise (RON) due to the ampliier) components. It has been shown [14] that the latter one can be handled by adding the variance o the RON to the measured images and the bacground and considering the resulting data as purely Poissonian. As a result o this, a statistical approach to the image reconstruction problem leads to the minimization o the sum o a given number p (equal to the number o images acquired by LBTI) o Kullbac- Leibler divergences min KL(, ω 1,..., ω p ) Ω ω 1,...,ωp Ωω p n 2 ( (g ) i (g ) i ln (ω + b ) i =1 ) + (ω + b ) i (g ) i, (1)
2 Progress In Electromagnetics Research Symposium Proceedings 2333 where: is the n n image to be recovered; g is the -th n n measured image, and b the corresponding n n bacground radiation; ω is the unnown n n point spread unction (PSF) related to g, which has to be estimated during the reconstruction procedure; Ω and Ω ω are the easible sets o and ω, accounting or physical constraints that object and PSFs have to satisy. Since we are assuming that the PSFs are unnown, problem (1) is a so-called blind deconvolution problem [7], which is highly ill-posed due to the ininite solutions available. Moreover, rom the optimization point o view the objective unction in (1) is nonconvex, thus leading to the presence o multiple stationary points. A possible way to address this ind o problems is to introduce suitable constraints on the unnown to reduce the set o easible solutions. In particular, in [11, 12] the ollowing easible sets have been considered { } { } n 2 Ω = R n2 0, i = 1 p n 2 n 2 (g b ) i, Ω ω = ω R n2 0 ω s, ω i = 1, p =1 where the upper bound s is estimated rom the nowledge o the SR. Moreover, in [12] an alternating minimization approach or the solution o (1) has been proposed based on the scaled gradient projection (SGP) method [3]. The aim o this paper is to exploit the same optimization approach to clash with a urther deect which characterizes intererometric images. As nown, the PSF can be described as the product o the Airy unction (representing the PSF o a single mirror) by the intererometric term (1 + cos θ), where θ depends on the current pixel and the orientation o the baseline (besides the other parameters o the acquisition system) 1. I we consider a (small) dierence in the two optical paths, i.e., there is a residual error when the system corrects or the so-called piston error, then the intererometric term becomes (1 + cos(θ + φ)). This phase error φ moves the position o the ringes in the direction o the baseline (and thereore in the direction orthogonal to the ringes themselves) with the result that the PSF is not symmetric with respect to the original PSF. The corresponding image is aected by the same eect. 2. THE CYCLIC BLOCK COORDINATE GRADIENT PROJECTION METHOD The cyclic bloc coordinate gradient projection (CBCGP) method is an alternating minimization algorithm originally proposed in [2] in the context o non-negative matrix actorization [6] and adapted in ollowing wors to astronomical blind deconvolution problems [4, 11, 12]. The main idea at the basis o this approach is a partial minimization over each bloc o variables in a Gauss-Seidel style perormed inexactly by means o the scaled gradient projection (SGP) method [3]. The SGP algorithm is a irst-order method which applies to a general constrained optimization problem min ϕ(x). (2) x Ω The main steps o SGP are reported in Algorithm 1, where D is the set o diagonal positive deinite matrices with diagonal elements bounded in a positive interval. I ϕ is equal to the Kullbac-Leibler divergence and the PSFs are normalized to unit volume, a requently used scaling matrix is the one borrowed rom the Richardson-Lucy method [8, 13] and deined as { 1 { [D ] ii = max µ, min µ, x ()}}, with µ > 1 (see, e.g., [3, 5, 10]). As concerns the choice o the steplength parameter α, we adopted a suitable alternation o the Barzilai and Borwein rules (see [3, 10] or urther details). The CBCGP scheme applied to problem (1) is reported in Algorithm 2. 1 For more details: θ(n 1, n 2 ) = A(n 1 cos α + n 2 sin α), where α is the hour angle o the baseline and A = 2π B/λ, being the pixel size o the image in rad/px, B the distance between the two mirrors and λ the observed wavelength.
3 2334 PIERS Proceedings, Prague, Czech Republic, July 6 9, 2015 Algorithm 2: Scaled gradient projection (SGP) method Choose the starting point x (0) Ω, set the parameters β, δ (0, 1), 0 < α min < α max. For = 0, 1, 2,... do the ollowing steps: Step 1. Choose the steplength parameter α [α min, α max ] and the scaling matrix D D; Step 2. Projection: y () = P Ω,D 1(x () α D ϕ(x () )); Step 3. Descent direction: d () = y () x () ; Step 4. Set λ = 1; Step 5. Bactracing loop: I ϕ(x () + λ d () ) ϕ(x () ) + βλ ϕ(x () ) T d () then go to Step 6; Else set λ = δλ and go to Step 5. i Step 6. Set x (+1) = x () + λ d (). Algorithm 3: Cyclic bloc coordinate gradient projection (CBCGP) method Choose the starting points (0) Ω, ω (0) 1,..., ω(0) p Ω ω, and two integers N, N ω 1. For h = 0, 1, 2,... do the ollowing steps: Step. Choose an integer 1 N (h) Algorithm 2 to problem starting rom the point (h). Step ω. For = 1,..., p do the ollowing steps: Choose an integer 1 N ω (h) Algorithm 2 to problem starting rom the point ω (h). N and compute (h+1) by applying N (h) iterations o min KL(, ω (h) 1,..., Ω ω(h) p ) (3) N ω and compute ω (h+1) by applying N (h) ω iterations o min KL( (h+1), ω (h+1) 1,..., ω (h+1) ω Ω 1, ω, ω(h) +1,..., ω(h) p ) (4) ω 3. NUMERICAL EXPERIMENTS We simulated several diraction limited, SR = 1, LBTI PSFs changing the values o the mentioned phase error. Each set o PSFs is a cube o three PSFs with baseline angles at 0, 60, and 120 and corresponds to a dierent value o φ (ranging rom 0, 0.2,..., 1.0). An horizontal cut o the 0 PSFs is shown in the let panel o Fig. 1. By choosing a model o star cluster (based on the nine brightest stars o the Pleiades, as described in [10, 11]), we used these PSFs or generating several multiple images. We adopted a pixel-size o 11 mas and we supposed to acquire images in M-band (λ = 4.8 µm). We applied our blind algorithm by using as initialization the autocorrelation o the ideal diraction limited PSF (i.e., the PSF with no phase error). We remar that this choice has a SR o about 0.35, as described in our previous papers, with positions o the ringes that do not correspond with those in the images. Thereore our goal is to veriy whether the algorithm is able to reconstruct both the object and the PSFs and what is the maximum phase error which gives satisying results. Following [11, 12], we chose N (h) = 50 and N (h) ω = 1 or all h,, and we perormed 1000 outer
4 Progress In Electromagnetics Research Symposium Proceedings 2335 (a) (b) (c) Figure 1: (a) The horizontal cut o the simulated PSFs or each o the phase error φ and or the baseline angle o 0. A positive value moves the position o the ringes towards let. (b) The input image corresponding to φ = 1.0 or the baseline angle o 0. (c) The 0 PSF used or the initialization. The images are shown in log scale with a zoom actor o 2. Table 1: Astrometric error (AE), magnitude average relative error (MARE) and root mean square error (RMSE) or dierent phase errors φ. For the three PSFs, the values in bracets denote the errors ater re-centering. The normalized value o the objective unction is also shown. φ AE (pixels) MARE RMSE 0 RMSE 60 RMSE 120 2/(pn 2 )KL 0.0 < % 0.7% 0.7% < % 0.6% 0.6% % (4.5%) 28.0% (16.5%) 26.6% (16.2%) % (15.9% ) 45.2% (18.3%) 20.0% (11.8%) % (20.0%) 64.1% (9.3%) 39.1% (9.2%) % (18.8%) 64.5% (9.1%) 39.0% (9.1%) (a) (b) Figure 2: (a) The horizontal cut o the reconstructed PSFs or each o the phase error φ. (b) The six reconstructed PSFs corresponding to φ = 0, 0.2,..., 1.0 (in order rom let to right, rom top to bottom) and or the baseline angle o 0. The images are shown in log scale with a zoom actor o 4. iterations. In order to evaluate the quality o the reconstructions we computed three errors: the astrometric error (AE), the magnitude average relative error (MARE), and the root mean square
5 2336 PIERS Proceedings, Prague, Czech Republic, July 6 9, 2015 error (RMSE). The three errors are deined as: AE = 1 q q (xi x i ) 2 + (y i ỹ i ) 2 ; MARE = 1 q q m i m i m i ; RMSE = ω ω 2 ω 2 where q is the number o stars, (x i, y i ) and ( x i, ỹ i ) are the reconstructed centroid and the original position o each star, m i and m i are the reconstructed and the true magnitudes, ω and ω are the reconstructed and true PSFs and 2 denotes the Euclidean norm. All the errors are reported in Table 1. As concerns RMSE, it is interesting to note that the values are quite high or all cases with φ > 0.2. Indeed, since each star is reconstructed with a small, but signiicant, astrometric error, this causes a small displacement o the reconstructed PSF with respect to the original one used in the computation o the RMSE. Thereore we re-centered each PSF by means o a sub-pixel-precision procedure and we computed again the RMSE (the new values are shown in parenthesis). O course, in case o real images this operation is not necessary. In conclusion, the proposed strategy seems to provide good reconstructions o the PSFs only i the measured images are aected by a low phase error, although the magnitude o the stars are in general well recovered. Future wor will include the addition o explicit regularization terms in the objective unction speciically aimed at preserving the eatures which characterize the LBTI PSFs. ACKNOWLEDGMENT The present paper has been supported by the Italian Ministry or University and Research (projects FIRB Futuro in Ricerca 2012, contract RBFR12M3AC, and PRIN 2012, contract 2012MTE38N) and the Italian Gruppo Nazionale per il Calcolo Scientiico (GNCS). REFERENCES 1. Bailey, V. P., P. M. Hinz, PA. T. uglisi, et al., Large binocular telescope intererometer adaptive optics: On-sy perormance and lessons learned, Proc. SPIE, Vol. 9148, , Bonettini, S., Inexact bloc coordinate descent methods with application to the nonnegative matrix actorization, IMA J. Numer. Anal., Vol. 31, No. 4, , Bonettini, S., R. Zanella, and L. Zanni, A scaled gradient projection method or constrained image deblurring, Inverse Probl., Vol. 25, No. 1, , Bonettini, S., A. Cornelio, and M. Prato, A new semiblind deconvolution approach or Fourierbased image restoration: An application in astronomy, SIAM J. Imaging Sci., Vol. 6, No. 3, , Bonettini, S. and M. Prato, Accelerated gradient methods or the X-ray imaging o solar lares, Inverse Probl., Vol. 30, No. 5, , Lee, D. D. and H. S. Seung, Learning the part o objects rom non-negative matrix actorization, Nature, Vol. 401, No. 6755, , Levin, A., Y. Weiss, F. Durand, and W. T. Freeman, Understanding and evaluating blind deconvolution algorithms, IEEE Conerence on Computer Vision and Pattern Recognition, 2009, CVPR 2009, , Paris, France, Jun Lucy, L. B., An iterative technique or the rectiication o observed distributions, Astronom. J., Vol. 79, No. 6, , Mahajan, V. N., Strehl ratio or primary aberrations in terms o their aberration variance, J. Opt. Soc. Am., Vol. 73, No. 6, , Prato, M., R. Cavicchioli, L. Zanni, P. Boccacci, and M. Bertero, Eicient deconvolution methods or astronomical imaging: Algorithms and IDL-GPU codes, Astron. Astrophys., Vol. 539, A133, Prato, M., A. La Camera, S. Bonettini, and M. Bertero, A convergent blind deconvolution method or post-adaptive-optics astronomical imaging, Inverse Probl., Vol. 29, No. 6, , Prato, M., A. La Camera, S. Bonettini, S. Rebegoldi, M. Bertero, and P. Boccacci, A blind deconvolution method or ground-based telescopes, New Astron., Vol. 40, 1 13, Richardson, W. H., Bayesian based iterative method o image restoration, J. Opt. Soc. Amer., Vol. 62, No. 1, 55 59, Snyder, D. L., A. M. Hammoud, and R. L. White, Image recovery rom data acquired with a charge-coupled-device camera, J. Opt. Soc. Am. A, Vol. 10, No. 5, , 1993.
Accelerated Gradient Methods for Constrained Image Deblurring
Accelerated Gradient Methods for Constrained Image Deblurring S Bonettini 1, R Zanella 2, L Zanni 2, M Bertero 3 1 Dipartimento di Matematica, Università di Ferrara, Via Saragat 1, Building B, I-44100
More informationScaled gradient projection methods in image deblurring and denoising
Scaled gradient projection methods in image deblurring and denoising Mario Bertero 1 Patrizia Boccacci 1 Silvia Bonettini 2 Riccardo Zanella 3 Luca Zanni 3 1 Dipartmento di Matematica, Università di Genova
More informationarxiv: v3 [math.na] 26 Feb 2015
arxiv:1406.6601v3 [math.na] 26 Feb 2015 New convergence results for the scaled gradient projection method S Bonettini 1 and M Prato 2 1 Dipartimento di Matematica e Informatica, Università di Ferrara,
More informationANALYSIS OF p-norm REGULARIZED SUBPROBLEM MINIMIZATION FOR SPARSE PHOTON-LIMITED IMAGE RECOVERY
ANALYSIS OF p-norm REGULARIZED SUBPROBLEM MINIMIZATION FOR SPARSE PHOTON-LIMITED IMAGE RECOVERY Aramayis Orkusyan, Lasith Adhikari, Joanna Valenzuela, and Roummel F. Marcia Department o Mathematics, Caliornia
More informationPhase Estimation in Differential-Interference-Contrast (DIC) Microscopy
Phase Estimation in Differential-Interference-Contrast (DIC) Microscopy Lola Bautista, Simone Rebegoldi, Laure Blanc-Féraud, Marco Prato, Luca Zanni, Arturo Plata To cite this version: Lola Bautista, Simone
More informationRegularized Multiplicative Algorithms for Nonnegative Matrix Factorization
Regularized Multiplicative Algorithms for Nonnegative Matrix Factorization Christine De Mol (joint work with Loïc Lecharlier) Université Libre de Bruxelles Dept Math. and ECARES MAHI 2013 Workshop Methodological
More informationNumerical Methods for Parameter Estimation in Poisson Data Inversion
DOI 10.1007/s10851-014-0553-9 Numerical Methods for Parameter Estimation in Poisson Data Inversion Luca Zanni Alessandro Benfenati Mario Bertero Valeria Ruggiero Received: 28 June 2014 / Accepted: 11 December
More informationIn many diverse fields physical data is collected or analysed as Fourier components.
1. Fourier Methods In many diverse ields physical data is collected or analysed as Fourier components. In this section we briely discuss the mathematics o Fourier series and Fourier transorms. 1. Fourier
More informationX-ray Diffraction. Interaction of Waves Reciprocal Lattice and Diffraction X-ray Scattering by Atoms The Integrated Intensity
X-ray Diraction Interaction o Waves Reciprocal Lattice and Diraction X-ray Scattering by Atoms The Integrated Intensity Basic Principles o Interaction o Waves Periodic waves characteristic: Frequency :
More informationProbabilistic Model of Error in Fixed-Point Arithmetic Gaussian Pyramid
Probabilistic Model o Error in Fixed-Point Arithmetic Gaussian Pyramid Antoine Méler John A. Ruiz-Hernandez James L. Crowley INRIA Grenoble - Rhône-Alpes 655 avenue de l Europe 38 334 Saint Ismier Cedex
More informationarxiv: v1 [math.na] 1 Jun 2015
VARIABLE METRIC INEXACT LINE SEARCH BASED METHODS FOR NONSMOOTH OPTIMIZATION S. BONETTINI, I. LORIS, F. PORTA, AND M. PRATO arxiv:1506.00385v1 [math.na] 1 Jun 2015 Abstract. We develop a new proximal gradient
More informationRestoration of interferometric images. IV. An algorithm for super-resolution of stellar systems
A&A 431, 747 755 (2005) DOI: 10.1051/0004-6361:20040366 c ESO 2005 Astronomy & Astrophysics Restoration of interferometric images IV. An algorithm for super-resolution of stellar systems B. Anconelli 1,
More informationTFY4102 Exam Fall 2015
FY40 Eam Fall 05 Short answer (4 points each) ) Bernoulli's equation relating luid low and pressure is based on a) conservation o momentum b) conservation o energy c) conservation o mass along the low
More informationFeasibility of a Multi-Pass Thomson Scattering System with Confocal Spherical Mirrors
Plasma and Fusion Research: Letters Volume 5, 044 200) Feasibility o a Multi-Pass Thomson Scattering System with Conocal Spherical Mirrors Junichi HIRATSUKA, Akira EJIRI, Yuichi TAKASE and Takashi YAMAGUCHI
More informationNumerical Methods - Lecture 2. Numerical Methods. Lecture 2. Analysis of errors in numerical methods
Numerical Methods - Lecture 1 Numerical Methods Lecture. Analysis o errors in numerical methods Numerical Methods - Lecture Why represent numbers in loating point ormat? Eample 1. How a number 56.78 can
More information1. Interference condition. 2. Dispersion A B. As shown in Figure 1, the path difference between interfering rays AB and A B is a(sin
asic equations or astronomical spectroscopy with a diraction grating Jeremy Allington-Smith, University o Durham, 3 Feb 000 (Copyright Jeremy Allington-Smith, 000). Intererence condition As shown in Figure,
More informationSIO 211B, Rudnick. We start with a definition of the Fourier transform! ĝ f of a time series! ( )
SIO B, Rudnick! XVIII.Wavelets The goal o a wavelet transorm is a description o a time series that is both requency and time selective. The wavelet transorm can be contrasted with the well-known and very
More informationITERATED SRINKAGE ALGORITHM FOR BASIS PURSUIT MINIMIZATION
ITERATED SRINKAGE ALGORITHM FOR BASIS PURSUIT MINIMIZATION Michael Elad The Computer Science Department The Technion Israel Institute o technology Haia 3000, Israel * SIAM Conerence on Imaging Science
More informationAnalytical expressions for field astigmatism in decentered two mirror telescopes and application to the collimation of the ESO VLT
ASTRONOMY & ASTROPHYSICS MAY II 000, PAGE 57 SUPPLEMENT SERIES Astron. Astrophys. Suppl. Ser. 44, 57 67 000) Analytical expressions or ield astigmatism in decentered two mirror telescopes and application
More informationSolving Continuous Linear Least-Squares Problems by Iterated Projection
Solving Continuous Linear Least-Squares Problems by Iterated Projection by Ral Juengling Department o Computer Science, Portland State University PO Box 75 Portland, OR 977 USA Email: juenglin@cs.pdx.edu
More informationarxiv: v1 [gr-qc] 18 Feb 2009 Detecting the Cosmological Stochastic Background of Gravitational Waves with FastICA
1 arxiv:0902.3144v1 [gr-qc] 18 Feb 2009 Detecting the Cosmological Stochastic Background o Gravitational Waves with FastICA Luca Izzo 1,2,3, Salvatore Capozziello 1 and MariaFelicia De Laurentis 1 1 Dipartimento
More informationAn IDL Based Image Deconvolution Software Package
An IDL Based Image Deconvolution Software Package F. Városi and W. B. Landsman Hughes STX Co., Code 685, NASA/GSFC, Greenbelt, MD 20771 Abstract. Using the Interactive Data Language (IDL), we have implemented
More informationPoint Spread Function of Symmetrical Optical System Apodised with Gaussian Filter
International Journal o Pure and Applied Physics. ISSN 973-776 Volume 4, Number (8), pp. 3-38 Research India Publications http://www.ripublication.com Point Spread Function o Symmetrical Optical System
More informationPole searching algorithm for Wide-field all-sky image analyzing monitoring system
Contrib. Astron. Obs. Skalnaté Pleso 47, 220 225, (2017) Pole searching algorithm for Wide-field all-sky image analyzing monitoring system J. Bednář, P. Skala and P. Páta Czech Technical University in
More informationPulling by Pushing, Slip with Infinite Friction, and Perfectly Rough Surfaces
1993 IEEE International Conerence on Robotics and Automation Pulling by Pushing, Slip with Ininite Friction, and Perectly Rough Suraces Kevin M. Lynch Matthew T. Mason The Robotics Institute and School
More informationScattered Data Approximation of Noisy Data via Iterated Moving Least Squares
Scattered Data Approximation o Noisy Data via Iterated Moving Least Squares Gregory E. Fasshauer and Jack G. Zhang Abstract. In this paper we ocus on two methods or multivariate approximation problems
More informationA Simple Explanation of the Sobolev Gradient Method
A Simple Explanation o the Sobolev Gradient Method R. J. Renka July 3, 2006 Abstract We have observed that the term Sobolev gradient is used more oten than it is understood. Also, the term is oten used
More informationVirtually everything you and I know about the Cosmos has been discovered via electromagnetic
Gravitational Waves Gravitational wave astronomy Virtually everything you and I know about the Cosmos has been discovered via electromagnetic observations. Some inormation has recently been gleaned rom
More informationSteepest descent method implementation on unconstrained optimization problem using C++ program
IOP Conerence Series: Materials Science and Engineering PAPER OPEN ACCESS Steepest descent method implementation on unconstrained optimization problem using C++ program o cite this article: H Napitupulu
More informationCHAPTER 8 ANALYSIS OF AVERAGE SQUARED DIFFERENCE SURFACES
CAPTER 8 ANALYSS O AVERAGE SQUARED DERENCE SURACES n Chapters 5, 6, and 7, the Spectral it algorithm was used to estimate both scatterer size and total attenuation rom the backscattered waveorms by minimizing
More informationarxiv: v1 [hep-lat] 26 Dec 2014
arxiv:42.7928v [hep-lat] 26 Dec 204 A perturbative sty o the chirally rotated Schrödinger Functional in QCD Stean Sint School o Mathematics, Trinity College, Dublin 2, Ireland E-mail: sint@maths.tcd.ie
More informationSPOC: An Innovative Beamforming Method
SPOC: An Innovative Beamorming Method Benjamin Shapo General Dynamics Ann Arbor, MI ben.shapo@gd-ais.com Roy Bethel The MITRE Corporation McLean, VA rbethel@mitre.org ABSTRACT The purpose o a radar or
More informationA NEW METHOD FOR ADAPTIVE OPTICS POINT SPREAD FUNCTION RECONSTRUCTION
Florence, Italy. May 2013 ISBN: 978-88-908876-0-4 DOI: 10.12839/AO4ELT3.13328 A NEW METHOD FOR ADAPTIVE OPTICS POINT SPREAD FUNCTION RECONSTRUCTION J. Exposito 1a, D. Gratadour 1, G. Rousset 1, Y. Clénet
More informationDealing with edge effects in least-squares image deconvolution problems
Astronomy & Astrophysics manuscript no. bc May 11, 05 (DOI: will be inserted by hand later) Dealing with edge effects in least-squares image deconvolution problems R. Vio 1 J. Bardsley 2, M. Donatelli
More informationReview of Prerequisite Skills for Unit # 2 (Derivatives) U2L2: Sec.2.1 The Derivative Function
UL1: Review o Prerequisite Skills or Unit # (Derivatives) Working with the properties o exponents Simpliying radical expressions Finding the slopes o parallel and perpendicular lines Simpliying rational
More informationRelating axial motion of optical elements to focal shift
Relating aial motion o optical elements to ocal shit Katie Schwertz and J. H. Burge College o Optical Sciences, University o Arizona, Tucson AZ 857, USA katie.schwertz@gmail.com ABSTRACT In this paper,
More informationSensor-fault Tolerance using Robust MPC with Set-based State Estimation and Active Fault Isolation
Sensor-ault Tolerance using Robust MPC with Set-based State Estimation and Active Fault Isolation Feng Xu, Sorin Olaru, Senior Member IEEE, Vicenç Puig, Carlos Ocampo-Martinez, Senior Member IEEE and Silviu-Iulian
More informationCOMP 408/508. Computer Vision Fall 2017 PCA for Recognition
COMP 408/508 Computer Vision Fall 07 PCA or Recognition Recall: Color Gradient by PCA v λ ( G G, ) x x x R R v, v : eigenvectors o D D with v ^v (, ) x x λ, λ : eigenvalues o D D with λ >λ v λ ( B B, )
More informationTwo-step self-tuning phase-shifting interferometry
Two-step sel-tuning phase-shiting intererometry J. Vargas, 1,* J. Antonio Quiroga, T. Belenguer, 1 M. Servín, 3 J. C. Estrada 3 1 Laboratorio de Instrumentación Espacial, Instituto Nacional de Técnica
More informationINPUT GROUND MOTION SELECTION FOR XIAOWAN HIGH ARCH DAM
3 th World Conerence on Earthquake Engineering Vancouver, B.C., Canada August -6, 24 Paper No. 2633 INPUT GROUND MOTION LECTION FOR XIAOWAN HIGH ARCH DAM CHEN HOUQUN, LI MIN 2, ZHANG BAIYAN 3 SUMMARY In
More informationOn the Girth of (3,L) Quasi-Cyclic LDPC Codes based on Complete Protographs
On the Girth o (3,L) Quasi-Cyclic LDPC Codes based on Complete Protographs Sudarsan V S Ranganathan, Dariush Divsalar and Richard D Wesel Department o Electrical Engineering, University o Caliornia, Los
More informationbetween electron energy levels. Using the spectrum of electron energy (1) and the law of energy-momentum conservation for photon absorption:
ENERGY MEASUREMENT O RELATIVISTIC ELECTRON BEAMS USING RESONANCE ABSORPTION O LASER LIGHT BY ELECTRONS IN A MAGNETIC IELD R.A. Melikian Yerevan Physics Institute, Yerevan ABSTRACT The possibility o a precise
More informationRelating axial motion of optical elements to focal shift
Relating aial motion o optical elements to ocal shit Katie Schwertz and J. H. Burge College o Optical Sciences, University o Arizona, Tucson AZ 857, USA katie.schwertz@gmail.com ABSTRACT In this paper,
More informationarxiv: v1 [cs.it] 12 Mar 2014
COMPRESSIVE SIGNAL PROCESSING WITH CIRCULANT SENSING MATRICES Diego Valsesia Enrico Magli Politecnico di Torino (Italy) Dipartimento di Elettronica e Telecomunicazioni arxiv:403.2835v [cs.it] 2 Mar 204
More informationGaussian Process Regression Models for Predicting Stock Trends
Gaussian Process Regression Models or Predicting Stock Trends M. Todd Farrell Andrew Correa December 5, 7 Introduction Historical stock price data is a massive amount o time-series data with little-to-no
More informationStrain and Stress Measurements with a Two-Dimensional Detector
Copyright ISSN (C) 97-, JCPDS-International Advances in X-ray Centre Analysis, or Volume Diraction 4 Data 999 5 Strain and Stress Measurements with a Two-Dimensional Detector Baoping Bob He and Kingsley
More informationEstimation and detection of a periodic signal
Estimation and detection o a periodic signal Daniel Aronsson, Erik Björnemo, Mathias Johansson Signals and Systems Group, Uppsala University, Sweden, e-mail: Daniel.Aronsson,Erik.Bjornemo,Mathias.Johansson}@Angstrom.uu.se
More informationAn efficient approximation method for the nonhomogeneous backward heat conduction problems
International Journal o Engineering and Applied Sciences IJEAS) ISSN: 394-3661 Volume-5 Issue-3 March 018 An eicient appromation method or the nonhomogeneous bacward heat conduction problems Jin Wen Xiuen
More informationYURI LEVIN AND ADI BEN-ISRAEL
Pp. 1447-1457 in Progress in Analysis, Vol. Heinrich G W Begehr. Robert P Gilbert and Man Wah Wong, Editors, World Scientiic, Singapore, 003, ISBN 981-38-967-9 AN INVERSE-FREE DIRECTIONAL NEWTON METHOD
More informationAPPENDIX 1 ERROR ESTIMATION
1 APPENDIX 1 ERROR ESTIMATION Measurements are always subject to some uncertainties no matter how modern and expensive equipment is used or how careully the measurements are perormed These uncertainties
More informationTLT-5200/5206 COMMUNICATION THEORY, Exercise 3, Fall TLT-5200/5206 COMMUNICATION THEORY, Exercise 3, Fall Problem 1.
TLT-5/56 COMMUNICATION THEORY, Exercise 3, Fall Problem. The "random walk" was modelled as a random sequence [ n] where W[i] are binary i.i.d. random variables with P[W[i] = s] = p (orward step with probability
More information3. Several Random Variables
. Several Random Variables. Two Random Variables. Conditional Probabilit--Revisited. Statistical Independence.4 Correlation between Random Variables. Densit unction o the Sum o Two Random Variables. Probabilit
More informationA Limited Memory, Quasi-Newton Preconditioner. for Nonnegatively Constrained Image. Reconstruction
A Limited Memory, Quasi-Newton Preconditioner for Nonnegatively Constrained Image Reconstruction Johnathan M. Bardsley Department of Mathematical Sciences, The University of Montana, Missoula, MT 59812-864
More informationAnalytical analysis of second-order Stokes wave in Brillouin ring fiber laser
University o Malaya From the SelectedWorks o Faisal Raiq Mahamd Adikan December, 2011 Analytical analysis o second-order Stokes wave in Brillouin ring iber laser Faisal Raiq Mahamd Adikan, University o
More informationRATIONAL FUNCTIONS. Finding Asymptotes..347 The Domain Finding Intercepts Graphing Rational Functions
RATIONAL FUNCTIONS Finding Asymptotes..347 The Domain....350 Finding Intercepts.....35 Graphing Rational Functions... 35 345 Objectives The ollowing is a list o objectives or this section o the workbook.
More informationOn the regularization properties of some spectral gradient methods
On the regularization properties of some spectral gradient methods Daniela di Serafino Department of Mathematics and Physics, Second University of Naples daniela.diserafino@unina2.it contributions from
More information9.1 The Square Root Function
Section 9.1 The Square Root Function 869 9.1 The Square Root Function In this section we turn our attention to the square root unction, the unction deined b the equation () =. (1) We begin the section
More informationRandom Error Analysis of Inertial Sensors output Based on Allan Variance Shaochen Li1, a, Xiaojing Du2,b and Junyi Zhai3,c
International Conerence on Civil, Transportation and Environment (ICCTE 06) Random Error Analysis o Inertial Sensors output Based on Allan Variance Shaochen Li, a, Xiaojing Du, and Junyi Zhai3,c School
More informationCHAPTER 1: INTRODUCTION. 1.1 Inverse Theory: What It Is and What It Does
Geosciences 567: CHAPTER (RR/GZ) CHAPTER : INTRODUCTION Inverse Theory: What It Is and What It Does Inverse theory, at least as I choose to deine it, is the ine art o estimating model parameters rom data
More informationMaximum Flow. Reading: CLRS Chapter 26. CSE 6331 Algorithms Steve Lai
Maximum Flow Reading: CLRS Chapter 26. CSE 6331 Algorithms Steve Lai Flow Network A low network G ( V, E) is a directed graph with a source node sv, a sink node tv, a capacity unction c. Each edge ( u,
More informationThe deconvolution of lunar brightness temperature based on maximum entropy method using Chang E-2 microwave data
Research in Astron. Astrophys. 201x Vol. X No. XX, 000 000 http://www.raa-journal.org http://www.iop.org/journals/raa Research in Astronomy and Astrophysics The deconvolution of lunar brightness temperature
More informationMSE405 Microstructure Characterization XRD-1 Lab X-ray diffraction in crystallography
X-ray diraction in crystallography I. Goals Crystallography is the science that studies the structure (and structure-derived properties) o crystals. Among its many tools, X-ray diraction (XRD) has had
More informationCISE-301: Numerical Methods Topic 1:
CISE-3: Numerical Methods Topic : Introduction to Numerical Methods and Taylor Series Lectures -4: KFUPM Term 9 Section 8 CISE3_Topic KFUPM - T9 - Section 8 Lecture Introduction to Numerical Methods What
More informationTESTING TIMED FINITE STATE MACHINES WITH GUARANTEED FAULT COVERAGE
TESTING TIMED FINITE STATE MACHINES WITH GUARANTEED FAULT COVERAGE Khaled El-Fakih 1, Nina Yevtushenko 2 *, Hacene Fouchal 3 1 American University o Sharjah, PO Box 26666, UAE kelakih@aus.edu 2 Tomsk State
More informationUNCERTAINTY EVALUATION OF SINUSOIDAL FORCE MEASUREMENT
XXI IMEKO World Congress Measurement in Research and Industry August 30 eptember 4, 05, Prague, Czech Republic UNCERTAINTY EVALUATION OF INUOIDAL FORCE MEAUREMENT Christian chlegel, Gabriela Kiekenap,Rol
More informationClicker questions. Clicker question 2. Clicker Question 1. Clicker question 2. Clicker question 1. the answers are in the lower right corner
licker questions the answers are in the lower right corner question wave on a string goes rom a thin string to a thick string. What picture best represents the wave some time ater hitting the boundary?
More informationExtreme Values of Functions
Extreme Values o Functions When we are using mathematics to model the physical world in which we live, we oten express observed physical quantities in terms o variables. Then, unctions are used to describe
More informationUltra Fast Calculation of Temperature Profiles of VLSI ICs in Thermal Packages Considering Parameter Variations
Ultra Fast Calculation o Temperature Proiles o VLSI ICs in Thermal Packages Considering Parameter Variations Je-Hyoung Park, Virginia Martín Hériz, Ali Shakouri, and Sung-Mo Kang Dept. o Electrical Engineering,
More informationSparsity-Aware Pseudo Affine Projection Algorithm for Active Noise Control
Sparsity-Aware Pseudo Aine Projection Algorithm or Active Noise Control Felix Albu *, Amelia Gully and Rodrigo de Lamare * Department o Electronics, Valahia University o argoviste, argoviste, Romania E-mail:
More informationA Brief Survey on Semi-supervised Learning with Graph Regularization
000 00 002 003 004 005 006 007 008 009 00 0 02 03 04 05 06 07 08 09 020 02 022 023 024 025 026 027 028 029 030 03 032 033 034 035 036 037 038 039 040 04 042 043 044 045 046 047 048 049 050 05 052 053 A
More informationDouble-slit interference of biphotons generated in spontaneous parametric downconversion from a thick crystal
INSTITUTE OF PHYSICS PUBLISHING JOURNAL OF OPTICS B: QUANTUM AND SEMICLASSICAL OPTICS J. Opt. B: Quantum Semiclass. Opt. 3 (2001 S50 S54 www.iop.org/journals/ob PII: S1464-4266(0115159-1 Double-slit intererence
More informationOn the FPA infrared camera transfer function calculation
On the FPA infrared camera transfer function calculation (1) CERTES, Université Paris XII Val de Marne, Créteil, France (2) LTM, Université de Bourgogne, Le Creusot, France by S. Datcu 1, L. Ibos 1,Y.
More informationy2 = 0. Show that u = e2xsin(2y) satisfies Laplace's equation.
Review 1 1) State the largest possible domain o deinition or the unction (, ) = 3 - ) Determine the largest set o points in the -plane on which (, ) = sin-1( - ) deines a continuous unction 3) Find the
More informationINVARIANT RELATIVE SATELLITE MOTION. Marco Sabatini Dip. Ingegneria Aerospaziale, Università di Roma La Sapienza
AC Workshop on Innovative Concepts. ESA-ESEC 8-9 January 8. INVARIAN RELAIVE SAELLIE MOION Marco Sabatini Dip. Ingegneria Aerospaziale, Università di Roma La Sapienza marcosabatini@hotmail.it Dario Izzo
More informationMODELLING PUBLIC TRANSPORT CORRIDORS WITH AGGREGATE AND DISAGGREGATE DEMAND
MODELLIG PUBLIC TRASPORT CORRIDORS WITH AGGREGATE AD DISAGGREGATE DEMAD Sergio Jara-Díaz Alejandro Tirachini and Cristián Cortés Universidad de Chile ITRODUCTIO Traditional microeconomic modeling o public
More informationELEG 3143 Probability & Stochastic Process Ch. 4 Multiple Random Variables
Department o Electrical Engineering University o Arkansas ELEG 3143 Probability & Stochastic Process Ch. 4 Multiple Random Variables Dr. Jingxian Wu wuj@uark.edu OUTLINE 2 Two discrete random variables
More informationPreconditioning of Radial Basis Function Interpolation Systems via Accelerated Iterated Approximate Moving Least Squares Approximation
Preconditioning o Radial Basis Function Interpolation Systems via Accelerated Iterated Approximate Moving Least Squares Approximation Gregory E. Fasshauer and Jack G. Zhang Abstract The standard approach
More informationBayesian Technique for Reducing Uncertainty in Fatigue Failure Model
9IDM- Bayesian Technique or Reducing Uncertainty in Fatigue Failure Model Sriram Pattabhiraman and Nam H. Kim University o Florida, Gainesville, FL, 36 Copyright 8 SAE International ABSTRACT In this paper,
More informationand ( x, y) in a domain D R a unique real number denoted x y and b) = x y = {(, ) + 36} that is all points inside and on
Mat 7 Calculus III Updated on 10/4/07 Dr. Firoz Chapter 14 Partial Derivatives Section 14.1 Functions o Several Variables Deinition: A unction o two variables is a rule that assigns to each ordered pair
More informationSupplementary material for Continuous-action planning for discounted infinite-horizon nonlinear optimal control with Lipschitz values
Supplementary material or Continuous-action planning or discounted ininite-horizon nonlinear optimal control with Lipschitz values List o main notations x, X, u, U state, state space, action, action space,
More informationLeast-Squares Spectral Analysis Theory Summary
Least-Squares Spectral Analysis Theory Summary Reerence: Mtamakaya, J. D. (2012). Assessment o Atmospheric Pressure Loading on the International GNSS REPRO1 Solutions Periodic Signatures. Ph.D. dissertation,
More informationSection 3.4: Concavity and the second Derivative Test. Find any points of inflection of the graph of a function.
Unit 3: Applications o Dierentiation Section 3.4: Concavity and the second Derivative Test Determine intervals on which a unction is concave upward or concave downward. Find any points o inlection o the
More informationSpectrum Management Methodology for WCDMA Systems Encompassing Uplink and Downlink
Spectrum Management Methodology or WCMA Systems Encompassing plink and ownlink J. asreddine, J. Pérez-Romero, O. Sallent, R. Agustí ept. o Signal Theory and Communications, niversitat Politècnica de Catalunya
More informationKalman filtering based probabilistic nowcasting of object oriented tracked convective storms
Kalman iltering based probabilistic nowcasting o object oriented traced convective storms Pea J. Rossi,3, V. Chandrasear,2, Vesa Hasu 3 Finnish Meteorological Institute, Finland, Eri Palménin Auio, pea.rossi@mi.i
More informationInterior-Point Methods as Inexact Newton Methods. Silvia Bonettini Università di Modena e Reggio Emilia Italy
InteriorPoint Methods as Inexact Newton Methods Silvia Bonettini Università di Modena e Reggio Emilia Italy Valeria Ruggiero Università di Ferrara Emanuele Galligani Università di Modena e Reggio Emilia
More informationPSFs for Laser Guide Star Tomographic Adaptive Optics Systems
PSFs or Laser Guide Star Tomograpic Adaptive Optics Systems Donald Gavel Center or Adaptive Optics UC Santa Cru PSF Reconstruction Worsop, HIA Victoria May 11, 4 UCO Lic Observatory NSF Center or Adaptive
More informationMidterm Summary Fall 08. Yao Wang Polytechnic University, Brooklyn, NY 11201
Midterm Summary Fall 8 Yao Wang Polytechnic University, Brooklyn, NY 2 Components in Digital Image Processing Output are images Input Image Color Color image image processing Image Image restoration Image
More information2. ETA EVALUATIONS USING WEBER FUNCTIONS. Introduction
. ETA EVALUATIONS USING WEBER FUNCTIONS Introduction So ar we have seen some o the methods or providing eta evaluations that appear in the literature and we have seen some o the interesting properties
More informationWhat is Image Deblurring?
What is Image Deblurring? When we use a camera, we want the recorded image to be a faithful representation of the scene that we see but every image is more or less blurry, depending on the circumstances.
More informationFluctuationlessness Theorem and its Application to Boundary Value Problems of ODEs
Fluctuationlessness Theorem and its Application to Boundary Value Problems o ODEs NEJLA ALTAY İstanbul Technical University Inormatics Institute Maslak, 34469, İstanbul TÜRKİYE TURKEY) nejla@be.itu.edu.tr
More informationDifferentiation. The main problem of differential calculus deals with finding the slope of the tangent line at a point on a curve.
Dierentiation The main problem o dierential calculus deals with inding the slope o the tangent line at a point on a curve. deinition() : The slope o a curve at a point p is the slope, i it eists, o the
More informationThu June 16 Lecture Notes: Lattice Exercises I
Thu June 6 ecture Notes: attice Exercises I T. Satogata: June USPAS Accelerator Physics Most o these notes ollow the treatment in the class text, Conte and MacKay, Chapter 6 on attice Exercises. The portions
More informationTikhonov Regularized Poisson Likelihood Estimation: Theoretical Justification and a Computational Method
Inverse Problems in Science and Engineering Vol. 00, No. 00, December 2006, 1 19 Tikhonov Regularized Poisson Likelihood Estimation: Theoretical Justification and a Computational Method Johnathan M. Bardsley
More informationDifferential Equaitons Equations
Welcome to Multivariable Calculus / Dierential Equaitons Equations The Attached Packet is or all students who are planning to take Multibariable Multivariable Calculus/ Dierential Equations in the all.
More informationLeveraging Machine Learning for High-Resolution Restoration of Satellite Imagery
Leveraging Machine Learning for High-Resolution Restoration of Satellite Imagery Daniel L. Pimentel-Alarcón, Ashish Tiwari Georgia State University, Atlanta, GA Douglas A. Hope Hope Scientific Renaissance
More informationGIMC-based Fault Detection and Its Application to Magnetic Suspension System
Proceedings o the 7th World Congress The International Federation o Automatic Control Seoul, Korea, Jul 6-, 28 GIC-based Fault Detection and Its Application to agnetic Suspension Sstem Yujiro Nakaso Toru
More informationCircuit Complexity / Counting Problems
Lecture 5 Circuit Complexity / Counting Problems April 2, 24 Lecturer: Paul eame Notes: William Pentney 5. Circuit Complexity and Uniorm Complexity We will conclude our look at the basic relationship between
More informationHyperparameter estimation in image restoration problems with partially known blurs
Hyperparameter estimation in image restoration problems with partially known blurs Nikolas P. Galatsanos Illinois Institute o Technology Department o Electrical and Computer Engineering 330 South Dearborn
More informationZERO-DISTANCE PULSE FRONTS OF STRETCHER AND ITS OPTICAL SYSTEM
ERODISTANCE PULSE RONTS O STRETCHER AND ITS OPTICAL SYSTEM Author: DOI: 10.12684/alt.1.70 Corresponding author: email: agitin@mbiberlin.de erodistance Pulse ronts o a Stretcher and its Optical System Max
More informationRESOLUTION MSC.362(92) (Adopted on 14 June 2013) REVISED RECOMMENDATION ON A STANDARD METHOD FOR EVALUATING CROSS-FLOODING ARRANGEMENTS
(Adopted on 4 June 203) (Adopted on 4 June 203) ANNEX 8 (Adopted on 4 June 203) MSC 92/26/Add. Annex 8, page THE MARITIME SAFETY COMMITTEE, RECALLING Article 28(b) o the Convention on the International
More information