Present Value in the Case of General Interest Rate and General Rate of Inflation
Početní základy finanční matematiky
sylabus 5. přednášky
(5.9)
24. 10. 2005
4 Úroková míra veličin, které nejsou v ustáleném stavu
4.1 Abstrakt
Ve finanční matematice se používají dva pojmy: jednoduché a složené úročení. Úroková míra za jednotku času je
v obou případech rovna
f(x + 1) - f(x)
f(x) ,
kde f je stavová funkce. Protože f nemusí být v bodě x + 1 vůbec definována (předpokládáme, že je definována
na nějakém pravém okolí x, x + h bodu x), počítá se v prvním případě úroková míra za jednotku času takto
 = 1 + f(x + h) - f(x)
f(x)
1
h
- 1
a ve druhém případě takto:
 = f(x + h) - f(x)
f(x) 
1
h
a říká se, že je f úročena složeně, pokud  nezávisí na h a že je úročena jednoduše, pokud  nezávisí na h. f je
úročena složeně, právě tehdy, když je v ustáleném stavu.
4.1.0.1 Pokud je f v ustáleném stavu (tj. roste exponenciálně), nezávisí  ani na x.
4.1.0.2 Pokud  nezávisí na x, je f nulová konstanta.
4.1.0.3 Proto je číslo  mnohem vhodnější míra než číslo  a v případech, které se dějí ze své a ne z naší vůle, jako je třeba
inflace používáme číslo .
4.1.0.4 Máme dvě různé míry růstu veličiny. Jejich nevýhodou je, že obě závisí nejen na bodě x, ve kterém růst měříme,
ale i na délce času h.
4.1.0.5 V makroekonomii se jako míra růstu veličiny f používá hodnota výrazu
 = f
(x)
f(x) = limh0
f(x + h) - f(x)
f(x) 
1
h
= limh0

4.1.0.6 Veličiny f, pro které je  konstantní jsou v ustáleném stavu. za povšimnutí stojí fakt, že se nepočítá limita , ale .
4.1.0.7 Ve finanční matematice také se používá pojem spojité úročení namísto složené úročení a v tom případě je úroková
míra jméno pro číslo , ne pro číslo . Tato neobratnost ve vyjařování (nejde o jiné úročení, stav kapitálu pokaždé
roste exponenciálně, nebo tvoří geometrickou posloupnopst, je-li čas diskrétní, ale o jiné číslo, které je nazváno
úrokovou mírou) je ve většině učebnic finanční matematiky znejasňována ještě tím, že se předvádějí různé početní
hříčky, kdy se z čísla  dostane . (známá je například tak, kdy se za h postupně dosazují prvky posloupnosti
(1/2n
)
n=0, číslúm (1/2)n
se říká frekvence připisování úroků)
4.1.0.8 Cílem této poznámky je:
4.1.0.9 Najít limitní tvar (pro h  0) míry , který lze použít namísto 
4.1.0.10 Najít vzorec pro budoucí hodnotu kapitálu, v němž bude mít pojem úroková míra tentýž význam, jaký má při
složeném úročení, totiž vzorec, do kterého bude možno dosadit dostatečně obecnou funkci, která bude vyjadřovat
závislost úrokové míry na čase a který přejde do známého tvaru xt = x0(1 + i)t
, pokud bude ůroková míra i
konstantní (poznamenejme, že vzorec xt = x0e it
používaný v učebnicích finanční matematiky tuto vlastnost
nemá).
4.1.0.11 Na případech, jejichž chování je autonomní (nezávislé na naší vůli), jako je inflace, ukázat, že volba tohoto
vzorce není libovolná, ale nutná: Totiž vzorec xt = x0e ln(1+) je limitním případem vzorce xt = x0(1 + i)t
tomto smyslu: pokud funkcionál, který přiřazuje úrokové míře budoucí hodnotu kapitálu má být spojitý, je budoucí
hodnota při úročení libovolnou úrokovou mírou vynucena známými budoucími hodnotami při úroočení (po částech)
konstantní úrokovou mírou, protože množina po částech konstantních funkcí je hustá v množině všech přípustných
funkcí.
41
Present Value in the Case of General Interest Rate and General Rate of Inflation
4.2 Odvození hlavní formule
1. Introduction:
4.2.0.1 Nech jsou aregované ceny nějaké funkce závislá na čase
F = CPI:R - R (1)
Míra inflace za dobu t1, t2 je relativní hodnota jejich přírůstku:
(t1, t2) = ( t1, t2 ) = F(t2) - F(t1)
F(t1) . (2)
Takto definováno je  zobrazení
: R2 - R (3)
4.2.0.2 1.1. Definition: (see. [2]) Steady state is a situation in which the various quantities splňují následující rovnici:
D(x)(t) = k  x(t) (4)
každá pro nějakou konstantu k. Tedy, každá funkce v ustáleném stavu závisí na dvou parametrech A a B a je tvaru
x  A  eBx
. (5)
4.2.0.3 Jediná quantity, kterou se nyní zabýváme jsou agregované ceny.
4.2.0.4 1.2. Lemma: Nech jsou v čase T1, T2 ceny v ustáleném stavu. Pak
: t1:t2: t1, t1 +   T1, T2  t2, t2 +   T1, T2 = ( t1, t1 +  ) = ( t2, t2 +  ) (6)
i. e. ( t, t + ) depends only on the length  of the interval ( t, t + but not on its beginning t. In a steady state
we can define
^: 0,T2 - T1 - R
x - ( T1, T1 + x ) (7)
Při tom platí:
^() = (1 + (t))

t - 1 (8)
takže známe-li hodnotu ^ v nějakém t  R můžeme rozšířit definiční obor ^ na celé R:
^: R - R
 - ^() = (1 + (t))

t - 1 (9)
Takto definované je ^ míra inflace agregovaných cen, které jsou v ustáleném stavu na R a které vyhovující rovnici
F() = F(0)  (1 +^(t))
t (10)
pro každé t  T0, T0 + t Tedy zejména, pokud zvolíme hodnotu t za jednotku času, máme:
^() = (1 +^(1))
- 1 = F(x + ) - F(x)
F(x) (11)
a naopak
^(1) = (1 +^())1
 - 1 = F(x + ) - F(x)
F(x)
1

. (12)
Takto můžeme definovat míru inflace za jednotku času pokud jsou agregované ceny v ustáleném stavu na nějakém,
libovolně malém, pravém okolí bodu x  R.
4.2.0.5 V obecném případě, i pokud není stav ustálený, můžeme normovat míru inflace naměřenou na intervalu t, t + 
s použitím vztahu (12). Definujeme:
42
Present Value in the Case of General Interest Rate and General Rate of Inflation
4.2.0.6 1.3. Definition: Míra inflace agregovaných cen F na intervalu t, t +  za jednotku času je:
~(t,) = F(t + ) - F(t)
F(t)
1/
. (13)
4.2.0.7 1.4. Note: ~ závisí jen na hodnotách F v krajních bodech intervalu t, t +  a je stejná, jako by byla hodnota
inflace v ustáleném stavu, ve kterém by agregované ceny v čase t byly F(t) a v čase t +  by byly F(t + ) tj.
podle (5) by byly splněny rovnice:
F(t) = AeBt
F(t + ) = AeB(t+)
(14)
které jsou ekvivaklentní s rovnicemi
B = ln F(t + )
F(t) -1
A = F(t + )
F(t)
- t

F(t).
(15)
4.2.0.8 Nás ovšem zajímá míra inflace  = lim0+ ~(-,). v jednom časovém okamžiku t. Pokud v tomto okamžiku
mají agregované ceny nespojitost prvního druhu je míra inflace  nenulová a je rovna
(t) = limt+ F() - limt- F()
limt- F() (16)
a míra inflace  za jednotku času bude rovna .
4.2.0.9 Pokud jsou v tomto okamžiku agregované ceny spojitě závislé na čase je míra inflace v tomto okamžiku 0. Pokud
jsou navíc v tomto okamžiku agregované ceny diferencovatelné podle času, existuje limita (15)   0 a platí:
lim0
A = e- D(F )(t)t
F (t) F(t)
lim0
B = D(F)(t)
F(t)
(17)
a míra inflace za jednotku času je
(t) = ^(1)AeB(t+1)
- AeBt
AeBt
= eB
- 1 = e
D(F )(t)
F (t) - 1 = eD(lnF)(t) (18)
a máme
ln((t) + 1) = D(lnF)(t) (19)
This is the main result of this part. Poslední rovnici můžeme psát také ve tvaru:
D(F)(t) = ln((t) + 1)  F(t) (20)
4.2.0.10 One can compare it with a equation
D(F)(t) = (t)  F(t) (21)
která se obvykle používá k definici míry růstu of different quantities v makroekonomii a jejíž speciálním případem
je (5). Podle toho, co jsme odvodili, by měla být makroekonomická míra růstu  rovna ln(1 + ).
4.2.0.11 Podotkněme ještě, že taylorův polynom stupmně 1v bodě 0funkce x  ln(1+x)je stejný jako taylorův polynom
identity, že tedy obě funkce mají dotyk řádu 1 (first order contact) a proto pro dostatečně malé hodnoty míry růstu
( nebo ) vychází obě míry ( a ln(1 + )) přibližně stejně.
4.2.0.12 2. Inverzní problém:
4.2.0.13 Nyní naopak předpokládejme, že známe míru inflace za jednotku času v každém bodě t  R.
4.2.0.14 (21) is differential equation for the state function f. The solution is
f(x) = e
x
0
ln((s)+1)ds
 C (22)
43
Present Value in the Case of General Interest Rate and General Rate of Inflation
We determine the value of C substituting constant function for : t  . It have to be fulfilled that:
f(x) = ( + 1)x
 f(0) (23)
but
f(x) = e
x
0
ln( +1)ds
 C = C  ( + 1)x
(24)
then:
C = f(0) (25)
And the state function has the form:
f:x  e
x
0
ln((s)+1)ds
f(0) (26)
where an (t) is interest rate per unit of time in time t.
4.2.0.15 2.5. Example
4.2.0.16 Let us suppose, that we can check the rate of inflation per unit of time in every moment (it can be the result of
measuring and some general theory). We know, that in time 0 the rate of unit time inflation was 0.1 and in time 1
the rate of unit time inflation was 0.2.
4.2.0.17 Předpokládejme, že uvnitř intervalu (0,1) se hodnoty inflace za jednotku času měnily podle jedné z následujících
možností (skokem v jednom y obou krajních bodů, lineárně, kvadraticky): in this four following ways:
1(u) := 0.1, if u < 1
0.2, if u  1 , 2(u) := u2
10 + 0.1, 3(u) := u
10 + 0.1, 4(u) := 0.1, if u  0
0.2, if u > 0 (27)
Míra inflace za dobu (0,1) byla v každém z těcto čtyř případů jiná. In general, we have
( 0,1 ) = e
1
0
(ln(1+(u)du
- 1 (28)
So, in our four cases we obtain:
 = 1 = 0.1:X(1) = e
1
0
ln(1.1) du
- 1 = 0.1 (29)
 = 2 = u 
1
10 u2 + 0.1:X(1) = e
1
0
ln(1.1+u2
/10) du
- 1 = 0.132945354. . . (30)
 = 3 = u 
1
10 u + 0.1:X(1) = e
1
0
ln(1.1+u/10) du
- 1 = 0.149637533. . . (31)
 = 4 = 0.2:X(1) = e
1
0
ln(1.2) du
- 1 = 0.2 (32)
4.2.0.18 3. Interest rate: Předpokládejme, že v každém okamžiku času t známe nominální hodnotu kapitálu f(t), jehož
množství roste nebo klesá pouze úročením. Potom v souladu s obvyklou terminologií nazýváme relative increment
of the capital during the time the interest rate xi:
( t1, t2 ) = x(t2) - x(t1)
x(t1) (33)
a tento vztah je stejný jako vztah (2), protože míra inflace je úroková míra, kterou se ůročí ceny.
4.2.0.19 Je celkem obvyklé, že banka vyhlašuje úrokovou sazbu p. a., i když ji mění častěji, než jednou za rok. Je
obvyklé, že střadatelé uloží peníze a chtějí je vybrat podle svých potřeb jindy, než po nějaké předem stanovené
periodě času a banky jim to umožní. Ve všech těchto případech se stanovuje úrok jako odměna za dočasné
přenechání práva disponovat s penězy zpravidla takto: Jsou-li Ii úrokové sazby p. a. konstantní na intervalech
Ii = (ti, ti+1); i = 0. . . n.
n-1
i=0
(1 + Ii)(ti+1-ti) (34)
(a úrok s touto mírou je připsán v okamžiku tn). Poznamenejme ještě, že jednotkou měření času je časový interval
na který je normovaná úroková míra (v našem případě je p. a., tedy jednotkou času je rok) a že banky obvykle
zaokrouhlují časové údaje na celé dny, což může být jev pouze dočasný.
44
Present Value in the Case of General Interest Rate and General Rate of Inflation
4.2.0.20 V obchodních vztazích je obvyklé, že se úroková míra stanovuje jako po částech konstantní funkce a její hodnota
se normuje, ale v obecných makroekonomických a ekonometrických úvahách, pokud například půjde o agregátní
úrokovou sazbu, to již po částech konstantní funkce být nemusí.
4.2.0.21 V učebnicích finanční matematiky bývá definováno tzv. spojité úročení, kde je hodnota stavové funkce x v závislosti
na úrokové míře  dána vztahem:
x(t) = x(0)e
t
0
(s) ds
(35)
namísto vztahu
x(t) = x(0)e
t
0
ln(1+(s))ds
(36)
jako v (22). Funkce 35 je řešením rovnice
D(x)(t) = (t)  x(t) (37)
namísto rovnice
D(x)(t) = ln(1 + (t))  x(t)
jako v (21). Funkce (35) ale není zobecněním složeného úročení (33) protože v ustáleném stavu nedává výsledek
x(t) = x(0)  (1 + )t
, (38)
proto je například nelze použít pro kvantitativní výpočty na datech v případě inflace a v příkladě (29) by nám vyšlo
e
1
0
0.1du - 1 = 0.105170918,
v příkladě (30) by nám vyšlo
e
1
0
1/10u2
+0.1du
- 1 = 0.142630812,
v příkladě (31) by nám vyšlo
e
1
0
1/10u+0.1du
- 1 = 0.161834243
a v příkladě (32) by nám vyšlo
e
1
0
0.2du
- 1 = 0.221402758
a první a poslední výsledek jsou očividně nesprávné. Lze tedy říci, že inflace se nechová jako tzv spojité úročení.
4.2.0.22 3.6. Theorem: The formula (38) je speciálním případem formule (22) pro konstantní úrokovou míru a formule
(34) je speciálním případem formule (22) pro po částech konstantní úrokovou míru.
4.2.0.23 Proof: Let us suppose, that (t) = I is constant:
e
t
0
ln(1+I) du
= e(tln(1+I)) = e(ln(1+I)t
) = (1 + I)t
(39)
q. e. d.
4.2.0.24 Now let us suppose, that  is piecewise constant and that it has value Ii in every point of the interval Ii =
(ti, ti+1); i = 0. . . n. Let A be the characteristic function of the set A, then (t) = (ti,ti+1)  Ii
= e
t
0
ln(1+(ti,ti+1)Ii du
= e
n-1
i=0
(ti+1 ln(1+Ii)-ti ln(1+Ii))
=
=
n-1
i=0
e(ti+1-ti)ln(1+Ii) =
n-1
i=0
eln(1+Ii)(ti+1-ti)
=
n-1
i=0
(1 + Ii)(ti+1-ti) (40)
q. e. d.
4.2.0.25 Funkcionál
  e
t
0
ln(1+(s)ds
je spojitý v topologii stejnoměrné konvergence. Důležitos věty 38 plyne z následujícího lemmatu:
4.2.0.26 3.7. Lemma: pro každou spojitou funkci  definovanou na uzavřeném intervalu existuje posloupnost piecewise
konstant funkcí i takových, že i  
45
Present Value in the Case of General Interest Rate and General Rate of Inflation
4.2.0.27 Proof: Předpokládejme, že f je spojitá. Podle předpokladu je Dom(f) kompaktní. Zvolíme nějaké  kladné.
Pro každé x  Dom(f) najdeme okolí O(x) takové, aby F(O(x))  O/2(f(x)). O(x) tvoří pokrytí Dom(f)
vybereme konečné podpokrytí . Definujeme  = minU (Diam(U)), kde (Diam(U)) je průměr množiny
U. Rozdělíme Dom(f) na n disjunktních podintervalů (J
i )i = 1n
délky . Pro každý interval J
i vybereme
nějaký bod xi, který leží uvnitř něj a označíme yi = f(xi). Definujeme: x  J
i = (x) = yi. Potom
x  Dom(f):|f(x) - (x)|  . A pro n =  1
2n
jsou splněny předpoklady lemma.
4.2.0.28 3.8. Corollary: (38) říká, že vzorec (34) je speciálním případem vorce (22). Pro funkci, která má jen konečně
mnho bodů nespojitosti můžeme provést aproximaci na každém intervalu, na kterém je spojitá jako v lemmatu
(40). A tedy vzorec (22) je limitním případem vzorce (34).
4.2.0.29 Následuje několik aplikací vzorce (22):
4.2.0.30 4. Boundary interest rate: If the interest rate is constant and positive, the state function is increasing and constant.
If it is positive, but very quickly decreasing, the state function will be concave. What interest rate makes the state
function afine (polynomial of degree 1)?
4.2.0.31 State function is
t  x(0)e
t
0
ln(1+(s))ds
(41)
its derivative is
t  x(0) ln(1 + (t))e
t
0
ln(1+(s))ds
(42)
and the second derivative is
t 
x(0)e
t
0
ln(1+(s))ds d
dt (t) + (ln(1 + (t)))2 + (ln(1 + (t)))2 (t)
1 + (t)
(43)
We are looking for the interest rate, which makes the second derivative equal to zero. If x(0) = 0 and  > 0,
second derivative is equal to zero for such a , which are the solutions of differential equation
d
dt
(t) + (ln(1 + (t)))
2 + (ln(1 + (t)))
2 (t) = 0 (44)
i. e. we solve the differential equation
d
dt
(t) = -((ln(1 + (t)))2 (1 + (t)). (45)
Its solution fulfilled the algebraical equation
- (ln(1 + (t)))
-1
+ t = C (46)
hence it is any of the functions
t  (t) = e( 1
t-C ) - 1 (47)
for all constant C. We express C using the initial condition
C = -
1
ln(1 + (0))
and we have
(t) = e
1
t+ 1
ln(1+(0))
- 1 = (1 + (0))
1
t ln(1+(0))+1
- 1 (48)
If we put this interest rate into the rule for interesting () we obtain:
x(t) = x(0) e
t
0
ln (1+(0))
1
s ln(1+(0))+1 ds
= x(0) (tln(1 + (0)) + 1) (49)
and this really is an affine function. We can conclude: If the interest rate has in the time 0 value (0) and if on the
dependence on the time grows more quickly (decrease slowly) than function
t  (1 + (0))
1
t  1
ln(1+(0))+1
- 1 (50)
46
Present Value in the Case of General Interest Rate and General Rate of Inflation
the state function is convex. If it decreases more quickly, the state function is concave.
4.2.0.32 5. Example with the real data of index of prices of non-regulated prices in Czech republic:.
4.2.0.33 We choose the year as a unit of time. We measured the index of non-regulated prices in every moment
P = 1993+ i
12 , where i is a natural number less or equal to 120.The values are relative. As a result of the mea-
surement, CPIis a functioncutting the points: (1993+ 1
12; 91.76),(1993+ 1
6; 92.97),(1993+ 1
4; 93.54), (1993+ 1
3;
93.94), (1993+ 5
12; 94.32), (1993+ 1
2; 94.62), (1993+ 7
12; 95.42), (1993+ 2
3; 96.10), (1993+ 3
4; 97.53), (1993+ 5
6;
98.62), (1993+11
12; 99.17), (1994; 100.0), (1994+ 1
12; 100.76), (1994+1
6; 101.12), (1994+1
4; 101.40), (1994+1
3;
101.88), (1994 + 5
12; 102.27), (1994 + 1
2; 103.45), (1994 + 7
12; 103.81), (1994 + 2
3; 104.62), (1994 + 3
4;
106.13), (1994+5
6; 107.45), (1994+11
12; 108.46), (1995; 109.23), (1995+ 1
12; 110.57), (1995+1
6; 111.62), (1995+1
4;
111.97), (1995 + 1
3; 112.72), (1995 + 5
12; 113.26), (1995 + 1
2; 114.13), (1995 + 7
12; 113.47), (1995 + 2
3;
113.41), (1995+3
4; 114.41), (1995+5
6; 115.20), (1995+11
12; 116.09), (1996; 116.82), (1996+ 1
12; 118.97), (1996+1
6;
119.62), (1996 + 1
4; 120.43), (1996 + 1
3; 121.15), (1996 + 5
12; 121.96), (1996 + 1
2; 123.06), (1996 + 7
12;
123.18), (1996+2
3; 122.67), (1996+3
4; 123.04), (1996+5
6; 123.76), (1996+11
12; 124.44), (1997; 125.23), (1997+ 1
12;
126.28), (1997 + 1
6; 126.71), (1997 + 1
4; 126.85), (1997 + 1
3; 127.45), (1997 + 5
12; 127.60), (1997 + 1
2;
129.42), (1997 + 7
12; 129.67), (1997 + 2
3; 130.77), (1997 + 3
4; 131.58), (1997 + 5
6; 132.33), (1997 + 11
12;
133.01), (1998; 133.75),(1998+ 1
12; 135.76), (1998+1
6; 136.71), (1998+1
4; 136.85), (1998+1
3; 137.12), (1998+ 5
12;
137.26), (1998 + 1
2; 137.81), (1998 + 7
12; 137.53), (1998 + 2
3; 137.12), (1998 + 3
4; 137.25), (1998 + 5
6;
136.84), (1998+11
12; 136.43), (1999; 136.02), (1999+ 1
12; 136.70), (1999+1
6; 136.57), (1999+1
4; 136.29), (1999+1
3;
136.84), (1999 + 5
12; 136.70), (1999 + 1
2; 136.98), (1999 + 7
12; 136.98), (1999 + 2
3; 137.11), (1999 + 3
4;
136.98), (1999+5
6; 136.98), (1999+11
12; 137.39), (2000; 138.21), (2000+ 1
12; 139.04), (2000+1
6; 139.32), (2000+1
4;
139.32), (2000 + 1
3; 139.32), (2000 + 5
12; 139.74), (2000 + 1
2; 140.71), (2000 + 7
12; 141.42), (2000 + 2
3;
141.70), (2000+3
4; 141.56), (2000+5
6; 141.98), (2000+11
12; 142.13), (2001; 142.41),(2001+ 1
12; 143.26), (2001+1
6;
143.26), (2001 + 1
4; 143.26), (2001 + 1
3; 143.84), (2001 + 5
12; 144.99), (2001 + 1
2; 146.87), (2001 + 7
12;
147.90), (2001+2
3; 147.46), (2001+3
4; 145.98), (2001+5
6; 145.84), (2001+11
12; 145.69), (2002; 145.98), (2002+ 1
12;
147.30), (2002+1
6; 147.30), (2002+1
4; 147.0), (2002+1
3; 147.30), (2002+ 5
12; 147.15), (2002+ 1
2; 146.71), (2002+
7
12; 147.15), (2002 + 2
3; 146.85), (2002 + 3
4; 145.82), (2002 + 5
6; 145.82), (2002 + 11
12; 145.68), (2003; 147.97).
Now we are going to approximate the values of CPI with the function. There are a lot of ways of how to do it.
We shall measure the accuracy of approximation as the sum of squares of differences between values of function
and measured values in points where the measurement is made. Trivial approximation is to put together points
by abscissas. We obtain piecewise affine function: and the accuracy is equal to 0. In this case the rate of the year
inflation  is the function:
(t) =



e14.52058082(14.52058082t-28847.75831)-1
- 1 t  1993.00000,1993.08333
e6.840273611(6.840273611t-13540.26521)-1
- 1 t  1993.08333,1993.16666
e4.799616031(4.799616031t-9472.894168)-1
- 1 t  1993.16666,1993.25000
e4.560182407(4.560182407t-8995.643826)-1
- 1 t  1993.25000,1993.33333
e3.600144006(3.600144006t-7081.966879)-1
- 1 t  1993.33333,1993.41666
e9.599232061(9.599232061t-19040.64915)-1
- 1 t  1993.41666,1993.50000
e8.160326413(8.160326413t-16172.19129)-1
- 1 t  1993.50000,1993.58333
e17.16068643(17.16068643t-34115.15901)-1
- 1 t  1993.58333,1993.66666
e13.07895368(13.07895368t-25977.54380)-1
- 1 t  1993.66666,1993.75000
e6.600264011(6.600264011t-13060.65643)-1
- 1 t  1993.75000,1993.83333
e9.960398416(9.960398416t-19760.20401)-1
- 1 t  1993.83333,1993.91666
e9.120364815(9.120364815t-18085.24781)-1
- 1 t  1993.91666,1993.99999
e4.320172807(4.320172807t-8513.664947)-1
- 1 t  1993.99999,1994.08332
e3.359731222(3.359731222t-6598.464123)-1
- 1 t  1994.08332,1994.16666
e5.760230409(5.760230409t-11385.45902)-1
- 1 t  1994.16666,1994.24999
e4.680187207(4.680187207t-9231.582863)-1
- 1 t  1994.24999,1994.33332
...
and the rate of inflation per the time interval 1993.5,1994 equals
CPI(1994)
CPI(1993.5) - 1 = e
1994
1993.5
ln(1+(z))dz
= 0.05596 (51)
47
Present Value in the Case of General Interest Rate and General Rate of Inflation
4.2.0.34 We can use more sophisticated approximation, forinstance: let us approximate a trend byfunction151.558469453+
53.8746490595(5/12)(x-1992)-108.487078769(3/4)(x-1992) which is the best approximationof the measured
values by means of pair of functions out of the set (i/12)(x-1992); (i/12)(x-1993); (x - 1992)(i/12); (x -
1993)(i/12); ln(x - 1992); ln(x - 1991)}24
i=1 and constant (accuracy (the sum of squares of distances between
values of the functions and the measured values) equal 237.324366060)Then let us approximate the rest by linear
combination of functions
{sin(x/i); cos(x/i);x sin(x/i);x cos(x/i)}5
i=1 (52)
Let us omit 15 of the functions from this set of functions, the absence of which brings the least loss of accuracy.
We obtain the following approximation of CPI: CPI(t) = 151.558469453+ 53.8746490595 (5/12)(x-1992) -
108.487078769(3/4)(x-1992)-0.08688205539+438.678532386sin(1/5x)-0.219746118617xsin(1/5x
) - 0.392448110715x cos(1/2 x ) +0.000156347532197 xsin(x ) +782.647022069 cos(1/2 x )
(accuracy (the sum of squares of distances between values of the functions and the measured values) equals
84.8709833772)Then the rate of inflation per year has in the time x the value
e53.8746490595 A
B
A = 5
12
x-1992
ln( 5
12) - 108.487078769 (3/4)x-1992
ln(3/4)+ 87.73570648 cos(1/5x)
-0.219746118617 sin(1/5x) - 0.04394922372x cos(1/5x)
-0.392448110715 cos(1/2x) + 0.1962240554x sin(1/2x)
+0.000156347532197 sin(x) + 0.000156347532197x cos(x)
-391.3235111 sin(1/2x)
B = 151.4715874+53.8746490595 5
12
x-1992
-108.487078769 (3/4)x-1992
+438.678532386 sin(1/5x)
-0.219746118617 x sin(1/5x) - 0.392448110715x cos(1/2x)
+0.000156347532197x sin(x) + 782.647022069 cos(1/2x)-1
and
CPI(1994)
CPI(1993.5) - 1 = e
1994
1993.5
ln(1+(z))dz
= 0.04737 (53)
4.2.0.35 But we can approximate the rate of inflation directly. If (tni )iJ , i < j  ti < tj are the moments, where CPI is
measured, then we have to fulfil the equation:
i  J:
ti+1
ti
ln(1 + (t)) dt = ln CPI(tn+1) - ln CPI(tn) (54)
and the accuracy of approximation should be measured as the measure of unreliableness of the equation.
4.2.0.36 More precisely: Let us choose the function:
: x - 0.04944210+ 0.009094682 cos(x) + 0.01215932 sin(2x) +
+0.01823419 sin(4x) + 0.05263253 sin(1/3x) - (55)
-0.009667420 cos(1/3x) - 0.07212055 cos(1/4x) + 0.02834145 sin(1/2x) .
The question is, how accurately this function corresponds to the measured values. The interesting thing is, that the
function does not approximate the measured values, but some other values, which depend on the measured ones.
And the attitude towards this two values is given by (54). So if we can compare the exactness of approximation
with the same measure as we did in the previous cases, we have to compute the number
120
i=1
e
ti+1
ti
ln(1+(t))dt
-
CPI(tn+1)
CPI(tn)
2
(56)
and this is the only relevant measure. In our case we obtain: 4.9914 We used the simplest function for the
approximation of rate of inflation we considered yet, and it gives us far the best approximation we have obtained!
48
4.2.0.37 References:
[1] Aczél, J.: Lectures on functional Equations and Their Applications, New York, Academic Press 1966
[2] Barro, Robert J., Sala-i-martin, Xavier: Economic Growth, MIT Press Cambridgr, Massachusetts, London,
England 1999, ISBN 0-262-02459-4
[3] Dieudonné, J.: Treatise on analysis, Vol. III, Pure and Applied Mathematics, Vol. 10-III, Academic Press. New
York ­ London, 1972
[4] Dupačová, J., Hurt, J., Šťepán, J.: Stochastic Modeling in economics and Finance, Kluwert Academic Publishers,
2002, ISBN 1-4020-0840-6.
[5] Studený V.: Functional Equation of the Rate of Inflation, e-print archive of Coronell University, 2003,
www.arxiv.org
[6] 6. Conclusion:
[7] Our sketch brings an open problem: what does the space of solution of (54) look like and what is the space of
functions, the rate of inflation lies in. How can it help us to understand what kind of measure is this one which we
called rate (of interest or inflation. . .). It shows us the possibility to approximate the measure of inflation except
of the index of prices and it suggests, that the solution should be better.
[8] 7. Apendix: Maple programs for computation. Let us suppose that in iota we have the function for approxi-
mation, in our case:
9] > iota:=.4944210e-1 + .9094682e-2 * cos(x*Pi) + .1215932e-1 * sin(2*x*Pi) + .1823419e-
1 * sin(4*x*Pi)+ .5263253e-1 * sin(1/3*x*Pi) - .9667420e-2 * cos(1/3*x*Pi) - .7212055e-
1 * cos(1/4*x*Pi)+ .2834145e-1 * sin(1/2*x*Pi);
[9] Let us suppose, that the values of CPI are in the file Values and the corresponding moments of time, in which
the values are measured, are in file Points. Then we can compute the preciseness of approximation using the
following program:
10] > Mistake := 0;
11] > for i to NumberOfPoints do
12] > A := evalf(int(ln(1+iota),x = Points 1] .. Points i+1]));
13] > A := Values 1]*exp(A);
14] > Mistake := Mistake+(A-Values i+1])^2
15] > end do:
16] > print(Sum('Delta'^2,i)=Mistake);
[1] References:
[1] Aczél, J.: Lectures on functional Equations and Their Applications, New York, Academic Press 1966
[2] Dieudonné, J.: Treatise on analysis, Vol. III, Pure and Applied Mathematics, Vol. 10-III, Academic Press. New
York ­ London, 1972
[3] Dupačová, J., Hurt, J., Šťepán, J.: Stochastic Modeling in economics and Finance, Kluwert Academic Publishers,
2002, ISBN 1-4020-0840-6.
[4] Studený V.: Functional Equation of the Rate of Inflation, e-print archive of Coronell University, 2003,
www.arxiv.org
49