Example Worksheet - Maple Help

All Products Maple MapleSim

Home : Support : Online Help : Science and Engineering : Scientific Error Analysis : Example Worksheet

Applications of the ScientificErrorAnalysis Package

The ScientificErrorAnalysis package provides representation and construction of numerical quantities that have a value and associated error, which is a measure of the degree of precision to which the quantity's value is known. Calculations of error analysis can be performed with these quantities.

Initialization

$restart$

>	$with (ScientificConstants) &colon;$ $with (ScientificErrorAnalysis) &colon;$

Specification of Quantities with Error

The same quantity-with-error can be entered three different ways.

>	$Quantity (20.0, 0.1)$

$Quantity (20.0, 0.1)$

(2.1)

>	$Quantity (20.0, 0.005,' relative')$

$Quantity (20.0, 0.1000)$

(2.2)

>	$Quantity (20.0, 1,' uld')$

$Quantity (20.0, 0.1)$

(2.3)

Other examples:

>	$Quantity (3.4&ExponentialE;15, 0.2&ExponentialE;15)$

$Quantity (3.4 10^{15}, 2. 10^{14})$

(2.4)

>	$Quantity (5.1&ExponentialE;−20, 0.1,' relative')$

$Quantity (5.1 10^{-20}, 5.1 10^{-21})$

(2.5)

>	$Quantity (4.2012883, 12,' uld')$

$Quantity (4.2012883, 0.0000012)$

(2.6)

Important: These quantities do not represent intervals in which possible values must be contained - they represent unknown values with central tendency. See any text on error analysis for the physical sciences or engineering.

A Note about the Display of Maple Floating-Point Numbers

If a floating-point number is not displayed in scientific notation, it may not be apparent whether trailing zeros are significant figures.

>	$Quantity (4.0&ExponentialE;2, 1,' uld')$

$Quantity (400., 10.)$

(2.1.1)

>	$SFloatMantissa (GetValue ()), SFloatMantissa (GetError ())$

$40, 1$

(2.1.2)

Simple Error Analysis - Propagation of Errors

Errors can be propagated through calculations by using the combine/errors command.

>	$q1 ≔ Quantity (10.0, 1.0)$

$q1 := ScientificErrorAnalysis:-Quantity (10.0, 1.0)$

(3.1)

>	$q2 ≔ Quantity (20.0, 3.0)$

$q2 := ScientificErrorAnalysis:-Quantity (20.0, 3.0)$

(3.2)

Error in the sum of the above two quantities:

>	$combine (q1 + q2,' errors')$

$ScientificErrorAnalysis:-Quantity (30.0, 3.162277660)$

(3.3)

Relative error in their sum:

>	$\frac{GetError ()}{GetValue ()}$

$0.1054092553$

(3.4)

Error in their difference:

>	$combine (q2 - q1,' errors')$

$ScientificErrorAnalysis:-Quantity (10.0, 3.162277660)$

(3.5)

Relative error in their difference:

>	$\frac{GetError ()}{GetValue ()}$

$0.3162277660$

(3.6)

As expected, the relative error in the difference is greater than that in the sum.

Error in the product of the above two quantities:

>	$combine (q1 q2,' errors')$

$ScientificErrorAnalysis:-Quantity (200.00, 36.05551275)$

(3.7)

Relative error in their product:

>	$\frac{GetError ()}{GetValue ()}$

$0.1802775638$

(3.8)

Error in their quotient:

>	$combine (\frac{q2}{q1},' errors')$

$ScientificErrorAnalysis:-Quantity (2.000000000, 0.3605551275)$

(3.9)

Relative error in their quotient:

>	$\frac{GetError ()}{GetValue ()}$

$0.1802775638$

(3.10)

As expected, the relative errors in the product and quotient are the same.

Error in a power:

>	$combine ({q1}^{2},' errors')$

$ScientificErrorAnalysis:-Quantity (100.00, 20.00000000)$

(3.11)

Error in differentiable functions:

>	$combine (\sin (q1),' errors')$

$ScientificErrorAnalysis:-Quantity (- 0.5440211109, 0.8390715291)$

(3.12)

>	$combine ({&ExponentialE;}^{\frac{q1}{q2}},' errors')$

$ScientificErrorAnalysis:-Quantity (1.648721271, 0.1486137271)$

(3.13)

Important: The above calculations are not interval arithmetic, and the quantities involved do not represent intervals in which values must be contained.

The Formulae of Error Propagation

The usual expressions for error propagation through simple functions are easily obtained.

>	$e1 ≔ Quantity (x, delx)$

$e1 := ScientificErrorAnalysis:-Quantity (x, delx)$

(4.1)

>	$e2 ≔ Quantity (y, dely)$

$e2 := ScientificErrorAnalysis:-Quantity (y, dely)$

(4.2)

Error in a sum:

>	$combine (e1 + e2,' errors')$

$ScientificErrorAnalysis:-Quantity (x + y, \sqrt{{dely}^{2} + {delx}^{2}})$

(4.3)

Error in a product:

>	$combine (e1 e2,' errors')$

$ScientificErrorAnalysis:-Quantity (x y, \sqrt{x^{2} {dely}^{2} + y^{2} {delx}^{2}})$

(4.4)

Relative error $^{2}$ in the product:

>	$expand (\frac{GetError {()}^{2}}{GetValue {()}^{2}})$

$\frac{{dely}^{2}}{y^{2}} + \frac{{delx}^{2}}{x^{2}}$

(4.5)

Error in a square:

>	$combine ({e1}^{2},' errors')$

$ScientificErrorAnalysis:-Quantity (x^{2}, 2 \sqrt{x^{2} {delx}^{2}})$

(4.6)

General Formula

The general formula for the error $σ_{y}$ in y, where y is a function of variables $x_{i}$ is:

$σ_{y}^{2} = \sum_{i = 1}^{N} {(\frac{\partial}{\partial x_{i}} y)}^{2} σ_{x_{i}}^{2}$

where $σ_{x_{i}}$ is the error in $x_{i}$ and the partials are evaluated at the central values of the $x_{i}$ .

With covariances $σ_{x_{i}, x_{j}}$ taken into account, the formula is:

$σ_{y}^{2} = \sum_{i = 1}^{N} {(\frac{\partial}{\partial x_{i}} y)}^{2} σ_{x_{i}}^{2} + 2 (\sum_{i = 1}^{N - 1} (\sum_{j = i + 1}^{N} (\frac{\partial}{\partial x_{i}} y) (\frac{\partial}{\partial x_{j}} y) σ_{x_{i}, x_{j}}))$

The covariance $σ_{x_{i}, x_{j}}$ may be expressed in terms of the correlation $r_{x_{i}, x_{j}}$ and errors $σ_{x_{i}}$ , $σ_{x_{j}}$ as follows:

$σ_{x_{i}, x_{j}} = r_{x_{i}, x_{j}} σ_{x_{i}} σ_{x_{j}}$

Working with Rounding Rules

Introduction

The initial default rounding rule is digits, for which no rounding is performed.

>	$combine (q1 + q2,' errors')$

$ScientificErrorAnalysis:-Quantity (30.0, 3.162277660)$

(5.1.1)

We can explicitly round a Quantity( ) object:

>	$ApplyRule (, round [2])$

$ScientificErrorAnalysis:-Quantity (30.0, 3.2)$

(5.1.2)

specify a rounding rule in combine/errors:

>	$combine (q1 + q2,' errors',' rule' = round [2])$

$ScientificErrorAnalysis:-Quantity (30.0, 3.2)$

(5.1.3)

>	$combine (q1 \cdot q2,' errors',' rule' = round [2])$

$ScientificErrorAnalysis:-Quantity (200., 36.)$

(5.1.4)

>	$combine ({q1}^{2},' errors',' rule' = round [2])$

$ScientificErrorAnalysis:-Quantity (100., 20.)$

(5.1.5)

or change the default rounding rule using the UseRule function.

>	$UseRule (round [2])$

>	$combine (q1 + q2,' errors')$

$ScientificErrorAnalysis:-Quantity (30.0, 3.2)$

(5.1.6)

>	$combine (q1 \cdot q2,' errors')$

$ScientificErrorAnalysis:-Quantity (200., 36.)$

(5.1.7)

>	$combine ({q1}^{2},' errors')$

$ScientificErrorAnalysis:-Quantity (100., 20.)$

(5.1.8)

>	$UseRule (digits)$

Another Note about Floating-Point Numbers

If a floating-point number is not displayed in scientific notation, it may not be apparent whether trailing zeros are significant figures.

>	$Quantity (23900.0, 1854.0)$

$ScientificErrorAnalysis:-Quantity (23900.0, 1854.0)$

(5.2.1)

>	$ApplyRule (%, round [2])$

$ScientificErrorAnalysis:-Quantity (23900., 1900.)$

(5.2.2)

>	$SFloatMantissa (GetValue ()), SFloatMantissa (GetError ())$

$239, 19$

(5.2.3)

Adding a New Rounding Rule

The rounding rule round3g[n] rounds the error to n figures for a significant figure of 3 or greater and n+1 figures for a significant figure of 1 or 2. Here, we add a similar rule round2g[n], which rounds the error to n figures for a significant figure of 2 or greater and n+1 figures for a significant figure of 1.

$AddRule (round2g = proc (x, y) local xp, yp, n &semi; if not type (procname,' indexed') then error `round2g` must be indexed &semi; end if &semi; n := op (procname) &semi; if not type (n,' posint') then error invalid index to `round2g` &semi; end if &semi; xp, yp := evalf (x), evalf (y) &semi; if iquo (SFloatMantissa (yp), 10 &Hat; ilog10 (SFloatMantissa (yp))) < 2 then n := n + 1 &semi; end if &semi; `if` (xp \neq 0, evalf [n + ilog10 (SFloatMantissa (xp)) - ilog10 (SFloatMantissa (yp)) + SFloatExponent (xp) - SFloatExponent (yp)] (xp), Float (0, - 1)), evalf [n] (yp) &semi; end proc)$

We can use the new rule round2g[n].

>	$HasRule (round2g)$

$true$

(5.3.1)

>	$Quantity (11.22, 2.3)$

$ScientificErrorAnalysis:-Quantity (11.22, 2.3)$

(5.3.2)

>	$ApplyRule (, round2g [1])$

$ScientificErrorAnalysis:-Quantity (11., 2.)$

(5.3.3)

>	$Quantity (11.22, 1.3)$

$ScientificErrorAnalysis:-Quantity (11.22, 1.3)$

(5.3.4)

>	$ApplyRule (, round2g [1])$

$ScientificErrorAnalysis:-Quantity (11.2, 1.3)$

(5.3.5)

Error Analysis with Physical Constants

A Simple Experiment

The period of oscillation T of a simple pendulum for small amplitude is given by:

$T = 2 π \sqrt{\frac{L}{g}}$

where g is the standard acceleration of gravity and L is the length of the string of the pendulum (or, more precisely, the distance from the point of suspension to the center of the pendulum bob).

A student constructs a pendulum to determine a value for g. The student measures L to be 153.2 centimeters (cm). It is decided that the uncertainty in this measurement is primarily due to the difficulty in measuring L to the center of the bob, and that this uncertainty is about plus or minus 2 millimeters (mm). (In other words, it is estimated that the measurement is probably within 2 millimeters of the correct value.) Thus, L is taken to be 153.2 +/- 0.2 centimeters. In standard units:

>	$Lq ≔ Quantity (1.532, 2,' uld')$

$Lq := ScientificErrorAnalysis:-Quantity (1.532, 0.002)$

(6.1.1)

Now, the student sets the pendulum swinging with a small amplitude (a few centimeters). Using a stopwatch, the student measures the time for 30 oscillations to be 75.3 seconds (s). The uncertainty in this measurement is due to the human reaction time in operating the stopwatch, and so the time measurement is taken to be 75.3 +/- 0.5 seconds.

>	$t30 ≔ Quantity (75.3, 0.5)$

$t30 := ScientificErrorAnalysis:-Quantity (75.3, 0.5)$

(6.1.2)

Because the number of oscillations (30) is exact, the period of a single oscillation is:

>	$Tq := combine (\frac{t30}{30},' errors', rule = round [2])$

$Tq := ScientificErrorAnalysis:-Quantity (2.510, 0.017)$

(6.1.3)

The previous formula rearranges to:

$g = \frac{4 π^{2} L}{T^{2}}$

and so, the student calculates the experimentally determined value of g to be:

>	$combine (\frac{4 {(3.1416)}^{2} Lq}{{Tq}^{2}},' errors', rule = {round}_{2})$

$ScientificErrorAnalysis:-Quantity (9.60, 0.13)$

(6.1.4)

The accepted value of g is:

>	$evalf (Constant (g))$

$9.80665$

(6.1.5)

Assuming the error in the experimental value represents a standard deviation, the accepted value thus lies farther than one standard deviation from the experimental central value. The student then must decide whether to repeat the experiment, or consider possibly overlooked sources of measurement error.

Another Example

A sample of a gas is measured to have a volume of 20.3 liters (L), pressure of 105 kilopascals (kPa), and temperature of 297 kelvins (K). If all the measurements have uncertainty of 1%, what is the value and uncertainty of the number of moles in the sample?

The ideal gas law is:

$P V = n R T$

where P is pressure, V is volume, T is temperature, and R is the ideal gas constant.

Constructing quantities-with-error for P, V, and T in standard units:

>	$Pq ≔ Quantity (105e3, 0.01,' relative')$

$Pq := ScientificErrorAnalysis:-Quantity (1.05 10^{5}, 1050.)$

(6.2.1)

>	$Vq ≔ Quantity (20.3&ExponentialE;−3, 0.01,' relative')$

$Vq := ScientificErrorAnalysis:-Quantity (0.0203, 0.000203)$

(6.2.2)

>	$Tq ≔ Quantity (297.0, 0.01,' relative')$

$Tq := ScientificErrorAnalysis:-Quantity (297.0, 2.970)$

(6.2.3)

Thus, the number of moles n is:

>	$combine (\frac{Pq Vq}{Constant (R) Tq},' errors',' rule' = round [2])$

$ScientificErrorAnalysis:-Quantity (0.863, 0.015)$

(6.2.4)

Thus, the sample has 0.863 moles (mol) of gas molecules, with a relative error of:

>	$\frac{GetError ()}{GetValue ()}$

$0.01738122827$

(6.2.5)

or about 2%.

Working at higher Digits

Some physical constants have their values determined to more than 10 significant figures, hence calculations involving these objects at the default setting of Digits may result in a loss of precision. If more precision is desired, Digits may be set to a higher value, 15 for example.

>	$GetConstant (R_{\infty})$

$Rydberg_constant, symbol = R_{\infty}, value = 1.0973731568549 10^{7}, uncertainty = 0.000083, units = \frac{1}{m}$

(6.3.1)

>	$combine (Constant (R_{\infty}) \cdot Constant (c),' errors',' rule = round [2]')$

$ScientificErrorAnalysis:-Quantity (3.289841961 10^{15}, 25000.)$

(6.3.2)

>	$Digits := 15 &colon;$

>	$combine (Constant (R_{\infty}) \cdot Constant (c),' errors',' rule = round [2]')$

$ScientificErrorAnalysis:-Quantity (3.289841960368 10^{15}, 25000.)$

(6.3.3)

>	$Digits ≔ 10 &colon;$

Error in a Derived Physical Constant

Let us examine the formula for the electron mass.

>	$GetConstant (m_{e})$

$electron_mass, symbol = m_{e}, derive = \frac{2 R_{\infty} h}{c α^{2}}$

(7.1)

As a demonstration, we manually construct the same expression using individual physical constant objects.

>	$e3 := \frac{2 Constant (R_{\infty}) Constant (h)}{Constant (c) {Constant (α)}^{2}}$

$e3 := \frac{2 Constant (R_{\infty}) Constant (h)}{Constant (c) {Constant (α)}^{2}}$

(7.2)

Propagating the errors gives us:

>	$combine (e3,' errors',' rule = round [2]')$

$ScientificErrorAnalysis:-Quantity (9.10938188 10^{-31}, 7.2 10^{-38})$

(7.3)

The above calculation is similar to the method ScientificConstants and ScientificErrorAnalysis use internally.

>	$(GetValue, GetError) (Constant (m_{e}))$

$9.109381882 10^{-31}, 7.179221838 10^{-38}$

(7.4)

>	$ApplyRule (Quantity (),' round [2]')$

$ScientificErrorAnalysis:-Quantity (9.10938188 10^{-31}, 7.2 10^{-38})$

(7.5)

Return to ScientificErrorAnalysis

Return to Index for Example Worksheets

Download Help Document

Maple

Maple 附加模块

Student Success Platform

提高存留率

MapleSim

MapleSim 附加模块

系统工程

项目服务

Maple T.A. and Möbius

教育

行业

汽车与航空航天

机器人

机器设计和工业自动化

其他

应用领域

产品价格

购买

院系/全校正版授权

Maplesoft精英维护升级计划

支持

产品培训

产品在线帮助

研讨会与活动

出版

资源中心

示例和应用

社区

关于 Maplesoft

媒体中心

用户社区

联系我们

Online Help

All Products Maple MapleSim

Maple

数学软件

Maple 附加模块

Student Success Platform

提高存留率

MapleSim

多学科系统级建模仿真

MapleSim 附加模块

系统工程

项目服务

Maple T.A. and Möbius

教育

行业

汽车与航空航天

机器人

机器设计和工业自动化

其他

应用领域

产品价格

购买

院系/全校正版授权

Maplesoft精英维护升级计划

支持

产品培训

产品在线帮助

研讨会与活动

出版

资源中心

示例和应用

社区

关于 Maplesoft

媒体中心

用户社区

联系我们

Online Help

All Products Maple MapleSim