arxiv:math/ v2 [math.na] 15 Nov 2005
|
|
- Loren Harper
- 5 years ago
- Views:
Transcription
1 arxiv:math/5247v2 [math.na] 5 Nov 25 A Least Squares Functional for Solving Inverse Sturm-Liouville Problems. Introduction Norbert Röhrl Fachbereich Mathematik, IADM, University of Stuttgart, 75 Stuttgart, Germany norbert.roehrl@mathematik.uni-stuttgart.de AMS classification scheme numbers: 65L9, 34A55 Abstract. We present a variational algorithm for solving the classical inverse Sturm-Liouville problem in one dimension when two spectra are given. All critical points of the least squares functional are at global minima, which suggests minimization by a (conjugate) gradient descent algorithm. Numerical examples show that the resulting algorithm works quite reliable without tuning for particular cases, even in the presence of noise. With the right choice of parameters the minimization procedure converges very quickly. The one dimensional inverse Sturm-Liouville problem and methods for its numerical solution are very well studied [8,, 2]. Still there is no general purpose algorithm, which efficiently satisfies all needs (see eg. [2]). We present a variational algorithm which proves to be very robust under noisy input, does not require any special tuning or additional input besides the partial spectra, and is quite efficient in the class of variational methods. In section 2 we will define the functional and prove the absence of true local minima, which could trap our minimization procedure. Numerical examples, the optimal choice of the weights and the performance under noisy input will be discussed in section 3. Section 4 is concerned with implementation details and compares this algorithm to related approaches, and section 5 finally contains the proof of the linear independence of the squares of eigenfunctions. 2. Definition and Properties of the Functional We consider the Sturm-Liouville equation u + q(x)u = λu on [, ] with q(x) L 2 ([, ], R) real, and separated boundary conditions (SL) u()cosα + u ()sin α =, u()cosβ + u ()sin β =. (αβ) For the asymptotics of the eigenvalues (λ n ) n= wrt. boundary conditions (αβ) there are three cases (see eg. [3, 4, 3, 9]): π 2 n 2 2sin(β α) sin α sin β + q ds + a n if sinαsin β λ n = π 2 (n + /2) 2 2 cos αcos β + sin(β α) + q ds + a n if sinαsin β = sin(β α) π 2 (n + ) 2 + (2.3) q ds + a n if sinα = sin β = where (a n ) l 2.
2 Least Squares Functional for Inverse Sturm-Liouville Problems 2 It is a classical result, that two spectra determine the potential uniquely. Fix three angles α, β, and γ with sin(β γ). Let λ q,i,n, i {, 2} denote the eigenvalues with respect to the boundary conditions (αβ) resp. (αγ). Theorem 2. (Borg [], Levinson [5]). Given potentials Q, q L [, ] with then Q = q. λ Q,i,n = λ q,i,n for all n N, i =, 2, Suppose we are given (partial) spectral data λ Q,i,n with (i, n) in I {, 2} N of an unknown potential Q. For a trial potential q and positive weights ω i,n, we define the functional G(q) = ω i,n (λ q,i,n λ Q,i,n ) 2. (2.4) If I is infinite and the sequence (ω i,n ) is summable, this converges because of the asymptotics of the eigenvalues (2.3) G(q) = ( ω i,n Q q + c i,n )2 <, where c i,n l 2, i =, 2. If moreover I = {, 2} N, the functional is zero if and only if q = Q (by theorem 2.). In the case of partial knowledge of the spectrum, a solution of G(q) = includes all information given about the unknown potential Q. But we only have a chance to get close to the original potential if there is not too much information in the higher eigenvalues, i.e. if the potential to be recovered is sufficiently smooth [7]. An important question is now, for which sequences (λ i,n ) there is a q with G(q) =, i.e. λ q,i,n = λ i,n. It is easy to see with the classical methods [3], that it is necessary for the eigenvalues to interlace λ,n < λ 2,n < λ,n+ or λ 2,n < λ,n < λ 2,n+. It is also known that if we choose either (i) sinα, sin β, sin γ (Levitan [6]) or (ii) α = β =, sinγ (Dahlberg Trubowitz [3]), then the interlacing property in connection with the correct asymptotics is sufficient for the existence of a potential q with G(q) =. An interesting feature of this functional is that all its critical points are at global minima. To prove this, we first compute the derivative of G. The derivative of λ q,i,n wrt. q in direction h is λ q,i,n [h] = hg 2 q,i,n dx, where g q,i,n denotes the normalized eigenfunction corresponding to λ q,i,n (see [9] for a nice proof). Thus the derivative of G is given by Ġ[h](q) = 2 ω i,n (λ q,i,n λ Q,i,n )gq,i,nh 2 dx. We note that if nω i,n is summable, then the gradient G(q) = 2 (n + )ω i,n (λ q,i,n λ Q,i,n ) n + g2 q,i,n (2.5) is in H because gq,i,n 2 H = O(n) [3] and λ q,i,n λ Q,i,n = O(). Theorem 5.2 below shows that the eigenfunctions gq,i,n 2 are linearly independent in H. immediately implies the essential convexity of the functional: This
3 Least Squares Functional for Inverse Sturm-Liouville Problems 3 Theorem 2.2. If I is finite or (nω i,n ) is summable, the functional G has no local minima at q with G(q) >, i.e. Ġ[h](q) = h H [, ] G(q) =. Thus a conjugate gradient algorithm will not get trapped in local minimas, as we will also observe in the examples. Finally, we want to address the question, which of the infinitely many q with G(q) = our algorithm will select. Let us first look at a related functional, which actually was also our first try. It is given by G(q) = ( )2 λ q,i,n λ Q,i,n + (Q q)dx, (2.6) with derivative G[h](q) = 2 ) (λ q,i,n λ Q,i,n (Q q)dx (g 2 q,i,n )h dx. In case the mean of Q is known, this functional works as well as the other. Its gradient flow leaves the mean of q constant by construction. The function gq,i,n 2 is the gradient of λ q,i,n q dx and the direction of strongest increase of λ q,i,n which leaves q dx fixed. Choosing α = β = and γ = π/2, the asymptotics of the squared eigenfunctions are given by g 2 q,i,n = { cos ( (2n + 2)πx ) + O(n ), i = cos ( (2n + )πx ) + O(n ), i = 2. (2.7) Hence the functions {gq,i,n 2 n N, i =, 2} {} are almost orthogonal for large n. Now, the derivative of our functional (2.4) can be written as Ġ[h](q) = 2 ω i,n (λ q,i,n λ Q,i,n ) ( (g 2 q,i,n ) + ) h dx, and the difference λ q+c,i,n (q + c)dx = λ q,i,n q dx is invariant under adding a constant function c to q. It follows that a gradient flow of our functional G leaves λ q,i,n q dx with (i, n) I almost invariant for n large enough. Since also a conjugate gradient descent is just an approximation of the gradient flow, in practice we do get some little changes in the higher eigenvalues. A similar argument holds for the case sin α, sin β, sin γ, but there the asymptotics for i =, 2 are equal. So we still have asymptotic orthogonality in n, but can not separate i = and i = Numerical Examples We use the standard Polak-Ribiere conjugate gradient descent algorithm [] to approximate the gradient flow and thus minimize the functional. To give the basic idea, we explain the simpler steepest descent: (i) choose initial potential q = q (ii) while G(q i ) too big do (a) compute the gradient G(q i ) (b) minimize the one dimensional function G(q i α G(q i )) wrt. α (c) set q i+ equal to the minimizing potential
4 Least Squares Functional for Inverse Sturm-Liouville Problems 4 G(q) = =.8295 λ = iter# = Figure. Plot using ω i,n =. The right two graphs show G(q) and 2 over the number of iterations. G(q) G(q) = = λ = iter# = 4 G(q) Figure 2. Plot using ω i,n = (n + ) 2. On the right side of the 2 plot we zoom in to show the convergence process in more detail. 2 This straight forward minimization scheme has a major disadvantage: consecutive gradients are always orthogonal. To avoid this, conjugate gradient descent computes the direction for the one dimensional minimization using the current and previous gradients. The first natural question is the optimal choice of the constants ω i,n. A short look at the leading coefficients of the asymptotics (2.3) may suggest something like ω i,n = (n + ) 2. But the leading term of λ q,i,n λ Q,i,n is (q Q)dx = O(), which corresponds to choosing ω i,n = and yields the fastest convergence in all examples. Figures and 2 were computed with ω i,n = resp. ω i,n = (n + ) 2. For all plots (except 2 and 3) we chose the lowest number of iterations after which the plot stabilizes, while usually already the second iteration reveals the global structure of the potential. By default we use 3 pairs of eigenvalues, q =, ω i,n =, and the boundary conditions α = β =, γ = π/2. We give the current value of the functional G(q) as well as the L 2 error 2 = q Q 2, and the maximum of the difference of the eigenvalues λ = max { λ q,i,n λ Q,i,n }. In addition if there is a visible difference, we also plot the original potential Q using a dashed line. We see that both G(q) and 2 converge much faster for ω i,n =. With ω i,n = (n + ) 2 on the other hand, q moves through smoother functions (as in figure 2). In contrast to the optical impression these are in general worse approximations wrt. G(q) and 2. But finally q, after around iterations, will also converge to the function shown in figure. From the good average convergence of G(q) in both cases, we can tell, that our algorithm also in practice does not get near a local minimum, as proved in theorem 2.2. The optimal choice of boundary conditions can already be guessed from the asymptotics of the squared eigenfunctions (2.7). Choosing α = β = and γ = π/2, the functions gq,i,n 2 are almost
5 Least Squares Functional for Inverse Sturm-Liouville Problems 5 G(q) = =.2837 λ = iter# = Figure 3. Plot with bounday conditions α = β = π/4, γ = π/4. The light graph shows the 62th iteration, where the iteration stabilizes. G(q) G(q) = =.8535 λ = iter# = 5 G(q) = =.362 λ = iter# = 5 Figure 4. Plots with added white noise of absolute value smaller or equal. resp... orthogonal for large n and i =, 2 (c.f. (2.7)). In contrast, if sinα, sin β, and sinγ are all non zero, gq,,n 2 and g2 q,2,n will get close for large n. Therefore, in the latter case, the algorithm will converge much slower. This can also clearly be observed in numerical examples, like in figure 3. But also there, after approximately 62 iterations, we finally will get a solution which is close to figure. Another important aspect for applications is stability against noise in the given spectral data. In figure 4 we computed two examples with random noise λ Q,i,n λ Q,i,n r, with r =. and r =., respectively. It is remarkable, that the convergence speed in G(q) is not significantly affected by the random noise; both examples reach G(q) 8 in 3 iterations. For testing the robustness, we also fed the functional with the following random sequence ,.6265, , , , 3.9, , 48.88, , 6.988, which is far from the generic asymptotics. Convergence of G(q) is slow but steady. In log scale the graph of G(q) over the number of iteration looks similar to those given above. From G(q) = initially, we get down to G(q) = in 45 iterations. Finally figure 5 shows the results for the other examples of [2, ].
6 Least Squares Functional for Inverse Sturm-Liouville Problems 6 G(q) = = G(q) = λ = iter# = 4 2 = λ = iter# = 3 Figure 5. Continuous examples. 4. Implementation and Comparison The implementation of the conjugate gradient minimization algorithm is straight forward (see eg. []), but our application crucially depends on a reliable eigenvalue problem solver. After some trial and error we settled for d2kdf from the Numerical Algorithms Group (NAG). For standard numerical routines we used the GNU Scientific Library. Another important point is the internal representation of the functions. We used cubic splines on 2 intervals of equal length. For 3 pairs of eigenvalues we already get satisfactory results at 8 intervals. But since the computation time depends only roughly linear on this number, we did not study its influence on the quality of the approximation. Because all high level routines, like the eigenvalue solver, only sample these functions on some points, the details of the representation are not that important. But since very many values are needed, this is the part in the algorithm, where speed optimizations are most valuable. Compared to the variational algorithm by Brown et al. [2], which uses other spectral data, our functional is more expensive to compute, because they only have to solve initial value problems, whereas the present algorithm requires solution of the eigenvalue problems. On the other hand they often need about - times as many iterations to get similar results. Another plus on our side is, that even in presence of noise, we know to get closer to our goal each step and do not have to regularize the process by limiting the number of iterations. The method of Rundell and Sacks [] uses the same spectral data as our algorithm and is without any doubt much faster, but needs the mean Q dx as additional input, which has to be guessed from the (partial) spectral data. Our method does this automatically. In the case of figure for example, the error of the mean is 2 3 % and using 5 pairs of eigenvalues we still only get an error of about.2%. This is probably better than an independent algorithm for computing the mean value from spectral data could do. 5. Linear Independence of Squared Eigenfunctions Define the Wronskian [f, g] = fg f g and the bilinear form Γ : H ([, ], R) 2 R (f, g) [f, g] dx which is bounded by Γ(f, g) f H g H, i.e. Γ(f, ) = f H. (We use the definition f H = f 2 L + f 2 2 L with distributional derivatives.) In particular Γ is 2 continuous on H. We have the following rules for the Wronskian:
7 Least Squares Functional for Inverse Sturm-Liouville Problems 7 (i) [fg, hj] = gj[f, h] + fh[g, j] = fj[g, h] + gh[f, j] (ii) If the functions f and f 2 fulfill the condition then [f, f 2 ](a) =. f i (a)cosα + f i(a)sin α = i =, 2, (iii) For two arbitrary solutions f and f 2 of the equation (SL) with eigenvalue parameters λ and λ 2 we have [f, f 2 ] = f f 2 f 2f = (λ λ 2 )f f 2. In this section we are only talking about a single q and therefore drop it from the index. Let s i,n and c i,n be the solutions of the differential equation (SL) for the eigenvalue parameter λ i,n and initial values s i,n () = sinβ, c i,n () = sinγ, s i,n () = cosβ, c i,n () = cosγ. It is well known and easy to see that the Wronskian of these solutions is constant and we can compute its value at : [s i,n, c i,n ] = sinβ cosγ + cosβ sinγ = sin(γ β). The normalized eigenfunctions are g,n = s,n / s,n 2 and g 2,n = c 2,n / c 2,n 2. Now we prove the central lemma, which builds on the ideas of a similar result for the Dirichlet case in the book of Pöschel and Trubowitz[9]. A result in the same spirit can also already be found in the paper of Borg []. Lemma 5.. Given three angles α, β, and γ with sin(β γ) and denote the L 2 -normalized eigenfunctions of the Sturm-Liouville equation with boundary conditions (αβ) and (αγ) by g i,n, i =, 2. With s i,n, c i,n as defined above we have the following relations for the squared eigenfunctions for all i, j {, 2} and m, n N : (i) Γ(g 2 i,n, g 2 i,m) = (ii) Γ(c i,n s i,n, g 2 j,m) = ( ) i sin(γ β)δ n,m δ i,j Proof. i) Γ(gi,n 2, g2 i,m ) = 2 g i,n g i,m [g i,n, g i,m ] dx If n = m, the term clearly vanishes. If n m, we get Γ(gi,n, 2 gi,m) 2 2 = λ i,n λ i,m [g i,n, g i,m ] [g i,n, g i,m ] dx = because g i,n and g i,m satisfy the same boundary conditions (rule ii). ii) λ i,n λ i,m [g i,n, g i,m ] 2 =, Γ(c i,n s i,n, gj,m 2 ) = ( ci,n g j,m [s i,n, g j,m ] + s i,n g j,m [c i,n, g j,m ] ) dx If i = j =, m = n, the first term vanishes and we are left with s,m g,m [c,m, g,m ] = and if i = j = 2, m = n, the second term vanishes and we get c 2,m g 2,m [s 2,m, g 2,m ] = g 2,m [c,m, s,m ] = sin(γ β), g 2 2,m[s 2,m, c 2,m ] = + sin(γ β).
8 Least Squares Functional for Inverse Sturm-Liouville Problems 8 If (i, n) (j, m) we compute Γ(c i,n s i,n, gj,m 2 ) = ( [ci,n, g j,m ] [s i,n, g j,m ] + [s i,n, g j,m ] [c i,n, g j,m ] ) dx = λ i,n λ j,m [s i,n, g j,m ][c i,n, g j,m ] λ i,n λ = j,m We note that λ i,n λ j,m in this case, because the eigenfunctions satisfy the same boundary conditions at and different boundary conditions at. Theorem 5.2. With the notations of the above lemma, the set of squared eigenfunctions { g 2 i,n (i, n) {, 2} N } is linearly independent in H. Proof. Suppose for some fixed (i, n) we have g 2 i,n = k N a k g 2 k in H, where a k R and g k = g jk,m k with (j k, m k ) (i, n). But this would imply ( ) i sin(γ β) = Γ(c i,n s i,n, gi,n (c 2 ) = Γ i,n s i,n, ) a k gk 2 = Γ(c i,n s i,n, a k gk 2 ) = k N k N Acknowledgments The author wishes to thank Ian Knowles and Rudi Weikard for drawing his attention to this topic and for the interesting discussions during his time at the University of Alabama at Birmingham. References [] G. Borg. Eine Umkehrung der Sturm-Liouvilleschen Eigenwertaufgabe. Bestimmung der Differentialgleichung durch die Eigenwerte. Acta Math., 78: 96, 946. [2] B. M. Brown, V. S. Samko, I. W. Knowles, and M. Marletta. Inverse spectral problem for the Sturm-Liouville equation. Inverse Problems, 9(): , 23. [3] B. E. J. Dahlberg and E. Trubowitz. The inverse Sturm-Liouville problem. III. Comm. Pure Appl. Math., 37(2): , 984. [4] E. L. Isaacson, H. P. McKean, and E. Trubowitz. The inverse Sturm-Liouville problem. II. Comm. Pure Appl. Math., 37():, 984. [5] N. Levinson. The inverse Sturm-Liouville problem. Mat. Tidsskr. B., 949:25 3, 949. [6] B. M. Levitan. Inverse Sturm-Liouville problems. VSP, Zeist, 987. Translated from the Russian by O. Efimov. [7] M. Marletta and R. Weikard. Weak stability for an inverse Sturm-Liouville problem with finite spectral data and complex potential. Inverse Problems, 2(4):275 29, 25. [8] J. R. McLaughlin. Analytical methods for recovering coefficients in differential equations from spectral data. SIAM Rev., 28():53 72, 986. [9] J. Pöschel and E. Trubowitz. Inverse spectral theory, volume 3 of Pure and Applied Mathematics. Academic Press Inc., Boston, MA, 987. [] W. H. Press, S. A. Teukolsky, W. T. Vetterling, and B. P. Flannery. Numerical recipes in C. Cambridge University Press, Cambridge, second edition, 992. [] W. Rundell and P. E. Sacks. Reconstruction techniques for classical inverse Sturm-Liouville problems. Math. Comp., 58(97):6 83, 992. [2] P. Schapotschnikow. Inverse Sturm-Liouville problems on trees. a variational approach. Master s thesis, School of Computer Science, Cardiff University, UK, 24. [3] E. C. Titchmarsh. Eigenfunction expansions associated with second-order differential equations. Part I. Second Edition. Clarendon Press, Oxford, 962.
Universität Stuttgart. Fachbereich Mathematik. A Least Squares Functional for Solving Inverse Sturm-Liouville Problems.
Universität Stuttgart Fachbereich Mathematik A Least Squares Functional for Solving Inverse Sturm-Liouville Problems Norbert Röhrl Preprint 25/4 Universität Stuttgart Fachbereich Mathematik A Least Squares
More informationINVERSE PROBLEMS FOR STURM-LIOUVILLE OPERATORS WITH BOUNDARY CONDITIONS DEPENDING ON A SPECTRAL PARAMETER
Electronic Journal of Differential Equations, Vol. 217 (217), No. 26, pp. 1 7. ISSN: 172-6691. URL: http://ejde.math.txstate.edu or http://ejde.math.unt.edu INVERSE PROBLEMS FOR STURM-LIOUVILLE OPERATORS
More informationOpen Society Foundation - Armenia. The inverse Sturm-Liouville problem with fixed boundary conditions. Chair of Differential Equations YSU, Armenia
Open Society Foundation - Armenia Physical and Mathematical Sciences 24, p. 8 M a t h e m a t i c s The inverse Sturm-Liouville problem with fixed boundary conditions YU. A. ASHRAFYAN, T. N. HARUTYUNYAN
More informationInverse Nodal Problems for Second Order Differential Operators with a Regular Singularity
International Journal of Difference Equations. ISSN 973-669 Volume 1 Number 6), pp. 41 47 c Research India Publications http://www.ripublication.com/ide.htm Inverse Nodal Problems for Second Order Differential
More informationA uniqueness result for one-dimensional inverse scattering
Math. Nachr. 285, No. 8 9, 941 948 (2012) / DOI 10.1002/mana.201100101 A uniqueness result for one-dimensional inverse scattering C. Bennewitz 1,B.M.Brown 2, and R. Weikard 3 1 Department of Mathematics,
More informationREPRESENTATION OF THE NORMING CONSTANTS BY TWO SPECTRA
Electronic Journal of Differential Equations, Vol. 0000, No. 59, pp. 0. ISSN: 07-669. URL: http://ejde.math.txstate.edu or http://ejde.math.unt.edu ftp ejde.math.txstate.edu REPRESENTATION OF THE NORMING
More informationIMRN Inverse Spectral Analysis with Partial Information on the Potential III. Updating Boundary Conditions Rafael del Rio Fritz Gesztesy
IMRN International Mathematics Research Notices 1997, No. 15 Inverse Spectral Analysis with Partial Information on the Potential, III. Updating Boundary Conditions Rafael del Rio, Fritz Gesztesy, and Barry
More informationAn Inverse Problem for the Matrix Schrödinger Equation
Journal of Mathematical Analysis and Applications 267, 564 575 (22) doi:1.16/jmaa.21.7792, available online at http://www.idealibrary.com on An Inverse Problem for the Matrix Schrödinger Equation Robert
More informationJournal of Mathematical Analysis and Applications
J. Math. Anal. Appl. 381 (2011) 506 512 Contents lists available at ScienceDirect Journal of Mathematical Analysis and Applications www.elsevier.com/locate/jmaa Inverse indefinite Sturm Liouville problems
More informationINVERSE NODAL PROBLEM FOR STURM-LIOUVILLE OPERATORS WITH COULOMB POTENTIAL M. Sat 1, E.S. Panakhov 2
International Journal of Pure and Applied Mathematics Volume 8 No. 2 212, 173-18 ISSN: 1311-88 printed version url: http://www.ipam.eu PA ipam.eu INVERSE NODAL PROBLEM FOR STURM-LIOUVILLE OPERATORS WITH
More informationResearch Article A Half-Inverse Problem for Impulsive Dirac Operator with Discontinuous Coefficient
Abstract and Applied Analysis Volume 213, Article ID 18189, 4 pages http://d.doi.org/1.1155/213/18189 Research Article A Half-Inverse Problem for Impulsive Dirac Operator with Discontinuous Coefficient
More informationAN INVERSE PROBLEMS FOR STURM-LIOUVILLE-TYPE DIFFERENTIAL EQUATION WITH A CONSTANT DELAY
SARAJEVO JOURNAL OF MATHEMATICS Vol.12 (24), No.1, (216), 83 88 DOI: 1.5644/SJM.12.1.6 AN INVERSE PROBLEMS FOR STURM-LIOUVILLE-TYPE DIFFERENTIAL EQUATION WITH A CONSTANT DELAY VLADIMIR VLADIČIĆ AND MILENKO
More informationarxiv:math-ph/ v1 27 Aug 2003
ON HALF-LINE SPECTRA FOR A CLASS OF NON-SELF-ADJOINT HILL OPERATORS KWANG C. SHIN arxiv:math-ph/3832v1 27 Aug 23 Abstract. In 198, Gasymov showed that non-self-adjoint Hill operators with complexvalued
More informationarxiv: v1 [math.ca] 27 Mar 2013
Boundary value problem with transmission conditions arxiv:133.6947v1 [math.ca] 27 Mar 213 K. Aydemir and O. Sh. Mukhtarov Department of Mathematics, Faculty of Science, Gaziosmanpaşa University, 625 Tokat,
More informationOn the existence of an eigenvalue below the essential spectrum
On the existence of an eigenvalue below the essential spectrum B.M.Brown, D.K.R.M c Cormack Department of Computer Science, Cardiff University of Wales, Cardiff, PO Box 916, Cardiff CF2 3XF, U.K. A. Zettl
More informationDesign and Optimization of Energy Systems Prof. C. Balaji Department of Mechanical Engineering Indian Institute of Technology, Madras
Design and Optimization of Energy Systems Prof. C. Balaji Department of Mechanical Engineering Indian Institute of Technology, Madras Lecture - 09 Newton-Raphson Method Contd We will continue with our
More informationSufficient conditions for functions to form Riesz bases in L 2 and applications to nonlinear boundary-value problems
Electronic Journal of Differential Equations, Vol. 200(200), No. 74, pp. 0. ISSN: 072-669. URL: http://ejde.math.swt.edu or http://ejde.math.unt.edu ftp ejde.math.swt.edu (login: ftp) Sufficient conditions
More informationConjugate Directions for Stochastic Gradient Descent
Conjugate Directions for Stochastic Gradient Descent Nicol N Schraudolph Thore Graepel Institute of Computational Science ETH Zürich, Switzerland {schraudo,graepel}@infethzch Abstract The method of conjugate
More informationarxiv:math/ v1 [math.sp] 18 Jan 2003
INVERSE SPECTRAL PROBLEMS FOR STURM-LIOUVILLE OPERATORS WITH SINGULAR POTENTIALS, II. RECONSTRUCTION BY TWO SPECTRA arxiv:math/31193v1 [math.sp] 18 Jan 23 ROSTYSLAV O. HRYNIV AND YAROSLAV V. MYKYTYUK Abstract.
More informationMulti-Robotic Systems
CHAPTER 9 Multi-Robotic Systems The topic of multi-robotic systems is quite popular now. It is believed that such systems can have the following benefits: Improved performance ( winning by numbers ) Distributed
More informationSTURM-LIOUVILLE PROBLEMS WITH RETARDED ARGUMENT AND A FINITE NUMBER OF TRANSMISSION CONDITIONS
Electronic Journal of Differential Equations, Vol. 217 (217), No. 31, pp. 1 8. ISSN: 172-6691. URL: http://ejde.math.txstate.edu or http://ejde.math.unt.edu STURM-LIOUVILLE PROBLEMS WITH RETARDED ARGUMENT
More informationCONSTRUCTION OF THE HALF-LINE POTENTIAL FROM THE JOST FUNCTION
CONSTRUCTION OF THE HALF-LINE POTENTIAL FROM THE JOST FUNCTION Tuncay Aktosun Department of Mathematics and Statistics Mississippi State University Mississippi State, MS 39762 Abstract: For the one-dimensional
More informationDeviation Measures and Normals of Convex Bodies
Beiträge zur Algebra und Geometrie Contributions to Algebra Geometry Volume 45 (2004), No. 1, 155-167. Deviation Measures Normals of Convex Bodies Dedicated to Professor August Florian on the occasion
More informationOn an inverse problem for Sturm-Liouville Equation
EUROPEAN JOURNAL OF PURE AND APPLIED MATHEMATICS Vol. 1, No. 3, 17, 535-543 ISSN 137-5543 www.ejpam.com Published by New York Business Global On an inverse problem for Sturm-Liouville Equation Döne Karahan
More informationUniqueness of the potential function for the vectorial Sturm-Liouville equation on a finite interval
RESEARCH Uniqueness of the potential function for the vectorial Sturm-Liouville equation on a finite interval Tsorng-Hwa Chang 1,2 and Chung-Tsun Shieh 1* Open Access * Correspondence: ctshieh@mail. tku.edu.tw
More informationGradient Method Based on Roots of A
Journal of Scientific Computing, Vol. 15, No. 4, 2000 Solving Ax Using a Modified Conjugate Gradient Method Based on Roots of A Paul F. Fischer 1 and Sigal Gottlieb 2 Received January 23, 2001; accepted
More informationZEROS OF THE JOST FUNCTION FOR A CLASS OF EXPONENTIALLY DECAYING POTENTIALS. 1. Introduction
Electronic Journal of Differential Equations, Vol. 20052005), No. 45, pp. 9. ISSN: 072-669. URL: http://ejde.math.tstate.edu or http://ejde.math.unt.edu ftp ejde.math.tstate.edu login: ftp) ZEROS OF THE
More information1 Expansion in orthogonal functions
Notes "orthogonal" January 9 Expansion in orthogonal functions To obtain a more useful form of the Green s function, we ll want to expand in orthogonal functions that are (relatively) easy to integrate.
More informationConjugate gradient algorithm for training neural networks
. Introduction Recall that in the steepest-descent neural network training algorithm, consecutive line-search directions are orthogonal, such that, where, gwt [ ( + ) ] denotes E[ w( t + ) ], the gradient
More informationAugmented Reality VU numerical optimization algorithms. Prof. Vincent Lepetit
Augmented Reality VU numerical optimization algorithms Prof. Vincent Lepetit P3P: can exploit only 3 correspondences; DLT: difficult to exploit the knowledge of the internal parameters correctly, and the
More informationEXISTENCE OF SOLUTIONS FOR BOUNDARY VALUE PROBLEMS ON INFINITE INTERVALS. We consider the second order nonlinear differential equation
Dynamic Systems and Applications 17 (28) 653-662 EXISTECE OF SOLUTIOS FOR BOUDARY VALUE PROBLEMS O IFIITE ITERVALS İSMAİL YASLA Department of Mathematics, Pamukkale University, Denizli 27, Turkey ABSTRACT.
More informationfür Mathematik in den Naturwissenschaften Leipzig
Max-Planck-Institut für Mathematik in den Naturwissenschaften Leipzig A Note on the Poincaré Inequality for Convex Domains by Mario Bebendorf Preprint no.: 24 23 ANoteonthePoincaré Inequality for Convex
More informationOn an Eigenvalue Problem Involving Legendre Functions by Nicholas J. Rose North Carolina State University
On an Eigenvalue Problem Involving Legendre Functions by Nicholas J. Rose North Carolina State University njrose@math.ncsu.edu 1. INTRODUCTION. The classical eigenvalue problem for the Legendre Polynomials
More informationInverse Problems, Information Content and Stratospheric Aerosols
Inverse Problems, Information Content and Stratospheric Aerosols May 22, 214 Table of contents 1 Inverse Problems 2 Aerosols 3 OSIRIS 4 Information Content 5 Non-Linear Techniques 6 Results Inverse Problems
More information2 Some estimates for the first eigenvalue of a Sturm-Liouville problem with a weighted integral condition
On some estimates for the first eigenvalue of a Sturm-Liouville problem with Dirichlet boundary conditions and a weighted integral condition S Ezhak, M Telnova Plekhanov Russian University of Economics,
More informationLine Search Methods for Unconstrained Optimisation
Line Search Methods for Unconstrained Optimisation Lecture 8, Numerical Linear Algebra and Optimisation Oxford University Computing Laboratory, MT 2007 Dr Raphael Hauser (hauser@comlab.ox.ac.uk) The Generic
More informationThe speed of Shor s R-algorithm
IMA Journal of Numerical Analysis 2008) 28, 711 720 doi:10.1093/imanum/drn008 Advance Access publication on September 12, 2008 The speed of Shor s R-algorithm J. V. BURKE Department of Mathematics, University
More informationIntroduction to gradient descent
6-1: Introduction to gradient descent Prof. J.C. Kao, UCLA Introduction to gradient descent Derivation and intuitions Hessian 6-2: Introduction to gradient descent Prof. J.C. Kao, UCLA Introduction Our
More informationRelevant sections from AMATH 351 Course Notes (Wainwright): 1.3 Relevant sections from AMATH 351 Course Notes (Poulin and Ingalls): 1.1.
Lecture 8 Qualitative Behaviour of Solutions to ODEs Relevant sections from AMATH 351 Course Notes (Wainwright): 1.3 Relevant sections from AMATH 351 Course Notes (Poulin and Ingalls): 1.1.1 The last few
More informationNonlinear equations and optimization
Notes for 2017-03-29 Nonlinear equations and optimization For the next month or so, we will be discussing methods for solving nonlinear systems of equations and multivariate optimization problems. We will
More informationLecture IX. Definition 1 A non-singular Sturm 1 -Liouville 2 problem consists of a second order linear differential equation of the form.
Lecture IX Abstract When solving PDEs it is often necessary to represent the solution in terms of a series of orthogonal functions. One way to obtain an orthogonal family of functions is by solving a particular
More informationComparison of Modern Stochastic Optimization Algorithms
Comparison of Modern Stochastic Optimization Algorithms George Papamakarios December 214 Abstract Gradient-based optimization methods are popular in machine learning applications. In large-scale problems,
More information( 2 + k 2 )Ψ = 0 (1) c sinhα sinφ. x = y = c sinβ. z =
Power series Schrodinger eigenfunctions for a particle on the torus. Mario Encinosa and Babak Etemadi Department of Physics Florida A & M University Tallahassee, Florida 3307 Abstract The eigenvalues and
More informationSteepest descent algorithm. Conjugate gradient training algorithm. Steepest descent algorithm. Remember previous examples
Conjugate gradient training algorithm Steepest descent algorithm So far: Heuristic improvements to gradient descent (momentum Steepest descent training algorithm Can we do better? Definitions: weight vector
More informationDETERMINATION OF AN UNKNOWN SOURCE TERM IN A SPACE-TIME FRACTIONAL DIFFUSION EQUATION
Journal of Fractional Calculus and Applications, Vol. 6(1) Jan. 2015, pp. 83-90. ISSN: 2090-5858. http://fcag-egypt.com/journals/jfca/ DETERMINATION OF AN UNKNOWN SOURCE TERM IN A SPACE-TIME FRACTIONAL
More informationSpectral Decimation for Families of Laplacians on the Sierpinski Gasket
Spectral Decimation for Families of Laplacians on the Sierpinski Gasket Seraphina Lee November, 7 Seraphina Lee Spectral Decimation November, 7 / 53 Outline Definitions: Sierpinski gasket, self-similarity,
More informationJOURNÉES ÉQUATIONS AUX DÉRIVÉES PARTIELLES
JOURNÉES ÉQUATIONS AUX DÉRIVÉES PARTIELLES HIROSHI ISOZAKI Some remarks on the multi-dimensional Borg-Levinson theorem Journées Équations aux dérivées partielles (1989), p. 1-6.
More informationProperties of the Scattering Transform on the Real Line
Journal of Mathematical Analysis and Applications 58, 3 43 (001 doi:10.1006/jmaa.000.7375, available online at http://www.idealibrary.com on Properties of the Scattering Transform on the Real Line Michael
More information8 Numerical methods for unconstrained problems
8 Numerical methods for unconstrained problems Optimization is one of the important fields in numerical computation, beside solving differential equations and linear systems. We can see that these fields
More informationGradient Descent. Sargur Srihari
Gradient Descent Sargur srihari@cedar.buffalo.edu 1 Topics Simple Gradient Descent/Ascent Difficulties with Simple Gradient Descent Line Search Brent s Method Conjugate Gradient Descent Weight vectors
More informationMinimization of Static! Cost Functions!
Minimization of Static Cost Functions Robert Stengel Optimal Control and Estimation, MAE 546, Princeton University, 2017 J = Static cost function with constant control parameter vector, u Conditions for
More information1 Introduction. Green s function notes 2018
Green s function notes 8 Introduction Back in the "formal" notes, we derived the potential in terms of the Green s function. Dirichlet problem: Equation (7) in "formal" notes is Φ () Z ( ) ( ) 3 Z Φ (
More informationTwo special equations: Bessel s and Legendre s equations. p Fourier-Bessel and Fourier-Legendre series. p
LECTURE 1 Table of Contents Two special equations: Bessel s and Legendre s equations. p. 259-268. Fourier-Bessel and Fourier-Legendre series. p. 453-460. Boundary value problems in other coordinate system.
More informationAsymptotic distributions of Neumann problem for Sturm-Liouville equation
Computational Methods for Differential Equations http://cmde.tabrizu.ac.ir Vol. 2, No., 24, pp. 9-25 Asymptotic distributions of Neumann problem for Sturm-Liouville equation H.R. Marasi Department of Mathematics,
More informationMATH 590: Meshfree Methods
MATH 590: Meshfree Methods Chapter 2 Part 3: Native Space for Positive Definite Kernels Greg Fasshauer Department of Applied Mathematics Illinois Institute of Technology Fall 2014 fasshauer@iit.edu MATH
More informationNonlinear Optimization
Nonlinear Optimization (Com S 477/577 Notes) Yan-Bin Jia Nov 7, 2017 1 Introduction Given a single function f that depends on one or more independent variable, we want to find the values of those variables
More informationInverse Sturm-Liouville problems using multiple spectra. Mihaela Cristina Drignei. A dissertation submitted to the graduate faculty
Inverse Sturm-Liouville problems using multiple spectra by Mihaela Cristina Drignei A dissertation submitted to the graduate faculty in partial fulfillment of the requirements for the degree of DOCTOR
More informationTHE RECONSTRUCTION OF GROUNDWATER PARAMETERS FROM HEAD DATA IN AN UNCONFINED AQUIFER
THE RECONSTRUCTION OF GROUNDWATER PARAMETERS FROM HEAD DATA IN AN UNCONFINED AQUIFER IAN KNOWLES AND AIMIN AN Abstract. The estimation of groundwater flow parameters from head measurements and other ancillary
More informationA new inversion method for dissipating electromagnetic problem
A new inversion method for dissipating electromagnetic problem Elena Cherkaev Abstract The paper suggests a new technique for solution of the inverse problem for Maxwell s equations in a dissipating medium.
More informationGradient Descent and Implementation Solving the Euler-Lagrange Equations in Practice
1 Lecture Notes, HCI, 4.1.211 Chapter 2 Gradient Descent and Implementation Solving the Euler-Lagrange Equations in Practice Bastian Goldlücke Computer Vision Group Technical University of Munich 2 Bastian
More informationAn a posteriori error estimate and a Comparison Theorem for the nonconforming P 1 element
Calcolo manuscript No. (will be inserted by the editor) An a posteriori error estimate and a Comparison Theorem for the nonconforming P 1 element Dietrich Braess Faculty of Mathematics, Ruhr-University
More informationIDENTIFICATION OF A POLYNOMIAL IN NONSEPARATED BOUNDARY CONDITIONS IN THE CASE OF A MULTIPLE ZERO EIGENVALUE
ISSN 2304-0122 Ufa Mathematical Journal. Vol. 7. No 1 (2015). P. 13-18. doi:10.13108/2015-7-1-13 UDC 517.984.54 IDENTIFICATION OF A POLYNOMIAL IN NONSEPARATED BOUNDARY CONDITIONS IN THE CASE OF A MULTIPLE
More informationNotes for CS542G (Iterative Solvers for Linear Systems)
Notes for CS542G (Iterative Solvers for Linear Systems) Robert Bridson November 20, 2007 1 The Basics We re now looking at efficient ways to solve the linear system of equations Ax = b where in this course,
More informationNotes for the Physics-Based Calculus workshop
Notes for the hysics-based Calculus workshop Adam Coffman June 9, 24 Trigonometric Functions We recall the following identities for trigonometric functions. Theorem.. For all x R, cos( x = cos(x and sin(
More informationKadath: a spectral solver and its application to black hole spacetimes
Kadath: a spectral solver and its application to black hole spacetimes Philippe Grandclément Laboratoire de l Univers et de ses Théories (LUTH) CNRS / Observatoire de Paris F-92195 Meudon, France philippe.grandclement@obspm.fr
More informationIn this worksheet we will use the eigenfunction expansion to solve nonhomogeneous equation.
Eigen Function Expansion and Applications. In this worksheet we will use the eigenfunction expansion to solve nonhomogeneous equation. a/ The theory. b/ Example: Solving the Euler equation in two ways.
More informationIterative Methods for Solving A x = b
Iterative Methods for Solving A x = b A good (free) online source for iterative methods for solving A x = b is given in the description of a set of iterative solvers called templates found at netlib: http
More informationFitting functions to data
1 Fitting functions to data 1.1 Exact fitting 1.1.1 Introduction Suppose we have a set of real-number data pairs x i, y i, i = 1, 2,, N. These can be considered to be a set of points in the xy-plane. They
More informationPhysics 403. Segev BenZvi. Numerical Methods, Maximum Likelihood, and Least Squares. Department of Physics and Astronomy University of Rochester
Physics 403 Numerical Methods, Maximum Likelihood, and Least Squares Segev BenZvi Department of Physics and Astronomy University of Rochester Table of Contents 1 Review of Last Class Quadratic Approximation
More informationMidterm Solution
18303 Midterm Solution Problem 1: SLP with mixed boundary conditions Consider the following regular) Sturm-Liouville eigenvalue problem consisting in finding scalars λ and functions v : [0, b] R b > 0),
More informationGENG2140, S2, 2012 Week 7: Curve fitting
GENG2140, S2, 2012 Week 7: Curve fitting Curve fitting is the process of constructing a curve, or mathematical function, f(x) that has the best fit to a series of data points Involves fitting lines and
More informationNearest Correlation Matrix
Nearest Correlation Matrix The NAG Library has a range of functionality in the area of computing the nearest correlation matrix. In this article we take a look at nearest correlation matrix problems, giving
More informationLaplace's equation: the potential between parallel plates
4/3/01 Electrostatics Laplace Solver 1 Laplace's equation: the potential between parallel plates Laplace's equation describing the electric potential in two dimensions is: ( x, y) 0 At right is the potential
More informationORTHOGONAL FUNCTIONS EXACT INVARIANT AND THE ADIABATIC LIMIT FOR TIME DEPENDENT HARMONIC OSCILLATORS
ORTHOGONAL FUNCTIONS EXACT INVARIANT AND THE ADIABATIC LIMIT FOR TIME DEPENDENT HARMONIC OSCILLATORS M. Fernández Guasti Depto. de Física, CBI. Universidad Autónoma Metropolitana - Iztapalapa, Av. San
More informationNormalization integrals of orthogonal Heun functions
Normalization integrals of orthogonal Heun functions Peter A. Becker a) Center for Earth Observing and Space Research, Institute for Computational Sciences and Informatics, and Department of Physics and
More informationChapter 4. Unconstrained optimization
Chapter 4. Unconstrained optimization Version: 28-10-2012 Material: (for details see) Chapter 11 in [FKS] (pp.251-276) A reference e.g. L.11.2 refers to the corresponding Lemma in the book [FKS] PDF-file
More informationOlga Boyko, Olga Martinyuk, and Vyacheslav Pivovarchik
Opuscula Math. 36, no. 3 016), 301 314 http://dx.doi.org/10.7494/opmath.016.36.3.301 Opuscula Mathematica HIGHER ORDER NEVANLINNA FUNCTIONS AND THE INVERSE THREE SPECTRA PROBLEM Olga Boyo, Olga Martinyu,
More informationNeural Network Training
Neural Network Training Sargur Srihari Topics in Network Training 0. Neural network parameters Probabilistic problem formulation Specifying the activation and error functions for Regression Binary classification
More informationOrdinary Differential Equations II
Ordinary Differential Equations II February 23 2017 Separation of variables Wave eq. (PDE) 2 u t (t, x) = 2 u 2 c2 (t, x), x2 c > 0 constant. Describes small vibrations in a homogeneous string. u(t, x)
More informationETNA Kent State University
Electronic Transactions on Numerical Analysis. Volume 37, pp. 166-172, 2010. Copyright 2010,. ISSN 1068-9613. ETNA A GRADIENT RECOVERY OPERATOR BASED ON AN OBLIQUE PROJECTION BISHNU P. LAMICHHANE Abstract.
More informationMath 115 Spring 11 Written Homework 10 Solutions
Math 5 Spring Written Homework 0 Solutions. For following its, state what indeterminate form the its are in and evaluate the its. (a) 3x 4x 4 x x 8 Solution: This is in indeterminate form 0. Algebraically,
More informationAPPENDIX A. Background Mathematics. A.1 Linear Algebra. Vector algebra. Let x denote the n-dimensional column vector with components x 1 x 2.
APPENDIX A Background Mathematics A. Linear Algebra A.. Vector algebra Let x denote the n-dimensional column vector with components 0 x x 2 B C @. A x n Definition 6 (scalar product). The scalar product
More informationSTURM-LIOUVILLE EIGENVALUE CHARACTERIZATIONS
Electronic Journal of Differential Equations, Vol. 2012 (2012), No. 123, pp. 1 13. ISSN: 1072-6691. URL: http://ejde.math.txstate.edu or http://ejde.math.unt.edu ftp ejde.math.txstate.edu STURM-LIOUVILLE
More informationarxiv: v1 [math.ca] 27 Mar 2013
Modified Expansion Theorem for Sturm-Liouville problem with transmission conditions arxiv:133.6898v1 [math.ca] 27 Mar 213 K.Aydemir and O. Sh. Mukhtarov Department of Mathematics, Faculty of Science, Gaziosmanpaşa
More informationInverse problem for Sturm Liouville operators with a transmission and parameter dependent boundary conditions
CJMS. 6()(17), 17-119 Caspian Journal of Mathematical Sciences (CJMS) University of Mazaran, Iran http://cjms.journals.umz.ac.ir ISSN: 1735-611 Inverse problem for Sturm Liouville operators with a transmission
More informationAN EIGENVALUE STUDY ON THE SUFFICIENT DESCENT PROPERTY OF A MODIFIED POLAK-RIBIÈRE-POLYAK CONJUGATE GRADIENT METHOD S.
Bull. Iranian Math. Soc. Vol. 40 (2014), No. 1, pp. 235 242 Online ISSN: 1735-8515 AN EIGENVALUE STUDY ON THE SUFFICIENT DESCENT PROPERTY OF A MODIFIED POLAK-RIBIÈRE-POLYAK CONJUGATE GRADIENT METHOD S.
More informationConjugate Gradient (CG) Method
Conjugate Gradient (CG) Method by K. Ozawa 1 Introduction In the series of this lecture, I will introduce the conjugate gradient method, which solves efficiently large scale sparse linear simultaneous
More informationOscillation theorems for the Dirac operator with spectral parameter in the boundary condition
Electronic Journal of Qualitative Theory of Differential Equations 216, No. 115, 1 1; doi: 1.14232/ejqtde.216.1.115 http://www.math.u-szeged.hu/ejqtde/ Oscillation theorems for the Dirac operator with
More information7.2 Steepest Descent and Preconditioning
7.2 Steepest Descent and Preconditioning Descent methods are a broad class of iterative methods for finding solutions of the linear system Ax = b for symmetric positive definite matrix A R n n. Consider
More information1 A complete Fourier series solution
Math 128 Notes 13 In this last set of notes I will try to tie up some loose ends. 1 A complete Fourier series solution First here is an example of the full solution of a pde by Fourier series. Consider
More informationNonlinear Optimization: What s important?
Nonlinear Optimization: What s important? Julian Hall 10th May 2012 Convexity: convex problems A local minimizer is a global minimizer A solution of f (x) = 0 (stationary point) is a minimizer A global
More informationA Multigrid Method for Two Dimensional Maxwell Interface Problems
A Multigrid Method for Two Dimensional Maxwell Interface Problems Susanne C. Brenner Department of Mathematics and Center for Computation & Technology Louisiana State University USA JSA 2013 Outline A
More informationLearning with Momentum, Conjugate Gradient Learning
Learning with Momentum, Conjugate Gradient Learning Introduction to Neural Networks : Lecture 8 John A. Bullinaria, 2004 1. Visualising Learning 2. Learning with Momentum 3. Learning with Line Searches
More informationOn Acceleration with Noise-Corrupted Gradients. + m k 1 (x). By the definition of Bregman divergence:
A Omitted Proofs from Section 3 Proof of Lemma 3 Let m x) = a i On Acceleration with Noise-Corrupted Gradients fxi ), u x i D ψ u, x 0 ) denote the function under the minimum in the lower bound By Proposition
More informationVasil Khalidov & Miles Hansard. C.M. Bishop s PRML: Chapter 5; Neural Networks
C.M. Bishop s PRML: Chapter 5; Neural Networks Introduction The aim is, as before, to find useful decompositions of the target variable; t(x) = y(x, w) + ɛ(x) (3.7) t(x n ) and x n are the observations,
More informationDegree Master of Science in Mathematical Modelling and Scientific Computing Mathematical Methods I Thursday, 12th January 2012, 9:30 a.m.- 11:30 a.m.
Degree Master of Science in Mathematical Modelling and Scientific Computing Mathematical Methods I Thursday, 12th January 2012, 9:30 a.m.- 11:30 a.m. Candidates should submit answers to a maximum of four
More informationA Uniqueness the Theorem for Singular Sturm-Liouville Problem
C.Ü. Fen-Edebiyat Fakültesi Fen Bilimleri Dergisi 3Cilt Sayı A Uniqueness the Theorem for Singular Sturm-Liouville Problem Etibar S. PANAKOV and ikmet KOYUNBAKAN Firat University Department of Mathematics
More informationIntroduction to Hilbert spaces. I.
Teoretisk Fysik, KTH September 5, 26 5A135 Fysikens matematiska metoder, ht6. Introduction to Hilbert spaces. I. My aim here is to give a down-to-earth introduction to the theory of Hilbert spaces which
More informationAn invitation to log geometry p.1
An invitation to log geometry James M c Kernan UCSB An invitation to log geometry p.1 An easy integral Everyone knows how to evaluate the following integral: 1 0 1 1 x 2 dx. An invitation to log geometry
More informationSpring 2014: Computational and Variational Methods for Inverse Problems CSE 397/GEO 391/ME 397/ORI 397 Assignment 4 (due 14 April 2014)
Spring 2014: Computational and Variational Methods for Inverse Problems CSE 397/GEO 391/ME 397/ORI 397 Assignment 4 (due 14 April 2014) The first problem in this assignment is a paper-and-pencil exercise
More information