This paper mainly proposes a new C-XSC （C- for eXtended Scientific Computing） software for the symmetric single step method and relaxation method for computing an enclosure for the solution set and compares the meth...This paper mainly proposes a new C-XSC （C- for eXtended Scientific Computing） software for the symmetric single step method and relaxation method for computing an enclosure for the solution set and compares the methods with others＇ and then makes some modifications and finally, examples illustrating the applicability of the proposed methods are given.展开更多
LSQR, a Lanczos bidiagonalization based Krylov subspace iterative method, and its mathematically equivalent conjugate gradient for least squares problems(CGLS) applied to normal equations system, are commonly used for...LSQR, a Lanczos bidiagonalization based Krylov subspace iterative method, and its mathematically equivalent conjugate gradient for least squares problems(CGLS) applied to normal equations system, are commonly used for large-scale discrete ill-posed problems. It is well known that LSQR and CGLS have regularizing effects, where the number of iterations plays the role of the regularization parameter. However, it has long been unknown whether the regularizing effects are good enough to find best possible regularized solutions. Here a best possible regularized solution means that it is at least as accurate as the best regularized solution obtained by the truncated singular value decomposition(TSVD) method. We establish bounds for the distance between the k-dimensional Krylov subspace and the k-dimensional dominant right singular space. They show that the Krylov subspace captures the dominant right singular space better for severely and moderately ill-posed problems than for mildly ill-posed problems. Our general conclusions are that LSQR has better regularizing effects for the first two kinds of problems than for the third kind, and a hybrid LSQR with additional regularization is generally needed for mildly ill-posed problems. Exploiting the established bounds, we derive an estimate for the accuracy of the rank k approximation generated by Lanczos bidiagonalization. Numerical experiments illustrate that the regularizing effects of LSQR are good enough to compute best possible regularized solutions for severely and moderately ill-posed problems, stronger than our theory, but they are not for mildly ill-posed problems and additional regularization is needed.展开更多
Parameter estimation of the 2 R-1 C model is usually performed using iterative methods that require high-performance processing units.Consequently,there is a strong motivation to develop less time-consuming and more p...Parameter estimation of the 2 R-1 C model is usually performed using iterative methods that require high-performance processing units.Consequently,there is a strong motivation to develop less time-consuming and more power-efficient parameter estimation methods.Such low-complexity algorithms would be suitable for implementation in portable microcontroller-based devices.In this study,we propose the quadratic interpolation non-iterative parameter estimation(QINIPE)method,based on quadratic interpolation of the imaginary part of the measured impedance,which enables more accurate estimation of the characteristic frequency.The 2 R-1 C model parameters are subsequently calculated from the real and imaginary parts of the measured impedance using a set of closed-form expressions.Comparative analysis conducted on the impedance data of the 2 R-1 C model obtained in both simulation and measurements shows that the proposed QINIPE method reduces the number of required measurement points by 80%in comparison with our previously reported non-iterative parameter estimation(NIPE)method,while keeping the relative estimation error to less than 1%for all estimated parameters.Both non-iterative methods are implemented on a microcontroller-based device;the estimation accuracy,RAM,flash memory usage,and execution time are monitored.Experiments show that the QINIPE method slightly increases the execution time by 0.576 ms(about 6.7%),and requires 24%(1.2 KB)more flash memory and just 2.4%(32 bytes)more RAM in comparison to the NIPE method.However,the impedance root mean square errors(RMSEs)of the QINIPE method are decreased to 42.8%(for the real part)and 64.5%(for the imaginary part)of the corresponding RMSEs obtained using the NIPE method.Moreover,we compared the QINIPE and the complex nonlinear least squares(CNLS)estimation of the 2 R-1 C model parameters.The results obtained show that although the estimation accuracy of the QINIPE is somewhat lower than the estimation accuracy of the CNLS,it is still satisfactory for many practical purposes and its execution time reduces to1/45–1/30.展开更多
文摘This paper mainly proposes a new C-XSC （C- for eXtended Scientific Computing） software for the symmetric single step method and relaxation method for computing an enclosure for the solution set and compares the methods with others＇ and then makes some modifications and finally, examples illustrating the applicability of the proposed methods are given.
基金supported by National Basic Research Program of China (Grant No. 2011CB302400)National Natural Science Foundation of China (Grant No. 11371219)
文摘LSQR, a Lanczos bidiagonalization based Krylov subspace iterative method, and its mathematically equivalent conjugate gradient for least squares problems(CGLS) applied to normal equations system, are commonly used for large-scale discrete ill-posed problems. It is well known that LSQR and CGLS have regularizing effects, where the number of iterations plays the role of the regularization parameter. However, it has long been unknown whether the regularizing effects are good enough to find best possible regularized solutions. Here a best possible regularized solution means that it is at least as accurate as the best regularized solution obtained by the truncated singular value decomposition(TSVD) method. We establish bounds for the distance between the k-dimensional Krylov subspace and the k-dimensional dominant right singular space. They show that the Krylov subspace captures the dominant right singular space better for severely and moderately ill-posed problems than for mildly ill-posed problems. Our general conclusions are that LSQR has better regularizing effects for the first two kinds of problems than for the third kind, and a hybrid LSQR with additional regularization is generally needed for mildly ill-posed problems. Exploiting the established bounds, we derive an estimate for the accuracy of the rank k approximation generated by Lanczos bidiagonalization. Numerical experiments illustrate that the regularizing effects of LSQR are good enough to compute best possible regularized solutions for severely and moderately ill-posed problems, stronger than our theory, but they are not for mildly ill-posed problems and additional regularization is needed.
基金Project supported by the Ministry of Science and Technology of the Republic of Srpska(No.19/6-020/961-143/18)the EU’s H2020 MSCA MEDLEM(No.690876).
文摘Parameter estimation of the 2 R-1 C model is usually performed using iterative methods that require high-performance processing units.Consequently,there is a strong motivation to develop less time-consuming and more power-efficient parameter estimation methods.Such low-complexity algorithms would be suitable for implementation in portable microcontroller-based devices.In this study,we propose the quadratic interpolation non-iterative parameter estimation(QINIPE)method,based on quadratic interpolation of the imaginary part of the measured impedance,which enables more accurate estimation of the characteristic frequency.The 2 R-1 C model parameters are subsequently calculated from the real and imaginary parts of the measured impedance using a set of closed-form expressions.Comparative analysis conducted on the impedance data of the 2 R-1 C model obtained in both simulation and measurements shows that the proposed QINIPE method reduces the number of required measurement points by 80%in comparison with our previously reported non-iterative parameter estimation(NIPE)method,while keeping the relative estimation error to less than 1%for all estimated parameters.Both non-iterative methods are implemented on a microcontroller-based device;the estimation accuracy,RAM,flash memory usage,and execution time are monitored.Experiments show that the QINIPE method slightly increases the execution time by 0.576 ms(about 6.7%),and requires 24%(1.2 KB)more flash memory and just 2.4%(32 bytes)more RAM in comparison to the NIPE method.However,the impedance root mean square errors(RMSEs)of the QINIPE method are decreased to 42.8%(for the real part)and 64.5%(for the imaginary part)of the corresponding RMSEs obtained using the NIPE method.Moreover,we compared the QINIPE and the complex nonlinear least squares(CNLS)estimation of the 2 R-1 C model parameters.The results obtained show that although the estimation accuracy of the QINIPE is somewhat lower than the estimation accuracy of the CNLS,it is still satisfactory for many practical purposes and its execution time reduces to1/45–1/30.
