Fel

För andra betydelser, se Fel (olika betydelser).

Fel, inom numerisk analys, anger differensen mellan ett värde och dess närmevärde. Fel spelar en stor roll inom kaos (icke-linjära ekvationer) där små variationer kan ge totalt olika resultat. Men det är även en central fråga när beräkningar ska utföras på datorer. Fel kan bero på flera olika saker.

Några inledande definitioner i ämnet:

• exakt värde: $a = 1,\!25 \pm 0,\!01$
• närmevärde: $\bar a = 1,\!25$
• absolut fel: $\Delta a = \bar a - a \qquad |\Delta a| \le 0,\!01$
• relativt fel: ${\Delta a \over a} \approx {\Delta a \over \bar a}, a \neq 0 \qquad {|\Delta a| \over |a|} \lesssim {0,\!01 \over 1,\!25} = 8*10^{-3}$

Anledningen att vi tar beloppet av vissa värden här ovanför är för att vi oftast bara är intresserade, eller bara vet, just beloppet av felet. Om absolutbeloppet av det absoluta felet $|\Delta a| \le 0,5 * 10^{-t}$ säger man att närmevärdet $\bar a$ har t st korrekta decimaler.

I ett närmevärde räknas alla siffror (förutom inledande nollor) på högre position än $10^{-t}$ som signifikanta siffror.

Felkällor

Fel i ett resultat kan komma från flera olika källor på grund av olika sorters approximationer. Den första approximationen som görs är oftast i den matematiska modellen, det vill säga att ett problem förenklas och idealiseras. Till exempel att man bortser från energiförluster, luftmotstånd och andra påverkningar.

Uppmätta värden är oftast inte exakta på grund av att mätinstrumenten har en ändlig precision. Det kan också vara så att felet kommer från tidigare beräkningar, från trunkering/diskretisering, det vill säga att man hugger av eller delar upp en oändlig serie, eller från avrundningar.

Bilden illustrerar hur feltyperna kommer in i beräkning av en integral:

  $R_T \$ pga metoden

  $R_{XX} \$ från fel i gränsen

  $R_{XF} \$ från fel i funktionsvärdena

$R_B \$ från avrundning av resultatet.

När man sedan ska göra beräkningar brukar man prata om tre typer av fel:

1. $R_X$ - Fel i indata är ofta oundvikliga och kan delas in i två undertyper:
   1. $R_{XF}$ - Fel som uppstår på grund av fel i använda funktionsvärden
   2. $R_{XX}$ - Fel som uppstår på grund av andra fel i indata
2. $R_T$ - Trunkeringsfel uppstår på grund av en använd approximativ metod, till exempel att en oändliga summa huggs av och ersätts med en ändlig
3. $R_B$ - Avrundningsfel uppstår då man räknar med ett ändligt antal siffror

Beloppet av det totala felet i ett värde blir $|R_{TOT}| = |R_{XX}| + |R_{XF}| + |R_T| + |R_B| \$

Avrundningsfel

Om $n=2,1212 \pm 0,2$ så är det onödigt att närmevärdet har fyra decimaler eftersom felet kan vara så stort som 0,2. Då kan vi välja att avrunda närmevärdet till 1 decimal. Vårt nya n blir således $n_1 = 2,1$ och vi får ett avrundningsfel $|R_B| = |n_1 - \bar n| = 0,0212 < 0,03$. För att skriva ut vårt nya n måste vi lägga till avrundningsfelet till den ursprungliga felgränsen så att $n = 2,1 \pm 0,23$.

Felet vid avrundning är alltid maximalt $0,5*10^{-t}$, där t är antalet decimaler efter avrundningen.

Felfortplantning av fel i indata

Om man har tal behäftade med fel, dvs närmevärden, och utför aritmetiska beräkningar på dem blir naturligtvis även svaret behäftat med fel. Den största möjliga storleken på felet kan beräknas med några olika metoder.

• För addition och subtraktion: det absoluta felet är summan av de absoluta felen i indata.
• För multiplikation och division: det relativa felet är summan av de relativa felen i indata.

Addition:

$y_1 = x_1 + x_2 \$ där $\Delta x_1 = \bar x_1 - x_1 och \Delta x_2 = \bar x_2 - x_2$ ger

$\Delta y_1 = \bar y_1 - y_1 = \bar x_1 + \bar x_2 - (x_1 + x_2) = \Delta x_1 + \Delta x_2$

$\Delta y_1 = \Delta x_1 + \Delta x_2 \Rightarrow |\Delta y_1| \le |\Delta x_1| + |\Delta x_2|$

Subtraktion:

$y_2 = x_1 - x_2 \$

$\Delta y_2 = \Delta x_1 - \Delta x_2 \Rightarrow |\Delta y_2| \le |\Delta x_1| + |\Delta x_2|$

Särskilt för subtraktion innebär detta att det relativa felet kan bli obegränsat stort om man subtraherar två tal som är ungefär lika stora (till exempel 100±2 - 100±2 = 0±4 vilket skulle ge det relativa felet 0/4). Detta fenomen kallas kancellation. Kancellation kan ibland minskas genom omskrivningar.

Funktion

En mer generellt tillämpbar metod är differentialkalkylens medelvärdessats. Som går att applicera på en deriverbar funktion f, till exempel f(x,y) = y * cos(x).

Här finns ett tal $\xi$ mellan $\bar x$ och x så att: $\Delta f = f(\bar x) - f(x) = f'(\xi)(\bar x - x)$ detsamma gäller för variabeln y.

När man praktiskt ska beräkna detta brukar man använda att $f'(\xi) \approx f'(\bar x)$ och sedan lägga till en säkerhetsmarginal, det vill säga avrunda felet uppåt.

Detta ger felfortplantningsformeln:

$\Delta f \approx {\delta f(\bar x, \bar y) \over \delta x} \Delta x + {\delta f(\bar x, \bar y) \over \delta y} \Delta y$

Och eftersom vi är intresserade av beloppet på felet kan vi skriva upp en generell Maximalfelsuppskattning:

$|\Delta f| \lesssim \sum_{k=1}^n |{\delta f \over \delta x_k} \Delta x_k|$

Lagring av tal i datorer

De flesta datorer arbetar med en fix informationslängd på tal, ett så kallat ord. Heltal kan lagras exakt i en dator, förutsatt att de inte är för stora, men reella tal kan få fel då de konverteras till basen 2, eller vid avkortningsfel på grund av att datorn jobbar med en fix ordlängd.

Ett reellt tal x lagras i en dator som ett flyttal i normaliserad flytande representation:

$x = m * \beta^{exp} \$ där $1 \le |m| < \beta,$ $\beta = 2$ (för dator) $\beta = 10$ (för räknedosa)

Totalt fel

Antag att funktionen f(x) beskriver ett fenomen exakt. Om $\hat{f}$ är närmevärdet till $f$ (approximationen $\hat{f}$ på det verkliga $f$) och $\hat{x}$ är närmevärdet till $x$ (mätvärde $\hat{x}$ till det verkliga värdet $x$) kommer totala felet att uppgå till

$\mbox{totalt fel} = \hat{f}(\hat{x}) - f(x).$

Ovanstånde kan även uttryckas med

$\hat{f}(\hat{x}) - f(x) = [\hat{f}(\hat{x}) - f(\hat{x})] + [f(\hat{x}) - f(x)],$

där

$e_b = \hat{f}(\hat{x}) - f(\hat{x})$

respektive

$e_d = f(\hat{x}) - f(x)$

är beräkningsfel och datafel.

(en. computational error respektive propagated data error.)

Beräkningsfelet är det fel som uppstår då datorer eller miniräknare evaluerar värdet. Datafel är det fel som uppstår som en följd av instrumentets fel och det exakta värdet.

Felmått

Fel kan mätas med två metoder: absolut fel och relativt fel.

Om $y$ är ett reellt tal och $\hat{y}$ är dess närmevärde, definieras:

• absolutfelet i $\hat{y}$ som $\Delta{y} = \hat{y} - y$ samt
• relativa felet i $\hat{y}$ som $\delta{y} = \frac{\Delta{y}}{y} = \frac{\hat{y} - y}{y}$ om $y \ne 0$

Relationen mellan absolutfel och relativt fel kan uttryckas med

$\hat{y} = y\cdot(1 + \delta{y})$

(men oftast finns inte det exakta värdet y att utgå från.)

För det mesta används ett mått på det absoluta eller relativa felet i värsta fall. Om ett värde anges till $\hat{x}$ och det absoluta felet till $\Delta{x}$ så ligger alltså det verkliga värdet $x$ någonstans mellan $\hat{x}-\Delta{x}$ och $\hat{x}+\Delta{x}$.

Källor

• Lars Eldén och Linde Wittmeyer-Koch (2009). Numeriska beräkningar -- analys och illustrationer med MATLAB (utgåva 4:7)