微分方程数值解作业 6
微分方程Problem 1
![](https://cdn.duanyll.com/img/20240525164212.png)
![](https://cdn.duanyll.com/img/20240525164402.png)
Question (a)
The equation to solve is
Use a rectangular grid
and the central difference scheme for first and second order derivatives
where
where
![](https://cdn.duanyll.com/img/20240525205836.png)
The boundary conditions can be written as
FDEs at the top row (where
Let
then the equation can be written as
where
where
where
and
It is easy to see that
Question (b)
With
and
To use the conjugate gradient method,
And derive a sufficient condition for
. , covered by the first condition.- The first row of
satifies when the first condition is satisfied. is irreducible, since the directed graph of is equivalent to the grid graph, which is connected.
The test suggests that
Using the same test, we can discover that
is a sufficient condition for
Question (c)
The symmetric condition requires
Question (d)
Combining the results from (b) and (c), we can find a sufficient condition for the utilization of the conjugate gradient method is
Problem 2
![](https://cdn.duanyll.com/img/20240525164256.png)
![](https://cdn.duanyll.com/img/20240525164420.png)
Question (a)
Equation (6.17)
for
The FDEs for
With the boundary condition
Let
we can express the above FDE system as
where
It is obvious that
Question (b)
Use second ordered forward difference at
and (6.17) at
That is
With
we express the above FDE system as
where
where
It is clear that
Problem 3
![](https://cdn.duanyll.com/img/20240525164443.png)
![](https://cdn.duanyll.com/img/20240525164454.png)
Question (a)
It is easy to see that the directed graph of
Question (b)
The Gershgorin-Taussky theorem states that the eigenvalues of a real symmetric matrix lies in the union of intervals
For the given matrix
Therefore
Problem 4
![](https://cdn.duanyll.com/img/20240525164524.png)
![](https://cdn.duanyll.com/img/20240525164542.png)
Question (a)
From equation (6.15) we can derive
Substracting the above equation with Equation (6.17)
and use the definition of error
The above equation has the same coefficient matrix as Equation (6.17). Same simplification can be applied to the boundary conditions, where the left side of the equaiton has the same coefficient matrix, and the constants on the right side are discarded, leaving only
Question (b)
The Schur inequality states that
That is
Therefore
Question (c)
From Section 6.2.2 we know
Therefore
That indicates that the error is bounded by the truncation error.