cats_v312.ps

最大似然估计算法

2689 y Fi(2)2575 2584 y Fc(!)2650 2580 y Fb(\024)p 2650
2592 37 4 v 2653 2633 a Fa(2)3445 2730 y Fn(\(7\))568
3052 y(The)38 b(step-v)n(ariable)g(white)g(noise)g(model)g(is)g(an)h
(ad-hoc)g(model)e(de)n(vised)h(so)g(that)568 3172 y(we)30
b(can)h(introduce)e(a)i(step)e(function)h(in)f(the)h(size)h(of)f(the)g
(white)f(nosie)h(component)568 3293 y(between)35 b(tw)o(o)g
(speci\002ed)h(times.)62 b(Y)-11 b(ou)36 b(can)g(let)f(the)g
(minimization)f(solv)o(e)g(for)i(the)568 3413 y(times)28
b(b)n(ut)h(this)f(is)h(v)o(ery)g(unstable)g(\(as)g(its)g(not)g(a)h
(continuous)d(function\))i(and)g(is)g(not)568 3534 y(recommended.)44
b(This)29 b(noise)f(model)h(should)f(AL)-7 b(W)-12 b(A)i(YS)29
b(be)h(used)f(in)g(conjunction)568 3654 y(with)36 b(a)g(white)h(noise)f
(or)g(v)n(ariable)g(white)g(noise)g(model.)66 b(The)36
b(model)g(should)g(be)568 3774 y(speci\002ed)25 b(as)g(follo)n(ws)568
4010 y Fd(--model)58 b(sw:b1997.0e1998.0)568 4245 y Fn(where)30
b Fk(b)g Fn(is)g(the)f(start)h(time)f(of)h(when)f(the)h(white)f(noise)g
(changes)h(and)g Fk(e)g Fn(is)f(the)h(end)568 4365 y(time.)g(Y)-11
b(ou)24 b(can)h(specify)f(start)g(or)h(end)g(times)e(outside)h(the)g
(range)h(of)g(the)f(time)g(series)568 4486 y(to)g(gi)n(v)o(e)g(a)h
(single)f(step)g(function)g(in)g(the)h(noise)f(\(b)n(ut)h(not)f
(both\).)568 4646 y(F)o(or)37 b(each)h(model)f(you)g(might)f(w)o(ant)h
(to)g(specify)g(the)h(sigmas)e(by)h(appending)g(the)568
4767 y(additional)23 b(parameter)i(as)g(follo)n(ws.)568
5002 y Fd(--model)58 b(wh:s0.0005/s0.0005/s0.0015)1918
5251 y Fn(7)p eop
%%Page: 8 8
8 7 bop 568 548 a Fn(The)25 b(/)h(seperation)f(allo)n(ws)g(you)g(to)g
(speci\002y)h(dif)n(ferent)f(sigmas)g(for)h(each)g(column)f(of)568
668 y(data.)30 b(The)22 b(sigma)f(\224option\224)g(is)h(almost)f
(meaningless)f(at)j(the)f(moment,)f(since)h(I)g(ha)n(v)o(e)568
789 y(virtually)32 b(eliminated)g(the)i(need)g(for)g(starting)e(v)n
(alues.)57 b(Ho)n(we)n(v)o(er)32 b(if)i(you)f(w)o(ant)g(to)568
909 y(fully)21 b(specify)h(a)g(model)f(and)h(just)g(compute)f(the)h
(weighted)f(least)h(squares)g(\(\226method)568 1029 y(W\))j(then)f(the)
h(sigmas)f(are)h(required)g(and)g(should)e(be)i(imput)f(as)h(sho)n(wn.)
324 1233 y Fo(--method)h Fk(<)p Fo(M/S/E/W)p Fk(>)e Fo(\(-B\))50
b Fn(Estimation)27 b(method.)41 b(Choose)28 b(from)h(Maximum)e(Lik)o
(e-)568 1353 y(lihood)17 b(Estimation)g(\(MLE\),)h(Empirical)g
(estimation,)g(Spectral)i(estimation)d(or)i(W)-8 b(eighted)568
1474 y(Least)24 b(Squares)i(Estimation.)i(Additional)23
b(parameter)j(should)d(be)i(M,E,S)g(or)g(W)-9 b(.)324
1677 y Fo(--columns)25 b Fk(<)p Fo(number)p Fk(>)i Fo(\(-C\))49
b Fn(Only)32 b(use)i(a)f(certain)h(column)e(from)h(the)g(data)g
(\002le.)57 b(The)568 1797 y(number)30 b(after)h(the)g(--columns)f
(de\002nes)h(which)f(columns)f(to)i(use.)48 b(The)31
b(number)f(is)568 1918 y(basically)g(the)g(decimal)g(representation)g
(of)h(a)g(binary)f(sequence;)j(1)e(for)g(use)f(and)h(0)568
2038 y(for)i(ignore.)55 b(So)33 b(if)g(there)g(are)g(three)h(columns)d
(and)i(you)g(wish)f(to)g(use)h(the)g(middle)568 2159
y(column)d(this)f(w)o(ould)i(be)g(in)f(binary)h(010,)g(which)g(is)f(2)h
(in)g(decimal)f(so)h(the)f(number)568 2279 y(w)o(ould)d(be)g(2.)40
b(If)28 b(you)f(w)o(ant)h(to)f(use)h(all)f(3)h(columns)e(then)h(the)h
(number)f(w)o(ould)g(be)h(7.)568 2399 y(This)22 b(is)h(essentially)f
(scalable,)i(if)f(in)g(the)h(future)f(there)h(is)f(a)g(\002le)h(with)f
(more)g(columns)568 2520 y(then)e(it)g(can)i(cope)e(with)g(this.)29
b(Of)22 b(course)g(if)g(you)f(w)o(ant)h(to)f(use)h(all)f(columns)g
(then)g(you)568 2640 y(dont)j(need)h(to)f(use)h(this)f(option.)324
2844 y Fo(--scale)h(scale)p 831 2844 30 4 v 35 w(factor)h(\(-S\))49
b Fn(.)54 b(Di)n(vide)31 b(the)h(data)g(by)g(scale)p
2412 2844 V 36 w(f)o(actor)h(before)f(estimating)f(the)568
2964 y(model.)324 3167 y Fo(--sinusoid)25 b Fk(<)p Fo(period)p
Fk(><)p Fo(time)p 1462 3167 V 36 w(\003ag)p Fk(><)p Fo(harmonics)p
Fk(>)g Fo(\(-A\))49 b Fn(Include)21 b(a)g(sinusoid)e(\(and)h
Fk(n)568 3288 y Fn(harmonics\))42 b(of)h(period)g Fk(<)p
Fn(period)p Fk(>)f Fn(to)g(the)h(linear)g(parameters)g(to)g(solv)o(e.)
84 b(Use)568 3408 y(time)p 752 3408 V 35 w(\003ag)25
b(to)f(describe)h(what)g(units)f(the)g(period)h(is)f(in)g(:)568
3611 y Fo(y)48 b Fn(years)568 3773 y Fo(d)h Fn(days)568
3935 y Fo(h)g Fn(hours)568 4097 y Fo(m)f Fn(minutes)568
4259 y Fo(s)g Fn(seconds)568 4462 y(T)-8 b(o)28 b(solv)o(e)f(for)i(an)g
(annual)f(sinusoid)f(use)i Fo(-)f(-sinusoid)h(1y)p Fn(.)42
b(T)-8 b(o)28 b(solv)o(e)g(for)g(an)h(annual)568 4583
y(sinusoid)23 b(and)h(a)i(semi-annual)d(use)i Fo(-)g(-sinusoid)g(1y1)p
Fn(.)324 4786 y Fo(--v)o(erbose)h(\(-V\))49 b Fn(v)o(erbose.)30
b(Output)24 b(more)g(information.)1918 5251 y(8)p eop
%%Page: 9 9
9 8 bop 324 548 a Fo(--speed)26 b Fk(<)p Fo(0/1/2/3)p
Fk(>)d Fo(\(-Z\))50 b Fn(Speed)32 b(up)g(the)g(computation.)51
b(At)32 b(e)n(v)o(ery)f(point)g(in)h(the)g(mini-)568
668 y(mization)27 b(algorithm)g(the)h(noise)f(parameters)i(are)g
(chosen,)g(the)f(co)o(v)n(ariance)g(matrix)568 789 y(is)38
b(calculated,)43 b(a)c(weighted)g(least)f(squares)h(is)g(performed)g
(on)g(the)g(data)g(and)g(the)568 909 y(residuals)c(calculated.)66
b(The)36 b(residuals,)i(along)e(with)g(the)g(noise)g(parameters,)j(are)
568 1029 y(used)21 b(to)h(calculate)h(the)e(MLE)h(v)n(alue)g(\(which)f
(is)h(the)g(objecti)n(v)o(e)e(function)i(in)f(the)h(min-)568
1150 y(imization\).)28 b(Ho)n(we)n(v)o(er)20 b(if)h(you)g(ha)n(v)o(e)g
(good)g(starting)g(parameters)g(then)g(the)g(residuals)568
1270 y(will)33 b(not)h(alter)g(much)g(between)h(dif)n(ferent)f(noise)g
(parameter)g(choices)h(and)f(so)g(the)568 1391 y(minimization)19
b(can)k(be)f(made)h(f)o(aster)f(by)g(using)f(the)i(same)f(residuals)f
(for)i(man)o(y)e(parts)568 1511 y(of)h(the)g(algorithm.)28
b(In)22 b(general)g(this)f(gi)n(v)o(es)g(results)g(that)g(are)i(a)f
(tenth)g(of)g(a)g(mm)f(or)h(less)568 1631 y(dif)n(ferent)i(than)g(the)g
(full)f(blo)n(wn)g(algorithm)g(for)i(a)f(considerable)g(increase)h(in)f
(speed.)568 1752 y(Note)d(this)f(w)o(as)i(more)f(of)h(an)f(issue)g(in)g
(the)g(pre)n(vious)g(v)o(ersion)f(of)h(the)h(softw)o(are,)g(when)568
1872 y(using)k(the)h(simple)o(x)f(algorithm.)36 b(This)27
b(ne)n(w)g(v)o(ersion)f(performs)h(at)g(the)h(same)f(speed)568
1993 y(as)j(the)h(pre)n(vious)e(e)n(v)o(en)h(without)f(this)h
(\224cheat\224.)48 b(The)31 b(higher)f(the)h(number)f(the)g(less)568
2113 y(times)23 b(the)i(residuals)f(are)i(estimated)d(in)i(the)g
(algorithm.)324 2309 y Fo(--help)h(\(-H\))49 b Fn(Bring)24
b(up)h(a)g(help)g(page.)324 2504 y Fo(--delta)g Fk(<)p
Fo(Delta)p Fk(>)g Fo(\(-D\))49 b Fn(This)27 b(sets)h(the)f(accurac)o(y)
i(of)f(the)f(starting)g(parameters.)40 b(Delta)28 b(is)568
2625 y(1.1)37 b(for)g(10\045)g(accurac)o(y)-6 b(,)41
b(1.01)36 b(for)i(1\045)f(...)68 b(This)37 b(is)f(used)h(in)g(the)h
(simple)o(x)d(algo-)568 2745 y(rithm)28 b(for)i(controlling)e(the)h
(size)g(of)h(the)f(initial)f(w)o(alks)h(through)g(parameter)g(space.)
568 2865 y Fo(Curr)n(ently)d(defunct.)324 3061 y(--tolerance)g
Fk(<)p Fo(tol1)p Fk(>)p Fo(/)p Fk(<)p Fo(tol2)p Fk(>)p
Fo(/)p Fk(<)p Fo(tol3)p Fk(>)c Fo(\(-T\))50 b Fn(T)-8
b(olerance)25 b(v)n(alues)e(for)i(the)f(simple)o(x)f(algo-)568
3181 y(rithm.)39 b Fo(tol1)27 b Fn(is)h(the)f(PRECISION)j(required)e
(for)g(all)g(v)n(alues)f(of)h(the)g(noise)f(parame-)568
3302 y(ters)f(at)h(solution.)35 b Fo(tol2)26 b Fn(is)h(the)f(PRECISION)
j(required)e(for)g(the)f(objecti)n(v)o(e)g(function)568
3422 y(\(the)31 b(MLE\))h(at)g(solution.)50 b Fo(tol3)31
b Fn(:)45 b(is)31 b(the)h(MAXIMUM)e(V)-13 b(ALUE)31 b(of)h(all)g(the)f
(noise)568 3542 y(parameters)25 b(\(OR-ed)g(with)f Fo(tol2)p
Fn(\).)31 b Fo(Curr)n(ently)26 b(defunct.)324 3738 y(--output)g
Fk(<)p Fo(output)p 1060 3738 30 4 v 37 w(\002le)p Fk(>)g
Fo(\(-O\))48 b Fn(All)30 b(results)f(are)i(written)f(to)g(this)f
(\002le.)47 b(If)31 b(this)e(option)g(is)568 3858 y(not)24
b(set)g(then)h(the)g(output)e(goes)i(to)f Fo(stdout)p
Fn(.)324 4054 y Fo(--\002letype)i Fk(<)p Fo(cats/psmsl)p
Fk(>)e Fo(\(-F\))49 b Fn(Filetype,)20 b(either)g Fo(cats)g
Fn(\(Def)o(ault\))g(or)g Fo(psmsl)p Fn(.)29 b(See)21
b(abo)o(v)o(e.)324 4250 y Fo(--psd\002le)26 b Fk(<)p
Fo(\002lename)p Fk(>)f Fo(\(-P\))49 b Fn(Create)23 b(the)f(po)n(wer)f
(spectrum)g(for)h(this)e(data)i(set)f(and)h(output)568
4370 y(it)k(to)h Fk(<)p Fn(\002lename)p Fk(>)p Fn(.)38
b(If)28 b(the)f(dataset)g(is)g(complete)f(\(no)h(gaps)g(or)g(missing)f
(data\))h(then)568 4490 y(the)d(psd)h(is)f(generated)h(from)g(the)f
(FFT)i(otherwise)e(scar)n(gle')-5 b(s)25 b(periodogram)f(is)g(used.)324
4686 y Fo(--notr)n(end)j(\(-E\))50 b Fn(Do)24 b(no)h(compute)f(a)h
(trend)g(to)f(the)h(data.)324 4882 y Fo(--co)o(v)p 539
4882 V 36 w(f)n(orm)f Fk(<)p Fo(1/2)p Fk(>)g Fo(\(-X\))50
b Fn(Specify)34 b(the)g(form)g(of)h(the)f(co)o(v)n(ariance)g(matrix.)58
b(In)34 b(the)g(pre-)568 5002 y(vious)h(v)o(ersions)h(the)g(co)o(v)n
(ariance)g(matrix)g(for)h(po)n(wer)n(-la)o(w)f(noise)g(w)o(as)g(formed)
h(as)1918 5251 y(9)p eop
%%Page: 10 10
10 9 bop 568 548 a Fn(described)28 b(in)h Fl(W)-5 b(illiams)27
b Fn([2003])h(i.e.)43 b(the)29 b(transformation)f(matrix)g(\(which)h
(is)f(used)568 668 y(to)d(produce)h(the)g(co)o(v)n(ariance)g(matrix\))f
(is)h(scaled)g(by)f(the)h(indi)n(vidual)e Fh(\001)p Fk(T)14
b Fn(s)26 b(\()p Fh(\001)p Fk(T)14 b(j)36 b Fh(=)568
789 y Fk(T)625 804 y Fg(j)686 789 y Fm(\000)26 b Fk(T)846
804 y Fg(j)t Fj(\000)p Fi(1)972 789 y Fn(\).)42 b(F)o(or)29
b(e)o(xtremely)e(une)n(v)o(enly)g(spaced)i(data,)g(then)f(this)g(is)g
(the)g(only)g(w)o(ay)568 909 y(to)e(do)h(this,)g(b)n(ut)f(it)h(can)g
(lead)g(to)g(uncertain)g(assumptions)d(when)j(trying)f(to)h(interpret)
568 1029 y(the)21 b(results,)g(in)f(particular)h(with)g(reference)h(to)
f(an)g(e)n(v)o(enly)f(spaced)h(dataset.)30 b(Alterna-)568
1150 y(ti)n(v)o(ely)-6 b(,)21 b(one)i(can)g(create)i(a)e(co)o(v)n
(ariance)g(matrix)f(assuming)g(a)i(certain)f(time-sampling)568
1270 y(and)30 b(then)f(remo)o(v)o(e)g(the)h(columns)e(and)i(ro)n(ws)f
(associated)h(with)f(missing)f(data.)46 b(The)568 1391
y(second)31 b(method)g(is,)i(I)f(belie)n(v)o(e,)g(generally)f(more)h
(stable)f(and)g(corresponds)h(to)f(the)568 1511 y(method)24
b(used)h(by)g Fl(Langbein)f Fn([2004].)32 b(It)25 b(is)g(also)g(the)g
(only)f(real)i(w)o(ay)f(of)g(computing)568 1631 y(the)32
b(co)o(v)n(ariance)g(matrix)g(for)h(certain)f(noise)g(models)g(such)g
(as)h(band-pass,)h(autore-)568 1752 y(gressi)n(v)o(e)22
b(and)i(mo)o(ving)d(a)n(v)o(erage.)31 b(The)24 b(def)o(ault)f(is)g(1,)h
(remo)o(ving)e(columns)h(and)g(ro)n(ws)568 1872 y(of)i(missing)d(data.)
324 2167 y Ff(3.1)119 b(Example)29 b(Command)h(Line)324
2335 y Fn(T)-8 b(o)33 b(estimate)f(the)h(amount)f(of)i(\003ick)o(er)f
(noise)g(and)g(white)g(noise)f(in)h(the)g(\002le)h(vyas.neu)e(as-)324
2433 y(suming)23 b(there)i(is)g(an)f(annual)h(signal)f(in)g(the)h
(series)g(we)g(can)g(use)g(the)f(command)324 2595 y Fe(cats)39
b(vyas.neu)g(--sinusoid)f(1y)i(--model)e(pl:k-1)h(--verbose)g(--output)
g(vyas.fn_mle)470 2757 y Fn(T)-8 b(o)37 b(estimate)g(the)g(spectral)g
(inde)o(x)g(plus)f(the)h(amplitudes)f(of)h(white)g(and)h(po)n(wer)n
(-la)o(w)324 2855 y(noise)22 b(in)h(the)g(east)g(and)g(up)g(component)f
(of)h(\002le)g(vyas.neu)g(assuming)e(an)j(annual)e(and)h(semi-)324
2952 y(annual)h(signal)g(and)h(using)f(the)g(f)o(ast)h(option)f(we)h
(can)g(use)g(the)f(command)324 3114 y Fe(cats)39 b(vyas.neu)g
(--sinusoid)f(1y1)h(--model)g(pl:)h(--columns)e(3)i(--verbose)e
(--speed)h(2)h(--output)f(vyas.pl_mle)324 3454 y Fp(4)143
b(T)-11 b(ypical)35 b(Output)324 3677 y Fn(Belo)n(w)28
b(is)g(the)g(output)f(for)i(the)f(North)g(component)f(of)i(the)f
(\002le)h(vyas.neu)e(using)h(the)g(com-)324 3798 y(mand)324
3984 y Fe(cats)39 b(--model)g(pl:k-1)g(--model)g(wh:)g(--columns)g(4)g
(--sinusoid)g(1y)g(--verbose)g(--output)f(vyas.fn_mle)g(vyas.neu)470
4213 y Fn(The)27 b(long)e(list)h(of)h(numbers)e(are)j(the)e(v)o(erte)o
(x)g(in)g(the)g(stages)g(of)h(the)g(simple)o(x)d(algorithm)324
4333 y(together)g(with)g(the)h(MLE)f(v)n(alue)g(and)h(the)g(noise)f
(amplitudes.)364 4520 y Fe(Sampling)38 b(frequency)h(1.15664e-05)f
(\(Hz\),)h(1.00)g(days)364 4599 y(Number)g(of)g(samples)g(1)h(period)f
(apart)g(=)g(565)h(of)f(587)364 4678 y(Number)g(of)g(points)g(in)h
(full)f(series)79 b(=)39 b(612)364 4757 y(Cats)g(Version)g(:)g(3.1.1)
364 4835 y(Cats)g(command)g(:)g(cats)h(--model)e(pl:k-1)h(--columns)g
(4)h(--sinusoid)e(1y)i(--verbose)e(--output)h(vyas.fn_mle)f(vyas.neu)
364 4914 y(Series[0])g(=)i(1)364 4993 y(Series[1])e(=)i(0)1893
5251 y Fn(10)p eop
%%Page: 11 11
11 10 bop 364 548 a Fe(Series[2])38 b(=)i(0)364 627 y(Number)f(of)g
(series)g(to)h(process)e(:)i(1)364 706 y(Data)f(from)g(file)g(:)h
(vyas.neu)364 785 y(cats)f(:)h(running)e(on)i(bilai)364
863 y(Linux)f(release)g(2.4.21-27.0.1.EL)e(\(version)i(#1)g(Mon)h(Dec)f
(20)h(18:47:51)e(EST)i(2004\))f(on)g(i686)364 942 y(userid)g(:)g(sdwil)
364 1100 y(Start)g(Time)g(:)h(Tue)f(Jan)g(25)h(14:21:23)f(2005)364
1258 y(work\(1\))78 b(=)40 b(19992.000000)364 1337 y(info)198
b(=)40 b(0)364 1416 y(Time)f(taken)g(to)h(create)e(covariance)h(matrix)
g(and)g(compute)g(eigen)g(value)g(and)h(vectors)e(:)i(8)g(seconds)364
1494 y(wh_only)e(=)i(0.00113253)e(\(3154.2463\),)h(cn_only)f(=)i
(0.00503155)e(\(3141.9072\))364 1573 y(Starting)g(a)i(one-dimensional)e
(minimisation)g(:)h(initial)g(angle)g(45.00)364 1652
y(angle)g(=)g(45.000000)g(mle)g(=)h(3176.04935824)e(radius)h(=)199
b(1.476166)39 b(wh)g(=)199 b(1.043807)39 b(cn)h(=)199
b(1.043807)364 1731 y(Next)39 b(choice)g(of)g(angle)g(=)h(27.811530)364
1810 y(angle)f(=)g(27.811530)g(mle)g(=)h(3184.98028986)e(radius)h(=)199
b(2.038823)39 b(wh)g(=)199 b(0.951243)39 b(cn)h(=)199
b(1.803313)364 1889 y(Next)39 b(choice)g(of)g(angle)g(=)h(17.188471)364
1968 y(angle)f(=)g(17.188471)g(mle)g(=)h(3181.49799071)e(radius)h(=)199
b(2.8