Source author record

Paul Van Dooren

Paul Van Dooren appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

10works
8topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

10 published item(s)

preprint2022arXiv

Root-max Problems, Hybrid Expansion-Contraction, and Quadratically Convergent Optimization of Passive Systems

We present quadratically convergent algorithms to compute the extremal value of a real parameter for which a given rational transfer function of a linear time-invariant system is passive. This problem is formulated for both continuous-time and discrete-time systems and is linked to the problem of finding a realization of a rational transfer function such that its passivity radius is maximized. Our new methods make use of the Hybrid Expansion-Contraction algorithm, which we extend and generalize to the setting of what we call root-max problems.

preprint2021arXiv

Linear system matrices of rational transfer functions

In this paper we derive new sufficient conditions for a linear system matrix $$S(λ):=\left[\begin{array}{ccc} T(λ) & -U(λ) \\ V(λ) & W(λ) \end{array}\right],$$ where $T(λ)$ is assumed regular, to be strongly irreducible. In particular, we introduce the notion of strong minimality, and the corresponding conditions are shown to be sufficient for a polynomial system matrix to be strongly minimal. A strongly irreducible or minimal system matrix has the same structural elements as the rational matrix $R(λ)= W(λ) + V(λ)T(λ)^{-1}U(λ)$, which is also known as the transfer function connected to the system matrix $S(λ)$. The pole structure, zero structure and null space structure of $R(λ)$ can be then computed with the staircase algorithm and the $QZ$ algorithm applied to pencils derived from $S(λ)$. We also show how to derive a strongly minimal system matrix from an arbitrary linear system matrix by applying to it a reduction procedure, that only uses unitary equivalence transformations. This implies that numerical errors performed during the reduction procedure remain bounded. Finally, we show how to perform diagonal scalings to an arbitrary pencil such that its row and column norms are all of the order of 1. Combined with the fact that we use unitary transformation in both the reduction procedure and the computation of the eigenstructure, this guarantees that we computed the exact eigenstructure of a perturbed linear system matrix, but where the perturbation is of the order of the machine precision.

preprint2020arXiv

Structured backward errors for eigenvalues of linear port-Hamiltonian descriptor systems

When computing the eigenstructure of matrix pencils associated with the passivity analysis of perturbed port-Hamiltonian descriptor system using a structured generalized eigenvalue method, one should make sure that the computed spectrum satisfies the symmetries that corresponds to this structure and the underlying physical system. We perform a backward error analysis and show that for matrix pencils associated with port-Hamiltonian descriptor systems and a given computed eigenstructure with the correct symmetry structure there always exists a nearby port-Hamiltonian descriptor system with exactly that eigenstructure. We also derive bounds for how near this system is and show that the stability radius of the system plays a role in that bound.

preprint2016arXiv

A framework for structured linearizations of matrix polynomials in various bases

We present a framework for the construction of linearizations for scalar and matrix polynomials based on dual bases which, in the case of orthogonal polynomials, can be described by the associated recurrence relations. The framework provides an extension of the classical linearization theory for polynomials expressed in non-monomial bases and allows to represent polynomials expressed in product families, that is as a linear combination of elements of the form $ϕ_i(λ) ψ_j(λ)$, where $\{ ϕ_i(λ) \}$ and $\{ ψ_j(λ) \}$ can either be polynomial bases or polynomial families which satisfy some mild assumptions. We show that this general construction can be used for many different purposes. Among them, we show how to linearize sums of polynomials and rational functions expressed in different bases. As an example, this allows to look for intersections of functions interpolated on different nodes without converting them to the same basis. We then provide some constructions for structured linearizations for $\star$-even and $\star$-palindromic matrix polynomials. The extensions of these constructions to $\star$-odd and $\star$-antipalindromic of odd degree is discussed and follows immediately from the previous results.

preprint2016arXiv

Robustness and Perturbations of Minimal Bases

Polynomial minimal bases of rational vector subspaces are a classical concept that plays an important role in control theory, linear systems theory, and coding theory. It is a common practice to arrange the vectors of any minimal basis as the rows of a polynomial matrix and to call such matrix simply a minimal basis. Very recently, minimal bases, as well as the closely related pairs of dual minimal bases, have been applied to a number of problems that include the solution of general inverse eigenstructure problems for polynomial matrices, the development of new classes of linearizations and $\ell$-ifications of polynomial matrices, and backward error analyses of complete polynomial eigenstructure problems solved via a wide class of strong linearizations. These new applications have revealed that although the algebraic properties of minimal bases are rather well understood, their robustness and the behavior of the corresponding dual minimal bases under perturbations have not yet been explored in the literature, as far as we know. Therefore, the main purpose of this paper is to study in detail when a minimal basis $M(λ)$ is robust under perturbations, i.e., when all the polynomial matrices in a neighborhood of $M(λ)$ are minimal bases, and, in this case, how perturbations of $M(λ)$ change its dual minimal bases. In order to study such problems, a new characterization of whether or not a polynomial matrix is a minimal basis in terms of a finite number of rank conditions is introduced and, based on it, we prove that polynomial matrices are generically minimal bases with very specific properties. In addition, some applications of the results of this paper are discussed.

preprint2016arXiv

Structured backward error analysis of linearized structured polynomial eigenvalue problems

We introduce a new class of structured matrix polynomials, namely, the class of M_A-structured matrix polynomials, to provide a common framework for many classes of structured matrix polynomials that are important in applications: the classes of (skew-)symmetric, (anti-)palindromic, and alternating matrix polynomials. Then, we introduce the families of M_A-structured strong block minimal bases pencils and of M_A-structured block Kronecker pencils,, and show that any M_A-structured odd-degree matrix polynomial can be strongly linearized via an M_A-structured block Kronecker pencil. Finally, for the classes of (skew-)symmetric, (anti-)palindromic, and alternating odd-degree matrix polynomials, the M_A-structured framework allows us to perform a global and structured backward stability analysis of complete structured polynomial eigenproblems, regular or singular, solved by applying to an M_A-structured block Kronecker pencil a structurally backward stable algorithm that computes its complete eigenstructure, like the palindromic-QR algorithm or the structured versions of the staircase algorithm. This analysis allows us to identify those M_A-structured block Kronecker pencils that yield a computed complete eigenstructure which is the exact one of a slightly perturbed structured matrix polynomial.These pencils include (modulo permutations) the well-known block-tridiagonal and block-antitridiagonal structure-preserving linearizations. Our analysis incorporates structure to the recent (unstructured) backward error analysis performed for block Kronecker linearizations by Dopico, Lawrence, Pérez and Van Dooren, and share with it its key features, namely, it is a rigorous analysis valid for finite perturbations, i.e., it is not a first order analysis, it provides precise bounds, and it is valid simultaneously for a large class of structure-preserving strong linearizations.

preprint2014arXiv

Low-rank Similarity Measure for Role Model Extraction

Computing meaningful clusters of nodes is crucial to analyze large networks. In this paper, we present a pairwise node similarity measure that allows to extract roles, i.e. group of nodes sharing similar flow patterns within a network. We propose a low rank iterative scheme to approximate the similarity measure for very large networks. Finally, we show that our low rank similarity score successfully extracts the different roles in random graphs and that its performances are similar to the pairwise similarity measure.

preprint2013arXiv

Fast community detection using local neighbourhood search

Communities play a crucial role to describe and analyse modern networks. However, the size of those networks has grown tremendously with the increase of computational power and data storage. While various methods have been developed to extract community structures, their computational cost or the difficulty to parallelize existing algorithms make partitioning real networks into communities a challenging problem. In this paper, we propose to alter an efficient algorithm, the Louvain method, such that communities are defined as the connected components of a tree-like assignment graph. Within this framework, we precisely describe the different steps of our algorithm and demonstrate its highly parallelizable nature. We then show that despite its simplicity, our algorithm has a partitioning quality similar to the original method on benchmark graphs and even outperforms other algorithms. We also show that, even on a single processor, our method is much faster and allows the analysis of very large networks.

preprint2012arXiv

Exploring the Mobility of Mobile Phone Users

Mobile phone datasets allow for the analysis of human behavior on an unprecedented scale. The social network, temporal dynamics and mobile behavior of mobile phone users have often been analyzed independently from each other using mobile phone datasets. In this article, we explore the connections between various features of human behavior extracted from a large mobile phone dataset. Our observations are based on the analysis of communication data of 100000 anonymized and randomly chosen individuals in a dataset of communications in Portugal. We show that clustering and principal component analysis allow for a significant dimension reduction with limited loss of information. The most important features are related to geographical location. In particular, we observe that most people spend most of their time at only a few locations. With the help of clustering methods, we then robustly identify home and office locations and compare the results with official census data. Finally, we analyze the geographic spread of users' frequent locations and show that commuting distances can be reasonably well explained by a gravity model.

preprint2011arXiv

Extracting spatial information from networks with low-order eigenvectors

We consider the problem of inferring meaningful spatial information in networks from incomplete information on the connection intensity between the nodes of the network. We consider two spatially distributed networks: a population migration flow network within the US, and a network of mobile phone calls between cities in Belgium. For both networks we use the eigenvectors of the Laplacian matrix constructed from the link intensities to obtain informative visualizations and capture natural geographical subdivisions. We observe that some low order eigenvectors localize very well and seem to reveal small geographically cohesive regions that match remarkably well with political and administrative boundaries. We discuss possible explanations for this observation by describing diffusion maps and localized eigenfunctions. In addition, we discuss a possible connection with the weighted graph cut problem, and provide numerical evidence supporting the idea that lower order eigenvectors point out local cuts in the network. However, we do not provide a formal and rigorous justification for our observations.