Skip to content
GitLab
Projects
Groups
Snippets
Help
Loading...
Help
What's new
7
Help
Support
Community forum
Keyboard shortcuts
?
Submit feedback
Contribute to GitLab
Sign in
Toggle navigation
Open sidebar
Erik Strand
nmm_2020_site
Commits
ece54d73
Commit
ece54d73
authored
Apr 30, 2020
by
Erik Strand
Browse files
Options
Browse Files
Download
Email Patches
Plain Diff
Add pade approximant exposition
parent
939034e2
Changes
1
Hide whitespace changes
Inline
Side-by-side
Showing
1 changed file
with
98 additions
and
0 deletions
+98
-0
_psets/10.md
_psets/10.md
+98
-0
No files found.
_psets/10.md
0 → 100644
View file @
ece54d73
---
title
:
Problem Set 10 (Functions)
---
## 1
{:.question}
Find the first five diagonal Pad
é
approximants [1/1], ..., [5/5] to $$e^x$$ around the
origin. Remember that the numerator and denominator can be multiplied by a constant to make the
numbers as convenient as possible. Evaluate the approximations at $$x = 1$$ and compare with the
correct value of $$e = 2.718281828459045$$. How is the error improving with the order? How does that
compare to the polynomial error?
For a fixed function $$f$$ and integers $$N
\g
eq 0$$ and $$M
\g
eq 0$$, the Pad
é
approximant
is the function
$$
[N/M]_ f(x) =
\f
rac{
\s
um_{n = 0}^N a_n x^n}{1 +
\s
um_{m = 1}^M b_m x^m}
$$
that matches as many terms as possible of the Taylor series of $$f$$. (We can define the approximant
about any point, but without loss of generality we'll only consider the origin.) There are $$N + M +
1$$ parameters in the formula above ($$a_0$$, $$
\l
dots$$, $$b_N$$, and $$b_1$$, $$
\l
dots$$,
$$b_M$$), so in general this means we can match $$f$$ up to order $$N + M + 1$$. (Though by
coincidence, or otherwise, we may end up matching some higher order terms as well.)
We can write the resulting system of equations explicitly by expanding the Taylor series of the
approximant to order $$L = N + M$$ (since this gives us $$N + M + 1$$ terms). To make the equations
simpler to write I will introduce a constant $$b_0 = 1$$.
$$
\f
rac{
\s
um_{n = 0}^N a_n x^n}{
\s
um_{m = 0}^M b_m x^m} =
\s
um_{l = 0}^L c_l x^l
$$
Then we multiply by the denominator.
$$
\b
egin{aligned}
\s
um_{n = 0}^N a_n x^n &=
\l
eft(
\s
um_{m = 0}^M b_m x^m
\r
ight)
\l
eft(
\s
um_{l = 0}^L c_l x^l
\r
ight)
\\
&=
\s
um_{m = 0}^M
\s
um_{l = 0}^L b_m c_l x^{m + l}
\e
nd{aligned}
$$
By setting the different powers of $$x$$ equal, this gives us $$L + 1 = N + M + 1$$ equations.
$$
\b
egin{cases}
a_n =
\s
um_{m = 0}^{
\m
in(n, M)} b_m c_{n - m} &
\t
ext{for } 0
\l
eq n
\l
eq N
\\
0 =
\s
um_{m = 0}^{
\m
in(n, M)} b_m c_{n - m} &
\t
ext{for } N < n
\l
eq L = N + M
\e
nd{cases}
$$
We know what $$c_0$$, $$
\l
dots$$, $$c_L$$ are, since they must match the Taylor series of $$f$$. (In
particular, $$c_l = f^{(l)} / l!$$.) So by solving this system of equations we can determine all of
the unknown parameters. In particular, the last $$M$$ equations (i.e. for $$N < n$$) can be solved
to determine $$b_1$$, $$
\l
dots$$, $$b_M$$. Then the first $$N + 1$$ equations immediately give
values for $$a_0$$, $$
\l
dots$$, $$a_N$$.
The first step can be written as a system of $$M + 1$$ equations.
$$
\b
egin{bmatrix}
1 & 0 & 0 &
\c
dots & 0 & 0
\\
c_{N + 1} & c_{N} & c_{N - 1} &
\c
dots & c_{N - M + 2} & c_{N - M + 1}
\\
c_{N + 2} & c_{N + 1} & c_{N} &
\c
dots & c_{N - M + 3} & c_{N - M + 2}
\\
\v
dots &
\v
dots &
\v
dots &
\d
dots &
\v
dots &
\v
dots
\\
c_{N + M - 1} & c_{N + M - 2} & c_{N + M - 3} &
\c
dots & c_{N} & c_{N - 1}
\\
c_{N + M} & c_{N + M - 1} & c_{N + M - 2} &
\c
dots & c_{N + 1} & c_{N}
\\
\e
nd{bmatrix}
\b
egin{bmatrix}
b_0
\\
b_1
\\
b_2
\\
\v
dots
\\
b_{M - 1}
\\
b_M
\e
nd{bmatrix}
=
\b
egin{bmatrix}
1
\\
0
\\
0
\\
\v
dots
\\
0
\\
0
\e
nd{bmatrix}
$$
Note that if $$N + 1 < M$$, then some upper diagonal chunk of the matrix will be zero, corresponding
to those entries where the subscript of $$c$$ would be negative.
Finally, there's no guarantee that this matrix will be nonsingular. (For example, consider the case
where $$N = M$$ and all derivatives of $$f$$ up to $$N + M$$ are zero. Then all rows of the matrix
except for the first will be zero.) In such a situation, one can try to throw out degenerate
equations, and pull new ones from higher order terms (i.e. consider $$L > N + M$$). However it's
possible that all additional rows will be degenerate; in this case, the approximant is legitimately
underdetermined. This means it can match the function exactly, and you can reduce $$M$$ until your
system is nonsingular.
Write
Preview
Markdown
is supported
0%
Try again
or
attach a new file
.
Attach a file
Cancel
You are about to add
0
people
to the discussion. Proceed with caution.
Finish editing this message first!
Cancel
Please
register
or
sign in
to comment