A necessary and sufficient condition for the existence of simultaneous (M,N)singular value decomposition of matrices is given.Some properties about the weighted partial ordering are discussed with the help of the deco...A necessary and sufficient condition for the existence of simultaneous (M,N)singular value decomposition of matrices is given.Some properties about the weighted partial ordering are discussed with the help of the decomposition.展开更多
In this paper we present a precise integration method based on high order multiple perturbation method and reduction method for solving a class of singular twopoint boundary value problems.Firstly,by employing the met...In this paper we present a precise integration method based on high order multiple perturbation method and reduction method for solving a class of singular twopoint boundary value problems.Firstly,by employing the method of variable coefficient dimensional expanding,the non-homogeneous ordinary differential equations(ODEs) are transformed into homogeneous ODEs.Then the interval is divided evenly,and the transfer matrix in every subinterval is worked out using the high order multiple perturbation method,and a set of algebraic equations is given in the form of matrix by the precise integration relation for each segment,which is worked out by the reduction method.Finally numerical examples are elaboratedd to validate the present method.展开更多
We have discussed and solved the boundary value problem with period 2aπ and the singular integral equation with kernel csc t-tv/a in solution having singularities of high order, where the smooth contour of integratio...We have discussed and solved the boundary value problem with period 2aπ and the singular integral equation with kernel csc t-tv/a in solution having singularities of high order, where the smooth contour of integration is in the strip 0<Rez<aπ.展开更多
This paper presents a high order multiplication perturbation method for sin- gularly perturbed two-point boundary value problems with the boundary layer at one end. By the theory of singular perturbations, the singula...This paper presents a high order multiplication perturbation method for sin- gularly perturbed two-point boundary value problems with the boundary layer at one end. By the theory of singular perturbations, the singularly perturbed two-point boundary value problems are first transformed into the singularly perturbed initial value problems. With the variable coefficient dimensional expanding, the non-homogeneous ordinary dif- ferential equations (ODEs) are transformed into the homogeneous ODEs, which are then solved by the high order multiplication perturbation method. Some linear and nonlinear numerical examples show that the proposed method has high precision.展开更多
A uniform high-order method is. presented for the numerical solution of a singular perturbation problem in conservative form. We firest replace the original second-order problem (1.1) by two equivalent first-order pro...A uniform high-order method is. presented for the numerical solution of a singular perturbation problem in conservative form. We firest replace the original second-order problem (1.1) by two equivalent first-order problems ( 1.4), i.e., the solution of (1.1) is a linear combination of the solutions of (1.4). Then we derive a uniformly O (hm+1) accurate scheme for the first-order problems (1.4), where m is an arbitrary nonnegative integer, so we can get a uniformly O (hm+1) accurate solution of the original problem (1.1) by relation (1.3). Some illustrative numerical results are also given.展开更多
Higher-order singular value decomposition (HOSVD) is an efficient way for data reduction and also eliciting intrinsic structure of multi-dimensional array data. It has been used in many applications, and some of the...Higher-order singular value decomposition (HOSVD) is an efficient way for data reduction and also eliciting intrinsic structure of multi-dimensional array data. It has been used in many applications, and some of them involve incomplete data. To obtain HOSVD of the data with missing values, one can first impute the missing entries through a certain tensor completion method and then perform HOSVD to the reconstructed data. However, the two-step procedure can be inefficient and does not make reliable decomposition. In this paper, we formulate an incomplete HOSVD problem and combine the two steps into solving a single optimization problem, which simultaneously achieves imputation of missing values and also tensor decomposition. We also present one algorithm for solving the problem based on block coordinate update (BCU). Global convergence of the algorithm is shown under mild assumptions and implies that of the popular higher-order orthogonality iteration (HOOI) method, and thus we, for the first time, give global convergence of HOOI. In addition, we compare the proposed method to state-of-the-art ones for solving incom- plete HOSVD and also low-rank tensor completion problems and demonstrate the superior performance of our method over other compared ones. Furthermore, we apply it to face recognition and MRI image reconstruction to show its practical performance.展开更多
The primary goal in the analysis of hierarchical distributed monitoring and control architectures is to study the spatiotemporal patterns of the interactions between areas or subsystems.In this paper,a novel conceptua...The primary goal in the analysis of hierarchical distributed monitoring and control architectures is to study the spatiotemporal patterns of the interactions between areas or subsystems.In this paper,a novel conceptual framework for distributed monitoring of power system oscillations using multiblock principal component analysis(MB-PCA)and higher-order singular value decomposition(HOSVD)is proposed to understand,characterize,and visualize the global behavior of the power system.The proposed framework can be used to evaluate the influence of a given area or utility on the oscillatory behavior,uncover low-dimensional structures from high-dimensional data,and analyze the effects of heterogeneous data on the modal characteristics and interpretation of power system.The metrics are then investigated to examine the relationships between the dynamic patterns and participation of individual data blocks in the global behavior of the system.Practical application of these techniques is demonstrated by case studies of two systems:a 14-machine test system and a 5449-bus 635-generator equivalent model of a large power system.展开更多
This paper proposes an extension of the algorithm in [1], as well as utilization of the wavelet transform in event detection, including High Impedance Fault (HIF). Techniques to analyze the abundant data of PMUs quick...This paper proposes an extension of the algorithm in [1], as well as utilization of the wavelet transform in event detection, including High Impedance Fault (HIF). Techniques to analyze the abundant data of PMUs quickly and effectively are paramount to increasing response time to events and unstable parameters. With the amount of data PMUs output, unstable parameters, tie line oscillations, and HIFs are often overlooked in the bulk of the data. This paper explores model-free techniques to attain stability information and determine events in real-time. When full system connectivity is unknown, many traditional methods requiring other bus measurements can be impossible or computationally extensive to apply. The traditional method of interest is analyzing the power flow Jacobian for singularities and system weak points, attained by applying singular value decomposition. This paper further develops upon the approach in [1] to expand the Discrete-Time Jacobian Eigenvalue Approximation (DDJEA), giving values to significant off-diagonal terms while establishing a generalized connectivity between correlated buses. Statistical linear models are applied over large data sets to prove significance to each term. Then the off diagonal terms are given time-varying weights to account for changes in topology or sensitivity to events using a reduced system model. The results of this novel method are compared to the present errors of the previous publication in order to quantify the degree of improvement that this novel method imposes. The effective bus eigenvalues are briefly compared to Prony analysis to check similarities. An additional application for biorthogonal wavelets is also introduced to detect event types, including the HIF, for PMU data.展开更多
A type of high-order integral observers for matrix second-order linear systems is proposed on the basis of generalized eigenstructure assignment via unified parametric approaches. Through establishing two general para...A type of high-order integral observers for matrix second-order linear systems is proposed on the basis of generalized eigenstructure assignment via unified parametric approaches. Through establishing two general parametric solutions to this type of generalized matrix second-order Sylvester matrix equations, two unified complete parametric methods for the proposed observer design problem are presented. Both methods give simple complete parametric expressions for the observer gain matrices. The first one mainly depends on a series of singular value decompositions, and is thus numerically simple and reliable; the second one utilizes the fight factorization of the system, and allows eigenvalues of the error system to be set undetermined and sought via certain optimization procedures. A spring-mass-dashpot system is utilized to illustrate the design procedure and show the effect of the proposed approach.展开更多
In practical applications, we often have to deal with high-order data, for example, a grayscale image and a video clip are intrinsically a 2nd-order tensor and a 3rd-order tensor, respectively. In order to satisty the...In practical applications, we often have to deal with high-order data, for example, a grayscale image and a video clip are intrinsically a 2nd-order tensor and a 3rd-order tensor, respectively. In order to satisty these high-order data, it is conventional to vectorize these data in advance, which often destroys the intrinsic structures of the data and includes the curse of dimensionality. For this reason, we consider the problem of high-order data representation and classification, and propose a tensor based fisher discriminant analysis (FDA), which is a generalized version of FDA, named as GFDA. Experimental results show our GFDA outperforms the existing methods, such as the 2-directional 2-dimensional principal component analysis ((2D)2pCA), 2-directional 2-dimensional linear discriminant analysis ((2D)2LDA), and multilinear discriminant analysis (MDA), in high-order data classification under a lower compression ratio.展开更多
为在大数据环境下处理高维矩阵和应用奇异值分解提供更高效的解决方案,从而加速数据分析和处理速度,通过研究随机投影以及Krylov子空间投影理论下关于高维矩阵求解特征值特征向量(奇异值奇异向量)问题,分别总结了6种高效计算方法并对其...为在大数据环境下处理高维矩阵和应用奇异值分解提供更高效的解决方案,从而加速数据分析和处理速度,通过研究随机投影以及Krylov子空间投影理论下关于高维矩阵求解特征值特征向量(奇异值奇异向量)问题,分别总结了6种高效计算方法并对其相关应用研究进行对比分析。结果表明,在谱聚类的应用上,通过降低核心步骤SVD(Singular Value Decomposition)的复杂度,使优化后的算法与原始谱聚类算法的精度相近,但大大缩短了运行时间,在1200维的数据下计算速度相较原算法快了10倍以上。同时,该方法应用于图像压缩领域,能有效地提高原有算法的运行效率,在精度不变的情况下,运行效率得到了1~5倍的提升。展开更多
基于高阶奇异值分解(High Order Singular Value Decomposition,HOSVD)降噪的信道预测算法对天线数较少引起的秩不足问题比较敏感,同时也难以应付较大多普勒频移的情况,从而引起信道估计性能和预测性能的急剧下降、损失信道容量.针对这...基于高阶奇异值分解(High Order Singular Value Decomposition,HOSVD)降噪的信道预测算法对天线数较少引起的秩不足问题比较敏感,同时也难以应付较大多普勒频移的情况,从而引起信道估计性能和预测性能的急剧下降、损失信道容量.针对这一问题,提出了一种改进的使用HOSVD降噪的信道预测算法.该算法先利用多输入多输出(Multiple-input Multiple-Output,MIMO)信道固有的空时相关性对采样得到的信道状态信息(Channel State Information,CSI)进行矩阵重排和数据平滑处理,随后基于信道的多维结构特性,使用HOSVD降低噪声的影响,继而重构信道矩阵,最后利用递归最小二乘滤波器对未来时刻的信道状态进行预测.仿真表明,所提算法的估计误差和预测误差性能均明显优于对比算法,这是因为所提算法通过矩阵重排和空时平滑,虚拟地增加了天线数,降低了秩缺失问题对估计和预测精度的影响,从而有效补偿了因误差所致的信道容量的损失.同时,对比天线数和多普勒频移对不同算法性能的影响可见,所提算法也能在大多普勒频移和天线数较少等不利条件下提供较好预测性能和信道容量,具有一定的优越性.展开更多
基金The Guangxi Science Foundation(0575032,06400161)the support program for 100 Young and Middle-aged Disciplinary Leaders in Guangxi Higher Education Institutions
文摘A necessary and sufficient condition for the existence of simultaneous (M,N)singular value decomposition of matrices is given.Some properties about the weighted partial ordering are discussed with the help of the decomposition.
基金supported by the National Natural Science Foundation of China (11132004 and 51078145)the Natural Science Foundation of Guangdong Province (9251064101000016)
文摘In this paper we present a precise integration method based on high order multiple perturbation method and reduction method for solving a class of singular twopoint boundary value problems.Firstly,by employing the method of variable coefficient dimensional expanding,the non-homogeneous ordinary differential equations(ODEs) are transformed into homogeneous ODEs.Then the interval is divided evenly,and the transfer matrix in every subinterval is worked out using the high order multiple perturbation method,and a set of algebraic equations is given in the form of matrix by the precise integration relation for each segment,which is worked out by the reduction method.Finally numerical examples are elaboratedd to validate the present method.
基金Supported by the National Natural Science Foundation of China(19971064 10161009)
文摘We have discussed and solved the boundary value problem with period 2aπ and the singular integral equation with kernel csc t-tv/a in solution having singularities of high order, where the smooth contour of integration is in the strip 0<Rez<aπ.
基金supported by the National Natural Science Foundation of China(Key Program)(Nos.11132004 and 51078145)
文摘This paper presents a high order multiplication perturbation method for sin- gularly perturbed two-point boundary value problems with the boundary layer at one end. By the theory of singular perturbations, the singularly perturbed two-point boundary value problems are first transformed into the singularly perturbed initial value problems. With the variable coefficient dimensional expanding, the non-homogeneous ordinary dif- ferential equations (ODEs) are transformed into the homogeneous ODEs, which are then solved by the high order multiplication perturbation method. Some linear and nonlinear numerical examples show that the proposed method has high precision.
文摘A uniform high-order method is. presented for the numerical solution of a singular perturbation problem in conservative form. We firest replace the original second-order problem (1.1) by two equivalent first-order problems ( 1.4), i.e., the solution of (1.1) is a linear combination of the solutions of (1.4). Then we derive a uniformly O (hm+1) accurate scheme for the first-order problems (1.4), where m is an arbitrary nonnegative integer, so we can get a uniformly O (hm+1) accurate solution of the original problem (1.1) by relation (1.3). Some illustrative numerical results are also given.
文摘Higher-order singular value decomposition (HOSVD) is an efficient way for data reduction and also eliciting intrinsic structure of multi-dimensional array data. It has been used in many applications, and some of them involve incomplete data. To obtain HOSVD of the data with missing values, one can first impute the missing entries through a certain tensor completion method and then perform HOSVD to the reconstructed data. However, the two-step procedure can be inefficient and does not make reliable decomposition. In this paper, we formulate an incomplete HOSVD problem and combine the two steps into solving a single optimization problem, which simultaneously achieves imputation of missing values and also tensor decomposition. We also present one algorithm for solving the problem based on block coordinate update (BCU). Global convergence of the algorithm is shown under mild assumptions and implies that of the popular higher-order orthogonality iteration (HOOI) method, and thus we, for the first time, give global convergence of HOOI. In addition, we compare the proposed method to state-of-the-art ones for solving incom- plete HOSVD and also low-rank tensor completion problems and demonstrate the superior performance of our method over other compared ones. Furthermore, we apply it to face recognition and MRI image reconstruction to show its practical performance.
文摘The primary goal in the analysis of hierarchical distributed monitoring and control architectures is to study the spatiotemporal patterns of the interactions between areas or subsystems.In this paper,a novel conceptual framework for distributed monitoring of power system oscillations using multiblock principal component analysis(MB-PCA)and higher-order singular value decomposition(HOSVD)is proposed to understand,characterize,and visualize the global behavior of the power system.The proposed framework can be used to evaluate the influence of a given area or utility on the oscillatory behavior,uncover low-dimensional structures from high-dimensional data,and analyze the effects of heterogeneous data on the modal characteristics and interpretation of power system.The metrics are then investigated to examine the relationships between the dynamic patterns and participation of individual data blocks in the global behavior of the system.Practical application of these techniques is demonstrated by case studies of two systems:a 14-machine test system and a 5449-bus 635-generator equivalent model of a large power system.
文摘This paper proposes an extension of the algorithm in [1], as well as utilization of the wavelet transform in event detection, including High Impedance Fault (HIF). Techniques to analyze the abundant data of PMUs quickly and effectively are paramount to increasing response time to events and unstable parameters. With the amount of data PMUs output, unstable parameters, tie line oscillations, and HIFs are often overlooked in the bulk of the data. This paper explores model-free techniques to attain stability information and determine events in real-time. When full system connectivity is unknown, many traditional methods requiring other bus measurements can be impossible or computationally extensive to apply. The traditional method of interest is analyzing the power flow Jacobian for singularities and system weak points, attained by applying singular value decomposition. This paper further develops upon the approach in [1] to expand the Discrete-Time Jacobian Eigenvalue Approximation (DDJEA), giving values to significant off-diagonal terms while establishing a generalized connectivity between correlated buses. Statistical linear models are applied over large data sets to prove significance to each term. Then the off diagonal terms are given time-varying weights to account for changes in topology or sensitivity to events using a reduced system model. The results of this novel method are compared to the present errors of the previous publication in order to quantify the degree of improvement that this novel method imposes. The effective bus eigenvalues are briefly compared to Prony analysis to check similarities. An additional application for biorthogonal wavelets is also introduced to detect event types, including the HIF, for PMU data.
基金This work was supported by the Chinese National Natural Science Foundation ( No. 69925308).
文摘A type of high-order integral observers for matrix second-order linear systems is proposed on the basis of generalized eigenstructure assignment via unified parametric approaches. Through establishing two general parametric solutions to this type of generalized matrix second-order Sylvester matrix equations, two unified complete parametric methods for the proposed observer design problem are presented. Both methods give simple complete parametric expressions for the observer gain matrices. The first one mainly depends on a series of singular value decompositions, and is thus numerically simple and reliable; the second one utilizes the fight factorization of the system, and allows eigenvalues of the error system to be set undetermined and sought via certain optimization procedures. A spring-mass-dashpot system is utilized to illustrate the design procedure and show the effect of the proposed approach.
文摘In practical applications, we often have to deal with high-order data, for example, a grayscale image and a video clip are intrinsically a 2nd-order tensor and a 3rd-order tensor, respectively. In order to satisty these high-order data, it is conventional to vectorize these data in advance, which often destroys the intrinsic structures of the data and includes the curse of dimensionality. For this reason, we consider the problem of high-order data representation and classification, and propose a tensor based fisher discriminant analysis (FDA), which is a generalized version of FDA, named as GFDA. Experimental results show our GFDA outperforms the existing methods, such as the 2-directional 2-dimensional principal component analysis ((2D)2pCA), 2-directional 2-dimensional linear discriminant analysis ((2D)2LDA), and multilinear discriminant analysis (MDA), in high-order data classification under a lower compression ratio.
文摘为在大数据环境下处理高维矩阵和应用奇异值分解提供更高效的解决方案,从而加速数据分析和处理速度,通过研究随机投影以及Krylov子空间投影理论下关于高维矩阵求解特征值特征向量(奇异值奇异向量)问题,分别总结了6种高效计算方法并对其相关应用研究进行对比分析。结果表明,在谱聚类的应用上,通过降低核心步骤SVD(Singular Value Decomposition)的复杂度,使优化后的算法与原始谱聚类算法的精度相近,但大大缩短了运行时间,在1200维的数据下计算速度相较原算法快了10倍以上。同时,该方法应用于图像压缩领域,能有效地提高原有算法的运行效率,在精度不变的情况下,运行效率得到了1~5倍的提升。
文摘基于高阶奇异值分解(High Order Singular Value Decomposition,HOSVD)降噪的信道预测算法对天线数较少引起的秩不足问题比较敏感,同时也难以应付较大多普勒频移的情况,从而引起信道估计性能和预测性能的急剧下降、损失信道容量.针对这一问题,提出了一种改进的使用HOSVD降噪的信道预测算法.该算法先利用多输入多输出(Multiple-input Multiple-Output,MIMO)信道固有的空时相关性对采样得到的信道状态信息(Channel State Information,CSI)进行矩阵重排和数据平滑处理,随后基于信道的多维结构特性,使用HOSVD降低噪声的影响,继而重构信道矩阵,最后利用递归最小二乘滤波器对未来时刻的信道状态进行预测.仿真表明,所提算法的估计误差和预测误差性能均明显优于对比算法,这是因为所提算法通过矩阵重排和空时平滑,虚拟地增加了天线数,降低了秩缺失问题对估计和预测精度的影响,从而有效补偿了因误差所致的信道容量的损失.同时,对比天线数和多普勒频移对不同算法性能的影响可见,所提算法也能在大多普勒频移和天线数较少等不利条件下提供较好预测性能和信道容量,具有一定的优越性.