Linear Operators Research Articles

Some norm estimates for the triangular projection on the space of bounded linear operators between two symmetric sequence spaces

Some norm estimates for the triangular projection on the space of bounded linear operators between two symmetric sequence spaces

An exponential spectral deferred correction method for multidimensional parabolic problems

We present some efficient algorithms based on an exponential time differencing spectral deferred correction (ETDSDC) method for multidimensional second and fourth-order parabolic problems with non-periodic boundary conditions including Dirichlet, Neumann, Robin boundary conditions. Similar to the Fourier spectral method for periodic problems, the key to the efficiency of our algorithms is to construct diagonal discrete linear operators via Legendre–Galerkin methods with Fourier-like basis functions. In combination with the ETDSDC scheme, the proposed methods are spectrally accurate in space and up to 10th-order accurate in time (as shown in this work). We demonstrate the high-order of convergence and efficiency of our algorithms in solving parabolic equations through a series of two-dimensional and three-dimensional examples including Ginzburg–Landau and Allen–Cahn equations.

Convolution equation and operators on the euclidean motion group

Let $G = \mathbb{R}^2\rtimes SO(2)$ be the Euclidean motion group, let g be the Lie algebra of G and let U(g) be the universal enveloping algebra of g. Then U(g) is an infinite dimensional, linear associative and non-commutative algebra consisting of invariant differential operators on G. The Dirac measure on G is represented by $\delta_G$, while the convolution product of functions or measures on G is represented by $\ast$. Among other notable results, it is demonstrated that for each u in U(g), there is a distribution E on G such that the convolution equation $u\ast E = \delta_G$ is solved by method of convolution. Further more, it is established that the (convolution) operator $A^\prime : C^\infty_c(G)\rightarrow C^\infty(G),$, which is defined as $A^\prime f = f\ast T^n\delta(t)$ extends to a bounded linear operator on $L^2(G)$, for $f\in C^\infty_c(G)$, the space of infinitely differentiable functions on G with compact support. Furthermore, we demonstrate that the left convolution operator LT denoted as $L_Tf = T\ast f$ commutes with left translation, for $T\in D^\prime(G)$.

A modified inertial proximal alternating direction method of multipliers with dual-relaxed term for structured nonconvex and nonsmooth problem

In this research, we introduce a novel optimization algorithm termed the dual-relaxed inertial alternating direction method of multipliers (DR-IADM), tailored for handling nonconvex and nonsmooth problems. These problems are characterized by an objective function that is a composite of three elements: a smooth composite function combined with a linear operator, a nonsmooth function, and a mixed function of two variables. To facilitate the iterative process, we adopt a straightforward parameter selection approach, integrate inertial components within each subproblem, and introduce two relaxed terms to refine the dual variable update step. Within a set of reasonable assumptions, we establish the boundedness of the sequence generated by our DR-IADM algorithm. Furthermore, leveraging the Kurdyka–Łojasiewicz (KŁ) property, we demonstrate the global convergence of the proposed method. To validate the practicality and efficacy of our algorithm, we present numerical experiments that corroborate its performance. In summary, our contribution lies in proposing DR-IADM for a specific class of optimization problems, proving its convergence properties, and supporting the theoretical claims with numerical evidence.

Bouligand–Newton type methods for non-smooth ill-posed problems

Abstract We consider Newton-type methods for solving nonlinear ill-posed inverse problems in Hilbert spaces where the forward operators are not necessarily G\^{a}teaux differentiable. Modifications are proposed with the non-existing Fr'{e}chet derivatives replaced by a family of bounded linear operators satisfying suitable properties. These bounded linear operators can be constructed by the Bouligand subderivatives which are defined as limits of Fr'{e}chet derivatives of the forward operator in differentiable points. The Bouligand subderivative mapping in general is not continuous unless the forward operator is G\^{a}teaux differentiable which
introduces challenges for convergence analysis of the corresponding Bouligand-Newton type methods. In this paper we will show that, under the discrepancy principle, these Bouligand-Newton type methods are iterative regularization methods of optimal order. Numerical results for an inverse problem arising from a non-smooth semi-linear elliptic equation are presented to test the performance of the methods.

Data-Driven Fault Detection and Isolation for Multirotor System Using Koopman Operator

This paper presents a data-driven fault detection and isolation (FDI) for a multirotor system using Koopman operator and Luenberger observer. Koopman operator is an infinite-dimensional linear operator that can transform nonlinear dynamical systems into linear ones. Using this transformation, our aim is to apply the linear fault detection method to the nonlinear system. Initially, a Koopman operator-based linear model is derived to represent the multirotor system, considering factors like non-diagonal inertial tensor, center of gravity variations, aerodynamic effects, and actuator dynamics. Various candidate lifting functions are evaluated for prediction performance and compared using the root mean square error to identify the most suitable one. Subsequently, a Koopman operator-based Luenberger observer is proposed using the lifted linear model to generate residuals for identifying faulty actuators. Simulation and experimental results demonstrate the effectiveness of the proposed observer in detecting actuator faults such as bias and loss of effectiveness, without the need for an explicitly defined fault dataset.

Sharp maximal function estimates for linear and multilinear pseudo-differential operators

In this paper, we study pointwise estimates for linear and multilinear pseudo-differential operators with exotic symbols in terms of the Fefferman-Stein sharp maximal function and Hardy-Littlewood type maximal function. Especially in the multilinear case, we use a multi-sublinear variant of the classical Hardy-Littlewood maximal function introduced by Lerner, Ombrosi, Pérez, Torres, and Trujillo-González [16], which provides more elaborate and natural weighted estimates in the multilinear setting.

Piecewise modified Prandtl-Ishlinskii model for precise nano-positioning

Piezoelectric ceramic, an indispensable micro-nano actuator, has advantages of small size, high sensitivity and strong controllability, which is widely applied in various fields like micro-nano processing, medicine, chemistry and materials. However, its nonlinear and asymmetric hysteresis features directly affect its positioning accuracy. The classical and generalized Prandtl-Ishlinskii models with linear symmetry operators are commonly used to deal with these issues, but they still face difficulties in asymmetric systems. Therefore, a piecewise modified Prandtl-Ishlinskii (PMPI) model is proposed to compensate nonlinear and asymmetric hysteresis in this paper. The new nonlinear operator and piecewise nonlinear envelope function are introduced into the PMPI model. Independent right and left parts are adopted to describe hysteresis increasing and decreasing stages, while the threshold and weight vectors correspond to them. Weight vectors are extended and a new diagonal matrix is constructed for avoiding excessive parameters and computation. Simulations and experiments are shown that the PMPI model has good compensation performance and effectiveness, which has better fitting, higher accuracy, fewer parameters and computation, and can be well applied to the precise nano-positioning systems with less nonlinear and asymmetric errors.

Parallel isogeometric boundary element analysis with T-splines on CUDA

We present a framework for parallel isogeometric boundary element analysis (BEA) of elastic solids on CUDA. To deal with traction discontinuities, we propose a BEA model that supports multiple nodes and semi-discontinuous elements. The multiplicity of a node is defined by the number of regions containing any element influenced by the node. A region is a group of connected elements delimited by a closed crease curve. The default shape function of an element is determined by a linear operator applied to a set of basis functions. A BEA model is supposed to be generated from a watertight boundary representation of a solid. In this paper, we employ bicubic analysis-suitable T-splines. In this case, the shape of an element is defined by its Bézier extraction operator applied to the tensor product of Bernstein polynomials of degree 3. We describe the data structures of the BEA model and the main algorithms of the analysis pipeline on CUDA. In particular, we describe two strategies for parallel assembling of the linear system of equations. We also introduce a novel approach for inside integration based on the subdivision of the singular region in triangles with constant aspect ratio. In the T-splines context, we extend the Bézier extraction to handle unstructured T-meshes with crease edges. Moreover, we propose a scheme for embedding the influence of linked tangency handles on the shape of an element directly into the Bézier extraction operator. Such an embedding enables the removal of the corresponding nodes from the BEA model and the application of an alternative collocation method we discuss in the paper. We present several experiments to evaluate the accuracy and efficiency of the proposed framework. The results demonstrate that the GPU can be advantageously employed for parallelizing T-spline-based isogeometric analysis using boundary elements, achieving a speedup of up to 29x compared to the sequential code on a current laptop. We make the BEA code available in a prototype in MATLAB with a graphical interface that allows users to apply boundary conditions and visualize analysis results on the boundary.

The Invariant Subspace Problem for Separable Hilbert Spaces

In this paper, we prove that every bounded linear operator on a separable Hilbert space has a non-trivial invariant subspace. This answers the well-known invariant subspace problem.

Globalization of distributed parameter self‐optimizing control

AbstractNumerous nonlinear distributed parameter systems (DPSs) operate within an extensive range due to process uncertainties. Their spatial distribution characteristic, combined with nonlinearity and uncertainty, poses challenges in optimal operation under two‐step real‐time optimization (RTO) and economic model predictive control (EMPC). Both methods necessitate substantial computational power for prompt online reoptimization. Recent local distributed parameter self‐optimizing control (DPSOC) achieves optimality without repetitive optimization. However, its effectiveness is confined to a narrow range around a nominal operation. Here, globalized DPSOC is introduced to surmount the limitation of the local DPSOC. A global loss functional concerning controlled variables (CVs) is formulated using linear operators and Fubini's theorem. Minimizing the loss with a numerical optimization procedure yields CVs exhibiting global optimality. Maintaining these CVs at constants ensures such a process has a minimal average loss in a large operating space. The effectiveness of the proposed method is substantiated through a transport reaction simulation.

Koopman analysis of the singularly perturbed van der Pol oscillator.

The Koopman operator framework holds promise for spectral analysis of nonlinear dynamical systems based on linear operators. Eigenvalues and eigenfunctions of the Koopman operator, the so-called Koopman eigenvalues and Koopman eigenfunctions, respectively, mirror global properties of the system's flow. In this paper, we perform the Koopman analysis of the singularly perturbed van der Pol system. First, we show the spectral signature depending on singular perturbation: how two Koopman principal eigenvalues are ordered and what distinct shapes emerge in their associated Koopman eigenfunctions. Second, we discuss the singular limit of the Koopman operator, which is derived through the concatenation of Koopman operators for the fast and slow subsystems. From the spectral properties of the Koopman operator for the singularly perturbed system and the singular limit, we suggest that the Koopman eigenfunctions inherit geometric properties of the singularly perturbed system. These results are applicable to general planar singularly perturbed systems with stable limit cycles.

Some Results about Acts over Monoid and Bounded Linear Operators

الأثر V بالنسبة إلى sinshT و خواصه قد تم دراسته في هذا البحث حيث تم دراسة علاقة الأثر المخلص والاثر المنتهى التولد والاثر المنفصل وربطها بالمؤثرات المتباينة حيث تم بهنة العلاقات التالية ان الاثر اذا وفقط اذا مقاس في حالة كون المؤثر هو عديم القوة وكذلك في حالة كون المؤثر شامل فان الاثر هو منتهي التولد اي ان الغضاء هو منتهي التولد وايضا تم برهن ان الاثر مخلص لكل مؤثر مقيد وك\لك قد تم التحقق من انه لاي مؤثر مقيد فان الاثر منفصل وفي حالة كون الاثر منفصل فان المؤثر سوف يكون متباين وايضا تم برهنة انه في حالة كون الفضاء يمتلك قاعدة فان الاثر سوف يكون كل عنصر فيه يمكن كتابته كتركيبة خطية ومن العلاقات المهمة التي قد تم برهانها انه اذا كان المؤثر مشابه لمؤثر اخر فان الاثر سوف يكون نوثيرين بوجد شرط على الفضاء .

General Upper Bounds for the Numerical Radii of Hilbert Space Operators

We present a collection upper bounds for the numerical radii of a certain 2 × 2 operator matrices. We use these bounds to improve on some known numerical radius inequalities for powers of Hilbert space operators. In particular, we show that if 𝐴 is a bounded linear operator on a complex Hilbert space, then 𝑤 2𝑟 (𝐴) ≤ 1+𝛼 8 ‖|𝐴| 2𝑟 +|𝐴 ∗ | 2𝑟‖+ 1+𝛼 4 𝑤(|𝐴| 𝑟 |𝐴 ∗ | 𝑟 )+ 1−𝛼 2 𝑤 𝑟 (𝐴 2 ) for every r ≥ 1 and α ∈ [0,1]. This substantially improves on the existing inequality 𝑤 2𝑟 (𝐴) ≤ 1 2 ‖|𝐴| 2𝑟 + |𝐴 ∗ | 2𝑟‖. Here 𝑤(. ) and ||. || denote the numerical radius and the usual operator norm, respectively.

Linear Operating Range and Extension of Three-Dimensional Active Zero-State PWM for Open-End Winding Drive

Linear Operating Range and Extension of Three-Dimensional Active Zero-State PWM for Open-End Winding Drive

The reconstruction problem for a multivalued linear operator's properties

The reconstruction problem for a multivalued linear operator's properties

Deformation of superintegrability in the Miwa-deformed Gaussian matrix model

We consider an arbitrary deformation of the Gaussian matrix model parametrized by Miwa variables za. One can look at it as a mixture of the Gaussian and logarithmic (Selberg) potentials, which are both superintegrable. The mixture is not, still one can find an explicit expression for an arbitrary Schur average as a linear transform of a polynomial made from the values of skew Schur functions at the Gaussian locus pk=δk,2. This linear operation includes multiplication with an exponential eza2/2 and a kind of Borel transform of the resulting product, which we call multiple and enhanced. The existence of such remarkable formulas appears intimately related to the theory of auxiliary K-polynomials, which appeared in superintegrable correlators at the Gaussian point (strict superintegrability). We also consider in great detail the generating function of correlators ⟨(TrX)k⟩ in this model, and discuss its integrable determinant representation. At last, we describe deformation of all results to the Gaussian β-ensemble. Published by the American Physical Society 2024

Contractivity of Möbius functions of operators

Let T be an injective bounded linear operator on a complex Hilbert space. We characterize the complex numbers λ,μ for which (I+λT)(I+μT)−1 is a contraction, the characterization being expressed in terms of the numerical range of the possibly unbounded operator T−1.When T=V, the Volterra operator on L2[0,1], this leads to a result of Khadkhuu, Zemánek and the second author, characterizing those λ,μ for which (I+λV)(I+μV)−1 is a contraction. Taking T=Vn, we further deduce that (I+λVn)(I+μVn)−1 is never a contraction if n≥2 and λ≠μ.

A fully implicit, asymptotic-preserving, semi-Lagrangian algorithm for the time dependent anisotropic heat transport equation

In this paper, we extend the operator-split asymptotic-preserving, semi-Lagrangian algorithm for time dependent anisotropic heat transport equation proposed in Chacón et al. (2014) [18] to use a fully implicit time integration with backward differentiation formulas. The proposed implicit method can deal with arbitrary heat-transport anisotropy ratios χ∥/χ⊥⋙1 (with χ∥, χ⊥ the parallel and perpendicular heat diffusivities, respectively) in complicated magnetic field topologies in an accurate and efficient manner. The implicit algorithm is second-order accurate temporally and demonstrates an accurate treatment at boundary layers (e.g., island separatrices), which was not ensured by the operator-split implementation. The condition number of the resulting algebraic system is independent of the anisotropy ratio, and is inverted with preconditioned GMRES. We propose a simple preconditioner that renders the finite-dimensional linear operator compact, resulting in mesh-independent convergence rates for topologically simple magnetic fields, and convergence rates scaling as ∼(NΔt)1/4 (with N the total mesh size and Δt the timestep) in topologically complex magnetic-field configurations. We demonstrate the accuracy and performance of the approach with test problems of varying complexity, including an analytically tractable boundary-layer problem in a straight magnetic field, and a topologically complex magnetic field featuring magnetic islands with extreme anisotropy ratios (χ∥/χ⊥=1010).

Linear and nonlinear pseudo-differential operators on p-adic fields

Linear and nonlinear pseudo-differential operators on p-adic fields