Nonhomogeneous linear system's particular solution and general solution structure

Introduction

A nonhomogeneous linear system is a system of linear differential equations where the right-hand side is not identically zero. Such systems are distinguished from homogeneous systems by the presence of external inputs that vary with time or space. The solution to a nonhomogeneous system consists of two parts: the general solution to the corresponding homogeneous system and a particular solution to the nonhomogeneous equation. This article explores the structure of particular solutions and their integration with homogeneous solutions to form the general solution of nonhomogeneous linear systems.

Definition and Properties of Nonhomogeneous Linear Systems

A nonhomogeneous linear system can be represented as:
$$ \mathbf{y}'(t) = A(t)\mathbf{y}(t) + \mathbf{f}(t), $$
where $ A(t) $ is a time-varying matrix, $ \mathbf{y}(t) $ is the state vector, and $ \mathbf{f}(t) $ is a nonhomogeneous term. The homogeneous counterpart of this system is:
$$ \mathbf{y}'(t) = A(t)\mathbf{y}(t). $$
The general solution to the nonhomogeneous system is the sum of the homogeneous solution and a particular solution to the nonhomogeneous equation. The homogeneous solution, $ \mathbf{y}_h(t) $, satisfies $ \mathbf{y}_h'(t) = A(t)\mathbf{y}_h(t) $, while the particular solution, $ \mathbf{y}_p(t) $, satisfies $ \mathbf{y}_p'(t) = A(t)\mathbf{y}_p(t) + \mathbf{f}(t) $.

The key property of nonhomogeneous systems is that the particular solution is unique up to the homogeneous solution. This uniqueness is guaranteed by the linearity of the system and the properties of linear operators. The structure of the general solution is thus:
$$ \mathbf{y}(t) = \mathbf{y}_h(t) + \mathbf{y}_p(t). $$
This decomposition ensures that the general solution captures both the transient behavior (due to the particular solution) and the steady-state behavior (due to the homogeneous solution).

Particular Solution: Methods and Techniques

Finding a particular solution to a nonhomogeneous system is often challenging due to the complexity of the nonhomogeneous term $ \mathbf{f}(t) $. Common methods include the method of undetermined coefficients, variation of parameters, and the method of annihilators.

Method of Undetermined Coefficients

This method is applicable when $ \mathbf{f}(t) $ is a polynomial, exponential, sine, or cosine function. For example, if $ \mathbf{f}(t) = e^{at} \cos(bt) $, the particular solution is assumed to be of the form:
$$ \mathbf{y}_p(t) = e^{at} \left( C \cos(bt) + D \sin(bt) \right). $$
The coefficients $ C $ and $ D $ are determined by substituting $ \mathbf{y}_p(t) $ into the nonhomogeneous equation and equating like terms. This method relies on the assumption that the particular solution does not overlap with the homogeneous solution, which is verified by ensuring that the assumed form of $ \mathbf{y}_p(t) $ does not coincide with the homogeneous solution.

Variation of Parameters

When the nonhomogeneous term is more complex, such as a product of polynomials and exponentials, the method of variation of parameters is employed. This method involves assuming a particular solution of the form:
$$ \mathbf{y}_p(t) = \mathbf{y}_1(t) \mathbf{y}_2(t) \int \frac{\mathbf{w}_2(t)}{\mathbf{w}(t)} dt + \mathbf{y}_2(t) \mathbf{y}_1(t) \int \frac{\mathbf{w}_1(t)}{\mathbf{w}(t)} dt, $$
where $ \mathbf{w}(t) = \det \begin{bmatrix} \mathbf{y}_1(t) & \mathbf{y}_2(t) \ \mathbf{y}_1'(t) & \mathbf{y}_2'(t) \end{bmatrix} $ is the Wronskian of the homogeneous solution vectors $ \mathbf{y}_1(t) $ and $ \mathbf{y}_2(t) $. This method requires solving a system of differential equations for the particular solution, ensuring that the particular solution is expressed in terms of the homogeneous solutions.

Method of Annihilators

For nonhomogeneous terms that are not expressible in terms of polynomials, exponentials, or trigonometric functions, the method of annihilators is used. This method involves finding a differential operator that "annihilates" the nonhomogeneous term, transforming the nonhomogeneous equation into a homogeneous equation. For example, if $ \mathbf{f}(t) = e^{at} $, the operator $ (D - a) $ annihilates $ e^{at} $, where $ D $ is the differential operator $ \frac{d}{dt} $. This technique is particularly useful for solving systems with nonhomogeneous terms that are not easily handled by the other methods.

Role of Homogeneous Solutions

The homogeneous solution plays a critical role in the structure of the general solution. It represents the system's behavior in the absence of the nonhomogeneous term $ \mathbf{f}(t) $, capturing the long-term trends or steady-state behavior. The homogeneous solution is determined by solving the homogeneous system:
$$ \mathbf{y}'(t) = A(t)\mathbf{y}(t). $$
The solution to this system depends on the eigenvalues and eigenvectors of the matrix $ A(t) $. If the system is stable, the homogeneous solution decays over time, while if it is unbounded, the solution may grow indefinitely.

The homogeneous solution is essential for ensuring that the general solution is complete. Without the homogeneous solution, the system would not account for the transient