{VERSION 6 0 "IBM INTEL NT" "6.0" } {USTYLETAB {CSTYLE "Maple Input" -1 0 "Courier" 0 1 255 0 0 1 0 1 0 0 1 0 0 0 0 1 }{CSTYLE "2D Math" -1 2 "Times" 0 1 0 0 0 0 0 0 2 0 0 0 0 0 0 1 }{CSTYLE "2D Comment" 2 18 "" 0 1 0 0 0 0 0 0 0 0 0 0 0 0 0 1 } {CSTYLE "2D Input" 2 19 "" 0 1 255 0 0 1 0 0 0 0 0 0 0 0 0 1 }{CSTYLE "2D Output" 2 20 "" 0 1 0 0 255 1 0 0 0 0 0 0 0 0 0 1 }{CSTYLE "" -1 256 "" 1 16 0 0 0 0 0 1 0 0 0 0 0 0 0 1 }{CSTYLE "" -1 257 "" 0 1 0 0 0 0 1 0 0 0 0 0 0 0 0 1 }{CSTYLE "" -1 258 "" 0 1 0 0 0 0 1 0 0 0 0 0 0 0 0 1 }{CSTYLE "" -1 259 "" 0 1 0 0 0 0 1 0 0 0 0 0 0 0 0 1 } {CSTYLE "" -1 260 "" 0 1 0 0 0 0 1 0 0 0 0 0 0 0 0 1 }{CSTYLE "" -1 261 "" 0 1 0 0 0 0 1 0 0 0 0 0 0 0 0 1 }{CSTYLE "" -1 262 "" 0 1 0 0 0 0 1 0 0 0 0 0 0 0 0 1 }{CSTYLE "" -1 263 "" 0 1 0 0 0 0 1 0 0 0 0 0 0 0 0 1 }{CSTYLE "" -1 264 "" 0 1 0 0 0 0 0 0 1 0 0 0 0 0 0 1 } {CSTYLE "" -1 265 "" 0 1 0 0 0 0 1 0 0 0 0 0 0 0 0 1 }{CSTYLE "" -1 266 "" 0 1 0 0 0 0 1 0 0 0 0 0 0 0 0 1 }{CSTYLE "" -1 267 "" 0 1 0 0 0 0 1 0 0 0 0 0 0 0 0 1 }{CSTYLE "" -1 268 "" 0 1 0 0 0 0 1 0 0 0 0 0 0 0 0 1 }{CSTYLE "" -1 269 "" 0 1 0 0 0 0 0 0 1 0 0 0 0 0 0 1 } {CSTYLE "" -1 270 "" 0 1 0 0 0 0 1 0 0 0 0 0 0 0 0 1 }{CSTYLE "" -1 271 "" 0 1 0 0 0 0 1 0 0 0 0 0 0 0 0 1 }{CSTYLE "" -1 272 "" 0 1 0 0 0 0 1 0 0 0 0 0 0 0 0 1 }{CSTYLE "" -1 273 "" 0 1 0 0 0 0 1 0 0 0 0 0 0 0 0 1 }{CSTYLE "" -1 274 "" 0 1 0 0 0 0 0 0 1 0 0 0 0 0 0 1 } {CSTYLE "" -1 275 "" 0 1 0 0 0 0 0 0 1 0 0 0 0 0 0 1 }{CSTYLE "" -1 276 "" 0 1 0 0 0 0 0 0 1 0 0 0 0 0 0 1 }{CSTYLE "" -1 277 "" 0 1 0 0 0 0 1 0 0 0 0 0 0 0 0 1 }{CSTYLE "" -1 278 "" 0 1 0 0 0 0 0 0 1 0 0 0 0 0 0 1 }{CSTYLE "" -1 279 "" 0 1 0 0 0 0 1 0 0 0 0 0 0 0 0 1 } {CSTYLE "" -1 280 "" 0 1 0 0 0 0 1 0 0 0 0 0 0 0 0 1 }{CSTYLE "" -1 281 "" 0 1 0 0 0 0 1 0 0 0 0 0 0 0 0 1 }{CSTYLE "" -1 282 "" 0 1 0 0 0 0 1 0 0 0 0 0 0 0 0 1 }{CSTYLE "" -1 283 "" 0 1 0 0 0 0 1 0 0 0 0 0 0 0 0 1 }{CSTYLE "" -1 284 "" 0 1 0 0 0 0 1 0 0 0 0 0 0 0 0 1 } {CSTYLE "" -1 285 "" 0 1 0 0 0 0 1 0 0 0 0 0 0 0 0 1 }{CSTYLE "" -1 286 "" 0 1 0 0 0 0 1 0 0 0 0 0 0 0 0 1 }{CSTYLE "" -1 287 "" 0 1 0 0 0 0 0 0 1 0 0 0 0 0 0 1 }{CSTYLE "" -1 288 "" 0 1 0 0 0 0 0 0 1 0 0 0 0 0 0 1 }{CSTYLE "" -1 289 "" 0 1 0 0 0 0 1 0 0 0 0 0 0 0 0 1 } {CSTYLE "" -1 290 "" 0 1 0 0 0 0 1 0 0 0 0 0 0 0 0 1 }{CSTYLE "" -1 291 "" 0 1 0 0 0 0 1 0 0 0 0 0 0 0 0 1 }{CSTYLE "" -1 292 "" 0 1 0 0 0 0 1 0 0 0 0 0 0 0 0 1 }{CSTYLE "" -1 293 "" 0 1 0 0 0 0 1 0 0 0 0 0 0 0 0 1 }{CSTYLE "" -1 294 "" 0 1 0 0 0 0 0 0 1 0 0 0 0 0 0 1 } {CSTYLE "" -1 295 "" 0 1 0 0 0 0 0 0 1 0 0 0 0 0 0 1 }{PSTYLE "Normal " -1 0 1 {CSTYLE "" -1 -1 "Times" 1 12 0 0 0 1 2 2 2 2 2 2 1 1 1 1 }1 1 0 0 0 0 1 0 1 0 2 2 0 1 }{PSTYLE "Text Output" -1 2 1 {CSTYLE "" -1 -1 "Courier" 1 10 0 0 255 1 0 0 0 0 0 1 3 0 3 0 }1 0 0 -1 -1 -1 0 0 0 0 0 0 -1 0 }{PSTYLE "Warning" 2 7 1 {CSTYLE "" -1 -1 "" 0 1 0 0 255 1 0 0 0 0 0 0 1 0 0 0 }0 0 0 -1 -1 -1 0 0 0 0 0 0 -1 0 }{PSTYLE "Maple O utput" -1 11 1 {CSTYLE "" -1 -1 "Times" 1 12 0 0 0 1 2 2 2 2 2 2 1 1 1 1 }3 3 0 0 0 0 1 0 1 0 2 2 0 1 }{PSTYLE "Maple Plot" -1 13 1 {CSTYLE "" -1 -1 "Times" 1 12 0 0 0 1 2 2 2 2 2 2 1 1 1 1 }3 1 0 0 0 0 1 0 1 0 2 2 0 1 }} {SECT 0 {EXCHG {PARA 0 "" 0 "" {TEXT 256 51 "Experimental uncertaintie s: the normal distribution" }}{PARA 0 "" 0 "" {TEXT -1 0 "" }}{PARA 0 "" 0 "" {TEXT -1 1435 "We consider the statistical analysis of random \+ experimental uncertainties in this worksheet. This will allow us to un derstand how to analyze repeated measurements of the same quantity, as suming that the uncertainties are caused by random measurement errors. These can include fluctuations in currents in a circuit (depending on someone else's loads in the power grid), variable reaction time when \+ a stopwatch is used, parallax problems when reading a length. Some of \+ these can also give rise to systematic errors, whereby all data values taken are consistently under- (or over-) estimated by the experimente r. These are much harder to track. The most typical systematic error t hat I observed in our laboratory over the years was the uncalibrated o ut-of-tune timebase of a time chart recorder (used in a Cavendish expe riment to track the motion of two lead balls via a photodiode pick-up \+ of a laser beam). This error was the easily spotted with a stopwatch, \+ but no student was able to pick up on it (they simply believed the rea ding on the knob). A 20% systematic error in the timebase resulted in \+ a wrong value of Newton's constant. Most students had a more accurate \+ value from the manual experiment (where they tracked the motion by mar kings on graph paper pasted onto a wall), and were surprised that the \+ automated experiment gave inferior results (and often inflated their u ncertainties to avoid conflict with the literature value)." }}{PARA 0 "" 0 "" {TEXT -1 0 "" }}{PARA 0 "" 0 "" {TEXT -1 386 "A detailed analy sis of experimental uncertainties is crucial in order to assess whethe r systematic errors are present in someone's experiment or not. We wil l not enter into the question here as to why the experimental uncertai nties can (usually) be assumed to follow a normal (Gaussian) distribut ion. The prerequisites for this situation are that there are many sour ces of small errors." }}{PARA 0 "" 0 "" {TEXT -1 0 "" }}{PARA 0 "" 0 " " {TEXT -1 412 "We are interested in a single quantity that is measure d repeatedly. For the measurement of linear (or other) relationships a n extension of the arguments presented here leads to the linear least \+ squares (or quadratic) fit. In this worksheet we will simulate random \+ errors by using normally distributed random numbers to test the theore tical ideas which are proposed for the analysis of experimental uncert ainties." }}{PARA 0 "" 0 "" {TEXT -1 0 "" }}{PARA 0 "" 0 "" {TEXT -1 425 "Every experimenter will intuitively form averages from repeated m easurements in order to arrive at an improved experimental result for \+ the quantity to be measured. Understanding experimental uncertainties \+ will not only lead to a proper justification, but also to a measure of the expected accuracy of the averaged result. For this purpose we hav e to adopt a model which tells us how the repeated measurements are di stributed." }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 8 "restart;" }}} {EXCHG {PARA 0 "" 0 "" {TEXT -1 109 "To illustrate the concept of aver age and standard deviation we follow the textbook approach (John R. Ta ylor: " }{TEXT 257 33 "An Introduction to Error Analysis" }{TEXT -1 418 ", chapters 4 and 5, University Science Books 1982, 2nd ed. 1998). Suppose that we have carried out a measurement for the spring constan t of a typical spring (out of a box containing a set of identically ma nufactured ones) repeatedly by timing the period of oscillation while \+ a well-known mass was attached. The series of 10 measurements is given as (in Newton/meter) the following list of values with 2-digit accura cy:" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 36 "kL:=[86,85,84,89,86, 88,88,85,83,85];" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%#kLG7,\"#')\"#&) \"#%)\"#*)F&\"#))F*F'\"#$)F'" }}}{EXCHG {PARA 0 "" 0 "" {TEXT -1 707 " The choice of 2-digit accuracy was reasonable. Our timing device gave \+ more digits of readout, but due to the substantial fluctuations in the second digit we decided to round the numbers to two digits. (Common m istakes made by students is to: (a) carry all 10 digits that their cal culator produced while dividing two numbers that were measured only to 2 or 3 figures; (b) to trust digitial devices that display much more \+ than the stated accuracy in the manual (digital voltmeters and ammeter s); (c) ignoring the fact that timing motion with digital stopwatches \+ does not eliminate human error, and reaction time; (d) measurements wi th calipers can be tricky due to bending, improper holding of the inst rument." }}{PARA 0 "" 0 "" {TEXT -1 0 "" }}{PARA 0 "" 0 "" {TEXT -1 400 "Suppose that all numbers would have been equal in the set of meas urements. Then there would be no deviation from the average. In that c ase we might ask ourselves whether more digits can be read from the in strument, and whether it is worthwhile to make an attempt at more prec ision. At this point we probably could proceed and carry out an uncert ainty estimate with the data set with 3-digit accuracy." }}{PARA 0 "" 0 "" {TEXT -1 0 "" }}{PARA 0 "" 0 "" {TEXT -1 50 "The average value is presumably our best estimate:" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 12 "N:=nops(kL);" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%\"NG\"#5" }} }{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 34 "k_a:=evalf(add(kL[i],i=1..N )/N,3);" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%$k_aG$\"$f)!\"\"" }}} {EXCHG {PARA 0 "" 0 "" {TEXT -1 141 "Before we proceed to calculate th e deviation we illustrate the data by forming a histogram (a frequency distribution or binning of the data)." }}{PARA 0 "" 0 "" {TEXT -1 189 "We choose a bin size and count how many measurements fall into ea ch bin. We could pick 10 bins around the calculated average value, or \+ simply cover the range between 80 and 90 with 20 bins." }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 8 "dk:=0.5;" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%#dkG$\"\"&!\"\"" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 22 "k_i:=i->80+(i-0.5)*dk;" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%$k_iGj+6#%\"iG6\"6$%)operatorG%&arrowGF(,&\"#!)\"\"\"*&,&9$F.$ \"\"&!\"\"F4F.%#dkGF.F.F(F(F(6#\"\"!" }}}{EXCHG {PARA 0 "" 0 "" {TEXT -1 173 "One way to proceed to bin the data is to subtract 80 from the \+ value, multiply the remainder by 2 (WHICH IS 1/DK!), and round it to a n integer (by adding 0.5 and truncating):" }}}{EXCHG {PARA 0 "> " 0 " " {MPLTEXT 1 0 36 "i_k:=k->trunc((k-80-0.5*dk)/dk+0.5);" }}{PARA 11 " " 1 "" {XPPMATH 20 "6#>%$i_kGj+6#%\"kG6\"6$%)operatorG%&arrowGF(-%&tru ncG6#,&*&,(9$\"\"\"\"#!)!\"\"*&$\"\"&F5F3%#dkGF3F5F3F9F5F3F7F3F(F(F(6# \"\"!" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 39 "j:=3; kL[j]; #PICK THE FIRST DATA POINT" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%\"jG\"\"$" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#\"#%)" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 41 "i_k(kL[j]); # it goes into bin number 12." }}{PARA 11 "" 1 "" {XPPMATH 20 "6#\"\")" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 48 "k_i(i_k(kL[j])),k_i(i_k(kL[j])+1); #YES, INDEED!" }}{PARA 11 " " 1 "" {XPPMATH 20 "6$$\"%v$)!\"#$\"%D%)F%" }}}{EXCHG {PARA 0 "" 0 "" {TEXT -1 169 "The defined functions are consistent, i.e., the index us ed in the binning process works out such that it assigns data values t o fall to the right of the bordering value." }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 60 "for i from 1 to 20 do: C[i]:=0: od: #INITIALIZE BI N COUNTERS" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 90 "for n from 1 \+ to N do: ict:=i_k(kL[n]); if ict>0 and ict<=20 then C[ict]:=C[ict]+1; \+ fi; od:" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 62 "with(plots): P1: =plot([k_a,t,t=0..2.5],color=red,thickness=3):" }}{PARA 7 "" 1 "" {TEXT -1 50 "Warning, the name changecoords has been redefined\n" }}} {EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 100 "P2:=plot([seq([k_i(ict)+0.5 *dk,t*C[ict],t=0..1],ict=1..20)],color=blue,thickness=2): display(P1,P 2);" }}{PARA 13 "" 1 "" {GLPLOT2D 564 304 304 {PLOTDATA 2 "69-%'CURVES G6%7S7$$\"3e++++++!f)!#;$\"\"!F,7$F($\"3WmmmT&)G\\a!#>7$F($\"3PL$ek`o! >5!#=7$F($\"3omm\"z>)G_:F47$F($\"3-nmT&QU!*3#F47$F($\"3HL$eRZXKi#F47$F ($\"3xm;z>,_=JF47$F($\"3v**\\7G$[8j$F47$F($\"35n;z%*frhTF47$F($\"3A+]i lFQ!p%F47$F($\"3@ML$3_\"=M_F47$F($\"3HnmTg(fJr&F47$F($\"3k++]7eP_iF47$ F($\"3Q++]Pf!Qz'F47$F($\"3@++](=ubJ(F47$F($\"37n;zW(*Q*y(F47$F($\"3#QL L3F-GN)F47$F($\"3=MLL$e'3I))F47$F($\"3?+]7.\"Feo7$F($\"3km;/Egw[7Feo7$F($\"3zm\"z%*f%)QI\"Feo7$F($ \"3/+voza'=N\"Feo7$F($\"3(om\"zWho.9Feo7$F($\"3-++]i>Ad9Feo7$F($\"32+] i:jf4:Feo7$F($\"39+DJ&>r-c\"Feo7$F($\"3++]P4q`;;Feo7$F($\"3;LL$eM%4n;F eo7$F($\"37++v$4v5s\"Feo7$F($\"3cm\"zWn*)*pFeo7$F($\"3CLL3-=!y(>Feo7$F($\"3))* \\7G8O;.#Feo7$F($\"3!pmm;*\\[$3#Feo7$F($\"3*pmT&Qz]O@Feo7$F($\"3iLekG= 4*=#Feo7$F($\"3F++]i4TPAFeo7$F($\"3qL$3F9!z#H#Feo7$F($\"3'pmm;%>KUBFeo 7$F($\"3/+DJqJ8&R#Feo7$F($\"3G+voa-oXCFeo7$F($\"3++++++++DFeo-%'COLOUR G6&%$RGBG$\"*++++\"!\")F+F+-%*THICKNESSG6#\"\"$-F$6%7S7$$\"3+++++++]!) F*F+FhuFhuFhuFhuFhuFhuFhuFhuFhuFhuFhuFhuFhuFhuFhuFhuFhuFhuFhuFhuFhuFhu FhuFhuFhuFhuFhuFhuFhuFhuFhuFhuFhuFhuFhuFhuFhuFhuFhuFhuFhuFhuFhuFhuFhuF huFhuFhu-F[u6&F]uF+F+F^u-Fbu6#\"\"#-F$6%7S7$$\"#\")F,F+FcvFcvFcvFcvFcv FcvFcvFcvFcvFcvFcvFcvFcvFcvFcvFcvFcvFcvFcvFcvFcvFcvFcvFcvFcvFcvFcvFcvF cvFcvFcvFcvFcvFcvFcvFcvFcvFcvFcvFcvFcvFcvFcvFcvFcvFcvFcvFcvF[vF]v-F$6% 7S7$$\"3+++++++]\")F*F+FivFivFivFivFivFivFivFivFivFivFivFivFivFivFivFi vFivFivFivFivFivFivFivFivFivFivFivFivFivFivFivFivFivFivFivFivFivFivFiv FivFivFivFivFivFivFivFivFivF[vF]v-F$6%7S7$$\"##)F,F+F_wF_wF_wF_wF_wF_w F_wF_wF_wF_wF_wF_wF_wF_wF_wF_wF_wF_wF_wF_wF_wF_wF_wF_wF_wF_wF_wF_wF_wF _wF_wF_wF_wF_wF_wF_wF_wF_wF_wF_wF_wF_wF_wF_wF_wF_wF_wF_wF[vF]v-F$6%7S7 $$\"3+++++++]#)F*F+FewFewFewFewFewFewFewFewFewFewFewFewFewFewFewFewFew FewFewFewFewFewFewFewFewFewFewFewFewFewFewFewFewFewFewFewFewFewFewFewF ewFewFewFewFewFewFewFewF[vF]v-F$6%7S7$$\"#$)F,F+7$F\\x$\"3emmm;arz@F07 $F\\x$\"3[LL$e9ui2%F07$F\\x$\"3nmmm\"z_\"4iF07$F\\x$\"3[mmmT&phN)F07$F \\x$\"3CLLe*=)H\\5F47$F\\x$\"3gmm\"z/3uC\"F47$F\\x$\"3%)***\\7LRDX\"F4 7$F\\x$\"3]mm\"zR'ok;F47$F\\x$\"3w***\\i5`h(=F47$F\\x$\"3WLLL3En$4#F47 $F\\x$\"3qmm;/RE&G#F47$F\\x$\"3\")*****\\K]4]#F47$F\\x$\"3$******\\PAv r#F47$F\\x$\"3)******\\nHi#HF47$F\\x$\"3jmm\"z*ev:JF47$F\\x$\"3?LLL347 TLF47$F\\x$\"3,LLLLY.KNF47$F\\x$\"3w***\\7o7Tv$F47$F\\x$\"3'GLLLQ*o]RF 47$F\\x$\"3A++D\"=lj;%F47$F\\x$\"31++vV&RY2aF47$F\\x$\"39mm;zXu9cF47$F\\x$\"3l******\\y))GeF47 $F\\x$\"3'*)***\\i_QQgF47$F\\x$\"3@***\\7y%3TiF47$F\\x$\"35****\\P![hY 'F47$F\\x$\"3kKLL$Qx$omF47$F\\x$\"3!)*****\\P+V)oF47$F\\x$\"3?mm\"zpe* zqF47$F\\x$\"3%)*****\\#\\'QH(F47$F\\x$\"3GKLe9S8&\\(F47$F\\x$\"3R*** \\i?=bq(F47$F\\x$\"3\"HLL$3s?6zF47$F\\x$\"3a***\\7`Wl7)F47$F\\x$\"3#pm mm'*RRL)F47$F\\x$\"3Qmm;a<.Y&)F47$F\\x$\"3=LLe9tOc()F47$F\\x$\"3u***** *\\Qk\\*)F47$F\\x$\"3CLL$3dg6<*F47$F\\x$\"3ImmmmxGp$*F47$F\\x$\"3A++D \"oK0e*F47$F\\x$\"3A++v=5s#y*F47$F\\x$\"\"\"F,F[vF]v-F$6%7S7$$\"3+++++ ++]$)F*F+FaalFaalFaalFaalFaalFaalFaalFaalFaalFaalFaalFaalFaalFaalFaalF aalFaalFaalFaalFaalFaalFaalFaalFaalFaalFaalFaalFaalFaalFaalFaalFaalFaa lFaalFaalFaalFaalFaalFaalFaalFaalFaalFaalFaalFaalFaalFaalFaalF[vF]v-F$ 6%7S7$$\"#%)F,F+7$FhalF_x7$FhalFbx7$FhalFex7$FhalFhx7$FhalF[y7$FhalF^y 7$FhalFay7$FhalFdy7$FhalFgy7$FhalFjy7$FhalF]z7$FhalF`z7$FhalFcz7$FhalF fz7$FhalFiz7$FhalF\\[l7$FhalF_[l7$FhalFb[l7$FhalFe[l7$FhalFh[l7$FhalF[ \\l7$FhalF^\\l7$FhalFa\\l7$FhalFd\\l7$FhalFg\\l7$FhalFj\\l7$FhalF]]l7$ FhalF`]l7$FhalFc]l7$FhalFf]l7$FhalFi]l7$FhalF\\^l7$FhalF_^l7$FhalFb^l7 $FhalFe^l7$FhalFh^l7$FhalF[_l7$FhalF^_l7$FhalFa_l7$FhalFd_l7$FhalFg_l7 $FhalFj_l7$FhalF]`l7$FhalF``l7$FhalFc`l7$FhalFf`l7$FhalFi`l7$FhalF\\al F[vF]v-F$6%7S7$$\"3+++++++]%)F*F+F]elF]elF]elF]elF]elF]elF]elF]elF]elF ]elF]elF]elF]elF]elF]elF]elF]elF]elF]elF]elF]elF]elF]elF]elF]elF]elF]e lF]elF]elF]elF]elF]elF]elF]elF]elF]elF]elF]elF]elF]elF]elF]elF]elF]elF ]elF]elF]elF]elF[vF]v-F$6%7S7$$\"#&)F,F+7$Fdel$\"3s******\\i9RlF07$Fde l$\"3/++vVA)GA\"F47$Fdel$\"3+++]Peui=F47$Fdel$\"3n****\\i3&o]#F47$Fdel $\"3%)***\\(oX*y9$F47$Fdel$\"3z***\\P9CAu$F47$Fdel$\"3E***\\P*zhdVF47$ Fdel$\"3]***\\P>fS*\\F47$Fdel$\"3s)**\\(=$f%GcF47$Fdel$\"3Q+++Dy,\"G'F 47$Fdel$\"33++]7&)\\\"Fe o7$Fdel$\"3w**\\P>:mk:Feo7$Fdel$\"3'***\\iv&QAi\"Feo7$Fdel$\"3%)***\\P PBWo\"Feo7$Fdel$\"3!******\\Nm'[Feo7$Fdel$\"3z*****\\@80+#Feo7$Fd el$\"31++]7,Hl?Feo7$Fdel$\"3()**\\P4w)R7#Feo7$Fdel$\"3;++]x%f\")=#Feo7 $Fdel$\"3!)**\\P/-a[AFeo7$Fdel$\"3g**\\(=Yb;J#Feo7$Fdel$\"3')****\\i@O tBFeo7$Fdel$\"3')**\\PfL'zV#Feo7$Fdel$\"3>+++!*>=+DFeo7$Fdel$\"3-++DE& 4Qc#Feo7$Fdel$\"3=+]P%>5pi#Feo7$Fdel$\"39+++bJ*[o#Feo7$Fdel$\"33++Dr\" [8v#Feo7$Fdel$\"3++++Ijy5GFeo7$Fdel$\"31+]P/)fT(GFeo7$Fdel$\"31+]i0j\" [$HFeo7$Fdel$FduF,F[vF]v-F$6%7S7$$\"3+++++++]&)F*F+Fh^mFh^mFh^mFh^mFh^ mFh^mFh^mFh^mFh^mFh^mFh^mFh^mFh^mFh^mFh^mFh^mFh^mFh^mFh^mFh^mFh^mFh^mF h^mFh^mFh^mFh^mFh^mFh^mFh^mFh^mFh^mFh^mFh^mFh^mFh^mFh^mFh^mFh^mFh^mFh^ mFh^mFh^mFh^mFh^mFh^mFh^mFh^mFh^mF[vF]v-F$6%7S7$$\"#')F,F+7$F__m$\"39L LLL3VfVF07$F__m$\"3'pmm;H[D:)F07$F__m$\"3LLLLe0$=C\"F47$F__m$\"3ILLL3R Br;F47$F__m$\"3Ymm;zjf)4#F47$F__m$\"3=LL$e4;[\\#F47$F__m$\"3p****\\i'y ]!HF47$F__m$\"3,LL$ezs$HLF47$F__m$\"3_****\\7iI_PF47$F__m$\"3#pmmm@Xt= %F47$F__m$\"3QLLL3y_qXF47$F__m$\"3i******\\1!>+&F47$F__m$\"3()******\\ Z/NaF47$F__m$\"3'*******\\$fC&eF47$F__m$\"3ELL$ez6:B'F47$F__m$\"3Smmm; =C#o'F47$F__m$\"3-mmmm#pS1(F47$F__m$\"3]****\\i`A3vF47$F__m$\"3slmmm(y 8!zF47$F__m$\"3V++]i.tK$)F47$F__m$\"39++](3zMu)F47$F__m$\"3#pmm;H_?<*F 47$F__m$\"3emm;zihl&*F47$F__m$\"39LLL3#G,***F47$F__m$\"3Feo7$F__m$\"3/++v.Uac>Feo7$F_ _m$F_vF,F[vF]v-F$6%7S7$$\"3+++++++]')F*F+FchmFchmFchmFchmFchmFchmFchmF chmFchmFchmFchmFchmFchmFchmFchmFchmFchmFchmFchmFchmFchmFchmFchmFchmFch mFchmFchmFchmFchmFchmFchmFchmFchmFchmFchmFchmFchmFchmFchmFchmFchmFchmF chmFchmFchmFchmFchmFchmF[vF]v-F$6%7S7$$\"#()F,F+FihmFihmFihmFihmFihmFi hmFihmFihmFihmFihmFihmFihmFihmFihmFihmFihmFihmFihmFihmFihmFihmFihmFihm FihmFihmFihmFihmFihmFihmFihmFihmFihmFihmFihmFihmFihmFihmFihmFihmFihmFi hmFihmFihmFihmFihmFihmFihmFihmF[vF]v-F$6%7S7$$\"3+++++++]()F*F+F_imF_i mF_imF_imF_imF_imF_imF_imF_imF_imF_imF_imF_imF_imF_imF_imF_imF_imF_imF _imF_imF_imF_imF_imF_imF_imF_imF_imF_imF_imF_imF_imF_imF_imF_imF_imF_i mF_imF_imF_imF_imF_imF_imF_imF_imF_imF_imF_imF[vF]v-F$6%7S7$$\"#))F,F+ 7$FfimFb_m7$FfimFe_m7$FfimFh_m7$FfimF[`m7$FfimF^`m7$FfimFa`m7$FfimFd`m 7$FfimFg`m7$FfimFj`m7$FfimF]am7$FfimF`am7$FfimFcam7$FfimFfam7$FfimFiam 7$FfimF\\bm7$FfimF_bm7$FfimFbbm7$FfimFebm7$FfimFhbm7$FfimF[cm7$FfimF^c m7$FfimFacm7$FfimFdcm7$FfimFgcm7$FfimFjcm7$FfimF]dm7$FfimF`dm7$FfimFcd m7$FfimFfdm7$FfimFidm7$FfimF\\em7$FfimF_em7$FfimFbem7$FfimFeem7$FfimFh em7$FfimF[fm7$FfimF^fm7$FfimFafm7$FfimFdfm7$FfimFgfm7$FfimFjfm7$FfimF] gm7$FfimF`gm7$FfimFcgm7$FfimFfgm7$FfimFigm7$FfimF\\hm7$FfimF_hmF[vF]v- F$6%7S7$$\"3+++++++]))F*F+F[]nF[]nF[]nF[]nF[]nF[]nF[]nF[]nF[]nF[]nF[]n F[]nF[]nF[]nF[]nF[]nF[]nF[]nF[]nF[]nF[]nF[]nF[]nF[]nF[]nF[]nF[]nF[]nF[ ]nF[]nF[]nF[]nF[]nF[]nF[]nF[]nF[]nF[]nF[]nF[]nF[]nF[]nF[]nF[]nF[]nF[]n F[]nF[]nF[vF]v-F$6%7S7$$\"#*)F,F+7$Fb]nF_x7$Fb]nFbx7$Fb]nFex7$Fb]nFhx7 $Fb]nF[y7$Fb]nF^y7$Fb]nFay7$Fb]nFdy7$Fb]nFgy7$Fb]nFjy7$Fb]nF]z7$Fb]nF` z7$Fb]nFcz7$Fb]nFfz7$Fb]nFiz7$Fb]nF\\[l7$Fb]nF_[l7$Fb]nFb[l7$Fb]nFe[l7 $Fb]nFh[l7$Fb]nF[\\l7$Fb]nF^\\l7$Fb]nFa\\l7$Fb]nFd\\l7$Fb]nFg\\l7$Fb]n Fj\\l7$Fb]nF]]l7$Fb]nF`]l7$Fb]nFc]l7$Fb]nFf]l7$Fb]nFi]l7$Fb]nF\\^l7$Fb ]nF_^l7$Fb]nFb^l7$Fb]nFe^l7$Fb]nFh^l7$Fb]nF[_l7$Fb]nF^_l7$Fb]nFa_l7$Fb ]nFd_l7$Fb]nFg_l7$Fb]nFj_l7$Fb]nF]`l7$Fb]nF``l7$Fb]nFc`l7$Fb]nFf`l7$Fb ]nFi`l7$Fb]nF\\alF[vF]v-F$6%7S7$$\"3+++++++]*)F*F+Fg`nFg`nFg`nFg`nFg`n Fg`nFg`nFg`nFg`nFg`nFg`nFg`nFg`nFg`nFg`nFg`nFg`nFg`nFg`nFg`nFg`nFg`nFg `nFg`nFg`nFg`nFg`nFg`nFg`nFg`nFg`nFg`nFg`nFg`nFg`nFg`nFg`nFg`nFg`nFg`n Fg`nFg`nFg`nFg`nFg`nFg`nFg`nFg`nF[vF]v-F$6%7S7$$\"#!*F,F+F]anF]anF]anF ]anF]anF]anF]anF]anF]anF]anF]anF]anF]anF]anF]anF]anF]anF]anF]anF]anF]a nF]anF]anF]anF]anF]anF]anF]anF]anF]anF]anF]anF]anF]anF]anF]anF]anF]anF ]anF]anF]anF]anF]anF]anF]anF]anF]anF]anF[vF]v-%+AXESLABELSG6%Q!6\"Fcan -%%FONTG6#%(DEFAULTG-%%VIEWG6$FhanFhan" 1 2 0 1 10 0 2 9 1 4 2 1.000000 45.000000 45.000000 0 0 "Curve 1" "Curve 2" "Curve 3" "Curve \+ 4" "Curve 5" "Curve 6" "Curve 7" "Curve 8" "Curve 9" "Curve 10" "Curve 11" "Curve 12" "Curve 13" "Curve 14" "Curve 15" "Curve 16" "Curve 17 " "Curve 18" "Curve 19" "Curve 20" "Curve 21" }}}}{EXCHG {PARA 0 "" 0 "" {TEXT -1 525 "Is this data set following a Gaussian distribution? A t first we might think that it does not (suspicious gap to the right o f the average). On second thought (and supported by testing done on th e binomial and Poisson distributions in a separate worksheet using a c hi-squared test of a hypothesis) we remain open-minded: the data set i s much to small to rule out that a much larger data set obtained with \+ the same measurement procedure will reach a Gaussian as its limiting d istribution (limit of infinitely many measurements!)" }}{PARA 0 "" 0 " " {TEXT -1 0 "" }}{PARA 0 "" 0 "" {TEXT -1 467 "The standard deviation of the data is a quantity that looks at the distances of the individu al measurements from the calculated average. The distances are added q uadratically in order not to cancel positive against negative errors. \+ We average the squared deviations, and then take the square root to ob tain a width (the deviation). When averaging the sum of the squared de viation we make a small adjustment to account for the fact that one de gree of freedom from the " }{TEXT 258 1 "N" }{TEXT -1 97 " available f rom the data set is lost, as the average was obtained from the data se t itself. Thus:" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 54 "sigma:=e valf(sqrt(add((kL[n]-k_a)^2,n=1..N)/(N-1)),3);" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%&sigmaG$\"$\">!\"#" }}}{EXCHG {PARA 0 "" 0 "" {TEXT -1 35 "This number represents the average " }{TEXT 262 35 "deviation o f individual data points" }{TEXT -1 276 " from the average value. It i s the typical uncertainty associated with an individual measurement. I f we make single measurements for other springs from the box of suppos edly identical springs, and find that the value is within sigma away f rom our calculated average (based on " }{TEXT 259 1 "N" }{TEXT -1 223 " measurements on the first spring), then we should not be concerned a t all about the question whether the springs are made identically to a sufficient degree, as they appear to be identical to within our measu ring accuracy." }}{PARA 0 "" 0 "" {TEXT -1 0 "" }}{PARA 0 "" 0 "" {TEXT -1 182 "We did promise, however, more. We promised to be able to assess from the distribution the accuracy of the determined average v alue, which is presumably the best value we have (after " }{TEXT 260 1 "N" }{TEXT -1 35 " measurements). The uncertainty or " }{TEXT 261 30 "standard deviation of the mean" }{TEXT -1 91 " - provided the data do have a Gaussian as its limiting distribution - can be estimated as :" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 25 "sig_avg:='sigma/sqrt(N )';" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%(sig_avgG*&%&sigmaG\"\"\"-%%s qrtG6#%\"NG!\"\"" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 15 "evalf(s ig_avg);" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#$\"+J.&*Rg!#5" }}}{EXCHG {PARA 0 "" 0 "" {TEXT -1 78 "This is smaller than the average uncertai nty of the individual data points by " }{TEXT 19 7 "sqrt(N)" }{TEXT -1 111 ", which is a good reason to take as many data as possible as t he reduction in this uncertainty with increasing " }{TEXT 263 1 "N" } {TEXT -1 9 " is slow." }}{PARA 0 "" 0 "" {TEXT -1 86 "Thus, the mean ( average) together with its deviation can be listed for our example as \+ " }}{PARA 0 "" 0 "" {TEXT -1 27 "(85.9 plusminus 0.6) [N/m]." }}} {EXCHG {PARA 0 "" 0 "" {TEXT 264 19 "Normal distribution" }}{PARA 0 " " 0 "" {TEXT -1 0 "" }}{PARA 0 "" 0 "" {TEXT -1 445 "To compare the hi stogram to a normal distribution it has to be properly normalized. We \+ have just determined the two parameters that control the shape of the \+ limiting Gaussian distribution, namely the position of the average (a \+ Gaussian is symmetric), and the width (the deviation). The Gaussian (o r normal) distribution is normalized by an integral to unit area. The \+ frequency distribution at present is normalized to the number of measu rements:" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 22 "add(C[ict],ict= 1..20);" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#\"#5" }}}{EXCHG {PARA 0 "" 0 "" {TEXT -1 41 "It is expedient to normalize it to unity:" }}} {EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 55 "for ict from 1 to 20 do: C[i ct]:=evalf(C[ict]/N,3); od:" }}}{EXCHG {PARA 0 "" 0 "" {TEXT -1 95 "Th e average value of the spring constant can now be obtained by a probab ilistic sample average:" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 49 " evalf(add((k_i(ict)+0.5*dk)*C[ict],ict=1..20),3);" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#$\"$f)!\"\"" }}}{EXCHG {PARA 0 "" 0 "" {TEXT -1 100 "Ot her averages can be obtained similarly, e.g., the deviation squared of the individual data points:" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 55 "evalf(add((k_i(ict)+0.5*dk-k_a)^2*C[ict],ict=1..20),3);" }} {PARA 11 "" 1 "" {XPPMATH 20 "6#$\"$G$!\"#" }}}{EXCHG {PARA 0 "> " 0 " " {MPLTEXT 1 0 17 "evalf(sqrt(%),3);" }}{PARA 11 "" 1 "" {XPPMATH 20 " 6#$\"$\"=!\"#" }}}{EXCHG {PARA 0 "" 0 "" {TEXT -1 83 "This is slightly less than the value obtained before, as it does not replace the 1/" } {TEXT 266 1 "N" }{TEXT -1 64 " factor in the formation of the average \+ squared deviation by 1/(" }{TEXT 265 1 "N" }{TEXT -1 4 "-1)." }}{PARA 0 "" 0 "" {TEXT -1 0 "" }}{PARA 0 "" 0 "" {TEXT -1 58 "To make compari son with the limiting distribution we plot:" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 227 "P1:=plot([k_a,t,t=0..0.25],color=red,thickness=2) : P4:=plot([[k_a-sig_avg,t,t=0..0.25],[k_a+sig_avg,t,t=0..0.25]],color =magenta,thickness=2): P5:=plot([[k_a-sigma,t,t=0..0.15],[k_a+sigma,t, t=0..0.15]],color=yellow,thickness=2):" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 90 "P3:=plot(exp(-(k-k_a)^2/(2*sigma^2))/(sqrt(2*Pi)*sigm a),k=80..90,color=green,thickness=3):" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 109 "P2:=plot([seq([k_i(ict)+0.5*dk,t*C[ict],t=0..1],ict= 1..20)],color=blue,thickness=3): display(P1,P2,P3,P4,P5);" }}{PARA 13 "" 1 "" {GLPLOT2D 608 216 216 {PLOTDATA 2 "6>-%'CURVESG6%7S7$$\"3e++++ ++!f)!#;$\"\"!F,7$F($\"3WmmmT&)G\\a!#?7$F($\"3PL$ek`o!>5!#>7$F($\"3omm \"z>)G_:F47$F($\"3immT&QU!*3#F47$F($\"34L$eRZXKi#F47$F($\"3\\m;z>,_=JF 47$F($\"3h**\\7G$[8j$F47$F($\"3Gm;z%*frhTF47$F($\"3S**\\ilFQ!p%F47$F($ \"3mLL$3_\"=M_F47$F($\"3tmmTg(fJr&F47$F($\"3_****\\7eP_iF47$F($\"3%)** **\\Pf!Qz'F47$F($\"3%*****\\(=ubJ(F47$F($\"3em;zW(*Q*y(F47$F($\"3+LL$3 F-GN)F47$F($\"3_KLL$e'3I))F47$F($\"3Q**\\7.\"Feo7$F($\"3km;/Egw[7Feo7$F($\"3Ym\"z%*f%)QI\"Feo7$F($ \"3))*\\(oza'=N\"Feo7$F($\"3`m;zWho.9Feo7$F($\"3!*****\\i>Ad9Feo7$F($ \"3u**\\i:jf4:Feo7$F($\"3!)*\\7`>r-c\"Feo7$F($\"3y**\\P4q`;;Feo7$F($\" 3;LL$eM%4n;Feo7$F($\"3&****\\P4v5s\"Feo7$F($\"3cm\"zWn*)*pFeo7$F($\"3CLL 3-=!y(>Feo7$F($\"3*)*\\7G8O;.#Feo7$F($\"3smmm\"*\\[$3#Feo7$F($\"3gm;aQ z]O@Feo7$F($\"3ILekG=4*=#Feo7$F($\"3%*****\\i4TPAFeo7$F($\"3KL$3F9!z#H #Feo7$F($\"3emmmT>KUBFeo7$F($\"31+DJqJ8&R#Feo7$F($\"31+voa-oXCFeo7$F($ \"3++++++++DFeo-%'COLOURG6&%$RGBG$\"*++++\"!\")F+F+-%*THICKNESSG6#\"\" #-F$6%7S7$$\"3+++++++]!)F*F+FhuFhuFhuFhuFhuFhuFhuFhuFhuFhuFhuFhuFhuFhu FhuFhuFhuFhuFhuFhuFhuFhuFhuFhuFhuFhuFhuFhuFhuFhuFhuFhuFhuFhuFhuFhuFhuF huFhuFhuFhuFhuFhuFhuFhuFhuFhuFhu-F[u6&F]uF+F+F^u-Fbu6#\"\"$-F$6%7S7$$ \"#\")F,F+FcvFcvFcvFcvFcvFcvFcvFcvFcvFcvFcvFcvFcvFcvFcvFcvFcvFcvFcvFcv FcvFcvFcvFcvFcvFcvFcvFcvFcvFcvFcvFcvFcvFcvFcvFcvFcvFcvFcvFcvFcvFcvFcvF cvFcvFcvFcvFcvF[vF]v-F$6%7S7$$\"3+++++++]\")F*F+FivFivFivFivFivFivFivF ivFivFivFivFivFivFivFivFivFivFivFivFivFivFivFivFivFivFivFivFivFivFivFi vFivFivFivFivFivFivFivFivFivFivFivFivFivFivFivFivFivF[vF]v-F$6%7S7$$\" ##)F,F+F_wF_wF_wF_wF_wF_wF_wF_wF_wF_wF_wF_wF_wF_wF_wF_wF_wF_wF_wF_wF_w F_wF_wF_wF_wF_wF_wF_wF_wF_wF_wF_wF_wF_wF_wF_wF_wF_wF_wF_wF_wF_wF_wF_wF _wF_wF_wF_wF[vF]v-F$6%7S7$$\"3+++++++]#)F*F+FewFewFewFewFewFewFewFewFe wFewFewFewFewFewFewFewFewFewFewFewFewFewFewFewFewFewFewFewFewFewFewFew FewFewFewFewFewFewFewFewFewFewFewFewFewFewFewFewF[vF]v-F$6%7S7$$\"#$)F ,F+7$F\\x$\"3[mmm;arz@F07$F\\x$\"3mLL$e9ui2%F07$F\\x$\"3-nmm\"z_\"4iF0 7$F\\x$\"3[mmmT&phN)F07$F\\x$\"3KLLe*=)H\\5F47$F\\x$\"3omm\"z/3uC\"F47 $F\\x$\"3#****\\7LRDX\"F47$F\\x$\"3emm\"zR'ok;F47$F\\x$\"3'****\\i5`h( =F47$F\\x$\"3gLLL3En$4#F47$F\\x$\"3wmm;/RE&G#F47$F\\x$\"3))*****\\K]4] #F47$F\\x$\"3:+++vB_Y2aF47$F\\x$\"3qmm;zXu9cF47$F \\x$\"3@+++]y))GeF47$F\\x$\"35****\\i_QQgF47$F\\x$\"3j***\\7y%3TiF47$F \\x$\"3o****\\P![hY'F47$F\\x$\"3[LLL$Qx$omF47$F\\x$\"3k+++v.I%)oF47$F \\x$\"3Amm\"zpe*zqF47$F\\x$\"37+++D\\'QH(F47$F\\x$\"3GKLe9S8&\\(F47$F \\x$\"3]++D1#=bq(F47$F\\x$\"3>LLL3s?6zF47$F\\x$\"3P++DJXaE\")F47$F\\x$ \"3[nmmm*RRL)F47$F\\x$\"3Smm;a<.Y&)F47$F\\x$\"3-MLe9tOc()F47$F\\x$\"3u ******\\Qk\\*)F47$F\\x$\"3!QLL3dg6<*F47$F\\x$\"3UnmmmxGp$*F47$F\\x$\"3 !3+]7oK0e*F47$F\\x$\"3M,+v=5s#y*F47$F\\x$\"3/+++++++5FeoF[vF]v-F$6%7S7 $$\"3+++++++]$)F*F+FaalFaalFaalFaalFaalFaalFaalFaalFaalFaalFaalFaalFaa lFaalFaalFaalFaalFaalFaalFaalFaalFaalFaalFaalFaalFaalFaalFaalFaalFaalF aalFaalFaalFaalFaalFaalFaalFaalFaalFaalFaalFaalFaalFaalFaalFaalFaalFaa lF[vF]v-F$6%7S7$$\"#%)F,F+7$FhalF_x7$FhalFbx7$FhalFex7$FhalFhx7$FhalF[ y7$FhalF^y7$FhalFay7$FhalFdy7$FhalFgy7$FhalFjy7$FhalF]z7$FhalF`z7$Fhal Fcz7$FhalFfz7$FhalFiz7$FhalF\\[l7$FhalF_[l7$FhalFb[l7$FhalFe[l7$FhalFh [l7$FhalF[\\l7$FhalF^\\l7$FhalFa\\l7$FhalFd\\l7$FhalFg\\l7$FhalFj\\l7$ FhalF]]l7$FhalF`]l7$FhalFc]l7$FhalFf]l7$FhalFi]l7$FhalF\\^l7$FhalF_^l7 $FhalFb^l7$FhalFe^l7$FhalFh^l7$FhalF[_l7$FhalF^_l7$FhalFa_l7$FhalFd_l7 $FhalFg_l7$FhalFj_l7$FhalF]`l7$FhalF``l7$FhalFc`l7$FhalFf`l7$FhalFi`l7 $FhalF\\alF[vF]v-F$6%7S7$$\"3+++++++]%)F*F+F]elF]elF]elF]elF]elF]elF]e lF]elF]elF]elF]elF]elF]elF]elF]elF]elF]elF]elF]elF]elF]elF]elF]elF]elF ]elF]elF]elF]elF]elF]elF]elF]elF]elF]elF]elF]elF]elF]elF]elF]elF]elF]e lF]elF]elF]elF]elF]elF]elF[vF]v-F$6%7S7$$\"#&)F,F+7$Fdel$\"3*)******\\ i9RlF07$Fdel$\"3,++vVA)GA\"F47$Fdel$\"3%*****\\Peui=F47$Fdel$\"3%***** \\i3&o]#F47$Fdel$\"3e***\\(oX*y9$F47$Fdel$\"3l***\\P9CAu$F47$Fdel$\"3S ***\\P*zhdVF47$Fdel$\"3Q***\\P>fS*\\F47$Fdel$\"3'))**\\(=$f%GcF47$Fdel $\"37+++Dy,\"G'F47$Fdel$\"3E****\\7&)\\\"Feo7$Fdel$\"3k**\\P>:mk:Feo7$Fdel$\"3!* **\\iv&QAi\"Feo7$Fdel$\"3u***\\PPBWo\"Feo7$Fdel$\"3y*****\\Nm'[F eo7$Fdel$\"3i*****\\@80+#Feo7$Fdel$\"3y****\\7,Hl?Feo7$Fdel$\"3')**\\P 4w)R7#Feo7$Fdel$\"3!*****\\x%f\")=#Feo7$Fdel$\"3o**\\P/-a[AFeo7$Fdel$ \"3g**\\(=Yb;J#Feo7$Fdel$\"3#)****\\i@OtBFeo7$Fdel$\"3q**\\PfL'zV#Feo7 $Fdel$\"3(********)>=+DFeo7$Fdel$\"3#****\\i_4Qc#Feo7$Fdel$\"31+]P%>5p i#Feo7$Fdel$\"3#******\\:$*[o#Feo7$Fdel$\"3')***\\7<[8v#Feo7$Fdel$\"3n ******Hjy5GFeo7$Fdel$\"3'***\\P/)fT(GFeo7$Fdel$\"3%)**\\i0j\"[$HFeo7$F del$\"3))**************HFeoF[vF]v-F$6%7S7$$\"3+++++++]&)F*F+Fi^mFi^mFi ^mFi^mFi^mFi^mFi^mFi^mFi^mFi^mFi^mFi^mFi^mFi^mFi^mFi^mFi^mFi^mFi^mFi^m Fi^mFi^mFi^mFi^mFi^mFi^mFi^mFi^mFi^mFi^mFi^mFi^mFi^mFi^mFi^mFi^mFi^mFi ^mFi^mFi^mFi^mFi^mFi^mFi^mFi^mFi^mFi^mFi^mF[vF]v-F$6%7S7$$\"#')F,F+7$F `_m$\"3)HLLL$3VfVF07$F`_m$\"3Inmm\"H[D:)F07$F`_m$\"3TLLLe0$=C\"F47$F`_ m$\"3ILLL3RBr;F47$F`_m$\"3imm;zjf)4#F47$F`_m$\"3MLL$e4;[\\#F47$F`_m$\" 3$)****\\i'y]!HF47$F`_m$\"3:LL$ezs$HLF47$F`_m$\"3%*****\\7iI_PF47$F`_m $\"3Anmm;_M(=%F47$F`_m$\"3aLLL3y_qXF47$F`_m$\"3v******\\1!>+&F47$F`_m$ \"3H+++]Z/NaF47$F`_m$\"3'*******\\$fC&eF47$F`_m$\"3ELL$ez6:B'F47$F`_m$ \"37mmm;=C#o'F47$F`_m$\"3ulmmm#pS1(F47$F`_m$\"3y****\\i`A3vF47$F`_m$\" 3$ommmwy8!zF47$F`_m$\"3b,+]i.tK$)F47$F`_m$\"39++](3zMu)F47$F`_m$\"3?nm m\"H_?<*F47$F`_m$\"3emm;zihl&*F47$F`_m$\"3ULLL3#G,***F47$F`_m$\"3GLLez w5V5Feo7$F`_m$\"3)****\\PQ#\\\"3\"Feo7$F`_m$\"3MLL$e\"*[H7\"Feo7$F`_m$ \"3/+++qvxl6Feo7$F`_m$\"3#)****\\_qn27Feo7$F`_m$\"3#****\\i&p@[7Feo7$F `_m$\"3#*****\\2'HKH\"Feo7$F`_m$\"3pmmmwanL8Feo7$F`_m$\"37+++v+'oP\"Fe o7$F`_m$\"3CLLeR<*fT\"Feo7$F`_m$\"3-+++&)Hxe9Feo7$F`_m$\"3Ymm\"H!o-*\\ \"Feo7$F`_m$\"35++DTO5T:Feo7$F`_m$\"3kmmmT9C#e\"Feo7$F`_m$\"33++D1*3`i \"Feo7$F`_m$\"3]LLL$*zym;Feo7$F`_m$\"3GLL$3N1#4Feo7$F`_m$\"3G++v.Uac>Feo7$F`_m$\"35++++++ +?FeoF[vF]v-F$6%7S7$$\"3+++++++]')F*F+FehmFehmFehmFehmFehmFehmFehmFehm FehmFehmFehmFehmFehmFehmFehmFehmFehmFehmFehmFehmFehmFehmFehmFehmFehmFe hmFehmFehmFehmFehmFehmFehmFehmFehmFehmFehmFehmFehmFehmFehmFehmFehmFehm FehmFehmFehmFehmFehmF[vF]v-F$6%7S7$$\"#()F,F+F[imF[imF[imF[imF[imF[imF [imF[imF[imF[imF[imF[imF[imF[imF[imF[imF[imF[imF[imF[imF[imF[imF[imF[i mF[imF[imF[imF[imF[imF[imF[imF[imF[imF[imF[imF[imF[imF[imF[imF[imF[imF [imF[imF[imF[imF[imF[imF[imF[vF]v-F$6%7S7$$\"3+++++++]()F*F+FaimFaimFa imFaimFaimFaimFaimFaimFaimFaimFaimFaimFaimFaimFaimFaimFaimFaimFaimFaim FaimFaimFaimFaimFaimFaimFaimFaimFaimFaimFaimFaimFaimFaimFaimFaimFaimFa imFaimFaimFaimFaimFaimFaimFaimFaimFaimFaimF[vF]v-F$6%7S7$$\"#))F,F+7$F himFc_m7$FhimFf_m7$FhimFi_m7$FhimF\\`m7$FhimF_`m7$FhimFb`m7$FhimFe`m7$ FhimFh`m7$FhimF[am7$FhimF^am7$FhimFaam7$FhimFdam7$FhimFgam7$FhimFjam7$ FhimF]bm7$FhimF`bm7$FhimFcbm7$FhimFfbm7$FhimFibm7$FhimF\\cm7$FhimF_cm7 $FhimFbcm7$FhimFecm7$FhimFhcm7$FhimF[dm7$FhimF^dm7$FhimFadm7$FhimFddm7 $FhimFgdm7$FhimFjdm7$FhimF]em7$FhimF`em7$FhimFcem7$FhimFfem7$FhimFiem7 $FhimF\\fm7$FhimF_fm7$FhimFbfm7$FhimFefm7$FhimFhfm7$FhimF[gm7$FhimF^gm 7$FhimFagm7$FhimFdgm7$FhimFggm7$FhimFjgm7$FhimF]hm7$FhimF`hmF[vF]v-F$6 %7S7$$\"3+++++++]))F*F+F]]nF]]nF]]nF]]nF]]nF]]nF]]nF]]nF]]nF]]nF]]nF]] nF]]nF]]nF]]nF]]nF]]nF]]nF]]nF]]nF]]nF]]nF]]nF]]nF]]nF]]nF]]nF]]nF]]nF ]]nF]]nF]]nF]]nF]]nF]]nF]]nF]]nF]]nF]]nF]]nF]]nF]]nF]]nF]]nF]]nF]]nF]] nF]]nF[vF]v-F$6%7S7$$\"#*)F,F+7$Fd]nF_x7$Fd]nFbx7$Fd]nFex7$Fd]nFhx7$Fd ]nF[y7$Fd]nF^y7$Fd]nFay7$Fd]nFdy7$Fd]nFgy7$Fd]nFjy7$Fd]nF]z7$Fd]nF`z7$ Fd]nFcz7$Fd]nFfz7$Fd]nFiz7$Fd]nF\\[l7$Fd]nF_[l7$Fd]nFb[l7$Fd]nFe[l7$Fd ]nFh[l7$Fd]nF[\\l7$Fd]nF^\\l7$Fd]nFa\\l7$Fd]nFd\\l7$Fd]nFg\\l7$Fd]nFj \\l7$Fd]nF]]l7$Fd]nF`]l7$Fd]nFc]l7$Fd]nFf]l7$Fd]nFi]l7$Fd]nF\\^l7$Fd]n F_^l7$Fd]nFb^l7$Fd]nFe^l7$Fd]nFh^l7$Fd]nF[_l7$Fd]nF^_l7$Fd]nFa_l7$Fd]n Fd_l7$Fd]nFg_l7$Fd]nFj_l7$Fd]nF]`l7$Fd]nF``l7$Fd]nFc`l7$Fd]nFf`l7$Fd]n Fi`l7$Fd]nF\\alF[vF]v-F$6%7S7$$\"3+++++++]*)F*F+Fi`nFi`nFi`nFi`nFi`nFi `nFi`nFi`nFi`nFi`nFi`nFi`nFi`nFi`nFi`nFi`nFi`nFi`nFi`nFi`nFi`nFi`nFi`n Fi`nFi`nFi`nFi`nFi`nFi`nFi`nFi`nFi`nFi`nFi`nFi`nFi`nFi`nFi`nFi`nFi`nFi `nFi`nFi`nFi`nFi`nFi`nFi`nFi`nF[vF]v-F$6%7S7$$\"#!*F,F+F_anF_anF_anF_a nF_anF_anF_anF_anF_anF_anF_anF_anF_anF_anF_anF_anF_anF_anF_anF_anF_anF _anF_anF_anF_anF_anF_anF_anF_anF_anF_anF_anF_anF_anF_anF_anF_anF_anF_a nF_anF_anF_anF_anF_anF_anF_anF_anF_anF[vF]v-F$6%7W7$$\"#!)F,$\"3'e()z \"[/dpsX5tf&=sF47$$\"3wLL$347TL)F*$ \"3Xv>#zS7Feo7$$\"32+]7=lj;%)F*$\"3 )Q\\WH,-NQ\"Feo7$$\"3g**\\PaR\"HgFeo7$$\"3?+](= >Y2a)F*$\"3Gc35@pR??Feo7$$\"3_Le*[K56b)F*$\"3g=sf2E&e/#Feo7$$\"3%om;zX u9c)F*$\"3kw+3_#Qb1#Feo7$$\"3wK$e9i\"=s&)F*$\"3gJUo!yL'z?Feo7$$\"35+++ &y))Ge)F*$\"3+tr%o(eD(3#Feo7$$\"3y+]ibOO$f)F*$\"3Q7%z'p\"z$)3#Feo7$$\" 3.++DE&QQg)F*$\"31o\\wmzA$3#Feo7$$\"3y*\\(=-N(Rh)F*$\"3kdAV^ZJs?Feo7$$ \"3_**\\7y%3Ti)F*$\"3mmXDBDmb?Feo7$$\"3a++v.[hY')F*$\"3S`(p)Q@$*)*>Feo 7$$\"3iLLLQx$om)F*$\"3kcs'z5Wj#>Feo7$$\"3_****\\P+V)o)F*$\"3y\\ERWW(*G =Feo7$$\"3im;zpe*zq)F*$\"3\"Feo7$$\"3\\+]7`Wl7))F*$ \"35N@Ck/ue5Feo7$$\"3Bmmm'*RRL))F*$\"3-J2Ob0)H&F47$F`io$\"3\\**\\(=-p6j&F47$F`io$\"3 u)****\\2Mg#fF47$F`io$\"3y**\\(=xZ&\\iF47$F`io$\"35+]i:$4wb'F47$F`io$ \"3O***\\(=#R!zoF47$F`io$\"3%***\\P4A@urF47$F`io$\"3-***\\i:'f#\\(F47$ F`io$\"3h**\\(of2L#yF47$F`io$\"3^**\\7yG>6\")F47$F`io$\"3/++voo6A%)F47 $F`io$\"3#*)****\\xJLu)F47$F`io$\"3+***\\P*ydd!*F47$F`io$\"35**\\(=1\"Feo7$F`io$\"3%****\\(Q(zS4\"Feo7$F `io$\"3%)*\\(=-,FC6Feo7$F`io$\"3%**\\P4tFe:\"Feo7$F`io$\"3w***\\73\"o' =\"Feo7$F`io$\"3&)*\\(oz;)*=7Feo7$F`io$\"3)******\\*44]7Feo7$F`io$\"3% ***\\7jZ!>G\"Feo7$F`io$\"3.+v=(4bMJ\"Feo7$F`io$\"3'*****\\xlWU8Feo7$F` io$\"3%***\\i&3ucP\"Feo7$F`io$\"3$)*****\\;$R09Feo7$F`io$\"3)**\\(=-*z qV\"Feo7$F`io$\"3#**\\7G:3uY\"Feo7$F`io$\"3%**************\\\"Feo-F[u6 &F]uF^uF^uF+Fau-F$6%7S7$$\"3C++++++\"y)F*F+7$FhbpFcio7$FhbpFfio7$FhbpF iio7$FhbpF\\jo7$FhbpF_jo7$FhbpFbjo7$FhbpFejo7$FhbpFhjo7$FhbpF[[p7$Fhbp F^[p7$FhbpFa[p7$FhbpFd[p7$FhbpFg[p7$FhbpFj[p7$FhbpF]\\p7$FhbpF`\\p7$Fh bpFc\\p7$FhbpFf\\p7$FhbpFi\\p7$FhbpF\\]p7$FhbpF_]p7$FhbpFb]p7$FhbpFe]p 7$FhbpFh]p7$FhbpF[^p7$FhbpF^^p7$FhbpFa^p7$FhbpFd^p7$FhbpFg^p7$FhbpFj^p 7$FhbpF]_p7$FhbpF`_p7$FhbpFc_p7$FhbpFf_p7$FhbpFi_p7$FhbpF\\`p7$FhbpF_` p7$FhbpFb`p7$FhbpFe`p7$FhbpFh`p7$FhbpF[ap7$FhbpF^ap7$FhbpFaap7$FhbpFda p7$FhbpFgap7$FhbpFjap7$FhbpF]bp7$FhbpF`bpFbbpFau-%+AXESLABELSG6%Q!6\"F ]fp-%%FONTG6#%(DEFAULTG-%%VIEWG6$;FfanF`anFbfp" 1 2 0 1 10 0 2 9 1 4 2 1.000000 45.000000 45.000000 0 0 "Curve 1" "Curve 2" "Curve 3" "Curv e 4" "Curve 5" "Curve 6" "Curve 7" "Curve 8" "Curve 9" "Curve 10" "Cur ve 11" "Curve 12" "Curve 13" "Curve 14" "Curve 15" "Curve 16" "Curve 1 7" "Curve 18" "Curve 19" "Curve 20" "Curve 21" "Curve 22" "Curve 23" " Curve 24" "Curve 25" "Curve 26" }}}}{EXCHG {PARA 0 "" 0 "" {TEXT -1 257 "We have added a few vertical lines to the comparison: the magenta -coloured lines indicate the standard deviation of the mean; the yello w coloured line indicate the standard deviation for the Gaussian distr ibution associated with the individual measurements." }}{PARA 0 "" 0 " " {TEXT -1 0 "" }}{PARA 0 "" 0 "" {TEXT -1 343 "It is the obligation o f the experimenter to demonstrate that the error in the data are indee d normally distributed, i.e., that the agreement with a Gaussian distr ibution improves for larger data sets. Once this is demonstrated, the \+ probability theory based on the normal distribution can be used to fur ther assess the confidence in our result." }}{PARA 0 "" 0 "" {TEXT -1 0 "" }}{PARA 0 "" 0 "" {TEXT -1 189 "First we demonstrate that the con tinuous Gaussian distribution results in the expected values for the t otal probability (area under the curve), the average value, and the de viation squared." }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 74 "int(exp (-(k-k_a)^2/(2*sigma^2))/(sqrt(2*Pi)*sigma),k=-infinity..infinity);" } }{PARA 11 "" 1 "" {XPPMATH 20 "6#$\"+++++5!\"*" }}}{EXCHG {PARA 0 "> \+ " 0 "" {MPLTEXT 1 0 76 "int(k*exp(-(k-k_a)^2/(2*sigma^2))/(sqrt(2*Pi)* sigma),k=-infinity..infinity);" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#$\"+ ,++!f)!\")" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 84 "int((k-k_a)^2 *exp(-(k-k_a)^2/(2*sigma^2))/(sqrt(2*Pi)*sigma),k=-infinity..infinity) ;" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#$\"+-+5[O!\"*" }}}{EXCHG {PARA 0 "" 0 "" {TEXT -1 85 "One can ask how much probability is contained wit hin the standard deviation interval:" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 75 "int(exp(-(k-k_a)^2/(2*sigma^2))/(sqrt(2*Pi)*sigma),k= k_a-sigma..k_a+sigma);" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#$\"+@\\*o#o! #5" }}}{EXCHG {PARA 0 "" 0 "" {TEXT -1 334 "We say that we have 68% co nfidence that a measurement will fall within this interval. Why do we \+ care about this statement (given that we have determined an accurate ' best' value by forming the average)? The answer is that what we learn \+ about individual measurements can be transferred later to assess the a ccuracy of the 'best' value." }}{PARA 0 "" 0 "" {TEXT -1 219 "So let u s continue with considering the individual measurements. A 68% confide nce sounds not so terribly high. We can carry out area calculations, i .e., confidence estimates for larger than one-sigma intervals. For the " }{XPPEDIT 18 0 "2*sigma;" "6#*&\"\"#\"\"\"%&sigmaGF%" }{TEXT -1 6 " - and " }{XPPEDIT 18 0 "3*sigma;" "6#*&\"\"$\"\"\"%&sigmaGF%" }{TEXT -1 19 "-intervals we have:" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 79 "int(exp(-(k-k_a)^2/(2*sigma^2))/(sqrt(2*Pi)*sigma),k=k_a-2*sigma.. k_a+2*sigma);" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#$\"+it*\\a*!#5" }}} {EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 79 "int(exp(-(k-k_a)^2/(2*sigma^ 2))/(sqrt(2*Pi)*sigma),k=k_a-3*sigma..k_a+3*sigma);" }}{PARA 11 "" 1 " " {XPPMATH 20 "6#$\"+T?+t**!#5" }}}{EXCHG {PARA 0 "" 0 "" {TEXT -1 217 "Therefore, once we have determined sigma for our distribution, we can say with a definite confidence that subsequent measurements (by o urselves, or someone else provided our systematics are identical) will fall into (" }{XPPEDIT 18 0 "t*sigma;" "6#*&%\"tG\"\"\"%&sigmaGF%" } {TEXT -1 17 ")-intervals, for " }{TEXT 267 1 "t" }{TEXT -1 9 "=1,2,... ." }}{PARA 0 "" 0 "" {TEXT -1 0 "" }}{PARA 0 "" 0 "" {TEXT -1 533 "Wha t relevance does this have for assessing our best value? One could per form a statistical analysis of the average value, i.e., run 10 dataset s with 10 measurements respectively, compute the averages, and then th e average of averages with known uncertainty. The previous assessment \+ of the standard deviation of the mean is a shortcut (but it is a guess as opposed to a measurement). It guesses that these measurements of t he averages ('best' values) will be distributed randomly according to \+ a Gaussian with a width estimated to be " }{TEXT 19 7 "sig_avg" } {TEXT -1 207 ". Therefore, we can consider the interval surrounding th e best estimate (red vertical line with magenta-coloured sidebars) to \+ be the 68% confidence interval for the true value. We can draw the int erval with " }{TEXT 19 9 "2*sig_avg" }{TEXT -1 218 " and be 95% confid ent that the true value lies within this interval. Taylor provides a p roof about the Gaussian distribution of the best values based on Gauss ian distributed individual data measurements in section 5.7." }}{PARA 0 "" 0 "" {TEXT -1 0 "" }}{PARA 0 "" 0 "" {TEXT -1 35 "The lessons to \+ be learned here are:" }}{PARA 0 "" 0 "" {TEXT -1 365 "a) the one-sigma interval is not a holy grail. it contains only a 68% confidence asses sment! One-sigma 'new' physics (propagated with big hoopla in Physical Review Letters or Physics Letters and suspected by the critical minds ) often goes away! For three-sigma events the chances are much less th at they are the result of hype rather than solid experimentation. By \+ " }{TEXT 270 1 "t" }{TEXT -1 72 "-sigma 'new' physics we mean unexpect ed results that lie outside of the " }{TEXT 271 1 "t" }{TEXT -1 51 "-s igma confidence interval of conventional physics." }}{PARA 0 "" 0 "" {TEXT -1 316 "b) the reader may wish to understand why the average rep resents the best estimate for the measured value. The argument involve s turning things around, i.e., assuming a limiting Gaussian distributi on with yet undetermined parameters to calculate the probability to ar rive at a certain measured data set (a product of " }{TEXT 268 1 "N" } {TEXT -1 298 " probabilities). A so-called principle of maximum likeli hood is invoked which results in a determination of the peak of the Ga ussian at the position given by the average of the measured data value s (the arguments are given in Taylor's book, and elaborated on in appl ied probability theory courses)." }}{PARA 0 "" 0 "" {TEXT -1 457 "c) t here are probably lessons to be learned also for error propagation. Gi ven that we are arriving at an understanding of uncertainties for one \+ variable, we can learn that normally distributed errors in two (or mor e) variables conspire to form errors in some final result (that may in volve some mathematical operation between the two basic variables). Ta ylor explains how addition of errors in quadrature follows from the no rmal distribution in section 5.6." }}{PARA 0 "" 0 "" {TEXT -1 0 "" }}} {EXCHG {PARA 0 "" 0 "" {TEXT 269 11 "Simulations" }}{PARA 0 "" 0 "" {TEXT -1 223 "We can verify our claims by generating virtual measureme nts in Maple, i.e., computer simulations of normally distributed measu rements. Then we apply the formulae and check whether our findings are consistent with the input." }}{PARA 0 "" 0 "" {TEXT -1 0 "" }}{PARA 0 "" 0 "" {TEXT -1 225 "The first task is to find a distribution of ra ndom numbers which follows a Gaussian. Normally pseudo-random number g enerators produce randoms that fill the [0,1] interval uniformly. Norm ally distributed randoms fill the range" }{TEXT 19 19 "-infinity..infi nity" }{TEXT -1 63 " in such a way that a histogram produces a Gaussia n of width 1." }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 12 "with(stats ):" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 35 "RNd:=[stats[random, n ormald](150)]:" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 23 "with(stat s[statplots]):" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 22 "histogram (RNd,area=1);" }}{PARA 13 "" 1 "" {GLPLOT2D 581 297 297 {PLOTDATA 2 "6 %-%)POLYGONSG6-7&7$$!+k$y_e#!\"*\"\"!7$F($\"+LLLL8!#67$$!++)fh:#F*F-7$ F1F+7&F37$F1$\"+)*******RF/7$$!+O7/F7$FAF+7&FC7$FA$\"+gmmm')F/7$$!*4T!)o)F*FF7$FIF+7&FK7$FI$ \"+mmmm;!#57$$!*XboR%F*FN7$FRF+7&FTFQ7$$!)\")pc5F*FN7$FWF+7&FY7$FW$\"+ LLLL6FP7$$\"*$e^&=%F*Ffn7$FinF+7&F[o7$Fin$\"+pmmm9FP7$$\"*Y,nZ)F*F^o7$ FaoF+7&Fco7$Fao$\"+&*******fF/7$$\"+5()yw7F*Ffo7$FioF+7&F[p7$Fio$\"+IL LL`F/7$$\"+us!fq\"F*F^p7$FapF+7&Fcp7$Fap$\"+jmmmmF/7$$\"+Qe-N@F*Ffp7$F ipF+-%&COLORG6&%$RGBG$\"\"(!\"\"F`q$\"#5Fbq-%%VIEWG6$;$!+%y3t0$F*$\"+e i02EF*%(DEFAULTG" 1 2 0 1 10 0 2 6 1 4 2 1.000000 45.000000 45.000000 0 0 "Curve 1" }}}}{EXCHG {PARA 0 "" 0 "" {TEXT -1 115 "Without the opt ion area=1 the histogram would be arranged in such a way that each bar would contains the same area." }}{PARA 0 "" 0 "" {TEXT -1 0 "" }} {PARA 0 "" 0 "" {TEXT -1 314 "We can use the numbers produced with a d eviation of 1 to generate errors by simply scaling the numbers with a \+ constant factor. The average comes out to be zero, and thus we can add the scaled numbers to some hypothetical measurement value. Suppose we have measured an angle of val=38 degrees with some uncertainty." }}} {EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 8 "val:=38;" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%$valG\"#Q" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 13 "N:=nops(RNd);" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%\"NG\"$]\"" }}} {EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 36 "MData:=[seq(val+0.5*RNd[i],i =1..N)]:" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 24 "histogram(MData ,area=1);" }}{PARA 13 "" 1 "" {GLPLOT2D 616 256 256 {PLOTDATA 2 "6%-%) POLYGONSG6-7&7$$\"+#3O2n$!\")\"\"!7$F($\"+XLLL8!#67$$\"+5?>#p$F*F-7$F1 F+7&F37$F1$\"+M+++SF/7$$\"+Qzk8PF*F67$F9F+7&F;7$F9$\"+Pjmm')F/7$$\"+nQ 5NPF*F>7$FAF+7&FC7$FA$\"+Tnmm')F/7$$\"+&zflv$F*FF7$FIF+7&FK7$FI$\"+\"o mmm\"!#57$$\"+Bd,yPF*FN7$FRF+7&FTFQ7$$\"+^;Z*z$F*FN7$FWF+7&FY7$FW$\"+V LLL6FP7$$\"+zv#4#QF*Ffn7$FinF+7&F[o7$Fin$\"+zmmm9FP7$$\"+2NQUQF*F^o7$F aoF+7&Fco7$Fao$\"+s(*****fF/7$$\"+O%RQ'QF*Ffo7$FioF+7&F[p7$Fio$\"+!QLL L&F/7$$\"+k`H&)QF*F^p7$FapF+7&Fcp7$Fap$\"+BnmmmF/7$$\"+#H^n!RF*Ffp7$Fi pF+-%&COLORG6&%$RGBG$\"\"(!\"\"F`q$\"#5Fbq-%%VIEWG6$;$\"+hX8ZOF*$\"+8G NIRF*%(DEFAULTG" 1 2 0 1 10 0 2 6 1 4 2 1.000000 45.000000 45.000000 0 0 "Curve 1" }}}}{EXCHG {PARA 0 "" 0 "" {TEXT -1 151 "Let us apply ou r expressions to recover the true value. In the example below it is im portant to carry more than 3 digits when carrying out the average!" }} }{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 37 "v_a:=evalf(add(MData[i],i=1 ..N)/N,7);" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%$v_aG$\"(Hjz$!\"&" }}} {EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 57 "sigma:=evalf(sqrt(add((MData [n]-v_a)^2,n=1..N)/(N-1)),5);" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%&si gmaG$\"&N:&!\"&" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 25 "sig_avg: ='sigma/sqrt(N)';" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%(sig_avgG*&%&si gmaG\"\"\"-%%sqrtG6#%\"NG!\"\"" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 24 "sig_avg:=evalf(sig_avg);" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>% (sig_avgG$\"+H^\"y?%!#6" }}}{EXCHG {PARA 0 "" 0 "" {TEXT -1 67 "We hav e a 68% confidence interval which can exclude the true value." }}} {EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 44 "[evalf(v_a-sig_avg,4),evalf( v_a+sig_avg,4)];" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#7$$\"%#z$!\"#$\"%+ QF&" }}}{EXCHG {PARA 0 "" 0 "" {TEXT -1 65 "For the 95% CI we almost c ertainly include the correct value for " }{TEXT 272 1 "N" }{TEXT -1 5 ">100:" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 48 "[evalf(v_a-2*sig_ avg,4),evalf(v_a+2*sig_avg,4)];" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#7$$ \"%)y$!\"#$\"%/QF&" }}}{EXCHG {PARA 0 "" 0 "" {TEXT -1 211 "If we don' t include the correct value in the two-sigma interval, the cause for t his discrepancy comes from the fact that the average of the sample tak en according to a normal distribution has a systematic drift:" }}} {EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 28 "evalf(add(RNd[i],i=1..N)/N); " }}{PARA 11 "" 1 "" {XPPMATH 20 "6#$!+$HehM(!#6" }}}{EXCHG {PARA 0 " " 0 "" {TEXT 274 11 "Exercise 1:" }}{PARA 0 "" 0 "" {TEXT -1 140 "Expl ore the question whether the true value lies within the 1-sigma and wi thin the 2-sigma confidence intervals around the 'best' answer by:" }} {PARA 0 "" 0 "" {TEXT -1 52 "a) changing the sample size (increase and decrease);" }}{PARA 0 "" 0 "" {TEXT -1 114 "b) changing the true valu e of the measured quantity (increase and decrease), while the error si ze is kept the same" }}{PARA 0 "" 0 "" {TEXT -1 91 "c) changing the er ror size, while keeping the original true value of the measured quanti ty." }}{PARA 0 "" 0 "" {TEXT -1 0 "" }}{PARA 0 "" 0 "" {TEXT -1 244 "K eep in mind that the simulation in Maple provides us with the luxury o f many repeated measurements. In the undergraduate laboratory setting \+ we often forego this luxury (unless experiments are automated), and pa y the price of reduced precision." }}{PARA 0 "" 0 "" {TEXT -1 0 "" }}} {EXCHG {PARA 0 "" 0 "" {TEXT -1 264 "We conclude the worksheet with an example of error propagation. Let us assume that the quantity measure d above represents an angle of refraction in a simple experiment with \+ a lightbeam at an air/glass interface. Let us assume that we know the \+ index of refraction (" }{TEXT 273 1 "n" }{TEXT -1 167 "=1.51 for glass ) to a higher degree of precision than the angle measurement. What is \+ the deviation in the incident angle. We first look at the transformed \+ data sample:" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 49 "fun:=x->eva lf(180/Pi*arcsin(sin(Pi/180*x)*1.51));" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%$funGj+6#%\"xG6\"6$%)operatorG%&arrowGF(-%&evalfG6#,$*(\"$!=\" \"\"%#PiG!\"\"-%'arcsinG6#*&-%$sinG6#,$*&#F2F1F2*&F3F29$F2F2F2F2$\"$^ \"!\"#F2F2F2F(F(F(6#\"\"!" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 23 "IAngle:=map(fun,MData):" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 10 "IAngle[1];" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#$\"+/pwNq!\")" }}} {EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 25 "histogram(IAngle,area=1);" } }{PARA 13 "" 1 "" {GLPLOT2D 579 243 243 {PLOTDATA 2 "6%-%)POLYGONSG6-7 &7$$\"+a**y\\k!\")\"\"!7$F($\"+LLLL8!#67$$\"+uT+>lF*F-7$F1F+7&F37$F1$ \"+mmmmYF/7$$\"+%R=#)e'F*F67$F9F+7&F;7$F9$\"+++++5!#57$$\"+9EVdmF*F>7$ FBF+7&FD7$FB$\"+mmmm')F/7$$\"+MokEnF*FG7$FJF+7&FL7$FJ$\"+++++?F@7$$\"+ a5'ez'F*FO7$FRF+7&FT7$FR$\"+mmmm;F@7$$\"+u_2loF*FW7$FZF+7&Ffn7$FZ$\"++ +++7F@7$$\"+%\\*GMpF*Fin7$F\\oF+7&F^o7$F\\o$\"+LLLL6F@7$$\"+9P].qF*Fao 7$FdoF+7&Ffo7$FdoF67$$\"+MzrsqF*F67$FjoF+7&F\\p7$Fjo$\"+KLLL`F/7$$\"+a @$>9(F*F_p7$FbpF+7&FdpFap7$$\"+uj96sF*F_p7$FgpF+-%&COLORG6&%$RGBG$\"\" (!\"\"F^q$\"#5F`q-%%VIEWG6$;$\"+7VltjF*$\"+;?G(G(F*%(DEFAULTG" 1 2 0 1 10 0 2 6 1 4 2 1.000000 46.000000 45.000000 0 0 "Curve 1" }}}} {EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 39 "IA_a:=evalf(add(IAngle[i],i= 1..N)/N,7);" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%%IA_aG$\"(8<$o!\"&" } }}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 61 "sigmaIA:=evalf(sqrt(add((I Angle[n]-IA_a)^2,n=1..N)/(N-1)),5);" }}{PARA 11 "" 1 "" {XPPMATH 20 "6 #>%(sigmaIAG$\"&]n\"!\"%" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 29 "sig_avgIA:='sigmaIA/sqrt(N)';" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%*s ig_avgIAG*&%(sigmaIAG\"\"\"-%%sqrtG6#%\"NG!\"\"" }}}{EXCHG {PARA 0 "> \+ " 0 "" {MPLTEXT 1 0 28 "sig_avgIA:=evalf(sig_avgIA);" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%*sig_avgIAG$\"+t " 0 "" {MPLTEXT 1 0 27 "IA_true:=evalf(fun( val),5);" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%(IA_trueG$\"&#Qo!\"$" }} }{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 50 "[evalf(IA_a-sig_avgIA,4),ev alf(IA_a+sig_avgIA,4)];" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#7$$\"%=o!\" #$\"%YoF&" }}}{EXCHG {PARA 0 "" 0 "" {TEXT -1 64 "We notice that the d eviation tripled. How can this be explained?" }}{PARA 0 "" 0 "" {TEXT -1 0 "" }}{PARA 0 "" 0 "" {TEXT -1 105 "We made the claim that the unc ertainty in the index of refraction is negligible. Thus we use Snell's law " }{TEXT 19 17 "sin(IA)=n*sin(RA)" }{TEXT -1 52 " to relate the u ncertainties in the reflected angle " }{TEXT 19 2 "RA" }{TEXT -1 42 " \+ to the uncertainty in the incident angle " }{TEXT 19 2 "IA" }{TEXT -1 57 ". Considering the differentials results in the statement:" }} {PARA 0 "" 0 "" {TEXT -1 5 "LHS: " }{TEXT 19 24 "d[sin(IA)]=cos(IA)*d[ IA]" }}{PARA 0 "" 0 "" {TEXT -1 5 "RHS: " }{TEXT 19 28 "n*d[sin(RA)]=n *cos(RA)*d[RA]" }}{PARA 0 "" 0 "" {TEXT -1 30 "We can solve the equati on for " }{TEXT 19 5 "d[IA]" }{TEXT -1 2 ": " }}{PARA 0 "" 0 "" {TEXT 19 33 "d[IA] = n*(cos(RA)/cos(IA))*d[RA]" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 51 "evalf(1.51*cos(val*Pi/180)/cos(Pi/180*IA_a)*sigma);" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#$\"+FYrf;!\"*" }}}{EXCHG {PARA 0 "" 0 "" {TEXT -1 26 "This calculated value for " }{TEXT 19 7 "sigmaIA" } {TEXT -1 111 " is very close indeed to the one measured above. Thus, i t is important to consider error propagation carefully!" }}{PARA 0 "" 0 "" {TEXT -1 0 "" }}{PARA 0 "" 0 "" {TEXT -1 282 "We can illustrate t he error propagation for this nonlinear function by a graph. This also serves to justify why a linearized treatment of the function in the n eighbourhood of the mean argument value is jsutified (the linearizatio n is implied by the use of first-oder differentials)." }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 40 "f_ex:=180/Pi*arcsin(1.51*sin(Pi/180 *x));" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%%f_exG,$*(\"$!=\"\"\"%#PiG! \"\"-%'arcsinG6#,$*&$\"$^\"!\"#F(-%$sinG6#,$*(F'F*F)F(%\"xGF(F(F(F(F(F (" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 44 "f_lin:=convert(taylor( f_ex,x=38,2),polynom);" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%&f_linG,&$ \"+!fdSV&!\")!\"\"*&$\"+;D\\HK!\"*\"\"\"%\"xGF.F." }}}{EXCHG {PARA 0 " > " 0 "" {MPLTEXT 1 0 6 "sigma;" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#$\" &N:&!\"&" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 49 "P1:=plot([f_ex, f_lin],x=37..39,color=[red,blue]):" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 184 "P2:=plot([[val-sigma,65],[val-sigma,subs(x=val-sigma ,f_lin)],[0,subs(x=val-sigma,f_lin)]],color=green): P3:=plot([[val,65] ,[val,subs(x=val,f_lin)],[0,subs(x=val,f_lin)]],color=yellow):" }} {PARA 0 "> " 0 "" {MPLTEXT 1 0 103 "P4:=plot([[val+sigma,65],[val+sigm a,subs(x=val+sigma,f_lin)],[0,subs(x=val+sigma,f_lin)]],color=brown): " }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 77 "display(P1,P2,P3,P4,sca ling=constrained,title=\"incident vs refracted angle\");" }}{PARA 13 " " 1 "" {GLPLOT2D 495 495 495 {PLOTDATA 2 "6+-%'CURVESG6$7S7$$\"#P\"\"! $\"3<&y:H&**=Ll!#;7$$\"3GLL$3VfVq$F-$\"3-%ej`$H\"ea'F-7$$\"3wm;H[D:3PF -$\"3kVX%)[*Rob'F-7$$\"3[LLe0$=Cr$F-$\"3g9V.B/HplF-7$$\"3@LL3RBr;PF-$ \"3%[T%e$3x=e'F-7$$\"3^m;zjf)4s$F-$\"3%p\">+6#eWf'F-7$$\"3VL$e4;[\\s$F -$\"3(\\STylrhg'F-7$$\"3K+]i'y]!HPF-$\"3.*fw]$3N=mF-7$$\"3PL$ezs$HLPF- $\"38;TZF:+JmF-7$$\"3u**\\7iI_PPF-$\"362!=@;oOk'F-7$$\"3Vmm;_M(=u$F-$ \"3g@**G)=dnl'F-7$$\"3oLL3y_qXPF-$\"3%RB8\"yqLomF-7$$\"3#)****\\1!>+v$ F-$\"3e%Q>%f9V\"o'F-7$$\"3!*****\\Z/NaPF-$\"3_$Rz!=Ek%p'F-7$$\"3K++]$f C&ePF-$\"3w_$\\3yMuq'F-7$$\"39L$ez6:Bw$F-$\"3KV.q:Q5>nF-7$$\"3imm;=C#o w$F-$\"3=P()=-c/LnF-7$$\"3kmmm#pS1x$F-$\"3x(=>yI8\\u'F-7$$\"3')**\\i`A 3vPF-$\"3s!zm\\V¥F-7$$\"3'pmmwy8!zPF-$\"3-$R&*o\"f7rnF-7$$\"3))**\\ i.tK$y$F-$\"3E^8g+Gt%y'F-7$$\"3O+](3zMuy$F-$\"3M\"*pT.lv(z'F-7$$\"3)om ;H_?53z^pPoF-7$$\"3jL$ezw5V!QF-$\"3%Q3#G%Ru>&oF-7$$\"3K+]PQ# \\\"3QF-$\"3_P:&[YrW'oF-7$$\"3PLLe\"*[H7QF-$\"31A\\([**R!yoF-7$$\"3J++ +dxd;QF-$\"3MC(f]\"p8#*oF-7$$\"3H++D0xw?QF-$\"3Ryo8Rv+1pF-7$$\"3x**\\i &p@[#QF-$\"3cpF-7$$\"3D++vgHKHQF-$\"3g@hC!H!eMpF-7$$\"3,nmmZvOL QF-$\"3=?pwpF-7$$\"3K++])Hxe%QF-$\"33Oi*['>(3*pF-7$$\"3sm;H!o-*\\QF -$\"33y&=`!*pZ+(F-7$$\"3o**\\7k.6aQF-$\"3#>rm)Q))Q>qF-7$$\"3wmm;WTAeQF -$\"3(GkJX*RxLqF-7$$\"3#)**\\i!*3`iQF-$\"3;MM/1T$*[qF-7$$\"3oLLL*zym'Q F-$\"3=K%)[%>MO1(F-7$$\"3SLL3N1#4(QF-$\"3-8'4:3q(yqF-7$$\"3Wm;HYt7vQF- $\"3!GPo!\\e)Q4(F-7$$\"3K+++xG**yQF-$\"3#*))f)e?qy5(F-7$$\"3ummT6KU$)Q F-$\"35U>z&*=,CrF-7$$\"3/LLLbdQ()QF-$\"3FEY#\\)[bQrF-7$$\"3#)**\\i`1h \"*QF-$\"3M()fT7M3fLa&=(F--%'COLOURG6&%$RGBG$\"*++++\"!\")$F*F*F_[l-F$6$7S7$F($\"3 \"3++#>Z1:lF-7$F/$\"3**o-XoM9HlF-7$F4$\"3%fE%*GJ$RTlF-7$F9$\"3ywF9W&p^ b'F-7$F>$\"3UyW>$4P!plF-7$FC$\"3[Ni?K(QGe'F-7$FH$\"3C,HI=Yj&f'F-7$FM$ \"3'3N*= ZH]mF-7$F[o$\"31E=rs&pEm'F-7$F`o$\"3_0&\\js+mn'F-7$Feo$\"3y+Dk#3*e!p'F -7$Fjo$\"32&>5fXpSq'F-7$F_p$\"3$e4!GE4J;nF-7$Fdp$\"3Gt-d=s'3t'F-7$Fip$ \"3lc?\\2$)>VnF-7$F^q$\"3%>Nm`HUvv'F-7$Fcq$\"37:.8b\"R-x'F-7$Fhq$\"3Y( RNEipTy'F-7$F]r$\"3IQ)o>sMuz'F-7$Fbr$\"3!=)H_haF6oF-7$Fgr$\"3Xl\\/\"e& )R#oF-7$F\\s$\"3SKp-i^pPoF-7$Fas$\"3_^'fOcN>&oF-7$Ffs$\"3-lJ9v=LkoF-7$ F[t$\"3)oS`6B?x(oF-7$F`t$\"3UnN'Hr^:*oF-7$Fet$\"3fr$z%>J30pF-7$Fjt$\"3 Qn`0^a<=pF-7$F_u$\"3a(3*3`ArKpF-7$Fdu$\"3kvnE;UxXpF-7$Fiu$\"3Ws\"z*\\1 sfpF-7$F^v$\"3_?-x\">eB(pF-7$Fcv$\"3E[+'HOuh)pF-7$Fhv$\"3Cy-L.V<**pF-7 $F]w$\"3]A(zn)Hw7qF-7$Fbw$\"3\"pvO'3%[g-(F-7$Fgw$\"31_A49q&*RqF-7$F\\x $\"3_oGAaEN`qF-7$Fax$\"3!=fJ5j^q1(F-7$Ffx$\"3)p(eRtrj!3(F-7$F[y$\"3+4) ))e(37$4(F-7$F`y$\"3?53t2'Gu5(F-7$Fey$\"3x0Sl8cA?rF-7$Fjy$\"3=1y1\"*)p Q8(F-7$F_z$\"3!3QBb?Hp9(F-7$Fdz$\"3]****RAK'4;(F--Fiz6&F[[lF_[lF_[lF\\ [l-F$6$7%7$$\"3?++++]Y[PF-$\"#lF*7$F\\el$\"3[******p?ermF-7$F_[lFael-F iz6&F[[lF_[lF\\[lF_[l-F$6$7%7$$\"#QF*F^el7$Fjel$\"34+++qR,QoF-7$F_[lF] fl-Fiz6&F[[lF\\[lF\\[lF_[l-F$6$7%7$$\"3!)********\\`^QF-F^el7$Fffl$\"3 q+++qeW/qF-7$F_[lFifl-Fiz6&F[[l$\")#)eqkF^[l$\"))eqk\"F^[lF`gl-%(SCALI NGG6#%,CONSTRAINEDG-%&TITLEG6#Q " 0 "" {MPLTEXT 1 0 13 "sigmaIA/IA_a;" }} {PARA 11 "" 1 "" {XPPMATH 20 "6#$\"+))3!=X#!#6" }}}{EXCHG {PARA 0 "> \+ " 0 "" {MPLTEXT 1 0 10 "sigma/val;" }}{PARA 11 "" 1 "" {XPPMATH 20 "6# $\"+6U=c8!#6" }}}{EXCHG {PARA 0 "" 0 "" {TEXT -1 75 "The relative devi ation doubled when going from refracted to incident angle." }}}{EXCHG {PARA 0 "" 0 "" {TEXT -1 0 "" }{TEXT 276 11 "Exercise 3:" }}{PARA 0 " " 0 "" {TEXT -1 150 "Change the value of the exactly known reflected a ngle and repeat the simulation. Compare the prediction of the simulati on on the relationship between " }{TEXT 19 5 "d[IA]" }{TEXT -1 5 " and " }{TEXT 19 5 "d[RA]" }{TEXT -1 148 ", and compare with the analytica l result based on error propagation analysis. What happens in the abov e example when the glass is replaced by water?" }}{PARA 0 "" 0 "" {TEXT -1 41 "The index of refraction for water equals " }{TEXT 277 1 " n" }{TEXT -1 6 "=1.33." }}{PARA 0 "" 0 "" {TEXT -1 0 "" }}{PARA 0 "" 0 "" {TEXT 275 11 "Exercise 4:" }}{PARA 0 "" 0 "" {TEXT -1 289 "Use so me other physics law that involves a nonlinearity to relate the uncert ainties in two variables (a measured variable with known uncertainty i s connected to a deduced variable). Perform a simulation, and confirm \+ your results by an analytical derivation (as done above for Snell's la w)." }}{PARA 0 "" 0 "" {TEXT -1 0 "" }}}{EXCHG {PARA 0 "" 0 "" {TEXT 287 17 "Systematic errors" }}{PARA 0 "" 0 "" {TEXT -1 0 "" }}{PARA 0 " " 0 "" {TEXT -1 817 "So far we have assumed that the only source of er rors were random, and that no systematic deviations were present, i.e. , that the experimental apparatus was optimized at the level of desire d precision. One possibility to check whether the measuring apparatus \+ is by comparison (e.g., identically built ammeters are compared in the ir readings of the same current, known objects are measured with diffe rent calipers or micrometer screws, etc.). In this way one can establi sh the presence of systematic errors in the equipment (e.g., one ammet er is measuring consistently high values, etc.). Often we do not have \+ the luxury of comparing equipment in this way and rely on tolerance sp ecifications of the manufacturer who was able and obliged to establish the systematic tolerances for the equipment as it left the factory." }}{PARA 0 "" 0 "" {TEXT -1 0 "" }}{PARA 0 "" 0 "" {TEXT -1 890 "An imp ortant question that arises in this context is how to deal with this s ystematic error. Suppose we have established a standard deviation inte rval for a measurement sequence. How should we add the systematic meas urement error to the statistical uncertainty? John R. Taylor points ou t (in chapter 4) that two situations can arise: either the error is kn own for the equipment never to exceed a certain tolerance (e.g., 1% of the reading; or 2% of the maximum reading on the scale). In this case one would simply add the systematic error to the statistical one. On \+ the other hand, it might be true that based on a sequence of compariso ns the systematic error is known to fall within a normal distribution, i.e., 70% of the instruments are better than 1%, but some of them are outside this interval. In this latter case the systematic error has t o be added in quadrature to the uncertainty." }}{PARA 0 "" 0 "" {TEXT -1 50 "D[x]_tot = sqrt(D[x]_systematic^2 + D[x]_random^2)" }}{PARA 0 " " 0 "" {TEXT -1 0 "" }}{PARA 0 "" 0 "" {TEXT -1 488 "A justification f or adding independent errors (and uncertainties) in quadrature is give n below. Usually it is the second method of adding systematic and rand om errors that we wish to use in a teaching laboratory. Even if equipm ent was certified to be within a certain tolerance range, we cannot be assured that after years of use (and possibly abuse) it was still per forming to norm. We have more confidence in the statement that with so me likelihood it still performed to specifications." }}{PARA 0 "" 0 " " {TEXT -1 0 "" }}}{EXCHG {PARA 0 "" 0 "" {TEXT 278 46 "Uncertainty in a function of several variables" }}{PARA 0 "" 0 "" {TEXT -1 0 "" }} {PARA 0 "" 0 "" {TEXT -1 188 "We have demonstrated in the previous sec tion how an error propagates for a non-linear function of one variable . Effectively, the function was linearized in the range around the mea n value " }{TEXT 19 3 "x_a" }{TEXT -1 15 " and the range " }{TEXT 19 22 "x_a-sigma .. x_a+sigma" }{TEXT -1 51 " was mapped into a correpson ding range surrounding " }{TEXT 19 6 "f(x_a)" }{TEXT -1 28 " using the first derivative " }{TEXT 280 1 "f" }{TEXT -1 2 "'(" }{TEXT 279 1 "x " }{TEXT -1 65 "). Depending on the magnitude of the derivative the de viation in " }{TEXT 281 1 "f" }{TEXT -1 8 " (i.e., " }{TEXT 19 7 "sigm a_f" }{TEXT -1 42 ") can increase or decrease as compared to " }{TEXT 19 5 "sigma" }{TEXT -1 79 ". Now we demonstrate for a function of two \+ variables that the uncertainties in " }{TEXT 286 1 "x" }{TEXT -1 5 " a nd " }{TEXT 285 1 "y" }{TEXT -1 61 " have to be added in quadrature to obtain the uncertainty in " }{TEXT 284 1 "f" }{TEXT -1 1 "(" }{TEXT 283 1 "x" }{TEXT -1 2 ", " }{TEXT 282 1 "y" }{TEXT -1 2 ")." }}{PARA 0 "" 0 "" {TEXT -1 0 "" }}{PARA 0 "" 0 "" {TEXT -1 75 "First we genera te a second sequence of normally distributed random numbers." }}} {EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 36 "RNd2:=[stats[random, normald ](150)]:" }}}{EXCHG {PARA 0 "" 0 "" {TEXT -1 0 "" }}{PARA 0 "" 0 "" {TEXT -1 287 "We pick Snell's law again, and consider the following pr oblem: suppose we measure corresponding incident and refracted angles \+ in order to determine the index of refraction of the optically dense m edium (we assume that for the purposes of the experiment air has an in dex of refraction of " }{TEXT 290 1 "n" }{TEXT -1 28 "=1, i.e., the va cuum value)." }}{PARA 0 "" 0 "" {TEXT -1 0 "" }}{PARA 0 "" 0 "" {TEXT -1 216 "We measure the incident angle with the same accuracy as the re fracted angle. The data set for the refracted angle is kept from the p revious section, and we produce a data set for the incident angles (co nsistent with " }{TEXT 291 1 "n" }{TEXT -1 7 "=1.51)." }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 43 "MDataIA:=[seq(IA_true+0.5*RNd2[i],i =1..N)]:" }}}{EXCHG {PARA 0 "" 0 "" {TEXT -1 220 "For the purposes of \+ our simulation we simply associate with each datapoint from the set of refracted angles a corresponding datapoint from the set of incident a ngles to obtain a set of values for the index of refraction:" }}} {EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 73 "MDatan:=[seq(evalf(sin(Pi/18 0*MDataIA[i])/sin(Pi/180*MData[i])),i=1..N)]:" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 25 "histogram(MDatan,area=1);" }}{PARA 13 "" 1 "" {GLPLOT2D 607 296 296 {PLOTDATA 2 "6%-%)POLYGONSG6-7&7$$\"+Al%zY\"!\"* \"\"!7$F($\"+'pmmm#!#67$$\"+OK3w9F*F-7$F1F+7&F37$F1$\"+qLLLLF/7$$\"+]* >U[\"F*F67$F9F+7&F;7$F9$\"+knmm')F/7$$\"+kmN#\\\"F*F>7$FAF+7&FC7$FA$\" +9+++7!#57$$\"+yL\\+:F*FF7$FJF+7&FL7$FJ$\"+^LLL:FH7$$\"+#4I'3:F*FO7$FR F+7&FT7$FR$\"+SJLLFH7$$\"+@N !\\_\"F*Fin7$F\\oF+7&F^o7$F\\o$\"+;MLLtF/7$$\"+N-/L:F*Fao7$FdoF+7&Ffo7 $Fdo$\"+YLLL6FH7$$\"+\\p " 0 "" {MPLTEXT 1 0 38 "n_a:=evalf(add(MDatan[i],i=1..N)/N,7);" } }{PARA 11 "" 1 "" {XPPMATH 20 "6#>%$n_aG$\"(d;^\"!\"'" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 60 "sigma_n:=evalf(sqrt(add((MDatan[n]- n_a)^2,n=1..N)/(N-1)),5);" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%(sigma_ nG$\"&Dw\"!\"'" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 29 "sig_avg_n :='sigma_n/sqrt(N)';" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%*sig_avg_nG* &%(sigma_nG\"\"\"-%%sqrtG6#%\"NG!\"\"" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 28 "sig_avg_n:=evalf(sig_avg_n);" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%*sig_avg_nG$\"+C_2R9!#7" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 30 "[n_a-sig_avg_n,n_a+sig_avg_n];" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#7$$\"+Dz@5:!\"*$\"+vg48:F&" }}}{EXCHG {PARA 0 "" 0 "" {TEXT -1 18 "The true value of " }{TEXT 292 1 "n" }{TEXT -1 55 "=1.51 \+ falls barely outside the 70% confidence interval." }}{PARA 0 "" 0 "" {TEXT -1 0 "" }}{PARA 0 "" 0 "" {TEXT -1 35 "How can we predict the un certainty " }{TEXT 19 7 "sigma_n" }{TEXT -1 67 " from the known uncert ainties in the refracted and incident angles?" }}{PARA 0 "" 0 "" {TEXT -1 233 "First we should verify that the uncertainty in the incid ent angle is of the same magnitude as for the refracted angle (note th at it independently given in this section, and not determined via Snel l's law as in the previous section!):" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 40 "IA_a:=evalf(add(MDataIA[i],i=1..N)/N,7);" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%%IA_aG$\"(Z,%o!\"&" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 63 "sigma_IA:=evalf(sqrt(add((MDataIA[n]-IA_a)^2,n=1.. N)/(N-1)),5);" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%)sigma_IAG$\"&i0&! \"&" }}}{EXCHG {PARA 0 "" 0 "" {TEXT -1 70 "We need to apply the addit ion in quadrature rule given for a function " }{TEXT 19 18 "q(x1, x2, \+ ..., xn)" }{TEXT -1 4 " as:" }}{PARA 0 "" 0 "" {TEXT 19 89 "D[q] = sqr t((diff(q, x1)*D[x1])^2 + (diff(q, x2)*D[x2])^2 + ... + (diff(q, xn)*D [xn] )^2)" }}{PARA 0 "" 0 "" {TEXT -1 60 "where the derivatives are ev aluated at the average value of " }{TEXT 19 17 "[x1, x2, ..., xn]" } {TEXT -1 1 "." }}{PARA 0 "" 0 "" {TEXT -1 25 "Our function is given as :" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 24 "RF_ind:=sin(x1)/sin(x2 );" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%'RF_indG*&-%$sinG6#%#x1G\"\"\" -F'6#%#x2G!\"\"" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 58 "Dn:=sqrt ((diff(RF_ind,x1)*Dx1)^2+(diff(RF_ind,x2)*Dx2)^2);" }}{PARA 11 "" 1 " " {XPPMATH 20 "6#>%#DnG*$,&*(-%$cosG6#%#x1G\"\"#-%$sinG6#%#x2G!\"#%$Dx 1GF,\"\"\"**-F.F*F,F-!\"%-F)F/F,%$Dx2GF,F3#F3F," }}}{EXCHG {PARA 0 "" 0 "" {TEXT -1 66 "Here x1 and x2 are the incident and refracted angles respectively." }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 82 "evalf(sub s(x1=Pi/180*IA_a,x2=Pi/180*v_a,Dx1=Pi/180*sigma_IA,Dx2=Pi/180*sigma,Dn ));" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#$\"+LMj?=!#6" }}}{EXCHG {PARA 0 "" 0 "" {TEXT -1 204 "It is important in the above expression to pro vide the conversion factors for the angles from degrees to radians not only in the arguments to the sine functions, but also for the uncerta inties themselves!" }}{PARA 0 "" 0 "" {TEXT -1 142 "The agreement of t he predicted uncertainty with the one observed in the simulation (whic h represents a statistical measurement!) is very good." }}{PARA 0 "" 0 "" {TEXT -1 0 "" }}}{EXCHG {PARA 0 "" 0 "" {TEXT 294 11 "Exercise 5: " }}{PARA 0 "" 0 "" {TEXT -1 135 "Repeat the above simulation with a l arger deviation in the incident and refracted angles. Pick different i ncident and refracted angles." }}{PARA 0 "" 0 "" {TEXT -1 0 "" }} {PARA 0 "" 0 "" {TEXT 295 11 "Exercise 6:" }}{PARA 0 "" 0 "" {TEXT -1 264 "Repeat the simulation with a very much reduced statistical sample as corresponds to teaching laboratory measurements (or as happens som etimes in life science research when circumstances do not permit large sample sizes), e.g., N=5. Are the conclusions still valid?" }}} {EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 0 "" }}}}{MARK "21 0 0" 0 } {VIEWOPTS 1 1 0 3 2 1804 1 1 1 1 }{PAGENUMBERS 0 1 2 33 1 1 }