Non-negatively constrained least squares and parameter choice by the residual periodogram for the inversion of electrochemical impedance spectroscopy data

The inverse problem associated with electrochemical impedance spectroscopy requiring the solution of a Fredholm integral equation of the first kind is considered. If the underlying physical model is not clearly determined, the inverse problem needs to be solved using a regularized linear least squares problem that is obtained from the discretization of the integral equation. For this system, it is shown that the model error can be made negligible by a change of variables and by extending the effective range of quadrature. This change of variables serves as a right preconditioner that significantly improves the condition of the system. Still, to obtain feasible solutions the additional constraint of non-negativity is required. Simulations with artificial, but realistic, data demonstrate that the use of non-negatively constrained least squares with a smoothing norm provides higher quality solutions than those obtained without the non-negative constraint. Using higher-order smoothing norms also reduces the error in the solutions. The L-curve and residual periodogram parameter choice criteria, which are used for parameter choice with regularized linear least squares, are successfully adapted to be used for the non-negatively constrained Tikhonov least squares problem. Although these results have been verified within the context of the analysis of electrochemical impedance spectroscopy, there is no reason to suppose that they would not be relevant within the broader framework of solving Fredholm integral equations for other applications.

[1]  Kendall E. Atkinson.  An introduction to numerical analysis , 1978 .

[2]  Rosemary A. Renaut,et al.  A Newton root-finding algorithm for estimating the regularization parameter for solving ill-conditioned least squares problems , 2009 .

[3]  Gene H. Golub,et al.  Matrix Computations, Third Edition , 1996 .

[4]  Per Christian Hansen,et al.  Regularization Tools version 4.0 for Matlab 7.3 , 2007, Numerical Algorithms.

[5]  P. Hansen,et al.  Exploiting Residual Information in the Parameter Choice for Discrete Ill-Posed Problems , 2006 .

[6]  Gene H. Golub,et al.  Matrix computations , 1983 .

[7]  Ellen Ivers-Tiffée,et al.  SOFC Modeling and Parameter Identification by Means of Impedance Spectroscopy , 2010 .

[8]  B. Rust.  Truncating the Singular Value Decomposition for Ill-Posed Problems , 1998 .

[9]  Ellen Ivers-Tiffée,et al.  Evaluation and Modeling of the Cell Resistance in Anode-Supported Solid Oxide Fuel Cells , 2008 .

[10]  B. Liu,et al.  Analysis of Impedance Spectra for Segmented-in-Series Tubular Solid Oxide Fuel Cells , 2010 .

[11]  C. Vogel.  Computational Methods for Inverse Problems , 1987 .

[12]  Stephen P. Boyd,et al.  Graph Implementations for Nonsmooth Convex Programs , 2008, Recent Advances in Learning and Control.

[13]  T. Matsui,et al.  Gas Transport Impedance in Segmented-in-Series Tubular Solid Oxide Fuel Cell , 2011 .

[14]  A. Matulis,et al.  Determination of the Distribution of the Relaxation Times from Dielectric Spectra , 2004 .

[15]  Begnaud Francis Hildebrand,et al.  Introduction to numerical analysis: 2nd edition , 1987 .

[16]  J. Navarro-Pedreño.  Numerical Methods for Least Squares Problems , 1996 .

[17]  J. Weese.  A reliable and fast method for the solution of Fredhol integral equations of the first kind based on Tikhonov regularization , 1992 .

[18]  Ellen Ivers-Tiffée,et al.  Combined Deconvolution and CNLS Fitting Approach Applied on the Impedance Response of Technical Ni ∕ 8YSZ Cermet Electrodes , 2008 .

[19]  Chris Chatfield,et al.  Introduction to Statistical Time Series. , 1976 .

[20]  J. Macdonald.  Exact and approximate nonlinear least‐squares inversion of dielectric relaxation spectra , 1995 .

[21]  P. Hansen.  Discrete Inverse Problems: Insight and Algorithms , 2010 .

[22]  H. Schichlein,et al.  Deconvolution of electrochemical impedance spectra for the identification of electrode reaction mechanisms in solid oxide fuel cells , 2002 .

[23]  E. Ivers-Tiffée,et al.  Impedance Study of Alternative ( La , Sr ) FeO3 − δ and ( La , Sr ) ( Co , Fe ) O3 − δ MIEC Cathode Compositions , 2010 .

[24]  Stefania Bellavia,et al.  An interior point Newton‐like method for non‐negative least‐squares problems with degenerate solution , 2006, Numer. Linear Algebra Appl..

[25]  Leigh C Ward,et al.  Determination of Cole parameters in multiple frequency bioelectrical impedance analysis using only the measurement of impedances , 2006, Physiological measurement.

[26]  Rosemary A. Renaut,et al.  Stability and error analysis of the polarization estimation inverse problem for microbial fuel cells , 2013 .

[27]  E. Barsoukov,et al.  Impedance spectroscopy : theory, experiment, and applications , 2005 .

[28]  P. Hansen.  Rank-Deficient and Discrete Ill-Posed Problems: Numerical Aspects of Linear Inversion , 1987 .

[29]  Dianne P. O'Leary,et al.  Residual periodograms for choosing regularization parameters for ill-posed problems , 2008 .

[30]  Matthias Hein,et al.  Non-negative least squares for high-dimensional linear models: consistency and sparse recovery without regularization , 2012, 1205.0953.

[31]  Jarom D. Hogue,et al.  Non-negatively constrained least squares and parameter choice by the residual periodogram for the inversion of dielectric relaxation spectra: Supplementary Materials , 2013, 1309.4506.

[32]  Charles L. Lawson,et al.  Solving least squares problems , 1976, Classics in applied mathematics.

[33]  F. Tietz,et al.  Time-Dependent Electrode Performance Changes in Intermediate Temperature Solid Oxide Fuel Cells , 2010 .

[34]  Kevin Barraclough,et al.  I and i , 2001, BMJ : British Medical Journal.

[35]  Inderjit S. Dhillon,et al.  A non-monotonic method for large-scale non-negative least squares , 2013, Optim. Methods Softw..

[36]  Edmond F. Bowden,et al.  The distribution of standard rate constants for electron transfer between thiol-modified gold electrodes and adsorbed cytochrome c , 1996 .