scholarly journals Scalability of k-Tridiagonal Matrix Singular Value Decomposition

Mathematics ◽  
2021 ◽  
Vol 9 (23) ◽  
pp. 3123
Author(s):  
Andrei Tănăsescu ◽  
Mihai Carabaş ◽  
Florin Pop ◽  
Pantelimon George Popescu

Singular value decomposition has recently seen a great theoretical improvement for k-tridiagonal matrices, obtaining a considerable speed up over all previous implementations, but at the cost of not ordering the singular values. We provide here a refinement of this method, proving that reordering singular values does not affect performance. We complement our refinement with a scalability study on a real physical cluster setup, offering surprising results. Thus, this method provides a major step up over standard industry implementations.

2012 ◽  
Vol 2012 ◽  
pp. 1-20 ◽  
Author(s):  
Muhammad Mohsin Riaz ◽  
Abdul Ghafoor

Singular value decomposition and information theoretic criterion-based image enhancement is proposed for through-wall imaging. The scheme is capable of discriminating target, clutter, and noise subspaces. Information theoretic criterion is used with conventional singular value decomposition to find number of target singular values. Furthermore, wavelet transform-based denoising is performed (to further suppress noise signals) by estimating noise variance. Proposed scheme works also for extracting multiple targets in heavy cluttered through-wall images. Simulation results are compared on the basis of mean square error, peak signal to noise ratio, and visual inspection.


2019 ◽  
Vol 22 (12) ◽  
pp. 2687-2698 ◽  
Author(s):  
Zhen Chen ◽  
Lifeng Qin ◽  
Shunbo Zhao ◽  
Tommy HT Chan ◽  
Andy Nguyen

This article introduces and evaluates the piecewise polynomial truncated singular value decomposition algorithm toward an effective use for moving force identification. Suffering from numerical non-uniqueness and noise disturbance, the moving force identification is known to be associated with ill-posedness. An important method for solving this problem is the truncated singular value decomposition algorithm, but the truncated small singular values removed by truncated singular value decomposition may contain some useful information. The piecewise polynomial truncated singular value decomposition algorithm extracts the useful responses from truncated small singular values and superposes it into the solution of truncated singular value decomposition, which can be useful in moving force identification. In this article, a comprehensive numerical simulation is set up to evaluate piecewise polynomial truncated singular value decomposition, and compare this technique against truncated singular value decomposition and singular value decomposition. Numerically simulated data are processed to validate the novel method, which show that regularization matrix [Formula: see text] and truncating point [Formula: see text] are the two most important governing factors affecting identification accuracy and ill-posedness immunity of piecewise polynomial truncated singular value decomposition.


Geophysics ◽  
1993 ◽  
Vol 58 (11) ◽  
pp. 1655-1661 ◽  
Author(s):  
Reinaldo J. Michelena

I perform singular value decomposition (SVD) on the matrices that result in tomographic velocity estimation from cross‐well traveltimes in isotropic and anisotropic media. The slowness model is parameterized in four ways: One‐dimensional (1-D) isotropic, 1-D anisotropic, two‐dimensional (2-D) isotropic, and 2-D anisotropic. The singular value distribution is different for the different parameterizations. One‐dimensional isotropic models can be resolved well but the resolution of the data is poor. One‐dimensional anisotropic models can also be resolved well except for some variations in the vertical component of the slowness that are not sensitive to the data. In 2-D isotropic models, “pure” lateral variations are not sensitive to the data, and when anisotropy is introduced, the result is that the horizontal and vertical component of the slowness cannot be estimated with the same spatial resolution because the null space is mostly related to horizontal and high frequency variations in the vertical component of the slowness. Since the distribution of singular values varies depending on the parametrization used, the effect of conventional regularization procedures in the final solution may also vary. When the model is isotropic, regularization translates into smoothness, and when the model is anisotropic regularization not only smooths but may also alter the anisotropy in the solution.


2009 ◽  
Vol 09 (03) ◽  
pp. 449-477 ◽  
Author(s):  
GAURAV BHATNAGAR ◽  
BALASUBRAMANIAN RAMAN

This paper presents a new robust reference watermarking scheme based on wavelet packet transform (WPT) and bidiagonal singular value decomposition (bSVD) for copyright protection and authenticity. A small gray scale logo is used as watermark instead of randomly generated Gaussian noise type watermark. A reference watermark is generated by original watermark and the process of embedding is done in wavelet packet domain by modifying the bidiagonal singular values. For the robustness and imperceptibly, watermark is embedded in the selected sub-bands, which are selected by taking into account the variance of the sub-bands, which serves as a measure of the watermark magnitude that could be imperceptibly embedded in each block. For this purpose, the variance is calculated in a small moving square window of size Sp× Sp(typically 3 × 3 or 5 × 5 window) centered at the pixel. A reliable watermark extraction is developed, in which the watermark bidiagonal singular values are extracted by considering the distortion caused by the attacks in neighboring bidiagonal singular values. Experimental evaluation demonstrates that the proposed scheme is able to withstand a variety of attacks and the superiority of the proposed method is carried out by the comparison which is made by us with the existing methods.


2008 ◽  
Vol 2008 ◽  
pp. 1-17 ◽  
Author(s):  
Eric J. Beh

The correspondence analysis of a two-way contingency table is now accepted as a very versatile tool for helping users to understand the structure of the association in their data. In cases where the variables consist of ordered categories, there are a number of approaches that can be employed and these generally involve an adaptation of singular value decomposition. Over the last few years, an alternative decomposition method has been used for cases where the row and column variables of a two-way contingency table have an ordinal structure. A version of this approach is also available for a two-way table where one variable has a nominal structure and the other variable has an ordinal structure. However, such an approach does not take into consideration the presence of the nominal variable. This paper explores an approach to correspondence analysis using an amalgamation of singular value decomposition and bivariate moment decomposition. A benefit of this technique is that it combines the classical technique with the ordinal analysis by determining the structure of the variables in terms of singular values and location, dispersion and higher-order moments.


Geophysics ◽  
2018 ◽  
Vol 83 (4) ◽  
pp. G25-G34 ◽  
Author(s):  
Saeed Vatankhah ◽  
Rosemary Anne Renaut ◽  
Vahid Ebrahimzadeh Ardestani

We develop a fast algorithm for solving the under-determined 3D linear gravity inverse problem based on randomized singular-value decomposition (RSVD). The algorithm combines an iteratively reweighted approach for [Formula: see text]-norm regularization with the RSVD methodology in which the large-scale linear system at each iteration is replaced with a much smaller linear system. Although the optimal choice for the low-rank approximation of the system matrix with [Formula: see text] rows is [Formula: see text], acceptable results are achievable with [Formula: see text]. In contrast to the use of the iterative LSQR algorithm for the solution of linear systems at each iteration, the singular values generated using RSVD yield a good approximation of the dominant singular values of the large-scale system matrix. Thus, the regularization parameter found for the small system at each iteration is dependent on the dominant singular values of the large-scale system matrix and appropriately regularizes the dominant singular space of the large-scale problem. The results achieved are comparable with those obtained using the LSQR algorithm for solving each linear system, but they are obtained at a reduced computational cost. The method has been tested on synthetic models along with real gravity data from the Morro do Engenho complex in central Brazil.


Sensors ◽  
2019 ◽  
Vol 19 (15) ◽  
pp. 3301 ◽  
Author(s):  
Liang Wu ◽  
Qian Xu ◽  
Janne Heikkilä ◽  
Zijun Zhao ◽  
Liwei Liu ◽  
...  

The navigation accuracy of a star sensor depends on the estimation accuracy of its optical parameters, and so, the parameters should be updated in real time to obtain the best performance. Current on-orbit calibration methods for star sensors mainly rely on the angular distance between stars, and few studies have been devoted to seeking new calibration references. In this paper, an on-orbit calibration method using singular values as the calibration reference is introduced and studied. Firstly, the camera model of the star sensor is presented. Then, on the basis of the invariance of the singular values under coordinate transformation, an on-orbit calibration method based on the singular-value decomposition (SVD) method is proposed. By means of observability analysis, an optimal model of the star combinations for calibration is explored. According to the physical interpretation of the singular-value decomposition of the star vector matrix, the singular-value selection for calibration is discussed. Finally, to demonstrate the performance of the SVD method, simulation calibrations are conducted by both the SVD method and the conventional angular distance-based method. The results show that the accuracy and convergence speed of both methods are similar; however, the computational cost of the SVD method is heavily reduced. Furthermore, a field experiment is conducted to verify the feasibility of the SVD method. Therefore, the SVD method performs well in the calibration of star sensors, and in particular, it is suitable for star sensors with limited computing resources.


2020 ◽  
Vol 28 (1) ◽  
pp. 105-118 ◽  
Author(s):  
Jodi L. Mead ◽  
James F. Ford

AbstractJoint inversion of multiple data types was studied as early as 1975 in [K. Vozoff and D. L. Jupp, Joint inversion of geophysical data, Geophys. J. Internat. 42 1975, 3, 977–991], where the authors used the singular value decomposition to determine the degree of ill-conditioning of joint inverse problems. The authors demonstrated in several examples that combining two physical models in a joint inversion, and by effectively stacking discrete linear models, improved the conditioning as compared to individual inversions. This work extends the notion of using the singular value decomposition to determine the conditioning of discrete joint inversion to using the singular value expansion to determine the well-posedness of joint operators. We provide a convergent technique for approximating the singular values of continuous joint operators. In the case of self-adjoint operators, we give an algebraic expression for the joint singular values in terms of the singular values of the individual operators. This expression allows us to show that while rare, there are situations where ill-posedness may be not improved through joint inversion and in fact can degrade the conditioning of an individual inversion. The expression also quantifies the benefits of including repeated measurements in an inversion. We give an example of joint inversion with two moderately ill-posed Green’s function solutions, and quantify the improvement over individual inversions. This work provides a framework in which to identify data types that are advantageous to combine in a joint inversion.


2012 ◽  
Vol 4 (3) ◽  
Author(s):  
Lei Cui ◽  
Jian S. Dai

With a new type of multifingered hands that raise a new philosophy in the construction and study of a multifingered hand, this paper is a follow-on study of the kinematics of the metamorphic multifingered hand based on finger constraint equations. The finger constraint equations lead to a comprehensive mathematical model of the hand with a reconfigurable palm which integrates all finger motions with the additional palm motion. Singular values of the partitioned Jacobian matrix in their analytical form are derived and applied to obtaining analytical solution to inverse kinematics of a complete robotic hand. The paper for the first time solves this integrated motion and the multifingered hand model with the singular value decomposition and extra degrees of freedom are examined with the singular value analysis to avoid the singularities. The work identifies finger displacement and velocity with effect from the articulated palm and presents a new way of analyzing a multifingered robotic hand.


Sign in / Sign up

Export Citation Format

Share Document