Recent advances in experimental neuroscience allow non-invasive studies of the white matter tracts in the human central nervous system, thus making available cutting-edge brain anatomical data describing these global ...
详细信息
Recent advances in experimental neuroscience allow non-invasive studies of the white matter tracts in the human central nervous system, thus making available cutting-edge brain anatomical data describing these global connectivity patterns. Through magnetic resonance imaging, this non-invasive technique is able to infer a snapshot of the cortical network within the living human brain. Here, we report on the initial success of a new weighted network communicability measure in distinguishing local and global differences between diseased patients and controls. This approach builds on recent advances in network science, where an underlying connectivity structure is used as a means to measure the ease with which information can flow between nodes. One advantage of our method is that it deals directly with the real-valued connectivity data, thereby avoiding the need to discretize the corresponding adjacency matrix, i.e. to round weights up to 1 or down to 0, depending upon some threshold value. Experimental results indicate that the new approach is able to extract biologically relevant features that are not immediately apparent from the raw connectivity data.
In this paper we consider the numerical solution of fractional differential equations by means of m-step recursions. The construction of such formulas can be obtained in many ways. Here we study a technique based on t...
详细信息
In this paper we consider the numerical solution of fractional differential equations by means of m-step recursions. The construction of such formulas can be obtained in many ways. Here we study a technique based on the rational approximation of the generating functions of fractional backward differentiation formulas (FBDFs). Accurate approximations lead to the definition of methods which simulate the underlying FBDF, with important computational advantages. Numerical experiments are presented.
The computation of matrix functions using quadrature formulas and rational approximations of very large structured matrices using tensor trains (TT), and quantized tensor trains (QTT) is considered here. The focus is ...
详细信息
The computation of matrix functions using quadrature formulas and rational approximations of very large structured matrices using tensor trains (TT), and quantized tensor trains (QTT) is considered here. The focus is on matrices with a small TT/QTT rank. Some analysis of the error produced by the use of the TT/QTT representation and the underlying approximation formula used is also provided. Promising experiments on exponential, power, Mittag-Leffler and logarithm function of multilevel Toeplitz matrices, that are among those which generate a low TT/QTT rank representation, are also provided, confirming that the proposed approach is feasible. (C) 2019 Elsevier B.V. All rights reserved.
A distance based on the exponential kernel of the adjacency matrix of a graph and representing how well two vertices connect to each other in a graph is defined and studied. This communicability cosine distance (CCD) ...
详细信息
A distance based on the exponential kernel of the adjacency matrix of a graph and representing how well two vertices connect to each other in a graph is defined and studied. This communicability cosine distance (CCD) is a Euclidean spherical distance accounting for the cosine of the angles spanned by the position vectors of the graph vertices in this space. The Euclidean distance matrix (EDM) of CCD is used to quantify the similarity between vertices in graphs and networks as well as to define a local vertex invariant-a closeness centrality measure, which discriminate very well vertices in small graphs. It allows to distinguish all nonidentical vertices, also characterizing all identity (asymmetric) graphs-those having only the identity automorphism-among all connected graphs of up to 9 vertices. It also characterizes several other classes of identity graphs. We also study real-world networks in term of both the discriminating power of the new centrality on their vertices as well as in ranking their vertices. We analyze some dictionary networks as well as the network of copurshasing of political books, remarking some of the main advantages of the new approaches studied here.
We propose a communication-driven mechanism for predicting triadic closure in complex networks. It is mathematically formulated on the basis of communicability distance functions that account for the quality of commun...
详细信息
We propose a communication-driven mechanism for predicting triadic closure in complex networks. It is mathematically formulated on the basis of communicability distance functions that account for the quality of communication between nodes in the network. We study 25 real-world networks and show that the proposed method correctly predicts 20% of triadic closures in these networks, in contrast to the 7.6% predicted by a random mechanism. We also show that the communication-driven method outperforms the random mechanism in explaining the clustering coefficient, average path length, and average communicability. The new method also displays some interesting features with regards to optimizing communication in networks.
We consider the problem of attaining either the maximal increase or reduction of the robustness of a complex network by means of a bounded modification of a subset of the edge weights. We propose two novel strategies ...
详细信息
We consider the problem of attaining either the maximal increase or reduction of the robustness of a complex network by means of a bounded modification of a subset of the edge weights. We propose two novel strategies combining Krylov subspace approximations with a greedy scheme and an interior point method employing either the Hessian or its approximation computed via the limited-memory Broyden-Fletcher-Goldfarb-Shanno algorithm (L-BFGS). The paper discusses the computational and modeling aspects of our methodology and illustrates the various optimization problems on networks that can be addressed within the proposed framework. Finally, in the numerical experiments we compare the performances of our algorithms with state-of-the-art techniques on synthetic and real-world networks.
We revisit the problem of computing submatrices of the Cramer-Rao bound (CRB), which lower bounds the variance of any unbiased estimator of a vector parameter theta. We explore iterative methods that avoid direct inve...
详细信息
We revisit the problem of computing submatrices of the Cramer-Rao bound (CRB), which lower bounds the variance of any unbiased estimator of a vector parameter theta. We explore iterative methods that avoid direct inversion of the Fisher information matrix, which can be computationally expensive when the dimension of theta is large. The computation of the bound is related to the quadratic matrix program, where there are highly efficient methods for solving it. We present several methods, and show that algorithms in prior work are special instances of existing optimization algorithms. Some of these methods converge to the bound monotonically, but in particular, algorithms converging nonmonotonically are much faster. We then extend the work to encompass the computation of the CRB when the Fisher information matrix is singular and when the parameter theta is subject to constraints. As an application, we consider the design of a data streaming algorithm for network measurement.
We propose a simple technique that, if combined with algorithms for computing functions of triangular matrices, can make them more efficient. Basically, such a technique consists in a specific scaling similarity trans...
详细信息
We propose a simple technique that, if combined with algorithms for computing functions of triangular matrices, can make them more efficient. Basically, such a technique consists in a specific scaling similarity transformation that reduces the departure from normality of a triangular matrix, thus decreasing its norm and in general its function condition number. It can easily be extended to non-triangular matrices, provided that it is combined with algorithms involving a prior Schur decomposition. Situations where the technique should be used or not will be discussed in detail. Special attention is devoted to particular algorithms like the inverse scaling and squaring to the matrix logarithm or inverse cosine and the scaling and squaring to the matrix exponential. The advantages of our proposal are supported by theoretical results and illustrated with numerical experiments, involving matrices of small, medium and large size.
In this paper we are interested in the approximation of fractional powers of self-adjoint positive operators. Starting from the integral representation of the operators, we apply the trapezoidal rule combined with a d...
详细信息
In this paper we are interested in the approximation of fractional powers of self-adjoint positive operators. Starting from the integral representation of the operators, we apply the trapezoidal rule combined with a double-exponential transform of the integrand function. In this work we show how to improve the existing error estimates for the scalar case and also extend the analysis to operators. We report some numerical experiments to show the reliability of the estimates obtained.
This paper studies the robust Ho filter design problem for discrete-time linear system with polytopic uncertainty. Differing from existing ideas in the literature to reduce filtering conservatism, the present article ...
详细信息
This paper studies the robust Ho filter design problem for discrete-time linear system with polytopic uncertainty. Differing from existing ideas in the literature to reduce filtering conservatism, the present article proposes a less conservative two-stage algorithm to handle this issue by stepwisely recovering and utilizing matrix variable information of analysis condition. This algorithm is formulated by virtue of the combined matrix inequality with a prefixed binary choice signal, which is composed of direct design condition with nonlinear matrix variable inequality and its corresponding linearized one: In general;a more relaxed matrix variable structure is derived based on improved scalar parameter approach from the perspective of eigenvalue theory, which is then applied to obtain less conservative filter design condition with linearized matrix variable inequality in the first stage of provided algorithm. With some solved variables in stage one, that nonlinear matrix variable inequality based design condition is transformed into linear version, and in the second stage it will be used to further optimize the filtering performance level via picking up the lost matrix variable information of stage one. Finally, the advantages of the given algorithm are clearly illustrated by two examples. (C) 2016 The Franklin Institute. Published by Elsevier Ltd. All rights reserved.
暂无评论