1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
|
%
% teil1.tex -- Beispiel-File für das Paper
%
% (c) 2020 Prof Dr Andreas Müller, Hochschule Rapperswil
%
\section{Padé-Approximation
\label{transfer:section:teil2}}
\rhead{}
\subsection{Idee
\label{transfer:pade:idee}}
Die Taylorapproximation ist für den Gebrauch als Ersatz des Tangenshyperbolicus als Transferfunktion nicht brauchbar. Die Padé-Approximation kann die grössten Probleme aber entschärfen und dies mit sehr begrenztem zusätzlichen Rechenaufwand. Dafür wird die Taylorapproximation in einen Bruch von zwei Polynom zerlegt.
\subsection{Definition
\label{transfer:pade:definition}}
Sei
\begin{equation}
R(x)=\frac{\sum_{j=0}^{m} a_{j} x^{j}}{1+\sum_{k=1}^{n} b_{k} x^{k}}=\frac{a_{0}+a_{1} x+a_{2} x^{2}+\cdots+a_{m} x^{m}}{1+b_{1} x+b_{2} x^{2}+\cdots+b_{n} x^{n}}
\end{equation}
und gilt
\begin{gather*}
f(0) =R(0) \\
f^{\prime}(0) =R^{\prime}(0) \\
f^{\prime \prime}(0) =R^{\prime \prime}(0) \\
\vdots \\
f^{(m+n)}(0) =R^{(m+n)}(0),
\end{gather*}
so ist $R(x)$ die Padé-Approximation von $f(x)$.
\subsection{Beispiel
\label{transfer:pade:beispiel}}
Sei $f(x) = \tanh (x)$ und $T_{5} \tanh(x ; a) = x-\frac{x^{3}}{3}+\frac{2 x^{5}}{15}$, dann gilt
$$
\begin{gathered}
[3 / 2]_{f}(x) = \frac{A_{0}+A_{1} x+A_{2} x^{2}+A_{3} x^{3}}{B_{0}+B_{1} x+B_{2} x^{2}}=x-\frac{x^{3}}{3}+\frac{2 x^{5}}{15}+O\left(x^{6}\right), B_{0} = 1,\\
\Downarrow \\
[3 / 2]_{f}(x) = \frac{15x+x^3}{15+6x^2}
\end{gathered}
$$
\begin{figure}
\centering
\begin{tikzpicture}
\begin{axis}[
xmin=-3.5, xmax=3.5,
ymin=-1.5, ymax=1.5,
axis lines=center,
axis on top=true,
domain=-3.5:3.5,
ylabel=$y$,
xlabel=$x$,
]
\addplot [mark=none,draw=red,thick] {tanh(\x)};
\node [right, red] at (axis cs: 1.4,0.7) {$\tanh(x)$};
\addplot [mark=none,draw=blue,ultra thick, samples=100, smooth] expression{x*(15+x^2)/(15+6*x^2)};
\node [right, blue] at (axis cs: -1.8,0.7) {$Padé$};
%% Add the asymptotes
\draw [blue, dotted, thick] (axis cs:-2.5,-1)-- (axis cs:0,-1);
\draw [blue, dotted, thick] (axis cs:+2.5,+1)-- (axis cs:0,+1);
\end{axis}
\end{tikzpicture}
\caption{$[3 / 2]_{f}(x)$
\label{motivation:figure:Pade32}}
\end{figure}
|
