b4m.ps

Interval Arithmetic Toolbox for Matlab

D<B512F814FF390FC01F80EC0FC0140715E0A515C0EC0F80EC3F00EBFFFE9038C00F80EC
07C0EC03E0A215F0A5EC07E0A2EC1FC0B6120014FC1C1C7E9B22>I<90381FE0209038FF
F8E03803F80F3807C003380F800148C7FC123E1560127E127C00FC1400A8007C1460127E
123E15C07E390F8001803907C003003803F80E3800FFFCEB1FE01B1C7D9B22>I<B512F8
14FF390FC01FC0EC07E0EC01F015F8140015FC157CA2157EA8157CA315F8A2EC01F0EC07
E0EC1FC0B6120014F81F1C7E9B25>I<EAFFFCA2EA0FC0B3A6EAFFFCA20E1C7E9B12>73
D<D8FFC0EB03FF6D5B000F15F0D80DF0130DA3D80CF81319A2017C1331A26D1361A26D13
C1A390380F8181A2903807C301A2EB03E6A2EB01FCA3EB00F8A23AFFC0701FFFA2281C7E
9B2D>77 D<EB3FE03801F07C3803C01E390F800F80391F0007C0001E1303003E14E0007E
14F0007C1301A200FC14F8A8007C14F0007E1303A2003E14E0003F13076C14C0390F800F
803903C01E003801F07C38003FE01D1C7D9B24>79 D<B512F014FE380FC03FEC0F806E7E
81A55D4A5A023EC7FCEBFFF0EBC07C80143F6E7EA6ED8180A2EC0FC13AFFFC07C300EC01
FE211C7E9B24>82 D<3807F820381FFEE0EA3C07EA7801EA700012F01460A26C130012FE
EAFFE0EA7FFE6C7E1480000F13C06C13E0EA007FEB03F01301130012C0A214E07E38F001
C0EAFC0338EFFF00EA83FC141C7D9B1B>I<007FB512E0A238781F810070138000601460
00E0147000C01430A400001400B03807FFFEA21C1C7E9B21>I<EA0FF8EA1C1E383E0F80
130714C0121C1200137FEA07E7EA1F07123C127C12F8A3EA780B383E13F8EA0FE115127F
9117>97 D<B4FCA2121FA9EB3F80EBE1E0EB8070EB0078143CA2143EA6143CA21478EB80
70381EC1E0381C3F00171D7F9C1B>I<EA03FCEA0E0EEA1C1F123C1278130E00F8C7FCA6
1278EB0180123C381C0300EA0E0EEA03F811127E9115>I<EB0FF0A21301A9EA01F9EA0F
07EA1C03EA3C011278A212F8A61278A2123CEA1C03380F0FFEEA03F9171D7E9C1B>I<EA
03FCEA0F07381C0380383C01C0127814E012F8A2B5FC00F8C7FCA3127814607E6C13C038
0F03803801FC0013127F9116>I<137F3801E3803803C7C0EA0787120FEB8380EB8000A5
EAFFF8A2EA0F80AEEA7FF0A2121D809C0F>I<3803F0F0380E1F38EA3C0F383807300078
1380A400381300EA3C0FEA1E1CEA33F00030C7FCA3EA3FFF14C06C13E014F0387801F838
F00078A300701370007813F0381E03C03807FF00151B7F9118>I<B4FCA2121FA9EB0FC0
EB31E0EB40F0EB80F8A21300AB38FFE7FFA2181D7F9C1B>I<121E123FA4121EC7FCA612
7FA2121FAEEAFFC0A20A1E7F9D0E>I<B4FCA2121FB3A7EAFFE0A20B1D7F9C0E>108
D<39FF0FC07E903831E18F3A1F40F20780D980FC13C0A2EB00F8AB3AFFE7FF3FF8A22512
7F9128>I<38FF0FC0EB31E0381F40F0EB80F8A21300AB38FFE7FFA218127F911B>I<EA01
FC380F0780381C01C0003C13E0387800F0A200F813F8A6007813F0A2383C01E0381E03C0
380F07803801FC0015127F9118>I<38FF3F80EBE1E0381F80F0EB0078147C143C143EA6
143C147C1478EB80F0EBC1E0EB3F0090C7FCA6EAFFE0A2171A7F911B>I<EAFE3E134738
1E8F80A3381F070090C7FCAAEAFFE0A211127F9114>114 D<EA1FD8EA3078EA601812E0
A2EAF000EAFF80EA7FE013F0EA1FF8EA07FCEA007CEAC01CA212E01318EAF830EACFC00E
127E9113>I<1203A45AA25AA2EA3FFC12FFEA1F00A9130CA4EA0F08EA0798EA03F00E1A
7F9913>I<38FF07F8A2EA1F00AC1301120F380786FFEA01F818127F911B>I<38FFC1FCA2
381F00601380000F13C0A23807C180A23803E300A213F7EA01F613FE6C5AA21378A21330
16127F9119>I<39FF8FF8FEA2391F03E030A3390F87F06013869038C6F8E03907CC78C0
A23903FC7D80EBF83D143F3901F01F00A20000131EEBE00EA21F127F9122>I<38FFC7FC
A2381F81C0380F83803807C700EA03EEEA01FC5B1200137C13FEEA01DF38039F80EA070F
380607C0380C03E038FF07FCA216127F9119>I E /Fo 3 110 df<ED7FC04B7EA25C5CA2
EC07F7A2EC0FE7EC1FC7A2EC3F87147F150714FEA2EB01FCEB03F8A2EB07F0EB0FE0A2EB
1FC0133F1480EB7F00A213FE485AA2485A485AA2485A121F5B485AA248C7FC12FEB812F8
17FCA46C16F8C83807E000AC020FB512F04A14F8A46E14F02E407DBF35>52
D<EA7FFC487EA4127F1200AD15FE913807FFC0023F13F04A13FC91B57E90B7FCDAFE0313
80ECF80002E0EB7FC04AEB3FE04A131F91C7EA0FF016074915F81603A317FC1601A81603
17F8A26D140717F0160F806EEB1FE0EE3FC06E137F9139F801FF80DAFC07130091B55A49
5C6E5B6E13E0D97E0F138090263C01FCC7FC2E3F7EBE35>98 D<267FC07C131F3BFFE3FF
80FFE001E701C17F01EF01E37F90B512E76C91B57E0007018713E1D9FE03138001FC4A7E
496C48137EA201F05BA301E05BB3A63C7FFE07FF81FFE0B54801C313F04A13C7A26E13C3
6C486C018113E0342C80AB35>109 D E end
TeXDict begin
%%PaperSize: a4

1 0 bop 894 283 a Fo(b4m)335 366 y Fn(A)16 b(free)f(in)o(terv)m(al)e
(arithmetic)g(to)q(olb)q(o)o(x)h(for)h(Matlab)g(based)g(on)g(BIAS)807
416 y(v)o(ersion)e(1.02.004)665 466 y(do)q(cumen)o(tation)f
Fm(\014)r Fn(-v)o(ersion)j(1.00)868 1327 y Fl(Jens)g(Zemk)o(e)594
1377 y(T)m(ec)o(hnisc)o(he)g(Univ)o(ersit\177)-21 b(at)14
b(Ham)o(burg-Harburg)742 1427 y(T)m(ec)o(hnisc)o(he)h(Informatik)d(I)q
(I)q(I)776 1476 y(Ei\031endorfer)j(Stra\031e)f(38)830
1526 y(21071)f(Ham)o(burg)690 1576 y Fk(email:)21 b(zemke@tu-harburg)o
(.de)429 2460 y(http://www.ti3.t)o(u-ha)o(rburg)o(.de/~)o(zemke)o
(/b4m/)o(index)o(.htm)o(l)p eop
2 1 bop 257 262 a Fj(Con)n(ten)n(ts)257 353 y Fn(1)39
b(In)o(tro)q(ducti)o(on)1080 b(3)257 444 y(2)39 b(The)15
b(to)q(olb)q(o)o(x)g(b4m)985 b(3)257 535 y(3)39 b(Installati)o(on)13
b(of)i(b4m)948 b(4)257 627 y(4)39 b(The)15 b(class)g(in)o(terv)m(al)978
b(6)257 718 y(5)39 b(Ov)o(erloaded)13 b(arithmetic)g(op)q(erators)663
b(7)257 809 y(6)39 b(Ov)o(erloaded)13 b(functions)g(with)i(one)g
(argumen)o(t)466 b(10)257 901 y(7)39 b(Ov)o(erloaded)13
b(comparison)h(op)q(erations)597 b(11)257 992 y(8)39
b(Miscellaneous)12 b(o)o(v)o(erloaded)i(op)q(erations)559
b(12)257 1083 y(9)39 b(In)o(terv)m(al)14 b(sp)q(eci\014c)h(op)q
(erations)757 b(14)257 1175 y(10)15 b(Real)g(op)q(erations)e(in)i(b4m)
832 b(17)257 1266 y(11)15 b(Do)q(cumen)o(tatio)o(n)e(and)i(help)804
b(18)257 1357 y(12)15 b(Do)o(wnload)f(BIAS)h(and)g(b4m)786
b(20)257 1449 y(A)27 b(Complete)14 b(list)g(of)h(M-\014les)g(used)f(in)
h(b4m)558 b(20)257 1540 y(B)29 b(Complete)14 b(list)g(of)h(Mex-\014les)
g(used)f(in)h(b4m)511 b(22)257 2032 y Fj(List)24 b(of)d(Figures)320
2123 y Fl(1)74 b(Setting)14 b(of)f(con)o(trols)h(in)g
Fk(BiasInt.h)35 b Fl(.)20 b(.)g(.)g(.)h(.)f(.)g(.)h(.)f(.)g(.)g(.)h(.)f
(.)g(.)h(.)f(.)g(.)g(.)66 b(5)320 2172 y(2)74 b(The)14
b(class)h Fk(interval)d Fl(in)h Fk(b4m)g Fl(.)21 b(.)f(.)g(.)h(.)f(.)g
(.)g(.)h(.)f(.)g(.)h(.)f(.)g(.)g(.)h(.)f(.)g(.)h(.)f(.)g(.)g(.)66
b(6)320 2222 y(3)74 b(The)14 b(arithmetic)f(op)q(erations)h(in)g
Fk(b4m)25 b Fl(.)20 b(.)g(.)g(.)h(.)f(.)g(.)h(.)f(.)g(.)g(.)h(.)f(.)g
(.)h(.)f(.)g(.)g(.)66 b(8)320 2272 y(4)74 b(Ov)o(erloaded)14
b(functions)g(in)g Fk(b4m)k Fl(.)i(.)g(.)h(.)f(.)g(.)g(.)h(.)f(.)g(.)h
(.)f(.)g(.)g(.)h(.)f(.)g(.)h(.)f(.)g(.)g(.)45 b(10)320
2322 y(5)74 b(Ov)o(erloaded)14 b(comparison)f(op)q(erators)i(in)e
Fk(b4m)20 b Fl(.)g(.)g(.)h(.)f(.)g(.)g(.)h(.)f(.)g(.)h(.)f(.)g(.)g(.)45
b(12)320 2372 y(6)74 b(Miscellaneous)14 b(o)o(v)o(erloaded)g(op)q
(erations)g(in)f Fk(b4m)36 b Fl(.)21 b(.)f(.)g(.)g(.)h(.)f(.)g(.)h(.)f
(.)g(.)g(.)45 b(13)320 2421 y(7)74 b(In)o(terv)n(al)14
b(sp)q(eci\014c)h(op)q(erations)f(in)g Fk(b4m)f Fl(.)20
b(.)g(.)g(.)h(.)f(.)g(.)h(.)f(.)g(.)g(.)h(.)f(.)g(.)h(.)f(.)g(.)g(.)45
b(15)320 2471 y(8)74 b(Real)13 b(op)q(erations)h(in)g
Fk(b4m)54 b Fl(.)20 b(.)h(.)f(.)g(.)h(.)f(.)g(.)g(.)h(.)f(.)g(.)h(.)f
(.)g(.)g(.)h(.)f(.)g(.)h(.)f(.)g(.)g(.)45 b(17)963 2628
y(2)p eop
3 2 bop 257 262 a Fj(1)67 b(In)n(tro)r(duction)257 352
y Fl(In)18 b(the)g(\014eld)f(of)g(n)o(umerical)f(computations)g(the)i
(fo)q(cus)g(has)g(c)o(hanged)g(during)f(the)h(last)257
402 y(decades.)30 b(F)m(or)17 b(most)f(problem)g(classes)j(standard)f
(algorithms)d(are)i(a)o(v)n(ailable)e(and)i(\(at)257
452 y(least\))c(in)e(the)i(case)g(of)e(small)f(sized)j(problems)e(a)h
(forw)o(ard/bac)o(kw)o(ard)f(error)i(analysis)e(has)257
502 y(b)q(een)k(carried)g(out.)i(These)e(standard)g(algorithms)c(deliv)
o(er)j(go)q(o)q(d)f(results)i(in)e(the)h(case)h(of)257
552 y(w)o(ell)f(p)q(osed)i(problems,)d(but)i(no)f(information)e
(whether)k(a)f(giv)o(en)f(problem)f(migh)o(t)g(b)q(e)i(ill)257
602 y(p)q(osed.)k(New)c(tasks)f(of)f(in)o(terest)i(are)320
682 y Fi(\017)20 b Fl(the)15 b(qualit)o(y)d(of)h(a)h(solution)f(of)g
(an)h(ill)e(p)q(osed)j(problem,)320 764 y Fi(\017)20
b Fl(the)15 b(solution)e(of)g(problems)g(with)g(uncertain)o(ties)i(in)f
(the)g(data)257 845 y(and/or)320 925 y Fi(\017)20 b Fl(to)15
b(pro)o(v)o(e)g(uniqueness)h(and)f(existence)i(of)e(solutions)f(\(e.g.)
21 b(in)15 b(the)h(\014eld)f(of)f(global)361 975 y(optimization\).)257
1056 y(F)m(or)g(this)h(tasks)g(in)e(most)h(cases)h(forw)o(ard)f(and)g
(bac)o(kw)o(ard)g(error)i(analysis)d(is)h(unsatisfac-)257
1106 y(tory)k(\(sometimes)d(ev)o(en)j(not)f(p)q(ossible\).)29
b(One)18 b(idea)f(is)g(to)g(use)h(p)q(erturbation)g(analysis)257
1156 y(`on)e(the)g(\015y')f(to)h(deriv)o(e)h(sharp)f(b)q(ounds)g(for)g
(the)g(p)q(erturbation)h(in)e(the)i(\014nite)f(precision)257
1205 y(solution)d(of)h(a)f(giv)o(en)h(problem.)i(This)e(can)g(b)q(e)h
(done)f(b)o(y)g(means)f(of)g(in)o(terv)n(al)g(arithmetic.)320
1255 y(In)h(con)o(trast)i(to)e(pure)i(\015oating)e(p)q(oin)o(t)g
(arithmetic)g(in)o(terv)n(al)g(arithmetic)f(can)i(b)q(e)h(im-)257
1305 y(plemen)o(ted)h(in)f(a)h(reliable)f(w)o(a)o(y)g(\(the)i(result)g
(of)e(an)g(in)o(terv)n(al)g(arithmetic)g(computation)257
1355 y(is)i(guaran)o(teed)h(to)f(include)h(the)g(true)g(solution\).)30
b(Readers)19 b(not)g(famili)o(ar)c(with)j(in)o(ter-)257
1405 y(v)n(al)c(arithmetic)g(are)h(referred)i(to)e(the)g(in)o(tro)q
(ductory)g(b)q(o)q(oks)h(of)e(Hansen)i(and)e(Neumaier)257
1454 y(\([1,)f(4]\).)320 1504 y(There)j(exist)h(man)o(y)c(implemen)o
(tations)g(of)i(in)o(terv)n(al)g(arithmetic.)22 b(Some)15
b(are)h(exten-)257 1554 y(sions)i(of)f(existing)g(languages)g(\(e.g.)29
b(libraries\),)18 b(but)f(ev)o(en)i(new)e(languages)g(including)257
1604 y(in)o(terv)n(al)c(arithmetic)g(ha)o(v)o(e)h(b)q(een)g(prop)q
(osed.)19 b(A)14 b(v)o(ery)g(fast)g(implemen)o(tatio)o(n)d(of)i(in)o
(terv)n(al)257 1654 y(arithmetic)f(is)h(giv)o(en)g(with)g(the)h(ANSI)f
(C)g(library)f(BIAS)i(b)o(y)f(Olaf)f(Kn)q(\177)-22 b(upp)q(el)14
b(\([3]\).)j(F)m(or)c(a)257 1704 y(complete)h(in)o(tro)q(duction)f(to)h
(BIAS)g(see)i([2)o(].)320 1753 y(On)j(the)h(other)f(hand,)h
Fh(Ma)m(tlab)g Fl(is)f(the)h(standard)f(to)q(ol)f(for)h(in)o(teractiv)o
(e)g(\015oating)257 1803 y(p)q(oin)o(t)14 b(computations,)e(but)i(no)g
(in)o(terv)n(al)f(arithmetic)f(coun)o(terpart)j(has)f(ev)o(olv)o(ed)g
(y)o(et.)257 1940 y Fj(2)67 b(The)22 b(to)r(olb)r(o)n(x)h(b4m)257
2031 y Fl(The)15 b(to)q(olb)q(o)o(x)e Fk(b4m)g Fl(has)h(b)q(een)h
(designed)320 2112 y Fi(\017)20 b Fl(to)14 b(use)g(the)g(ANSI)g(C)g(in)
o(terv)n(al)f(arithmetic)f(library)h(BIAS)h(in)g(an)f(in)o(teractiv)o
(e)h(w)o(a)o(y)257 2192 y(and)320 2273 y Fi(\017)20 b
Fl(to)11 b(add)g(in)o(terv)n(al)f(arithmetic)g(op)q(erations)i(and)e
(inclusion)h(algorithms)e(to)i(the)g(stan-)361 2323 y(dard)j
(\015oating)f(p)q(oin)o(t)g(en)o(vironmen)o(t)g Fh(Ma)m(tlab)p
Fl(.)257 2403 y Fk(b4m)g Fl(consists)h(of)f(a)g(new)h(datat)o(yp)q(e)f
Fk(interval)f Fl(and)h(of)f(arithmetic)h(op)q(erations)g(b)q(et)o(w)o
(een)257 2453 y(argumen)o(ts)j(of)g(t)o(yp)q(e)h Fk(double)f
Fl(and)g Fk(interval)p Fl(.)25 b(It)16 b(also)g(con)o(tains)h
(functions)g(handling)257 2503 y(data)g(of)g(t)o(yp)q(e)h
Fk(interval)p Fl(.)27 b(All)16 b(these)j(elemen)o(tary)d(op)q(erations)
i(are)g(p)q(erformed)f(using)963 2628 y(3)p eop
4 3 bop 257 262 a Fl(BIAS)13 b(\(Basic)g(In)o(terv)n(al)f(Arithmetic)g
(Subroutines\))h(as)f(prop)q(osed)h(and)g(implem)o(en)o(ted)e(b)o(y)257
311 y(Kn)q(\177)-22 b(upp)q(el)16 b(in)e(1993.)20 b(BIAS)c(do)q(es)f
(not)g(supp)q(ort)h(the)g(use)f(of)g(sparse)h(and)e(complex)g(op)q(er-)
257 361 y(ations,)f(so)h(up)g(to)g(no)o(w)g(no)f(op)q(erations)i(in)o
(v)o(olving)c(argumen)o(ts)i(of)h(t)o(yp)q(e)g Fk(sparse)f
Fl(and/or)257 411 y Fk(complex)18 b Fl(are)i(p)q(ossible.)36
b(Nev)o(ertheless,)23 b(it)c(is)g(p)q(ossible)h(to)f(extend)i
Fk(b4m)e Fl(to)g(p)q(erform)257 461 y(computations)13
b(in)g(the)i(complex)d(\014eld)i(and)g(to)g(p)q(erform)f(sparse)i
(computations.)320 511 y Fk(b4m)d Fl(tries)j(to)e(com)o(bine)f(the)i
(e\016ciency)h(of)e(BIAS)h(and)f(the)h(ease)h(of)d(use)j(of)h
Fh(Ma)m(tlab)p Fl(.)257 560 y(All)c(in)o(terv)n(al)f(data)h(is)g
(stored)h(in)e(BIAS)i(format,)d(so)i(the)h(computing)d(time)h(splits)h
(roughly)257 610 y(in)o(to)17 b(the)i(time)d Fh(Ma)m(tlab)j
Fl(needs)g(for)e(the)h(in)o(terpretation)g(of)f(the)i(co)q(de)f(and)g
(the)g(time)257 660 y(BIAS)d(needs)g(to)f(p)q(erform)f(the)h(op)q
(erations.)320 710 y(The)19 b(in)o(terface)g(b)q(et)o(w)o(een)h(BIAS)f
(and)f Fh(Ma)m(tlab)i Fl(is)e(implemen)o(ted)e(in)i(the)h(form)e(of)257
760 y(Mex-\014les.)i(Those)13 b(are)g(stored)h(hidden)f(from)e(the)j
(user)g(in)e(priv)n(ate)g(sub)q(directories.)20 b(The)257
809 y(user)d(in)o(terface)e(consists)i(of)d(ordinary)h(M-\014les)g
(with)g(inline)g(help,)g(do)q(cumen)o(tation)f(and)257
859 y(error)h(handling.)i(Ev)o(ery)d(M-\014le)f(calls)h(appropriate)f
(Mex-\014le\(s\).)19 b(It)14 b(is)g(tried)g(to)f(ac)o(hiev)o(e)257