% Vlhkost dřeva
% Robert Mařík
% jaro 2017
# Úvod
Posláním tohoto dokumentu je jednak ukázka toho, jak může vypadat
projekt do Aplikované matematiky a jednak vysvětlení významu parciální
derivace. Protože se jedná i o učební materiál, je delší než
požadovaná stránka a půl. Protože budu pracovat se vzorečkem který je
jednoduchý na zápis do počítače, přidávám do dokumentu i vzorce se
složitějším formátováním, jako například definici parciální derivace.
# Teorie parciální derivace
Je-li $f(x,y)$ funkce dvou proměnných, je parciální derivace podle $x$
definována vztahem
$$\frac{\partial f}{\partial x}=\lim_{\Delta x\to 0}\frac{f(x+\Delta x,y)-f(x,y)}{\Delta x}$$
a analogicky je definována parciální derivace podle $y$. Fyzikálně se
u parciální derivace podle $x$ jedná o rychlost s jakou se mění funkce
$f$ při změnách funkce $x$. Pokud se veličina $x$ změní o $\Delta x$,
změní se funkce $f$ přibližně o
$\frac{\partial f}{\partial x} \Delta x$. Analogické tvrzení platí pro
proměnnou $y$.
# Vlhkost dřeva

Nejjednodušší a nejrychlejší metodou měření vlhkosti dřeva je použití
měřiče určeného k tomuto účelu. Tento přístroj měří elektrické
vlastnosti dřeva a na základě korelací mezi těmito veličinami a
vlhkostí dřeva zobrazuje přímo výslednou vlhkost. Tato metoda je
rychlá, ale není příliš přesná, protože korelace mezi vlhkostí a
elektrickými vlastnostmi dřeva není moc těsná. Přesnější ale nepoměrně
složitější metodou je stanovení hmotnosti dřeva v syrovém stavu,
hmotnosti dřeva po vysušení v peci a výpočet podílu vody ve
dřevě. Přesná metodika měření (při jaké teplotě se dřevo suší, jak má
vypadat měřený vzorek, jak se ním má před měřením a po měření
zacházet, jak se pozná že vzorek je dostatečně vysušený) nás teď
nezbude zajímat, podíváme se pouze na vzorec, do nějž se dosazují
naměřené hodnoty.
Budeme pracovat s definicí vlhkosti dřeva ve tvaru
$$M=\frac{m-m_0}{m_0}$$
kde $m$ je syrového hmotnost dřeva, $m_0$ hmotnost sušiny a $M$
vlhkost vyjádřená desetinným číslem (nikoliv v procentech). Pro snazší
derivování si funkci upravíme na tvar
$$M=m m_0^{-1}-1$$
a derivováním dostávme
$$\begin{aligned}\frac{\partial M}{\partial m}&=m_0^{-1}=\frac 1{m_0},\\
\frac{\partial M}{\partial m_0}&=m(-1)m_0^{-2}=-\frac m{m_0^2}.
\end{aligned}$$
Změna parametru $m$ o $\Delta m$ tedy vyvolá změnu $M$ přibližně o
$\frac 1{m_0}\Delta m$. Podobně, změna parametru $m_0$ o $\Delta m_0$
tedy vyvolá změnu $M$ přibližně o hodnotu $-\frac m{m_0^2}\Delta
m_0$. Všimněte si, že zvětšení hodnoty $m$ způsobí zvětšení hodnoty
$M$ (derivace podle $m$ je kladná), kdežto zvětšení hodnoty $m_0$
způsobí zmenšení hodnoty $M$ (derivace podle $m_0$ je
záporná). Hodnoty $\Delta m$ a $\Delta m_0$ často reprezentují
nepřesnosti vzniklé při měření a vektor
$$\left(\frac 1{m_0}\Delta
m,-\frac m{m_0^2}\Delta m_0\right)$$
charakterizuje vliv chyby ve měření vstupních dat na spolehlivost
výsledku. Zpravidla se za chybu výsledku považuje délka tohoto
vektoru, tj.
$$\sqrt{\left(\frac 1{m_0}\Delta m\right)^2+\left(-\frac m{m_0^2}\Delta m_0\right)^2}$$
# Příklad
Při měření vlhkosti dřeva byla naměřena hmotnost syrového vzorku $$ m=36\, \mathrm{g}\pm 1\,\mathrm{g} $$ a hmotnost sušiny $$ m_0=22\, \mathrm{g}\pm 1\,\mathrm{g} .$$ Určete vlhkost dřeva.
*Řešení:* Ze vzorce pro vlhkost dostáváme
$$M=\frac{36-22}{22}=0.6364.$$
Parciální derivace v bodě $(m,m_0)=(36,22)$ jsou
$$\begin{aligned}
\frac{\partial M}{\partial m}(36,22)&=\frac 1{22}=0.04545,\\
\frac{\partial M}{\partial m_0}(36,22)&=-\frac {36}{22^2}=-0.07438.
\end{aligned}
$$
Číslo $0.07438$ je zhruba o polovinu větší než $0.04545$.
Vidíme tedy, že parciální derivace podle $m_0$ je zhruba o polovinu větší
než parciální derivace podle $m$ a nepřesnost měření v $m_0$ se proto
projeví na spolehlivosti výsledku o polovinu větší měrou. Navíc, v
praxi bude nepřesnost měření $m_0$ větší než nepřesnost měření $m$,
protože kromě chyby způsobené vážením se na nepřesnosti měření
hmotnosti $m_0$ může podepsat i proces sušení. Dřevo nemuselo být
dostatečně vysušené, nebo naopak mohlo v peci začít hořet pokud
experimentátor spěchal a nedodržel technologický postup.
Znaménko derivace nerozhoduje, protože při stanovení chyby měření
nevíme, zda skutečná hodnota je větší nebo menší. Délka vektoru
charakterizujícího chybu je po zaokrouhlení na jednu platnou cifru
směrem nahoru rovna
$$\sqrt{(0.04545)^2+(-0.07438)^2}=0.09$$
a výsledná vlhkost dřeva je
$$M=0.64\pm 0.09,$$
tj. v procentech
$$M=64\%\pm 9\%.$$
# Odkazy
* PDF verze
* markdown verze
* html verze
* html slidy
* příkazy pro generování všech verzí