next up previous


Weighted residual methods for solving ODE and PDE
Below we discuss some approximation methods to solve arbitrary linear differential equations. Both, ordinal (ODE) and partial (PDE) differential equations, are considered. With limited scope we consider as example the following one-dimensional differential equation
equation648
where u(x) is the unknown function and g(x) is a known function. The capital L denotes a linear differential operator, which specifies the actual form of the differential equation (73) (for example tex2html_wrap_inline4666 or tex2html_wrap_inline4668). The boundary conditions are given by
equation652

Multiplying (73) by an arbitrary weight function v(x) and integrating over the interval [a,b] one obtains
equation654
Evidently (73) and (75) are equivalent, because v(x) is an arbitrary function. Now we seek a numerical solution to (75), (74) in the form
equation656
Here tex2html_wrap_inline4676, ..., tex2html_wrap_inline4678 are functions of x and a1, ..., an are unknown coefficients.

In vector form (76) becomes
equation658
where
equation662

In (75) we may substitute u by u* to obtain
equation669
However, substituting u(x) by its approximation u*(x) in (73), generally it appears that (73) is not satisfied exactly, e.g.
equation671
Here e(x) is a measure for the error.

It follows from (79)-(80) that
equation673
Obviously, the residual, e(x), depends on the unknown parameters given by vector tex2html_wrap_inline4698. Therefore the coefficients a1, ..., an must be determined so, that expression (81) is satisfied.

Generally
equation676
where V1, ..., Vn are known functions of x and c1, .., cn are certain parameters. In terms of vector notation (82) reads
equation678
where
equation681
Evidently
equation688
and therefore (see (81))
equation692
Relation (86) holds for arbitrary cT- matrices, i.e.
equation694
or
equation698
Now, we have n equations (88) to determine coefficients a1,...,an. Inserting (77) in (80) yields
equation700
and the condition (86) can be rewritten as
equation705
Introducing the matrix K and the vector tex2html_wrap_inline4722 as
equation711
we can write (90) in compact form
equation715
Finally, we have n linear equations (92) for determing n coefficients a1,...,an.




next up previous