added the paper on bandwidth tuning

This commit is contained in:
vasiloglou
2007-07-31 14:25:51 +00:00
parent e800d3272c
commit 408c73bf19
28 changed files with 190926 additions and 480 deletions
File diff suppressed because it is too large Load Diff
File diff suppressed because it is too large Load Diff
+373
View File
@@ -0,0 +1,373 @@
%!PS-Adobe-3.0 EPSF-3.0
%%Creator: MATLAB, The Mathworks, Inc.
%%Title: C:\Nick\MLSP2006_figures\fig1.1c.eps
%%CreationDate: 04/14/2006 21:13:49
%%DocumentNeededFonts: Helvetica
%%DocumentProcessColors: Cyan Magenta Yellow Black
%%LanguageLevel: 2
%%Pages: 1
%%BoundingBox: 63 215 550 589
%%EndComments
%%BeginProlog
% MathWorks dictionary
/MathWorks 160 dict begin
% definition operators
/bdef {bind def} bind def
/ldef {load def} bind def
/xdef {exch def} bdef
/xstore {exch store} bdef
% operator abbreviations
/c /clip ldef
/cc /concat ldef
/cp /closepath ldef
/gr /grestore ldef
/gs /gsave ldef
/mt /moveto ldef
/np /newpath ldef
/cm /currentmatrix ldef
/sm /setmatrix ldef
/rm /rmoveto ldef
/rl /rlineto ldef
/s {show newpath} bdef
/sc {setcmykcolor} bdef
/sr /setrgbcolor ldef
/sg /setgray ldef
/w /setlinewidth ldef
/j /setlinejoin ldef
/cap /setlinecap ldef
/rc {rectclip} bdef
/rf {rectfill} bdef
% page state control
/pgsv () def
/bpage {/pgsv save def} bdef
/epage {pgsv restore} bdef
/bplot /gsave ldef
/eplot {stroke grestore} bdef
% orientation switch
/portraitMode 0 def /landscapeMode 1 def /rotateMode 2 def
% coordinate system mappings
/dpi2point 0 def
% font control
/FontSize 0 def
/FMS {/FontSize xstore findfont [FontSize 0 0 FontSize neg 0 0]
makefont setfont} bdef
/ISOLatin1Encoding where {pop /WindowsLatin1Encoding 256 array bdef
ISOLatin1Encoding WindowsLatin1Encoding copy pop
/.notdef/.notdef/quotesinglbase/florin/quotedblbase/ellipsis/dagger
/daggerdbl/circumflex/perthousand/Scaron/guilsinglleft/OE/.notdef/.notdef
/.notdef/.notdef/quoteleft/quoteright/quotedblleft/quotedblright/bullet
/endash/emdash/tilde/trademark/scaron/guilsinglright/oe/.notdef/.notdef
/Ydieresis WindowsLatin1Encoding 128 32 getinterval astore pop}
{/WindowsLatin1Encoding StandardEncoding bdef} ifelse
/reencode {exch dup where {pop load} {pop StandardEncoding} ifelse
exch dup 3 1 roll findfont dup length dict begin
{ 1 index /FID ne {def}{pop pop} ifelse } forall
/Encoding exch def currentdict end definefont pop} bdef
/isroman {findfont /CharStrings get /Agrave known} bdef
/FMSR {3 1 roll 1 index dup isroman {reencode} {pop pop} ifelse
exch FMS} bdef
/csm {1 dpi2point div -1 dpi2point div scale neg translate
dup landscapeMode eq {pop -90 rotate}
{rotateMode eq {90 rotate} if} ifelse} bdef
% line types: solid, dotted, dashed, dotdash
/SO { [] 0 setdash } bdef
/DO { [.5 dpi2point mul 4 dpi2point mul] 0 setdash } bdef
/DA { [6 dpi2point mul] 0 setdash } bdef
/DD { [.5 dpi2point mul 4 dpi2point mul 6 dpi2point mul 4
dpi2point mul] 0 setdash } bdef
% macros for lines and objects
/L {lineto stroke} bdef
/MP {3 1 roll moveto 1 sub {rlineto} repeat} bdef
/AP {{rlineto} repeat} bdef
/PDlw -1 def
/W {/PDlw currentlinewidth def setlinewidth} def
/PP {closepath eofill} bdef
/DP {closepath stroke} bdef
/MR {4 -2 roll moveto dup 0 exch rlineto exch 0 rlineto
neg 0 exch rlineto closepath} bdef
/FR {MR stroke} bdef
/PR {MR fill} bdef
/L1i {{currentfile picstr readhexstring pop} image} bdef
/tMatrix matrix def
/MakeOval {newpath tMatrix currentmatrix pop translate scale
0 0 1 0 360 arc tMatrix setmatrix} bdef
/FO {MakeOval stroke} bdef
/PO {MakeOval fill} bdef
/PD {currentlinewidth 2 div 0 360 arc fill
PDlw -1 eq not {PDlw w /PDlw -1 def} if} def
/FA {newpath tMatrix currentmatrix pop translate scale
0 0 1 5 -2 roll arc tMatrix setmatrix stroke} bdef
/PA {newpath tMatrix currentmatrix pop translate 0 0 moveto scale
0 0 1 5 -2 roll arc closepath tMatrix setmatrix fill} bdef
/FAn {newpath tMatrix currentmatrix pop translate scale
0 0 1 5 -2 roll arcn tMatrix setmatrix stroke} bdef
/PAn {newpath tMatrix currentmatrix pop translate 0 0 moveto scale
0 0 1 5 -2 roll arcn closepath tMatrix setmatrix fill} bdef
/vradius 0 def /hradius 0 def /lry 0 def
/lrx 0 def /uly 0 def /ulx 0 def /rad 0 def
/MRR {/vradius xdef /hradius xdef /lry xdef /lrx xdef /uly xdef
/ulx xdef newpath tMatrix currentmatrix pop ulx hradius add uly
vradius add translate hradius vradius scale 0 0 1 180 270 arc
tMatrix setmatrix lrx hradius sub uly vradius add translate
hradius vradius scale 0 0 1 270 360 arc tMatrix setmatrix
lrx hradius sub lry vradius sub translate hradius vradius scale
0 0 1 0 90 arc tMatrix setmatrix ulx hradius add lry vradius sub
translate hradius vradius scale 0 0 1 90 180 arc tMatrix setmatrix
closepath} bdef
/FRR {MRR stroke } bdef
/PRR {MRR fill } bdef
/MlrRR {/lry xdef /lrx xdef /uly xdef /ulx xdef /rad lry uly sub 2 div def
newpath tMatrix currentmatrix pop ulx rad add uly rad add translate
rad rad scale 0 0 1 90 270 arc tMatrix setmatrix lrx rad sub lry rad
sub translate rad rad scale 0 0 1 270 90 arc tMatrix setmatrix
closepath} bdef
/FlrRR {MlrRR stroke } bdef
/PlrRR {MlrRR fill } bdef
/MtbRR {/lry xdef /lrx xdef /uly xdef /ulx xdef /rad lrx ulx sub 2 div def
newpath tMatrix currentmatrix pop ulx rad add uly rad add translate
rad rad scale 0 0 1 180 360 arc tMatrix setmatrix lrx rad sub lry rad
sub translate rad rad scale 0 0 1 0 180 arc tMatrix setmatrix
closepath} bdef
/FtbRR {MtbRR stroke } bdef
/PtbRR {MtbRR fill } bdef
/stri 6 array def /dtri 6 array def
/smat 6 array def /dmat 6 array def
/tmat1 6 array def /tmat2 6 array def /dif 3 array def
/asub {/ind2 exch def /ind1 exch def dup dup
ind1 get exch ind2 get sub exch } bdef
/tri_to_matrix {
2 0 asub 3 1 asub 4 0 asub 5 1 asub
dup 0 get exch 1 get 7 -1 roll astore } bdef
/compute_transform {
dmat dtri tri_to_matrix tmat1 invertmatrix
smat stri tri_to_matrix tmat2 concatmatrix } bdef
/ds {stri astore pop} bdef
/dt {dtri astore pop} bdef
/db {2 copy /cols xdef /rows xdef mul dup string
currentfile
3 index 0 eq {/ASCIIHexDecode filter}
{/ASCII85Decode filter 3 index 2 eq {/RunLengthDecode filter} if }
ifelse exch readstring pop
/bmap xdef pop pop} bdef
/it {gs np dtri aload pop moveto lineto lineto cp c
cols rows 8 compute_transform
{bmap} image gr}bdef
/il {newpath moveto lineto stroke}bdef
currentdict end def
%%EndProlog
%%BeginSetup
MathWorks begin
0 cap
end
%%EndSetup
%%Page: 1 1
%%BeginPageSetup
%%PageBoundingBox: 63 215 550 589
MathWorks begin
bpage
%%EndPageSetup
%%BeginObject: obj1
bplot
/dpi2point 12 def
portraitMode 0204 7344 csm
561 274 5843 4486 rc
85 dict begin %Colortable dictionary
/c0 { 0.000000 0.000000 0.000000 sr} bdef
/c1 { 1.000000 1.000000 1.000000 sr} bdef
/c2 { 0.900000 0.000000 0.000000 sr} bdef
/c3 { 0.000000 0.820000 0.000000 sr} bdef
/c4 { 0.000000 0.000000 0.800000 sr} bdef
/c5 { 0.910000 0.820000 0.320000 sr} bdef
/c6 { 1.000000 0.260000 0.820000 sr} bdef
/c7 { 0.000000 0.820000 0.820000 sr} bdef
c0
1 j
1 sg
0 0 6913 5185 rf
6 w
0 4225 5356 0 0 -4225 899 4614 4 MP
PP
-5356 0 0 4225 5356 0 0 -4225 899 4614 5 MP stroke
4 w
DO
SO
6 w
0 sg
899 389 mt 6255 389 L
899 4614 mt 6255 4614 L
6255 4614 mt 6255 389 L
899 4614 mt 899 389 L
899 4614 mt 6255 4614 L
899 4614 mt 899 389 L
899 4614 mt 899 4560 L
899 389 mt 899 442 L
%%IncludeResource: font Helvetica
/Helvetica /WindowsLatin1Encoding 120 FMSR
866 4759 mt
(0) s
1791 4614 mt 1791 4560 L
1791 389 mt 1791 442 L
1691 4759 mt
(100) s
2684 4614 mt 2684 4560 L
2684 389 mt 2684 442 L
2584 4759 mt
(200) s
3577 4614 mt 3577 4560 L
3577 389 mt 3577 442 L
3477 4759 mt
(300) s
4469 4614 mt 4469 4560 L
4469 389 mt 4469 442 L
4369 4759 mt
(400) s
5362 4614 mt 5362 4560 L
5362 389 mt 5362 442 L
5262 4759 mt
(500) s
6255 4614 mt 6255 4560 L
6255 389 mt 6255 442 L
6155 4759 mt
(600) s
899 4614 mt 952 4614 L
6255 4614 mt 6201 4614 L
561 4658 mt
(-0.08) s
899 4085 mt 952 4085 L
6255 4085 mt 6201 4085 L
561 4129 mt
(-0.06) s
899 3557 mt 952 3557 L
6255 3557 mt 6201 3557 L
561 3601 mt
(-0.04) s
899 3029 mt 952 3029 L
6255 3029 mt 6201 3029 L
561 3073 mt
(-0.02) s
899 2501 mt 952 2501 L
6255 2501 mt 6201 2501 L
798 2545 mt
(0) s
899 1973 mt 952 1973 L
6255 1973 mt 6201 1973 L
631 2017 mt
(0.02) s
899 1445 mt 952 1445 L
6255 1445 mt 6201 1445 L
631 1489 mt
(0.04) s
899 917 mt 952 917 L
6255 917 mt 6201 917 L
631 961 mt
(0.06) s
899 389 mt 952 389 L
6255 389 mt 6201 389 L
631 433 mt
(0.08) s
899 389 mt 6255 389 L
899 4614 mt 6255 4614 L
6255 4614 mt 6255 389 L
899 4614 mt 899 389 L
gs 899 389 5357 4226 rc
9 -653 9 436 9 -298 9 -233 9 -48 9 1122 9 -675 9 558
9 103 9 -746 9 -397 9 -35 9 1195 8 -1245 9 894 9 481
9 -1165 9 192 9 -138 9 39 9 80 9 527 9 -560 9 -517
9 964 9 -120 9 -408 8 656 9 -105 9 92 9 401 9 -900
9 820 9 -996 9 -57 9 0 9 -232 9 142 9 -86 9 -77
8 405 9 -185 9 1270 9 -1070 9 944 9 -1259 9 -313 9 944
9 -266 9 389 9 455 9 -514 9 618 9 -631 8 -38 9 -656
9 -64 9 330 9 94 9 -347 9 1238 9 -1303 9 -11 9 165
9 1274 9 -1400 9 306 9 679 8 -1021 9 9 9 99 9 1339
9 57 9 -1227 9 -185 9 552 9 324 9 -319 9 673 9 -817
9 -335 8 1336 9 -436 9 -680 9 147 9 -332 9 799 9 514
9 -1822 9 1234 9 -183 9 175 9 -306 9 154 9 -633 8 463
9 62 9 782 9 26 9 32 9 -144 9 -206 9 -774 9 -273
9 -200 9 683 9 -27 9 -299 9 1202 8 35 9 -413 9 -1148
9 1446 9 -1199 9 823 9 84 9 -17 9 -134 9 -45 9 182
9 333 9 -1478 8 1174 9 -1037 9 154 9 -66 9 805 9 -934
9 792 9 -37 9 -845 9 598 9 114 9 -558 9 860 9 -805
8 466 9 -805 9 950 9 -787 9 669 9 -498 9 1267 9 -261
9 -1243 9 639 9 563 9 371 9 -571 8 242 9 -617 9 201
9 649 9 -1347 9 114 9 680 9 129 9 -522 9 251 9 793
9 -68 9 -1344 9 1342 8 -887 9 739 9 -579 9 -650 9 -123
9 1157 9 -302 9 -898 9 188 9 199 9 814 9 -619 9 736
9 -862 8 1121 9 -278 9 -989 9 611 9 -462 9 600 9 237
9 -1029 9 1313 9 13 9 -501 9 -969 9 22 8 1165 9 -972
9 1288 9 -1386 9 1112 9 -929 9 361 9 -176 9 -245 9 1066
9 -1306 9 667 9 -376 9 239 8 67 9 -281 9 153 9 -178
9 136 9 -334 9 1397 9 -637 9 640 9 -417 9 -1099 9 304
9 -145 9 -31 8 88 9 123 9 -295 9 -178 9 792 9 720
9 177 9 -1167 9 546 9 -46 9 -884 9 834 9 -265 8 -378
9 -221 9 1286 9 -1095 9 20 9 1113 9 -176 9 -856 9 -267
9 1375 9 -875 9 648 9 -593 9 -41 8 99 9 910 9 -1259
9 -32 9 131 9 172 9 -235 9 -42 9 1215 9 -1468 9 935
9 400 9 -679 9 -935 8 544 9 760 9 -845 9 1332 9 -645
9 516 9 -1445 9 209 9 729 9 566 9 -538 9 -109 9 -889
8 1520 9 -1575 9 205 9 969 9 -533 9 -637 9 1291 9 -681
9 -6 9 474 9 365 9 -1354 9 30 9 841 8 605 9 -1269
9 1347 9 -308 9 -596 9 671 9 -792 9 38 9 580 9 -1013
9 302 9 168 9 -427 3585 2056 300 MP stroke
8 -1224 9 51 9 -713 9 354 9 -291 9 385 9 -27 9 -395
9 -137 9 802 9 -444 9 -344 9 471 9 -389 8 323 9 507
9 297 9 -231 9 40 9 11 9 -718 9 663 9 -676 9 1142
9 -6 9 -1038 9 -389 9 231 8 -190 9 348 9 492 9 -948
9 1291 9 -160 9 -994 9 262 9 -280 9 647 9 -734 9 27
9 -100 8 611 9 -149 9 901 9 -129 9 -86 9 -101 9 109
9 -398 9 563 9 -456 9 -163 9 -739 9 1075 9 -932 8 990
9 -751 9 234 9 -49 9 473 9 -933 9 640 9 -460 9 242
9 258 9 568 9 -1168 9 786 9 24 8 -978 9 1121 9 331
9 -231 9 -1054 9 249 9 643 9 -1018 9 -12 9 1539 9 -1333
9 571 9 124 8 -832 9 500 9 -515 9 -18 9 926 9 -103
9 80 9 -101 9 -218 9 100 9 479 9 -1071 9 55 9 -209
8 1007 9 -47 9 -825 9 1231 9 -250 9 -881 9 831 9 -303
9 -22 9 -613 9 1203 9 -720 9 660 9 -273 8 -820 9 221
9 641 9 -810 9 -253 9 1269 9 -1160 9 793 9 38 9 -663
9 1101 9 -113 9 -1223 8 -22 9 478 9 -493 9 249 9 -251
9 1441 9 -1213 9 -116 9 1448 9 -1001 9 284 9 252 9 -1039
9 1402 8 -762 9 437 9 -604 9 -507 9 353 9 514 9 295
9 -671 9 464 9 -444 9 541 9 186 9 -775 8 1031 9 -772
9 277 9 131 9 -726 9 72 9 668 9 -356 9 206 9 -949
9 3 9 86 9 188 9 -53 8 215 9 589 9 -809 9 481
9 -555 9 65 9 553 9 18 9 264 9 -265 9 -15 9 54
9 -240 9 -93 8 943 9 -863 9 -330 9 615 9 -525 9 -277
9 86 9 277 9 63 9 -13 9 -265 9 950 9 -402 8 452
9 -53 9 -24 9 -630 9 -262 9 278 9 860 9 -1408 9 1021
9 -974 9 605 9 332 9 -806 9 1206 8 -1090 9 1056 9 -1083
9 704 9 -140 9 429 9 -531 9 -584 9 1053 9 -506 9 303
9 -537 9 131 9 439 8 -654 9 507 9 -89 9 112 9 -815
9 1069 9 -315 9 367 9 -1192 9 1119 9 -443 9 339 9 -843
8 738 9 185 9 -1136 9 1253 9 -250 9 177 9 82 9 -567
9 564 9 111 9 -1324 9 775 9 -463 9 821 8 -281 9 148
9 -79 9 -175 9 178 9 -54 9 -792 9 1118 9 -423 9 190
9 8 9 -498 9 -147 9 807 8 84 9 124 9 -105 9 -220
9 54 9 -636 9 -289 9 1020 9 -1075 9 239 9 858 9 -874
9 342 8 -603 9 1018 9 218 9 -576 9 121 9 -711 9 804
9 -686 9 131 9 270 9 -470 9 837 9 -22 9 120 8 -723
9 -201 9 1044 9 -1093 9 126 9 -136 9 433 9 369 9 -391
9 -349 9 475 9 107 916 3360 300 MP stroke
9 323 907 3037 2 MP stroke
gr
end %%Color Dict
eplot
%%EndObject
epage
end
showpage
%%Trailer
%%EOF
File diff suppressed because it is too large Load Diff
File diff suppressed because it is too large Load Diff
+381
View File
@@ -0,0 +1,381 @@
%!PS-Adobe-3.0 EPSF-3.0
%%Creator: MATLAB, The Mathworks, Inc.
%%Title: C:\Nick\MLSP2006_figures\fig1c.eps
%%CreationDate: 04/14/2006 21:07:13
%%DocumentNeededFonts: Helvetica
%%DocumentProcessColors: Cyan Magenta Yellow Black
%%LanguageLevel: 2
%%Pages: 1
%%BoundingBox: 63 215 550 589
%%EndComments
%%BeginProlog
% MathWorks dictionary
/MathWorks 160 dict begin
% definition operators
/bdef {bind def} bind def
/ldef {load def} bind def
/xdef {exch def} bdef
/xstore {exch store} bdef
% operator abbreviations
/c /clip ldef
/cc /concat ldef
/cp /closepath ldef
/gr /grestore ldef
/gs /gsave ldef
/mt /moveto ldef
/np /newpath ldef
/cm /currentmatrix ldef
/sm /setmatrix ldef
/rm /rmoveto ldef
/rl /rlineto ldef
/s {show newpath} bdef
/sc {setcmykcolor} bdef
/sr /setrgbcolor ldef
/sg /setgray ldef
/w /setlinewidth ldef
/j /setlinejoin ldef
/cap /setlinecap ldef
/rc {rectclip} bdef
/rf {rectfill} bdef
% page state control
/pgsv () def
/bpage {/pgsv save def} bdef
/epage {pgsv restore} bdef
/bplot /gsave ldef
/eplot {stroke grestore} bdef
% orientation switch
/portraitMode 0 def /landscapeMode 1 def /rotateMode 2 def
% coordinate system mappings
/dpi2point 0 def
% font control
/FontSize 0 def
/FMS {/FontSize xstore findfont [FontSize 0 0 FontSize neg 0 0]
makefont setfont} bdef
/ISOLatin1Encoding where {pop /WindowsLatin1Encoding 256 array bdef
ISOLatin1Encoding WindowsLatin1Encoding copy pop
/.notdef/.notdef/quotesinglbase/florin/quotedblbase/ellipsis/dagger
/daggerdbl/circumflex/perthousand/Scaron/guilsinglleft/OE/.notdef/.notdef
/.notdef/.notdef/quoteleft/quoteright/quotedblleft/quotedblright/bullet
/endash/emdash/tilde/trademark/scaron/guilsinglright/oe/.notdef/.notdef
/Ydieresis WindowsLatin1Encoding 128 32 getinterval astore pop}
{/WindowsLatin1Encoding StandardEncoding bdef} ifelse
/reencode {exch dup where {pop load} {pop StandardEncoding} ifelse
exch dup 3 1 roll findfont dup length dict begin
{ 1 index /FID ne {def}{pop pop} ifelse } forall
/Encoding exch def currentdict end definefont pop} bdef
/isroman {findfont /CharStrings get /Agrave known} bdef
/FMSR {3 1 roll 1 index dup isroman {reencode} {pop pop} ifelse
exch FMS} bdef
/csm {1 dpi2point div -1 dpi2point div scale neg translate
dup landscapeMode eq {pop -90 rotate}
{rotateMode eq {90 rotate} if} ifelse} bdef
% line types: solid, dotted, dashed, dotdash
/SO { [] 0 setdash } bdef
/DO { [.5 dpi2point mul 4 dpi2point mul] 0 setdash } bdef
/DA { [6 dpi2point mul] 0 setdash } bdef
/DD { [.5 dpi2point mul 4 dpi2point mul 6 dpi2point mul 4
dpi2point mul] 0 setdash } bdef
% macros for lines and objects
/L {lineto stroke} bdef
/MP {3 1 roll moveto 1 sub {rlineto} repeat} bdef
/AP {{rlineto} repeat} bdef
/PDlw -1 def
/W {/PDlw currentlinewidth def setlinewidth} def
/PP {closepath eofill} bdef
/DP {closepath stroke} bdef
/MR {4 -2 roll moveto dup 0 exch rlineto exch 0 rlineto
neg 0 exch rlineto closepath} bdef
/FR {MR stroke} bdef
/PR {MR fill} bdef
/L1i {{currentfile picstr readhexstring pop} image} bdef
/tMatrix matrix def
/MakeOval {newpath tMatrix currentmatrix pop translate scale
0 0 1 0 360 arc tMatrix setmatrix} bdef
/FO {MakeOval stroke} bdef
/PO {MakeOval fill} bdef
/PD {currentlinewidth 2 div 0 360 arc fill
PDlw -1 eq not {PDlw w /PDlw -1 def} if} def
/FA {newpath tMatrix currentmatrix pop translate scale
0 0 1 5 -2 roll arc tMatrix setmatrix stroke} bdef
/PA {newpath tMatrix currentmatrix pop translate 0 0 moveto scale
0 0 1 5 -2 roll arc closepath tMatrix setmatrix fill} bdef
/FAn {newpath tMatrix currentmatrix pop translate scale
0 0 1 5 -2 roll arcn tMatrix setmatrix stroke} bdef
/PAn {newpath tMatrix currentmatrix pop translate 0 0 moveto scale
0 0 1 5 -2 roll arcn closepath tMatrix setmatrix fill} bdef
/vradius 0 def /hradius 0 def /lry 0 def
/lrx 0 def /uly 0 def /ulx 0 def /rad 0 def
/MRR {/vradius xdef /hradius xdef /lry xdef /lrx xdef /uly xdef
/ulx xdef newpath tMatrix currentmatrix pop ulx hradius add uly
vradius add translate hradius vradius scale 0 0 1 180 270 arc
tMatrix setmatrix lrx hradius sub uly vradius add translate
hradius vradius scale 0 0 1 270 360 arc tMatrix setmatrix
lrx hradius sub lry vradius sub translate hradius vradius scale
0 0 1 0 90 arc tMatrix setmatrix ulx hradius add lry vradius sub
translate hradius vradius scale 0 0 1 90 180 arc tMatrix setmatrix
closepath} bdef
/FRR {MRR stroke } bdef
/PRR {MRR fill } bdef
/MlrRR {/lry xdef /lrx xdef /uly xdef /ulx xdef /rad lry uly sub 2 div def
newpath tMatrix currentmatrix pop ulx rad add uly rad add translate
rad rad scale 0 0 1 90 270 arc tMatrix setmatrix lrx rad sub lry rad
sub translate rad rad scale 0 0 1 270 90 arc tMatrix setmatrix
closepath} bdef
/FlrRR {MlrRR stroke } bdef
/PlrRR {MlrRR fill } bdef
/MtbRR {/lry xdef /lrx xdef /uly xdef /ulx xdef /rad lrx ulx sub 2 div def
newpath tMatrix currentmatrix pop ulx rad add uly rad add translate
rad rad scale 0 0 1 180 360 arc tMatrix setmatrix lrx rad sub lry rad
sub translate rad rad scale 0 0 1 0 180 arc tMatrix setmatrix
closepath} bdef
/FtbRR {MtbRR stroke } bdef
/PtbRR {MtbRR fill } bdef
/stri 6 array def /dtri 6 array def
/smat 6 array def /dmat 6 array def
/tmat1 6 array def /tmat2 6 array def /dif 3 array def
/asub {/ind2 exch def /ind1 exch def dup dup
ind1 get exch ind2 get sub exch } bdef
/tri_to_matrix {
2 0 asub 3 1 asub 4 0 asub 5 1 asub
dup 0 get exch 1 get 7 -1 roll astore } bdef
/compute_transform {
dmat dtri tri_to_matrix tmat1 invertmatrix
smat stri tri_to_matrix tmat2 concatmatrix } bdef
/ds {stri astore pop} bdef
/dt {dtri astore pop} bdef
/db {2 copy /cols xdef /rows xdef mul dup string
currentfile
3 index 0 eq {/ASCIIHexDecode filter}
{/ASCII85Decode filter 3 index 2 eq {/RunLengthDecode filter} if }
ifelse exch readstring pop
/bmap xdef pop pop} bdef
/it {gs np dtri aload pop moveto lineto lineto cp c
cols rows 8 compute_transform
{bmap} image gr}bdef
/il {newpath moveto lineto stroke}bdef
currentdict end def
%%EndProlog
%%BeginSetup
MathWorks begin
0 cap
end
%%EndSetup
%%Page: 1 1
%%BeginPageSetup
%%PageBoundingBox: 63 215 550 589
MathWorks begin
bpage
%%EndPageSetup
%%BeginObject: obj1
bplot
/dpi2point 12 def
portraitMode 0204 7344 csm
561 274 5843 4486 rc
85 dict begin %Colortable dictionary
/c0 { 0.000000 0.000000 0.000000 sr} bdef
/c1 { 1.000000 1.000000 1.000000 sr} bdef
/c2 { 0.900000 0.000000 0.000000 sr} bdef
/c3 { 0.000000 0.820000 0.000000 sr} bdef
/c4 { 0.000000 0.000000 0.800000 sr} bdef
/c5 { 0.910000 0.820000 0.320000 sr} bdef
/c6 { 1.000000 0.260000 0.820000 sr} bdef
/c7 { 0.000000 0.820000 0.820000 sr} bdef
c0
1 j
1 sg
0 0 6913 5185 rf
6 w
0 4225 5356 0 0 -4225 899 4614 4 MP
PP
-5356 0 0 4225 5356 0 0 -4225 899 4614 5 MP stroke
4 w
DO
SO
6 w
0 sg
899 389 mt 6255 389 L
899 4614 mt 6255 4614 L
6255 4614 mt 6255 389 L
899 4614 mt 899 389 L
899 4614 mt 6255 4614 L
899 4614 mt 899 389 L
899 4614 mt 899 4560 L
899 389 mt 899 442 L
%%IncludeResource: font Helvetica
/Helvetica /WindowsLatin1Encoding 120 FMSR
866 4759 mt
(0) s
1791 4614 mt 1791 4560 L
1791 389 mt 1791 442 L
1691 4759 mt
(100) s
2684 4614 mt 2684 4560 L
2684 389 mt 2684 442 L
2584 4759 mt
(200) s
3577 4614 mt 3577 4560 L
3577 389 mt 3577 442 L
3477 4759 mt
(300) s
4469 4614 mt 4469 4560 L
4469 389 mt 4469 442 L
4369 4759 mt
(400) s
5362 4614 mt 5362 4560 L
5362 389 mt 5362 442 L
5262 4759 mt
(500) s
6255 4614 mt 6255 4560 L
6255 389 mt 6255 442 L
6155 4759 mt
(600) s
899 4614 mt 952 4614 L
6255 4614 mt 6201 4614 L
561 4658 mt
(-0.05) s
899 4191 mt 952 4191 L
6255 4191 mt 6201 4191 L
561 4235 mt
(-0.04) s
899 3769 mt 952 3769 L
6255 3769 mt 6201 3769 L
561 3813 mt
(-0.03) s
899 3346 mt 952 3346 L
6255 3346 mt 6201 3346 L
561 3390 mt
(-0.02) s
899 2924 mt 952 2924 L
6255 2924 mt 6201 2924 L
561 2968 mt
(-0.01) s
899 2501 mt 952 2501 L
6255 2501 mt 6201 2501 L
798 2545 mt
(0) s
899 2079 mt 952 2079 L
6255 2079 mt 6201 2079 L
631 2123 mt
(0.01) s
899 1656 mt 952 1656 L
6255 1656 mt 6201 1656 L
631 1700 mt
(0.02) s
899 1234 mt 952 1234 L
6255 1234 mt 6201 1234 L
631 1278 mt
(0.03) s
899 811 mt 952 811 L
6255 811 mt 6201 811 L
631 855 mt
(0.04) s
899 389 mt 952 389 L
6255 389 mt 6201 389 L
631 433 mt
(0.05) s
899 389 mt 6255 389 L
899 4614 mt 6255 4614 L
6255 4614 mt 6255 389 L
899 4614 mt 899 389 L
gs 899 389 5357 4226 rc
9 56 9 -62 9 33 9 -40 9 124 9 -156 9 71 9 -49
9 0 9 1 9 81 9 56 9 -143 8 152 9 -137 9 -31
9 139 9 -40 9 17 9 -31 9 6 9 -62 9 49 9 143
9 -193 9 23 9 26 8 -33 9 -15 9 1 9 -20 9 62
9 -64 9 68 9 -50 9 62 9 74 9 -78 9 23 9 50
8 -85 9 -48 9 -20 9 21 9 -26 9 107 9 374 9 -434
9 38 9 -55 9 -18 9 26 9 -38 9 23 8 1 9 142
9 -4 9 -114 9 39 9 -13 9 -59 9 149 9 3 9 -61
9 -101 9 111 9 5 9 -91 8 105 9 -11 9 -98 9 -21
9 0 9 117 9 -96 9 42 9 -45 9 22 9 -45 9 28
9 98 8 -122 9 27 9 -5 9 48 9 25 9 -71 9 -24
9 809 9 -765 9 -9 9 -4 9 34 9 -19 9 104 8 -104
9 4 9 -48 9 -3 9 1 9 6 9 11 9 52 9 90
9 -14 9 -89 9 -35 9 75 9 -99 8 4 9 20 9 145
9 -170 9 137 9 -107 9 -8 9 3 9 8 9 -3 9 -9
9 -4 9 150 8 -147 9 98 9 -90 9 -2 9 2 9 139
9 -120 9 -12 9 136 9 -144 9 28 9 95 9 -118 9 55
8 -64 9 164 9 -136 9 131 9 -131 9 23 9 -73 9 19
9 186 9 -182 9 -4 9 -19 9 29 8 -17 9 50 9 -17
9 -27 9 91 9 -13 9 -63 9 12 9 32 9 -24 9 -53
9 5 9 155 9 -160 8 65 9 -43 9 22 9 122 9 2
9 -145 9 14 9 132 9 -2 9 -45 9 -91 9 -6 9 -5
9 5 8 -24 9 20 9 4 9 11 9 70 9 -74 9 -4
9 104 9 -132 9 11 9 17 9 111 9 -17 8 -103 9 106
9 -121 9 113 9 -108 9 125 9 -109 9 53 9 68 9 -129
9 159 9 -120 9 23 9 -4 8 -9 9 37 9 -78 9 90
9 -88 9 134 9 -158 9 25 9 -27 9 28 9 104 9 -58
9 -48 9 113 8 -115 9 11 9 104 9 139 9 -223 9 -50
9 -3 9 100 9 -75 9 -7 9 159 9 -151 9 40 8 56
9 63 9 -164 9 132 9 -21 9 -126 9 13 9 108 9 31
9 -147 9 11 9 -2 9 2 9 -3 8 5 9 -26 9 105
9 5 9 -4 9 -85 9 50 9 54 9 -122 9 158 9 -115
9 -35 9 13 9 947 8 -835 9 -117 9 127 9 -140 9 40
9 -27 9 146 9 -85 9 -51 9 -30 9 35 9 -14 9 202
8 -224 9 248 9 -81 9 -135 9 46 9 216 9 -264 9 53
9 -59 9 10 9 -10 9 159 9 -15 9 -139 8 -23 9 68
9 -73 9 28 9 39 9 -48 9 8 9 51 9 -55 9 118
9 -118 9 40 9 40 3585 4218 300 MP stroke
8 3438 9 -2 9 42 9 -54 9 17 9 -5 9 2 9 40
9 -35 9 -6 9 52 9 -9 9 -42 9 41 8 -40 9 -6
9 -60 9 -54 9 120 9 8 9 -16 9 16 9 30 9 -513
9 111 9 348 9 55 9 12 8 -48 9 -16 9 -56 9 111
9 -250 9 153 9 56 9 -13 9 16 9 -35 9 71 9 1
9 -35 8 -20 9 0 9 -215 9 174 9 48 9 13 9 -10
9 3 9 -42 9 47 9 -24 9 71 9 -42 9 34 8 -99
9 92 9 -7 9 -36 9 -103 9 115 9 -7 9 5 9 -3
9 -3 9 -71 9 78 9 -5 9 -8 8 15 9 -53 9 -236
9 230 9 92 9 -44 9 -107 9 163 9 -6 9 -628 9 585
9 -54 9 14 8 42 9 -17 9 54 9 10 9 -103 9 60
9 1 9 -44 9 59 9 -48 9 2 9 62 9 -37 9 3
8 -2 9 -7 9 7 9 -76 9 19 9 92 9 -31 9 -56
9 45 9 43 9 -218 9 177 9 -94 9 90 8 52 9 -47
9 7 9 -7 9 51 9 -119 9 119 9 -41 9 -9 9 28
9 -124 9 43 9 103 8 -2 9 0 9 1 9 -44 9 45
9 -297 9 286 9 9 9 -497 9 431 9 -15 9 -17 9 55
9 -314 8 331 9 -119 9 81 9 66 9 -11 9 -48 9 -114
9 157 9 -70 9 23 9 -37 9 -3 9 49 8 -328 9 304
9 -19 9 52 9 -13 9 -6 9 -100 9 106 9 16 9 35
9 7 9 -48 9 -10 9 55 8 -60 9 20 9 -5 9 -26
9 72 9 -47 9 9 9 -40 9 -84 9 83 9 38 9 -12
9 32 9 -45 8 -333 9 331 9 58 9 -39 9 -10 9 53
9 -42 9 32 9 -33 9 7 9 32 9 -88 9 -3 8 -93
9 13 9 5 9 130 9 0 9 -2 9 -261 9 266 9 3
9 -2 9 -30 9 -15 9 78 9 -113 8 120 9 -116 9 69
9 -49 9 8 9 -18 9 31 9 31 9 -134 9 130 9 9
9 41 9 -64 9 -98 8 159 9 -38 9 -13 9 4 9 52
9 -43 9 -35 9 -110 9 188 9 -42 9 -28 9 -5 9 63
8 -48 9 -52 9 112 9 -195 9 104 9 -10 9 38 9 -2
9 -68 9 -130 9 263 9 -84 9 32 9 -21 8 20 9 -109
9 115 9 -5 9 -18 9 -54 9 113 9 -179 9 149 9 -45
9 -13 9 53 9 -4 9 -47 8 7 9 -33 9 9 9 62
9 9 9 24 9 19 9 -62 9 15 9 34 9 -79 9 66
9 -53 8 77 9 -170 9 -50 9 197 9 -60 9 38 9 -31
9 29 9 42 9 -68 9 28 9 -90 9 82 9 -38 8 41
9 7 9 -132 9 166 9 -37 9 30 9 -34 9 -3 9 2
9 6 9 -26 9 -9 916 810 300 MP stroke
9 39 907 771 2 MP stroke
gr
end %%Color Dict
eplot
%%EndObject
epage
end
showpage
%%Trailer
%%EOF
File diff suppressed because it is too large Load Diff
File diff suppressed because it is too large Load Diff
+394
View File
@@ -0,0 +1,394 @@
%!PS-Adobe-3.0 EPSF-3.0
%%Creator: MATLAB, The Mathworks, Inc.
%%Title: C:\Nick\MLSP2006_figures\fig2c.eps
%%CreationDate: 04/14/2006 19:11:41
%%DocumentNeededFonts: Helvetica
%%DocumentProcessColors: Cyan Magenta Yellow Black
%%LanguageLevel: 2
%%Pages: 1
%%BoundingBox: 69 215 550 589
%%EndComments
%%BeginProlog
% MathWorks dictionary
/MathWorks 160 dict begin
% definition operators
/bdef {bind def} bind def
/ldef {load def} bind def
/xdef {exch def} bdef
/xstore {exch store} bdef
% operator abbreviations
/c /clip ldef
/cc /concat ldef
/cp /closepath ldef
/gr /grestore ldef
/gs /gsave ldef
/mt /moveto ldef
/np /newpath ldef
/cm /currentmatrix ldef
/sm /setmatrix ldef
/rm /rmoveto ldef
/rl /rlineto ldef
/s {show newpath} bdef
/sc {setcmykcolor} bdef
/sr /setrgbcolor ldef
/sg /setgray ldef
/w /setlinewidth ldef
/j /setlinejoin ldef
/cap /setlinecap ldef
/rc {rectclip} bdef
/rf {rectfill} bdef
% page state control
/pgsv () def
/bpage {/pgsv save def} bdef
/epage {pgsv restore} bdef
/bplot /gsave ldef
/eplot {stroke grestore} bdef
% orientation switch
/portraitMode 0 def /landscapeMode 1 def /rotateMode 2 def
% coordinate system mappings
/dpi2point 0 def
% font control
/FontSize 0 def
/FMS {/FontSize xstore findfont [FontSize 0 0 FontSize neg 0 0]
makefont setfont} bdef
/ISOLatin1Encoding where {pop /WindowsLatin1Encoding 256 array bdef
ISOLatin1Encoding WindowsLatin1Encoding copy pop
/.notdef/.notdef/quotesinglbase/florin/quotedblbase/ellipsis/dagger
/daggerdbl/circumflex/perthousand/Scaron/guilsinglleft/OE/.notdef/.notdef
/.notdef/.notdef/quoteleft/quoteright/quotedblleft/quotedblright/bullet
/endash/emdash/tilde/trademark/scaron/guilsinglright/oe/.notdef/.notdef
/Ydieresis WindowsLatin1Encoding 128 32 getinterval astore pop}
{/WindowsLatin1Encoding StandardEncoding bdef} ifelse
/reencode {exch dup where {pop load} {pop StandardEncoding} ifelse
exch dup 3 1 roll findfont dup length dict begin
{ 1 index /FID ne {def}{pop pop} ifelse } forall
/Encoding exch def currentdict end definefont pop} bdef
/isroman {findfont /CharStrings get /Agrave known} bdef
/FMSR {3 1 roll 1 index dup isroman {reencode} {pop pop} ifelse
exch FMS} bdef
/csm {1 dpi2point div -1 dpi2point div scale neg translate
dup landscapeMode eq {pop -90 rotate}
{rotateMode eq {90 rotate} if} ifelse} bdef
% line types: solid, dotted, dashed, dotdash
/SO { [] 0 setdash } bdef
/DO { [.5 dpi2point mul 4 dpi2point mul] 0 setdash } bdef
/DA { [6 dpi2point mul] 0 setdash } bdef
/DD { [.5 dpi2point mul 4 dpi2point mul 6 dpi2point mul 4
dpi2point mul] 0 setdash } bdef
% macros for lines and objects
/L {lineto stroke} bdef
/MP {3 1 roll moveto 1 sub {rlineto} repeat} bdef
/AP {{rlineto} repeat} bdef
/PDlw -1 def
/W {/PDlw currentlinewidth def setlinewidth} def
/PP {closepath eofill} bdef
/DP {closepath stroke} bdef
/MR {4 -2 roll moveto dup 0 exch rlineto exch 0 rlineto
neg 0 exch rlineto closepath} bdef
/FR {MR stroke} bdef
/PR {MR fill} bdef
/L1i {{currentfile picstr readhexstring pop} image} bdef
/tMatrix matrix def
/MakeOval {newpath tMatrix currentmatrix pop translate scale
0 0 1 0 360 arc tMatrix setmatrix} bdef
/FO {MakeOval stroke} bdef
/PO {MakeOval fill} bdef
/PD {currentlinewidth 2 div 0 360 arc fill
PDlw -1 eq not {PDlw w /PDlw -1 def} if} def
/FA {newpath tMatrix currentmatrix pop translate scale
0 0 1 5 -2 roll arc tMatrix setmatrix stroke} bdef
/PA {newpath tMatrix currentmatrix pop translate 0 0 moveto scale
0 0 1 5 -2 roll arc closepath tMatrix setmatrix fill} bdef
/FAn {newpath tMatrix currentmatrix pop translate scale
0 0 1 5 -2 roll arcn tMatrix setmatrix stroke} bdef
/PAn {newpath tMatrix currentmatrix pop translate 0 0 moveto scale
0 0 1 5 -2 roll arcn closepath tMatrix setmatrix fill} bdef
/vradius 0 def /hradius 0 def /lry 0 def
/lrx 0 def /uly 0 def /ulx 0 def /rad 0 def
/MRR {/vradius xdef /hradius xdef /lry xdef /lrx xdef /uly xdef
/ulx xdef newpath tMatrix currentmatrix pop ulx hradius add uly
vradius add translate hradius vradius scale 0 0 1 180 270 arc
tMatrix setmatrix lrx hradius sub uly vradius add translate
hradius vradius scale 0 0 1 270 360 arc tMatrix setmatrix
lrx hradius sub lry vradius sub translate hradius vradius scale
0 0 1 0 90 arc tMatrix setmatrix ulx hradius add lry vradius sub
translate hradius vradius scale 0 0 1 90 180 arc tMatrix setmatrix
closepath} bdef
/FRR {MRR stroke } bdef
/PRR {MRR fill } bdef
/MlrRR {/lry xdef /lrx xdef /uly xdef /ulx xdef /rad lry uly sub 2 div def
newpath tMatrix currentmatrix pop ulx rad add uly rad add translate
rad rad scale 0 0 1 90 270 arc tMatrix setmatrix lrx rad sub lry rad
sub translate rad rad scale 0 0 1 270 90 arc tMatrix setmatrix
closepath} bdef
/FlrRR {MlrRR stroke } bdef
/PlrRR {MlrRR fill } bdef
/MtbRR {/lry xdef /lrx xdef /uly xdef /ulx xdef /rad lrx ulx sub 2 div def
newpath tMatrix currentmatrix pop ulx rad add uly rad add translate
rad rad scale 0 0 1 180 360 arc tMatrix setmatrix lrx rad sub lry rad
sub translate rad rad scale 0 0 1 0 180 arc tMatrix setmatrix
closepath} bdef
/FtbRR {MtbRR stroke } bdef
/PtbRR {MtbRR fill } bdef
/stri 6 array def /dtri 6 array def
/smat 6 array def /dmat 6 array def
/tmat1 6 array def /tmat2 6 array def /dif 3 array def
/asub {/ind2 exch def /ind1 exch def dup dup
ind1 get exch ind2 get sub exch } bdef
/tri_to_matrix {
2 0 asub 3 1 asub 4 0 asub 5 1 asub
dup 0 get exch 1 get 7 -1 roll astore } bdef
/compute_transform {
dmat dtri tri_to_matrix tmat1 invertmatrix
smat stri tri_to_matrix tmat2 concatmatrix } bdef
/ds {stri astore pop} bdef
/dt {dtri astore pop} bdef
/db {2 copy /cols xdef /rows xdef mul dup string
currentfile
3 index 0 eq {/ASCIIHexDecode filter}
{/ASCII85Decode filter 3 index 2 eq {/RunLengthDecode filter} if }
ifelse exch readstring pop
/bmap xdef pop pop} bdef
/it {gs np dtri aload pop moveto lineto lineto cp c
cols rows 8 compute_transform
{bmap} image gr}bdef
/il {newpath moveto lineto stroke}bdef
currentdict end def
%%EndProlog
%%BeginSetup
MathWorks begin
0 cap
end
%%EndSetup
%%Page: 1 1
%%BeginPageSetup
%%PageBoundingBox: 69 215 550 589
MathWorks begin
bpage
%%EndPageSetup
%%BeginObject: obj1
bplot
/dpi2point 12 def
portraitMode 0204 7344 csm
628 274 5776 4486 rc
85 dict begin %Colortable dictionary
/c0 { 0.000000 0.000000 0.000000 sr} bdef
/c1 { 1.000000 1.000000 1.000000 sr} bdef
/c2 { 0.900000 0.000000 0.000000 sr} bdef
/c3 { 0.000000 0.820000 0.000000 sr} bdef
/c4 { 0.000000 0.000000 0.800000 sr} bdef
/c5 { 0.910000 0.820000 0.320000 sr} bdef
/c6 { 1.000000 0.260000 0.820000 sr} bdef
/c7 { 0.000000 0.820000 0.820000 sr} bdef
c0
1 j
1 sg
0 0 6913 5185 rf
6 w
0 4225 5356 0 0 -4225 899 4614 4 MP
PP
-5356 0 0 4225 5356 0 0 -4225 899 4614 5 MP stroke
4 w
DO
SO
6 w
0 sg
899 389 mt 6255 389 L
899 4614 mt 6255 4614 L
6255 4614 mt 6255 389 L
899 4614 mt 899 389 L
899 4614 mt 6255 4614 L
899 4614 mt 899 389 L
1337 4614 mt 1337 4560 L
1337 389 mt 1337 442 L
%%IncludeResource: font Helvetica
/Helvetica /WindowsLatin1Encoding 120 FMSR
1271 4759 mt
(50) s
1784 4614 mt 1784 4560 L
1784 389 mt 1784 442 L
1684 4759 mt
(100) s
2231 4614 mt 2231 4560 L
2231 389 mt 2231 442 L
2131 4759 mt
(150) s
2678 4614 mt 2678 4560 L
2678 389 mt 2678 442 L
2578 4759 mt
(200) s
3125 4614 mt 3125 4560 L
3125 389 mt 3125 442 L
3025 4759 mt
(250) s
3572 4614 mt 3572 4560 L
3572 389 mt 3572 442 L
3472 4759 mt
(300) s
4019 4614 mt 4019 4560 L
4019 389 mt 4019 442 L
3919 4759 mt
(350) s
4466 4614 mt 4466 4560 L
4466 389 mt 4466 442 L
4366 4759 mt
(400) s
4913 4614 mt 4913 4560 L
4913 389 mt 4913 442 L
4813 4759 mt
(450) s
5360 4614 mt 5360 4560 L
5360 389 mt 5360 442 L
5260 4759 mt
(500) s
5807 4614 mt 5807 4560 L
5807 389 mt 5807 442 L
5707 4759 mt
(550) s
6255 4614 mt 6255 4560 L
6255 389 mt 6255 442 L
6155 4759 mt
(600) s
899 4614 mt 952 4614 L
6255 4614 mt 6201 4614 L
628 4658 mt
(-0.4) s
899 4144 mt 952 4144 L
6255 4144 mt 6201 4144 L
628 4188 mt
(-0.2) s
899 3675 mt 952 3675 L
6255 3675 mt 6201 3675 L
798 3719 mt
(0) s
899 3205 mt 952 3205 L
6255 3205 mt 6201 3205 L
698 3249 mt
(0.2) s
899 2736 mt 952 2736 L
6255 2736 mt 6201 2736 L
698 2780 mt
(0.4) s
899 2266 mt 952 2266 L
6255 2266 mt 6201 2266 L
698 2310 mt
(0.6) s
899 1797 mt 952 1797 L
6255 1797 mt 6201 1797 L
698 1841 mt
(0.8) s
899 1327 mt 952 1327 L
6255 1327 mt 6201 1327 L
798 1371 mt
(1) s
899 858 mt 952 858 L
6255 858 mt 6201 858 L
698 902 mt
(1.2) s
899 389 mt 952 389 L
6255 389 mt 6201 389 L
698 433 mt
(1.4) s
899 389 mt 6255 389 L
899 4614 mt 6255 4614 L
6255 4614 mt 6255 389 L
899 4614 mt 899 389 L
gs 899 389 5357 4226 rc
9 0 9 0 9 0 9 0 9 0 9 0 9 0 9 0
9 0 9 0 9 0 9 0 9 0 9 0 8 0 9 0
9 0 9 0 9 0 9 0 9 0 9 0 9 0 9 0
9 0 9 0 9 0 9 0 9 0 9 0 9 0 8 0
9 0 9 0 9 0 9 0 9 0 9 0 9 0 9 0
9 0 9 0 9 0 9 0 9 0 9 0 9 0 9 0
8 0 9 0 9 0 9 0 9 0 9 0 9 0 9 0
9 0 9 0 9 0 9 0 9 0 9 0 9 0 9 0
9 0 8 0 9 0 9 0 9 0 9 0 9 0 9 0
9 0 9 0 9 0 9 0 9 0 9 0 9 0 9 0
9 0 9 0 8 0 9 0 9 0 9 0 9 0 9 0
9 0 9 0 9 0 9 0 9 0 9 0 9 0 9 0
9 0 9 0 9 0 8 0 9 0 9 0 9 0 9 0
9 0 9 0 9 0 9 0 9 0 9 0 9 0 9 0
9 0 9 0 9 0 9 0 8 0 9 0 9 0 9 0
9 0 9 0 9 0 9 0 9 0 9 0 9 0 9 0
9 0 9 0 9 0 9 0 9 0 9 0 8 0 9 0
9 0 9 0 9 0 9 0 9 0 9 0 9 0 9 0
9 0 9 0 9 0 9 0 9 0 9 0 9 0 8 0
9 0 9 0 9 0 9 0 9 0 9 0 9 0 9 0
9 0 9 0 9 0 9 0 9 0 9 0 9 0 9 0
8 0 9 0 9 0 9 0 9 0 9 0 9 0 9 0
9 0 9 0 9 0 9 0 9 0 9 0 9 0 9 0
9 0 8 0 9 0 9 0 9 0 9 0 9 0 9 0
9 0 9 0 9 0 9 0 9 0 9 0 9 0 9 0
9 0 9 0 8 0 9 0 9 0 9 0 9 0 9 0
9 0 9 0 9 0 9 0 9 0 9 0 9 0 9 0
9 0 9 0 9 0 8 0 9 0 9 0 9 0 9 0
9 0 9 0 9 0 9 0 9 0 9 0 9 0 9 0
9 0 9 0 9 0 9 0 8 0 9 0 9 0 9 0
9 0 9 0 9 0 9 0 9 0 9 0 9 0 9 0
9 0 9 0 9 0 9 0 9 0 8 0 9 0 9 0
9 0 9 0 9 0 9 0 9 0 9 0 9 0 9 0
9 0 9 0 9 0 9 0 9 0 9 0 8 0 9 0
9 0 9 0 9 0 9 0 9 0 9 0 9 0 9 0
9 0 9 0 9 0 9 0 9 0 9 0 9 0 9 0
8 0 9 2348 9 0 9 0 9 0 9 0 9 0 9 0
9 0 9 0 9 0 3402 1327 300 MP stroke
9 0 9 0 9 0 9 0 9 0 9 0 8 0 9 0
9 0 9 0 9 0 9 0 9 0 9 0 9 0 9 0
9 0 9 0 9 0 9 0 9 0 9 0 9 0 8 0
9 0 9 0 9 0 9 0 9 0 9 0 9 0 9 0
9 0 9 0 9 0 9 0 9 0 9 0 9 0 9 0
8 0 9 0 9 0 9 0 9 0 9 0 9 0 9 0
9 0 9 0 9 0 9 0 9 0 9 0 9 0 9 0
9 0 8 0 9 0 9 0 9 0 9 0 9 0 9 0
9 0 9 0 9 0 9 0 9 0 9 0 9 0 9 0
9 0 9 0 8 0 9 0 9 0 9 0 9 0 9 0
9 0 9 0 9 0 9 0 9 0 9 0 9 0 9 0
9 0 9 0 9 0 8 0 9 0 9 0 9 0 9 0
9 0 9 0 9 0 9 0 9 0 9 0 9 0 9 0
9 0 9 0 9 0 9 0 8 0 9 0 9 0 9 0
9 0 9 0 9 0 9 0 9 0 9 0 9 0 9 0
9 0 9 0 9 0 9 0 9 0 8 0 9 0 9 0
9 0 9 0 9 0 9 0 9 0 9 0 9 0 9 0
9 0 9 0 9 0 9 0 9 0 9 0 9 0 8 0
9 0 9 0 9 0 9 0 9 0 9 0 9 0 9 0
9 0 9 0 9 0 9 0 9 0 9 0 9 0 9 0
8 0 9 0 9 0 9 0 9 0 9 0 9 0 9 0
9 0 9 0 9 0 9 0 9 0 9 0 9 0 9 0
9 0 8 0 9 0 9 0 9 0 9 0 9 0 9 0
9 0 9 0 9 0 9 0 9 0 9 0 9 0 9 0
9 0 9 0 8 0 9 0 9 0 9 0 9 0 9 0
9 0 9 0 9 0 9 0 9 0 9 0 9 0 9 0
9 0 9 0 9 0 8 0 9 0 9 0 9 0 9 0
9 0 9 0 9 0 9 0 9 0 9 0 9 0 9 0
9 0 9 0 9 0 9 0 8 0 9 0 9 0 9 0
9 0 9 0 9 0 9 0 9 0 9 0 9 0 9 0
9 0 9 0 9 0 9 0 9 0 8 0 9 0 9 0
9 0 9 0 9 0 9 0 9 0 9 0 9 0 9 0
9 0 9 0 9 0 9 0 9 0 9 0 8 0 9 0
9 0 9 0 9 0 9 0 9 0 9 0 9 0 9 0
9 0 9 0 9 0 9 0 9 0 9 0 9 0 8 0
899 1327 281 MP stroke
gr
end %%Color Dict
eplot
%%EndObject
epage
end
showpage
%%Trailer
%%EOF
File diff suppressed because it is too large Load Diff
File diff suppressed because it is too large Load Diff
+394
View File
@@ -0,0 +1,394 @@
%!PS-Adobe-3.0 EPSF-3.0
%%Creator: MATLAB, The Mathworks, Inc.
%%Title: C:\Nick\MLSP2006_figures\fig3c.eps
%%CreationDate: 04/14/2006 19:17:24
%%DocumentNeededFonts: Helvetica
%%DocumentProcessColors: Cyan Magenta Yellow Black
%%LanguageLevel: 2
%%Pages: 1
%%BoundingBox: 69 215 550 589
%%EndComments
%%BeginProlog
% MathWorks dictionary
/MathWorks 160 dict begin
% definition operators
/bdef {bind def} bind def
/ldef {load def} bind def
/xdef {exch def} bdef
/xstore {exch store} bdef
% operator abbreviations
/c /clip ldef
/cc /concat ldef
/cp /closepath ldef
/gr /grestore ldef
/gs /gsave ldef
/mt /moveto ldef
/np /newpath ldef
/cm /currentmatrix ldef
/sm /setmatrix ldef
/rm /rmoveto ldef
/rl /rlineto ldef
/s {show newpath} bdef
/sc {setcmykcolor} bdef
/sr /setrgbcolor ldef
/sg /setgray ldef
/w /setlinewidth ldef
/j /setlinejoin ldef
/cap /setlinecap ldef
/rc {rectclip} bdef
/rf {rectfill} bdef
% page state control
/pgsv () def
/bpage {/pgsv save def} bdef
/epage {pgsv restore} bdef
/bplot /gsave ldef
/eplot {stroke grestore} bdef
% orientation switch
/portraitMode 0 def /landscapeMode 1 def /rotateMode 2 def
% coordinate system mappings
/dpi2point 0 def
% font control
/FontSize 0 def
/FMS {/FontSize xstore findfont [FontSize 0 0 FontSize neg 0 0]
makefont setfont} bdef
/ISOLatin1Encoding where {pop /WindowsLatin1Encoding 256 array bdef
ISOLatin1Encoding WindowsLatin1Encoding copy pop
/.notdef/.notdef/quotesinglbase/florin/quotedblbase/ellipsis/dagger
/daggerdbl/circumflex/perthousand/Scaron/guilsinglleft/OE/.notdef/.notdef
/.notdef/.notdef/quoteleft/quoteright/quotedblleft/quotedblright/bullet
/endash/emdash/tilde/trademark/scaron/guilsinglright/oe/.notdef/.notdef
/Ydieresis WindowsLatin1Encoding 128 32 getinterval astore pop}
{/WindowsLatin1Encoding StandardEncoding bdef} ifelse
/reencode {exch dup where {pop load} {pop StandardEncoding} ifelse
exch dup 3 1 roll findfont dup length dict begin
{ 1 index /FID ne {def}{pop pop} ifelse } forall
/Encoding exch def currentdict end definefont pop} bdef
/isroman {findfont /CharStrings get /Agrave known} bdef
/FMSR {3 1 roll 1 index dup isroman {reencode} {pop pop} ifelse
exch FMS} bdef
/csm {1 dpi2point div -1 dpi2point div scale neg translate
dup landscapeMode eq {pop -90 rotate}
{rotateMode eq {90 rotate} if} ifelse} bdef
% line types: solid, dotted, dashed, dotdash
/SO { [] 0 setdash } bdef
/DO { [.5 dpi2point mul 4 dpi2point mul] 0 setdash } bdef
/DA { [6 dpi2point mul] 0 setdash } bdef
/DD { [.5 dpi2point mul 4 dpi2point mul 6 dpi2point mul 4
dpi2point mul] 0 setdash } bdef
% macros for lines and objects
/L {lineto stroke} bdef
/MP {3 1 roll moveto 1 sub {rlineto} repeat} bdef
/AP {{rlineto} repeat} bdef
/PDlw -1 def
/W {/PDlw currentlinewidth def setlinewidth} def
/PP {closepath eofill} bdef
/DP {closepath stroke} bdef
/MR {4 -2 roll moveto dup 0 exch rlineto exch 0 rlineto
neg 0 exch rlineto closepath} bdef
/FR {MR stroke} bdef
/PR {MR fill} bdef
/L1i {{currentfile picstr readhexstring pop} image} bdef
/tMatrix matrix def
/MakeOval {newpath tMatrix currentmatrix pop translate scale
0 0 1 0 360 arc tMatrix setmatrix} bdef
/FO {MakeOval stroke} bdef
/PO {MakeOval fill} bdef
/PD {currentlinewidth 2 div 0 360 arc fill
PDlw -1 eq not {PDlw w /PDlw -1 def} if} def
/FA {newpath tMatrix currentmatrix pop translate scale
0 0 1 5 -2 roll arc tMatrix setmatrix stroke} bdef
/PA {newpath tMatrix currentmatrix pop translate 0 0 moveto scale
0 0 1 5 -2 roll arc closepath tMatrix setmatrix fill} bdef
/FAn {newpath tMatrix currentmatrix pop translate scale
0 0 1 5 -2 roll arcn tMatrix setmatrix stroke} bdef
/PAn {newpath tMatrix currentmatrix pop translate 0 0 moveto scale
0 0 1 5 -2 roll arcn closepath tMatrix setmatrix fill} bdef
/vradius 0 def /hradius 0 def /lry 0 def
/lrx 0 def /uly 0 def /ulx 0 def /rad 0 def
/MRR {/vradius xdef /hradius xdef /lry xdef /lrx xdef /uly xdef
/ulx xdef newpath tMatrix currentmatrix pop ulx hradius add uly
vradius add translate hradius vradius scale 0 0 1 180 270 arc
tMatrix setmatrix lrx hradius sub uly vradius add translate
hradius vradius scale 0 0 1 270 360 arc tMatrix setmatrix
lrx hradius sub lry vradius sub translate hradius vradius scale
0 0 1 0 90 arc tMatrix setmatrix ulx hradius add lry vradius sub
translate hradius vradius scale 0 0 1 90 180 arc tMatrix setmatrix
closepath} bdef
/FRR {MRR stroke } bdef
/PRR {MRR fill } bdef
/MlrRR {/lry xdef /lrx xdef /uly xdef /ulx xdef /rad lry uly sub 2 div def
newpath tMatrix currentmatrix pop ulx rad add uly rad add translate
rad rad scale 0 0 1 90 270 arc tMatrix setmatrix lrx rad sub lry rad
sub translate rad rad scale 0 0 1 270 90 arc tMatrix setmatrix
closepath} bdef
/FlrRR {MlrRR stroke } bdef
/PlrRR {MlrRR fill } bdef
/MtbRR {/lry xdef /lrx xdef /uly xdef /ulx xdef /rad lrx ulx sub 2 div def
newpath tMatrix currentmatrix pop ulx rad add uly rad add translate
rad rad scale 0 0 1 180 360 arc tMatrix setmatrix lrx rad sub lry rad
sub translate rad rad scale 0 0 1 0 180 arc tMatrix setmatrix
closepath} bdef
/FtbRR {MtbRR stroke } bdef
/PtbRR {MtbRR fill } bdef
/stri 6 array def /dtri 6 array def
/smat 6 array def /dmat 6 array def
/tmat1 6 array def /tmat2 6 array def /dif 3 array def
/asub {/ind2 exch def /ind1 exch def dup dup
ind1 get exch ind2 get sub exch } bdef
/tri_to_matrix {
2 0 asub 3 1 asub 4 0 asub 5 1 asub
dup 0 get exch 1 get 7 -1 roll astore } bdef
/compute_transform {
dmat dtri tri_to_matrix tmat1 invertmatrix
smat stri tri_to_matrix tmat2 concatmatrix } bdef
/ds {stri astore pop} bdef
/dt {dtri astore pop} bdef
/db {2 copy /cols xdef /rows xdef mul dup string
currentfile
3 index 0 eq {/ASCIIHexDecode filter}
{/ASCII85Decode filter 3 index 2 eq {/RunLengthDecode filter} if }
ifelse exch readstring pop
/bmap xdef pop pop} bdef
/it {gs np dtri aload pop moveto lineto lineto cp c
cols rows 8 compute_transform
{bmap} image gr}bdef
/il {newpath moveto lineto stroke}bdef
currentdict end def
%%EndProlog
%%BeginSetup
MathWorks begin
0 cap
end
%%EndSetup
%%Page: 1 1
%%BeginPageSetup
%%PageBoundingBox: 69 215 550 589
MathWorks begin
bpage
%%EndPageSetup
%%BeginObject: obj1
bplot
/dpi2point 12 def
portraitMode 0204 7344 csm
628 274 5776 4486 rc
85 dict begin %Colortable dictionary
/c0 { 0.000000 0.000000 0.000000 sr} bdef
/c1 { 1.000000 1.000000 1.000000 sr} bdef
/c2 { 0.900000 0.000000 0.000000 sr} bdef
/c3 { 0.000000 0.820000 0.000000 sr} bdef
/c4 { 0.000000 0.000000 0.800000 sr} bdef
/c5 { 0.910000 0.820000 0.320000 sr} bdef
/c6 { 1.000000 0.260000 0.820000 sr} bdef
/c7 { 0.000000 0.820000 0.820000 sr} bdef
c0
1 j
1 sg
0 0 6913 5185 rf
6 w
0 4225 5356 0 0 -4225 899 4614 4 MP
PP
-5356 0 0 4225 5356 0 0 -4225 899 4614 5 MP stroke
4 w
DO
SO
6 w
0 sg
899 389 mt 6255 389 L
899 4614 mt 6255 4614 L
6255 4614 mt 6255 389 L
899 4614 mt 899 389 L
899 4614 mt 6255 4614 L
899 4614 mt 899 389 L
1337 4614 mt 1337 4560 L
1337 389 mt 1337 442 L
%%IncludeResource: font Helvetica
/Helvetica /WindowsLatin1Encoding 120 FMSR
1271 4759 mt
(50) s
1784 4614 mt 1784 4560 L
1784 389 mt 1784 442 L
1684 4759 mt
(100) s
2231 4614 mt 2231 4560 L
2231 389 mt 2231 442 L
2131 4759 mt
(150) s
2678 4614 mt 2678 4560 L
2678 389 mt 2678 442 L
2578 4759 mt
(200) s
3125 4614 mt 3125 4560 L
3125 389 mt 3125 442 L
3025 4759 mt
(250) s
3572 4614 mt 3572 4560 L
3572 389 mt 3572 442 L
3472 4759 mt
(300) s
4019 4614 mt 4019 4560 L
4019 389 mt 4019 442 L
3919 4759 mt
(350) s
4466 4614 mt 4466 4560 L
4466 389 mt 4466 442 L
4366 4759 mt
(400) s
4913 4614 mt 4913 4560 L
4913 389 mt 4913 442 L
4813 4759 mt
(450) s
5360 4614 mt 5360 4560 L
5360 389 mt 5360 442 L
5260 4759 mt
(500) s
5807 4614 mt 5807 4560 L
5807 389 mt 5807 442 L
5707 4759 mt
(550) s
6255 4614 mt 6255 4560 L
6255 389 mt 6255 442 L
6155 4759 mt
(600) s
899 4614 mt 952 4614 L
6255 4614 mt 6201 4614 L
628 4658 mt
(-0.4) s
899 4144 mt 952 4144 L
6255 4144 mt 6201 4144 L
628 4188 mt
(-0.2) s
899 3675 mt 952 3675 L
6255 3675 mt 6201 3675 L
798 3719 mt
(0) s
899 3205 mt 952 3205 L
6255 3205 mt 6201 3205 L
698 3249 mt
(0.2) s
899 2736 mt 952 2736 L
6255 2736 mt 6201 2736 L
698 2780 mt
(0.4) s
899 2266 mt 952 2266 L
6255 2266 mt 6201 2266 L
698 2310 mt
(0.6) s
899 1797 mt 952 1797 L
6255 1797 mt 6201 1797 L
698 1841 mt
(0.8) s
899 1327 mt 952 1327 L
6255 1327 mt 6201 1327 L
798 1371 mt
(1) s
899 858 mt 952 858 L
6255 858 mt 6201 858 L
698 902 mt
(1.2) s
899 389 mt 952 389 L
6255 389 mt 6201 389 L
698 433 mt
(1.4) s
899 389 mt 6255 389 L
899 4614 mt 6255 4614 L
6255 4614 mt 6255 389 L
899 4614 mt 899 389 L
gs 899 389 5357 4226 rc
9 -49 9 -807 9 825 9 -817 9 808 9 -831 9 -51 9 629
9 242 9 -44 9 58 9 -15 9 37 9 -408 8 405 9 -33
9 12 9 15 9 -10 9 -677 9 -268 9 374 9 549 9 -338
9 345 9 30 9 19 9 -763 9 721 9 -11 9 29 8 -28
9 30 9 -46 9 93 9 -123 9 125 9 -72 9 2 9 1
9 19 9 -1051 9 26 9 1032 9 -840 9 -57 9 857 9 19
8 -521 9 541 9 2 9 18 9 -4 9 -31 9 52 9 -56
9 -910 9 684 9 226 9 45 9 -26 9 -29 9 -811 9 806
9 7 8 -902 9 944 9 -882 9 828 9 7 9 4 9 -339
9 343 9 -4 9 -225 9 181 9 58 9 -936 9 915 9 10
9 4 9 -767 8 229 9 184 9 -412 9 759 9 -4 9 -1
9 -294 9 -685 9 1014 9 -621 9 60 9 513 9 24 9 -5
9 -71 9 2 9 63 8 -2 9 -5 9 -559 9 580 9 -35
9 -338 9 366 9 14 9 -7 9 -7 9 -939 9 550 9 151
9 241 9 -693 9 -43 9 727 8 -4 9 4 9 55 9 -66
9 18 9 -628 9 -53 9 125 9 551 9 2 9 10 9 53
9 -553 9 564 9 -70 9 4 9 -855 9 773 8 -761 9 820
9 5 9 -959 9 450 9 -473 9 987 9 10 9 -2 9 -1018
9 82 9 712 9 -215 9 -126 9 -306 9 -37 9 765 8 -794
9 258 9 -137 9 -122 9 728 9 -750 9 31 9 177 9 548
9 -171 9 -141 9 -381 9 957 9 -992 9 300 9 705 9 -1104
8 854 9 42 9 -7 9 221 9 -688 9 233 9 -573 9 550
9 -451 9 451 9 183 9 10 9 -183 9 131 9 -746 9 581
9 -236 8 -308 9 -49 9 752 9 -613 9 941 9 -534 9 -446
9 855 9 -10 9 -274 9 48 9 236 9 -533 9 77 9 -450
9 667 9 196 8 -647 9 416 9 -657 9 335 9 -153 9 -129
9 -62 9 168 9 -176 9 93 9 927 9 -868 9 794 9 -238
9 211 9 -230 9 -369 8 238 9 400 9 -157 9 -177 9 -404
9 236 9 167 9 -134 9 -291 9 467 9 -40 9 -548 9 56
9 684 9 186 9 -913 9 263 8 -137 9 -221 9 14 9 98
9 -69 9 886 9 -895 9 897 9 -900 9 412 9 327 9 305
9 -1082 9 690 9 -405 9 57 9 525 8 -770 9 904 9 -291
9 -386 9 -122 9 -97 9 255 9 584 9 -552 9 -60 9 200
9 -502 9 37 9 -34 9 759 9 -689 9 587 8 155 9 -826
9 187 9 449 9 -362 9 623 9 -92 9 42 9 -814 9 512
9 382 9 -673 9 -223 9 111 9 202 9 -266 9 636 9 -125
8 -389 9 1302 9 421 9 -604 9 -81 9 -37 9 1087 9 -1049
9 -102 9 204 9 -108 3402 2412 300 MP stroke
9 32 9 -12 9 690 9 -756 9 52 9 260 8 -285 9 551
9 460 9 -1007 9 40 9 -35 9 274 9 -263 9 70 9 870
9 -968 9 135 9 923 9 -316 9 -712 9 -86 9 115 8 -50
9 19 9 27 9 978 9 -1035 9 69 9 -96 9 89 9 48
9 60 9 -168 9 56 9 -3 9 779 9 -748 9 -149 9 112
8 294 9 -122 9 -181 9 1056 9 -1058 9 -88 9 88 9 -25
9 -14 9 -18 9 49 9 -23 9 -82 9 54 9 52 9 1083
9 -1045 8 -3 9 1012 9 9 9 -589 9 -495 9 50 9 -27
9 -34 9 23 9 59 9 -22 9 -62 9 283 9 -345 9 96
9 1002 9 -966 8 -69 9 86 9 -104 9 59 9 -64 9 13
9 10 9 1115 9 -809 9 -256 9 355 9 -436 9 1150 9 -1052
9 -73 9 27 9 -19 8 -51 9 481 9 -404 9 36 9 -130
9 1110 9 -381 9 -669 9 -21 9 58 9 27 9 -62 9 926
9 -940 9 7 9 92 9 410 8 -158 9 -323 9 295 9 -252
9 -29 9 14 9 -51 9 53 9 -81 9 105 9 91 9 -173
9 367 9 -277 9 -9 9 -114 9 331 8 -239 9 200 9 311
9 -488 9 1035 9 -1070 9 34 9 1047 9 -1163 9 89 9 97
9 866 9 -528 9 272 9 -511 9 -41 9 824 9 -632 8 574
9 -133 9 -540 9 205 9 192 9 -310 9 663 9 -48 9 -758
9 855 9 -552 9 312 9 -505 9 776 9 -99 9 -674 9 664
8 -54 9 -434 9 332 9 -71 9 108 9 -209 9 319 9 95
9 32 9 -411 9 -93 9 422 9 -543 9 442 9 -465 9 617
9 -419 8 425 9 -654 9 222 9 -444 9 623 9 124 9 -823
9 769 9 -577 9 281 9 445 9 -20 9 52 9 -352 9 -577
9 111 9 172 8 533 9 -617 9 -98 9 425 9 -567 9 450
9 125 9 -350 9 124 9 235 9 -378 9 22 9 755 9 -934
9 772 9 154 9 -802 8 227 9 -207 9 794 9 -176 9 41
9 -351 9 412 9 -429 9 -267 9 696 9 -646 9 505 9 199
9 -481 9 -305 9 776 9 -24 8 2 9 -302 9 -401 9 -116
9 782 9 -13 9 -828 9 130 9 385 9 342 9 -237 9 86
9 -142 9 24 9 -63 9 -468 9 427 8 345 9 -5 9 -20
9 82 9 -538 9 357 9 70 9 -254 9 4 9 -116 9 498
9 -590 9 -2 9 -298 9 824 9 -614 9 587 8 94 9 -672
9 475 9 -323 9 453 9 21 9 74 9 -631 9 525 9 36
9 -406 9 -166 9 454 9 -65 9 201 9 24 9 -718 8 435
899 1638 281 MP stroke
gr
end %%Color Dict
eplot
%%EndObject
epage
end
showpage
%%Trailer
%%EOF
File diff suppressed because it is too large Load Diff
File diff suppressed because it is too large Load Diff
File diff suppressed because it is too large Load Diff
File diff suppressed because it is too large Load Diff
File diff suppressed because it is too large Load Diff
+426
View File
@@ -0,0 +1,426 @@
%!PS-Adobe-3.0 EPSF-3.0
%%Creator: MATLAB, The Mathworks, Inc.
%%Title: C:\Nick\MLSP2006_figures\fig4f.eps
%%CreationDate: 04/11/2006 20:53:56
%%DocumentNeededFonts: Helvetica
%%DocumentProcessColors: Cyan Magenta Yellow Black
%%LanguageLevel: 2
%%Pages: 1
%%BoundingBox: 72 215 550 583
%%EndComments
%%BeginProlog
% MathWorks dictionary
/MathWorks 160 dict begin
% definition operators
/bdef {bind def} bind def
/ldef {load def} bind def
/xdef {exch def} bdef
/xstore {exch store} bdef
% operator abbreviations
/c /clip ldef
/cc /concat ldef
/cp /closepath ldef
/gr /grestore ldef
/gs /gsave ldef
/mt /moveto ldef
/np /newpath ldef
/cm /currentmatrix ldef
/sm /setmatrix ldef
/rm /rmoveto ldef
/rl /rlineto ldef
/s {show newpath} bdef
/sc {setcmykcolor} bdef
/sr /setrgbcolor ldef
/sg /setgray ldef
/w /setlinewidth ldef
/j /setlinejoin ldef
/cap /setlinecap ldef
/rc {rectclip} bdef
/rf {rectfill} bdef
% page state control
/pgsv () def
/bpage {/pgsv save def} bdef
/epage {pgsv restore} bdef
/bplot /gsave ldef
/eplot {stroke grestore} bdef
% orientation switch
/portraitMode 0 def /landscapeMode 1 def /rotateMode 2 def
% coordinate system mappings
/dpi2point 0 def
% font control
/FontSize 0 def
/FMS {/FontSize xstore findfont [FontSize 0 0 FontSize neg 0 0]
makefont setfont} bdef
/ISOLatin1Encoding where {pop /WindowsLatin1Encoding 256 array bdef
ISOLatin1Encoding WindowsLatin1Encoding copy pop
/.notdef/.notdef/quotesinglbase/florin/quotedblbase/ellipsis/dagger
/daggerdbl/circumflex/perthousand/Scaron/guilsinglleft/OE/.notdef/.notdef
/.notdef/.notdef/quoteleft/quoteright/quotedblleft/quotedblright/bullet
/endash/emdash/tilde/trademark/scaron/guilsinglright/oe/.notdef/.notdef
/Ydieresis WindowsLatin1Encoding 128 32 getinterval astore pop}
{/WindowsLatin1Encoding StandardEncoding bdef} ifelse
/reencode {exch dup where {pop load} {pop StandardEncoding} ifelse
exch dup 3 1 roll findfont dup length dict begin
{ 1 index /FID ne {def}{pop pop} ifelse } forall
/Encoding exch def currentdict end definefont pop} bdef
/isroman {findfont /CharStrings get /Agrave known} bdef
/FMSR {3 1 roll 1 index dup isroman {reencode} {pop pop} ifelse
exch FMS} bdef
/csm {1 dpi2point div -1 dpi2point div scale neg translate
dup landscapeMode eq {pop -90 rotate}
{rotateMode eq {90 rotate} if} ifelse} bdef
% line types: solid, dotted, dashed, dotdash
/SO { [] 0 setdash } bdef
/DO { [.5 dpi2point mul 4 dpi2point mul] 0 setdash } bdef
/DA { [6 dpi2point mul] 0 setdash } bdef
/DD { [.5 dpi2point mul 4 dpi2point mul 6 dpi2point mul 4
dpi2point mul] 0 setdash } bdef
% macros for lines and objects
/L {lineto stroke} bdef
/MP {3 1 roll moveto 1 sub {rlineto} repeat} bdef
/AP {{rlineto} repeat} bdef
/PDlw -1 def
/W {/PDlw currentlinewidth def setlinewidth} def
/PP {closepath eofill} bdef
/DP {closepath stroke} bdef
/MR {4 -2 roll moveto dup 0 exch rlineto exch 0 rlineto
neg 0 exch rlineto closepath} bdef
/FR {MR stroke} bdef
/PR {MR fill} bdef
/L1i {{currentfile picstr readhexstring pop} image} bdef
/tMatrix matrix def
/MakeOval {newpath tMatrix currentmatrix pop translate scale
0 0 1 0 360 arc tMatrix setmatrix} bdef
/FO {MakeOval stroke} bdef
/PO {MakeOval fill} bdef
/PD {currentlinewidth 2 div 0 360 arc fill
PDlw -1 eq not {PDlw w /PDlw -1 def} if} def
/FA {newpath tMatrix currentmatrix pop translate scale
0 0 1 5 -2 roll arc tMatrix setmatrix stroke} bdef
/PA {newpath tMatrix currentmatrix pop translate 0 0 moveto scale
0 0 1 5 -2 roll arc closepath tMatrix setmatrix fill} bdef
/FAn {newpath tMatrix currentmatrix pop translate scale
0 0 1 5 -2 roll arcn tMatrix setmatrix stroke} bdef
/PAn {newpath tMatrix currentmatrix pop translate 0 0 moveto scale
0 0 1 5 -2 roll arcn closepath tMatrix setmatrix fill} bdef
/vradius 0 def /hradius 0 def /lry 0 def
/lrx 0 def /uly 0 def /ulx 0 def /rad 0 def
/MRR {/vradius xdef /hradius xdef /lry xdef /lrx xdef /uly xdef
/ulx xdef newpath tMatrix currentmatrix pop ulx hradius add uly
vradius add translate hradius vradius scale 0 0 1 180 270 arc
tMatrix setmatrix lrx hradius sub uly vradius add translate
hradius vradius scale 0 0 1 270 360 arc tMatrix setmatrix
lrx hradius sub lry vradius sub translate hradius vradius scale
0 0 1 0 90 arc tMatrix setmatrix ulx hradius add lry vradius sub
translate hradius vradius scale 0 0 1 90 180 arc tMatrix setmatrix
closepath} bdef
/FRR {MRR stroke } bdef
/PRR {MRR fill } bdef
/MlrRR {/lry xdef /lrx xdef /uly xdef /ulx xdef /rad lry uly sub 2 div def
newpath tMatrix currentmatrix pop ulx rad add uly rad add translate
rad rad scale 0 0 1 90 270 arc tMatrix setmatrix lrx rad sub lry rad
sub translate rad rad scale 0 0 1 270 90 arc tMatrix setmatrix
closepath} bdef
/FlrRR {MlrRR stroke } bdef
/PlrRR {MlrRR fill } bdef
/MtbRR {/lry xdef /lrx xdef /uly xdef /ulx xdef /rad lrx ulx sub 2 div def
newpath tMatrix currentmatrix pop ulx rad add uly rad add translate
rad rad scale 0 0 1 180 360 arc tMatrix setmatrix lrx rad sub lry rad
sub translate rad rad scale 0 0 1 0 180 arc tMatrix setmatrix
closepath} bdef
/FtbRR {MtbRR stroke } bdef
/PtbRR {MtbRR fill } bdef
/stri 6 array def /dtri 6 array def
/smat 6 array def /dmat 6 array def
/tmat1 6 array def /tmat2 6 array def /dif 3 array def
/asub {/ind2 exch def /ind1 exch def dup dup
ind1 get exch ind2 get sub exch } bdef
/tri_to_matrix {
2 0 asub 3 1 asub 4 0 asub 5 1 asub
dup 0 get exch 1 get 7 -1 roll astore } bdef
/compute_transform {
dmat dtri tri_to_matrix tmat1 invertmatrix
smat stri tri_to_matrix tmat2 concatmatrix } bdef
/ds {stri astore pop} bdef
/dt {dtri astore pop} bdef
/db {2 copy /cols xdef /rows xdef mul dup string
currentfile
3 index 0 eq {/ASCIIHexDecode filter}
{/ASCII85Decode filter 3 index 2 eq {/RunLengthDecode filter} if }
ifelse exch readstring pop
/bmap xdef pop pop} bdef
/it {gs np dtri aload pop moveto lineto lineto cp c
cols rows 8 compute_transform
{bmap} image gr}bdef
/il {newpath moveto lineto stroke}bdef
currentdict end def
%%EndProlog
%%BeginSetup
MathWorks begin
0 cap
end
%%EndSetup
%%Page: 1 1
%%BeginPageSetup
%%PageBoundingBox: 72 215 550 583
MathWorks begin
bpage
%%EndPageSetup
%%BeginObject: obj1
bplot
/dpi2point 12 def
portraitMode 0204 7344 csm
664 341 5735 4419 rc
85 dict begin %Colortable dictionary
/c0 { 0.000000 0.000000 0.000000 sr} bdef
/c1 { 1.000000 1.000000 1.000000 sr} bdef
/c2 { 0.900000 0.000000 0.000000 sr} bdef
/c3 { 0.000000 0.820000 0.000000 sr} bdef
/c4 { 0.000000 0.000000 0.800000 sr} bdef
/c5 { 0.910000 0.820000 0.320000 sr} bdef
/c6 { 1.000000 0.260000 0.820000 sr} bdef
/c7 { 0.000000 0.820000 0.820000 sr} bdef
c0
1 j
1 sg
0 0 6913 5185 rf
6 w
0 -4225 5356 0 0 4225 899 389 4 MP
PP
-5356 0 0 -4225 5356 0 0 4225 899 389 5 MP stroke
gs 899 389 5357 4226 rc
/mwscm { [/Indexed /DeviceRGB 63
<
000000 040404 080808 0c0c0c 101010 141414 181818 1c1c1c 202020 242424
282828 2c2c2c 303030 343434 383838 3c3c3c 404040 444444 484848 4c4c4c
505050 555555 595959 5d5d5d 616161 656565 696969 6d6d6d 717171 757575
797979 7d7d7d 818181 858585 898989 8d8d8d 919191 959595 999999 9d9d9d
a1a1a1 a5a5a5 aaaaaa aeaeae b2b2b2 b6b6b6 bababa bebebe c2c2c2 c6c6c6
cacaca cecece d2d2d2 d6d6d6 dadada dedede e2e2e2 e6e6e6 eaeaea eeeeee
f2f2f2 f6f6f6 fafafa ffffff
>
] setcolorspace } bdef mwscm
gs np 899 389 mt 0 4225 rl 5356 0 rl 0 -4225 rl cp c np
[5356 0 0 4225 899 389] cc
<< % Image dictionary
/ImageType 1
/Width 600 /Height 600
/BitsPerComponent 8
/Decode [0 255]
/ImageMatrix [600.000000 0 0 600.000000 0 0]
/DataSource currentfile
/ASCII85Decode filter
/RunLengthDecode filter
>>
image
!!%R7JNn]MJNp+u!!%R7JNn]MJNp+u!!%R7JNn]MJNp+u!!%R7JNn]MJNp+u!!%R7JNn]MJNp+u!!%R7
JNn]MJNp+u!!%R7JNn]MJNp+u!!%R7JNn]MJNp+u!!%R7JNn]MJNp+u!!%R7JNn]MJNp+u!!%R7JNn]M
JNp+u!!%R7JNn]MJNp+u!!%R7JNn]MJNp+u!!%R7JNn]MJNp+u!!%R7JNn]MJNp+u!!%R7JNn]MJNp+u
!!%R7JNn]MJNp+u!!%R7JNn]MJNp+u!!%R7JNn]MJNp+u!!%R7JNn]MJNp+u!!%R7JNn]MJNp+u!!%R7
JNn]MJNp+u!!%R7JNn]MJNp+u!!%R7JNn]MJNp+u!!%R7JNn]MJNp+u!!%R7JNn]MJNp+u!!%R7JNn]M
JNp+u!!%R7JNn]MJNp+u!!%R7JNn]MJNp+u!!%R7JNn]MJNp+u!!%R7JNn]MJNp+u!!%R7JNn]MJNp+u
!!%R7JNn]MJNp+u!!%R7JNn]MJNp+u!!%R7JNn]MJNp+u!!%R7JNn]MJNp+u!!rQg!&oHnJNn]MJNp8$
!!rQg!&oHnJNn]MJNp8$!"Aik!&K0jJNn]MJNp8$!#,>r!%iadJNn]MJNpD(!'1$C!#ki$!%WUbJNn]M
JNpD(!'1$C!#GPu!&/sgJNn]MJNpD(!&X[>!"]&n!'GfsJNn]MJNpD(!%n17!#'oJJNn]MJNp8$!%e+6
!#to%!&oHnJNn]MJNpD(!&"78!"jcHJNn]MJNp8$!'C0E!!Ij;JNn]MJNp+u!#>Jt!%NOaJNn]MJNpD(
!'1$C!$V>+!&=I;!&smA!'>`rJNn]MJNpP,!"]&n!&/sgJNn]MJNp8$!!%R7JNn]MJNp8$!&=I;!$hJ-
!&FO<!'#NoJNn]MJNpP,!&jg@!%\%5!&=I;!&=I;!'5ZqJNn]MJNp\0!&4C:!"XWFJNn]MJNp+u!!%R7
JNn]MJNp8$!&OU=!%%V/!&FO<!&]<lJNn]MJNp\0!':*D!&FO<!%e+6!&FO<!&T6kJNn]MJNpD(!#ki$
!%!1\JNn]MJNp8$!#>Jt!%NOaJNn]MJNpP,!':*D!&X[>!%@h2!&"78!'5ZqJNn]MJNph4!'U<G!&OU=
!&FO<!%n17!&=I;!'PltJNn]MJNpP,!%%V/!#gDQJNn]MJNp8$!%Rt4!#:&LJNn]MJNpD(!&smA!&"78
!%In3!&/sgJNn]MJNp\0!':*D!&X[>!&4C:!%e+6!&T6kJNn]MJNp8$!!%R7JNn]MJNp+u!!%R7JNn]M
JNpD(!''sB!%\%5!$M8*!'>`rJNn]MJNpP,!':*D!%Rt4!#gDQJNn]MJNp+u!!%R7JNn]MJNp+u!!%R7
JNn]MJNp8$!''sB!"Jol!'5ZqJNn]MJNp8$!!%R7JNn]MJNp+u!!%R7JNn]MJNp+u!!%R7JNn]MJNp8$
!'C0E!#bc#!%WUbJNn]MJNp8$!!%R7JNn]MJNp+u!!%R7JNn]MJNp+u!!%R7JNn]MJNp8$!%In3!#C,M
JNn]MJNp+u!!%R7JNn]MJNp+u!!%R7JNn]MJNp+u!!%R7JNn]MJNp+u!!%R7JNn]MJNp+u!!%R7JNn]M
JNp+u!!%R7JNn]MJNp+u!!%R7JNn]MJNp+u!!%R7JNn]MJNp+u!!%R7JNn]MJNp+u!!%R7JNn]MJNp+u
!!%R7JNn]MJNp+u!!%R7JNn]MJNp+u!!%R7JNn]MJNp+u!!%R7JNn]MJNp+u!!%R7JNn]MJNp+u!!%R7
JNn]MJNp+u!!%R7JNn]MJNp+u!!%R7JNn]MJNp+u!!%R7JNn]MJNp+u!!%R7JNn]MJNp+u!!%R7JNn]M
JNp+u!!%R7JNn]MJNp+u!!%R7JNn]MJNp+u!!%R7JNn]MJNp+u!!%R7JNn]MJNp+u!!%R7JNn]MJNp+u
!!%R7JNn]MJNp+u!!%R7JNn]MJNp+u!!%R7JNn]MJNp+u!!%R7JNn]MJNp+u!!%R7JNn]MJNp+u!!%R7
JNn]MJNp+u!"Sum!&9$hJNn]MJNp8$!!`Ee!',TpJNn]MJNp8$!!%R7JNn]MJNp+u!!%R7JNn]MJNp8$
!&=I;!"OQEJNn]MJNp8$!'1$C!!\!=JNn]MJNp+u!!%R7JNn]MJNp+u!!%R7JNn]MJNp+u!!%R7JNn]M
JNp+u!!%R7JNn]MJNp+u!!%R7JNn]MJNp+u!!%R7JNn]MJNp+u!!%R7JNn]MJNp+u!!%R7JNn]MJNp+u
!!%R7JNn]MJNp+u!!%R7JNn]MJNp+u!!%R7JNn]MJNp+u!!%R7JNn]MJNp+u!!%R7JNn]MJNp+u!!%R7
JNn]MJNp+u!!%R7JNn]MJNp+u!!%R7JNn]MJNp+u!!iKf!'#NoJNn]MJNp8$!!%R7JNn]MJNp+u!!%R7
JNn]MJNp+u!!%R7JNn]MJNp8$!''sB!!e'>JNn]MJNp+u!!W?d!'5ZqJNn]MJNp8$!!%R7JNn]MJNp+u
!!%R7JNn]MJNp+u!!rQg!&oHnJNn]MJNpD(!'C0E!"]&n!&]<lJNn]MJNp8$!"Sum!&9$hJNn]MJNp8$
!"&Wh!&fBmJNn]MJNpD(!&smA!"8cj!'GfsJNn]MJNpD(!&X[>!"4?BJNn]MJNp8$!&=I;!"OQEJNn]M
JNp8$!&jg@!""3@JNn]MJNp8$!'L6F!!@d:JNn]MJNp+u!!%R7JNn]MJNp+u!!%R7JNn]MJNp+u!!%R7
JNn]MJNp+u!!%R7JNn]MJNp+u!!%R7JNn]MJNp+u!!%R7JNn]MJNp+u!!%R7JNn]MJNp+u!!%R7JNn]M
JNp+u!!%R7JNn]MJNp+u!!%R7JNn]MJNp+u!!%R7JNn]MJNp+u!!%R7JNn]MJNp+u!!%R7JNn]MJNp+u
!!%R7JNn]MJNp+u!!%R7JNn]MJNp+u!!%R7JNn]MJNp+u!!%R7JNn]MJNp+u!!%R7JNn]MJNp+u!!%R7
JNn]MJNp+u!!%R7JNn]MJNp+u!!%R7JNn]MJNp+u!!%R7JNn]MJNp+u!!%R7JNn]MJNp+u!!%R7JNn]M
JNp+u!!%R7JNn]MJNp+u!!%R7JNn]MJNp+u!!%R7JNn]MJNp+u!!%R7JNn]MJNp+u!!%R7JNn]MJNp+u
!!%R7JNn]MJNp+u!!%R7JNn]MJNp+u!!%R7JNn]MJNp+u!!%R7JNn]MJNp+u!!%R7JNn]MJNp+u!!%R7
JNn]MJNp+u!!%R7JNn]MJNp+u!!%R7JNn]MJNp+u!!%R7JNn]MJNp+u!!%R7JNn]MJNp+u!!%R7JNn]M
JNp+u!!%R7JNn]MJNp+u!!%R7JNn]MJNp+u!!%R7JNn]MJNp+u!!%R7JNn]MJNp+u!!%R7JNn]MJNp+u
!!%R7JNn]MJNp+u!!%R7JNn]MJNp+u!!%R7JNn]MJNp+u!!%R7JNn]MJNp+u!!%R7JNn]MJNp+u!!%R7
JNn]MJNp+u!!%R7JNn]MJNp+u!!%R7JNn]MJNp+u!!%R7JNn]MJNp+u!!%R7JNn]MJNp+u!!%R7JNn]M
JNp+u!!%R7JNn]MJNp+u!!%R7JNn]MJNp+u!!%R7JNn]MJNp+u!!%R7JNn]MJNp+u!!%R7JNn]MJNp+u
!!%R7JNn]MJNp+u!!%R7JNn]MJNp+u!!%R7JNn]MJNp+u!!%R7JNn]MJNp+u!!%R7JNn]MJNp+u!!<-a
!'PltJNn]MJNp8$!!%R7JNn]MJNp+u!!%R7JNn]MJNp+u!#Y]"!%3=^JNn]MJNpD(!'U<G!#GPu!%WUb
JNn]MJNp8$!!%R7JNn]MJNp+u!!%R7JNn]MJNp8$!&"78!$hJ-!&OI9!'>`rJNn]MJNp\0!&"78!$2&'
!&sa=!'GfsJNn]MJNpP,!!%R7JNn]MJNp+u!"o&l!''sB!&]<lJNn]MJNp\0!&FO<!$_D,!&+=9!'>`r
JNn]MJNpP,!&jg@!$(u&!%`[cJNn]MJNp8$!!%R7JNn]MJNp+u!!%R7JNn]MJNp8$!%@h2!#L2NJNn]M
JNp8$!%In3!#C,MJNn]MJNp+u!!%R7JNn]MJNpD(!&sa=!!n-?JNn]MJNpP,!':*D!'9s@!#ki$!&/sg
JNn]MJNp\0!'Bm=!#5Ds!&/sgJNn]MJNp8$!!%R7JNn]MJNpP,!&XC6!"4?BJNn]MJNp8$!%n17!"siI
JNn]MJNp8$!&"78!"jcHJNn]MJNp+u!!%R7JNn]MJNp+u!!%R7JNn]MJNp+u!!%R7JNn]MJNp+u!!%R7
JNn]MJNp+u!!%R7JNn]MJNp+u!!%R7JNn]MJNp+u!!%R7JNn]MJNp+u!!%R7JNn]MJNp+u!!%R7JNn]M
JNp+u!!%R7JNn]MJNp+u!!%R7JNn]MJNp+u!!%R7JNn]MJNp+u!!%R7JNn]MJNp+u!!%R7JNn]MJNp+u
!!%R7JNn]MJNp+u!!%R7JNn]MJNp+u!!%R7JNn]MJNp+u!!%R7JNn]MJNp+u!!%R7JNn]MJNp+u!!%R7
JNn]MJNp+u!!%R7JNn]MJNp+u!!%R7JNn]MJNp+u!!%R7JNn]MJNp+u!!%R7JNn]MJNp+u!!%R7JNn]M
JNp+u!!%R7JNn]MJNp+u!!%R7JNn]MJNp+u!!%R7JNn]MJNp+u!!%R7JNn]MJNp+u!!%R7JNn]MJNp+u
!!%R7JNn]MJNp+u!!%R7JNn]MJNp+u!!%R7JNn]MJNp+u!!%R7JNn]MJNp+u!!%R7JNn]MJNp+u!!%R7
JNn]MJNp+u!!%R7JNn]MJNp+u!!%R7JNn]MJNp+u!!%R7JNn]MJNp+u!!%R7JNn]MJNp+u!!%R7JNn]M
JNp+u!!%R7JNn]MJNp+u!!%R7JNn]MJNp+u!!%R7JNn]MJNp+u!!%R7JNn]MJNp+u!!%R7JNn]MJNp+u
!!%R7JNn]MJNp+u!!%R7JNn]MJNp+u!!%R7JNn]MJNp+u!!%R7JNn]MJNp+u!!%R7JNn]MJNp+u!!%R7
JNn]MJNp+u!!%R7JNn]MJNp+u!!%R7JNn]MJNp+u!!%R7JNn]MJNp+u!!%R7JNn]MJNp+u!!%R7JNn]M
JNp+u!!%R7JNn]MJNp+u!!%R7JNn]MJNp+u!!%R7JNn]MJNp+u!!%R7JNn]MJNp+u!!%R7JNn]MJNp+u
!!%R7JNn]MJNp+u!!%R7JNn]MJNp+u!!%R7JNn]MJNp+u!!%R7JNn]MJNp+u!!%R7JNn]MJNp+u!!%R7
JNn]MJNp+u!!%R7JNn]MJNp+u!!%R7JNn]MJNp+u!!%R7JNn]MJNp+u!!%R7JNn]MJNp+u!!%R7JNn]M
JNp+u!!%R7JNn]MJNp+u!!%R7JNn]MJNp+u!!%R7JNn]MJNp+u!!%R7JNn]MJNp+u!!%R7JNn]MJNp+u
!!%R7JNn]MJNp+u!!%R7JNn]MJNp+u!!%R7JNn]MJNp+u!!%R7JNn]MJNp+u!!%R7JNn]MJNp+u!!%R7
JNn]MJNp+u!!%R7JNn]MJNp+u!!%R7JNn]MJNp+u!!%R7JNn]MJNp+u!!rQg!&oHnJNn]MJNp8$!!%R7
JNn]MJNp+u!!%R7JNn]MJNp+u!!%R7JNn]MJNp8$!&smA!!n-?JNn]MJNp+u!!%R7JNn]MJNp+u!!%R7
JNn]MJNp+u!!%R7JNn]MJNp+u!!%R7JNn]MJNp+u!!%R7JNn]MJNp+u!!%R7JNn]MJNp+u!!%R7JNn]M
JNp+u!!%R7JNn]MJNp+u!!%R7JNn]MJNp+u!!%R7JNn]MJNp+u!!%R7JNn]MJNp+u!!%R7JNn]MJNp+u
!!%R7JNn]MJNp+u!!%R7JNn]MJNp+u!!%R7JNn]MJNp+u!!%R7JNn]MJNp+u!!%R7JNn]MJNp+u!!%R7
JNn]MJNp+u!!%R7JNn]MJNp+u!!%R7JNn]MJNp+u!!%R7JNn]MJNp+u!!%R7JNn]MJNp+u!!%R7JNn]M
JNp+u!!%R7JNn]MJNp+u!!%R7JNn]MJNp+u!!%R7JNn]MJNp+u!!%R7JNn]MJNp+u!!%R7JNn]MJNp+u
!!%R7JNn]MJNp+u!!%R7JNn]MJNp+u!!%R7JNn]MJNp+u!!%R7JNn]MJNp+u!!%R7JNn]MJNp+u!!%R7
JNn]MJNp+u!!%R7JNn]MJNp+u!!%R7JNn]MJNp+u!!%R7JNn]MJNp+u!!%R7JNn]MJNp+u!!%R7JNn]M
JNp+u!!%R7JNn]MJNp+u!!%R7JNn]MJNp+u!!%R7JNn]MJNp+u!!%R7JNn]MJNp+u!!%R7JNn]MJNp+u
!!%R7JNn]MJNp+u!!%R7JNn]MJNp+u!!%R7JNn]MJNp+u!!%R7JNn]MJNp+u!<GB85Cl%M5Cl&"563qK
JNn]MJNn]MWBU@_JNn]MJNn]MWBU@_JNn]MJNn]MWBU@_JNn]MJNn]MWBU@_JNn]MJNn]MWBU@_JNn]M
JNn]MWBU@_JNn]MJNn]MWBU@_JNn]MJNn]MWBU@_JNn]MJNn]MWBU@_JNn]MJNn]MWBU@_JNn]MJNn]M
WBU@_JNn]MJNn]MWBU@_JNn]MJNn]MWBU@_JNn]MJNn]MWBU@_JNn]MJNn]MWBU@_JNn]MJNn]MWBU@_
JNn]MJNn]MWBU@_JNn]MJNn]MWBU@_JNn]MJNn]MWBU@_JNn]MJNn]MWBU@_JNn]MJNn]MWBU@_JNn]M
JNn]MWBU@_JNn]MJNn]MWBU@_JNn]MJNn]MWBU@_JNn]MJNn]MWBU@_JNn]MJNn]MWBU@_JNn]MJNn]M
WBU@_JNn]MJNn]MWBU@_JNn]MJNn]MWBU@_JNn]MJNn]MWBU@_JNn]MJNn]MWBU@_JNn]MJNn]MWBU@_
JNn]MJNn]MWBU@_JNn]MJNn]MWBU@_JNn]MJNn]MWBU@_JNn]MJNn]MWBU@_JNn]MJNn]MWBU@_JNn]M
JNn]MWBU@_JNn]MJNn]MWBU@_JNn]MJNn]MWBU@_JNn]MJNn]MWBU@_JNn]MJNn]MWBU@_JNn]MJNn]M
WBU@_JNn]MJNn]MWBU@_JNn]MJNn]MWBU@_JNn]MJNn]MWBU@_JNn]MJNn]MWBU@_JNn]MJNn]MWBU@_
JNn]MJNn]MWBU@_JNn]MJNn]MWBU@_JNn]MJNn]MWBU@_JNn]MJNn]MWBU@_JNn]MJNn]MWBU@_JNn]M
JNn]MWBU@_JNn]MJNn]MWBU@_JNn]MJNn]MWBU@_JNn]MJNn]MWBU@_JNn]MJNn]MWBU@_JNn]MJNn]M
WBU@_JNn]MJNn]MWBU@_JNn]MJNn]MWBU@_JNn]MJNn]MWBU@_JNn]MJNn]MWBU@_JNn]MJNn]MWBU@_
JNn]MJNn]MWBU@_JNn]MJNn]MWBU@_JNn]MJNn]MWBU@_JNn]MJNn]MWBU@_JNn]MJNn]MWBU@_JNn]M
JNn]MWBU@_JNn]MJNn]MWBU@_JNn]MJNn]MWBU@_JNn]MJNn]MWBU@_JNn]MJNn]MWBU@_JNn]MJNn]M
WBU@_JNn]MJNn]MWBU@_JNn]MJNn]MWBU@_JNn]MJNn]MWBU@_JNn]MJNn]MWBU@_JNn]MJNn]MWBU@_
JNn]MJNn]MWBU@_JNn]MJNn]MWBU@_JNn]MJNn]MWBU@_JNn]MJNn]MWBU@_JNn]MJNn]MWBU@_JNn]M
JNn]MWBU@_JNn]MJNn]MWBU@_JNn]MJNn]MWBU@_JNn]MJNn]MWBU@_JNn]MJNn]MWBU@_JNn]MJNn]M
WBU@_JNn]MJNn]MWBU@_JNn]MJNn]MWBU@_JNn]MJNn]MWBU@_JNn]MJNn]MWBU@_JNn]MJNn]MWBU@_
JNn]MJNn]MWBU@_JNn]MJNn]MWBU@_JNn]MJNn]MWBU@_JNn]MJNn]MWBU@_JNn]MJNn]MWBU@_JNn]M
JNn]MWBU@_JNn]MJNn]MWBU@_JNn]MJNn]MWBU@_JNn]MJNn]MWBU@_JNn]MJNn]MWBU@_JNn]MJNn]M
WBU@_JNn]MJNn]MWBU@_JNn]MJNn]MWBU@_JNn]MJNn]MWBU@_JNn]MJNn]MWBU@_JNn]MJNn]MWBU@_
JNn]MJNn]MWBU@_JNn]MJNn]MWBU@_JNn]MJNn]MWBU@_JNn]MJNn]MWBU@_JNn]MJNn]MWBU@_JNn]M
JNn]MWBU@_JNn]MJNn]MWBU@_JNn]MJNn]MWBU@_JNn]MJNn]MWBU@_JNn]MJNn]MWBU@_JNn]MJNn]M
WBU@_JNn]MJNn]MWBU@_JNn]MJNn]MWBU@_JNn]MJNn]MWBU@_JNn]MJNn]MWBU@_JNn]MJNn]MWBU@_
JNn]MJNn]MWBU@_JNn]MJNn]MWBU@_JNn]MJNn]MWBU@_JNn]MJNn]MWBU@_JNn]MJNn]MWBU@_JNn]M
JNn]MWBU@_JNn]MJNn]MWBU@_JNn]MJNn]MWBU@_JNn]MJNn]MWBU@_JNn]MJNn]MWBU@_JNn]MJNn]M
WBU@_JNn]MJNn]MWBU@_JNn]MJNn]MWBU@_JNn]MJNn]MWBU@_JNn]MJNn]MWBU@_JNn]MJNn]MWBU@_
JNn]MJNn]MWBU@_JNn]MJNn]MWBU@_JNn]MJNn]MWBU@_JNn]MJNn]MWBU@_JNn]MJNn]MWBU@_JNn]M
JNn]MWBU@_JNn]MJNn]MWBU@_JNn]MJNn]MWBU@_JNn]MJNn]MWBU@_JNn]MJNn]MWBU@_JNn]MJNn]M
WBU@_JNn]MJNn]MWBU@_JNn]MJNn]MWBU@_JNn]MJNn]MWBU@_JNn]MJNn]MWBU@_JNn]MJNn]MWBU@_
JNn]MJNn]MWBU@_JNn]MJNn]MWBU@_JNn]MJNn]MWBU@_JNn]MJNn]MWBU@_JNn]MJNn]MWBU@_JNn]M
JNn]MWBU@_JNn]MJNn]MWBU@_JNn]MJNn]MWBU@_JNn]MJNn]MWBU@_JNn]MJNn]MWBU@_JNn]MJNn]M
WBU@_JNn]MJNn]MWBU@_JNn]MJNn]MWBU@_JNn]MJNn]MWBU@_JNn]MJNn]MWBU@_JNn]MJNn]MWBU@_
JNn]MJNn]MWBU@_JNn]MJNn]MWBU@_JNn]MJNn]MWBU@_JNn]MJNn]MWBU@_JNn]MJNn]MWBU@_JNn]M
JNn]MWBU@_JNn]MJNn]MWBU@_JNn]MJNn]MWBU@_JNn]MJNn]MWBU@_JNn]MJNn]MWBU@_JNn]MJNn]M
WBU@_JNn]MJNn]MWBU@_JNn]MJNn]MWBU@_JNn]MJNn]MWBU@_JNn]MJNn]MWBU@_JNn]MJNn]MWBU@_
JNn]MJNn]MWBU@_JNn]MJNn]MWBU@_JNn]MJNn]MWBU@_JNn]MJNn]MWBU@_JNn]MJNn]MWBU@_JNn]M
JNn]MWBU@_JNn]MJNn]MWBU@_JNn]MJNn]MWBU@_JNn]MJNn]MWBU@_JNn]MJNn]MWBU@_JNn]MJNn]M
WBU@_JNn]MJNn]MWBU@_JNn]MJNn]MWBU@_JNn]MJNn]MWBU@_JNn]MJNn]MWBU@_JNn]MJNn]MWBU@_
JNn]MJNn]MWBU@_JNn]MJNn]MWBU@_JNn]MJNn]MWBU@_JNn]MJNn]MWBU@_JNn]MJNn]MWBU@_JNn]M
JNn]MWBU@_JNn]MJNn]MWBU@_JNn]MJNn]MWBU@_J,
~>
gr
gr
4 w
DO
SO
6 w
0 sg
899 4614 mt 6255 4614 L
899 389 mt 6255 389 L
6255 389 mt 6255 4614 L
899 389 mt 899 4614 L
899 4614 mt 6255 4614 L
899 389 mt 899 4614 L
1787 4614 mt 1787 4560 L
1787 389 mt 1787 442 L
%%IncludeResource: font Helvetica
/Helvetica /WindowsLatin1Encoding 120 FMSR
1687 4759 mt
(100) s
2679 4614 mt 2679 4560 L
2679 389 mt 2679 442 L
2579 4759 mt
(200) s
3572 4614 mt 3572 4560 L
3572 389 mt 3572 442 L
3472 4759 mt
(300) s
4465 4614 mt 4465 4560 L
4465 389 mt 4465 442 L
4365 4759 mt
(400) s
5357 4614 mt 5357 4560 L
5357 389 mt 5357 442 L
5257 4759 mt
(500) s
6250 4614 mt 6250 4560 L
6250 389 mt 6250 442 L
6150 4759 mt
(600) s
899 1089 mt 952 1089 L
6255 1089 mt 6201 1089 L
664 1133 mt
(100) s
899 1793 mt 952 1793 L
6255 1793 mt 6201 1793 L
664 1837 mt
(200) s
899 2497 mt 952 2497 L
6255 2497 mt 6201 2497 L
664 2541 mt
(300) s
899 3202 mt 952 3202 L
6255 3202 mt 6201 3202 L
664 3246 mt
(400) s
899 3906 mt 952 3906 L
6255 3906 mt 6201 3906 L
664 3950 mt
(500) s
899 4610 mt 952 4610 L
6255 4610 mt 6201 4610 L
664 4654 mt
(600) s
899 4614 mt 6255 4614 L
899 389 mt 6255 389 L
6255 389 mt 6255 4614 L
899 389 mt 899 4614 L
end %%Color Dict
eplot
%%EndObject
epage
end
showpage
%%Trailer
%%EOF
+52 -58
View File
@@ -34,79 +34,73 @@
\@writefile{toc}{\contentsline {subsubsection}{\numberline {3.1.3}Ball trees}{5}{subsubsection.16}}
\@writefile{toc}{\contentsline {subsubsection}{\numberline {3.1.4}Nearest Neighbor Algorithm}{6}{subsubsection.18}}
\@writefile{toc}{\contentsline {subsubsection}{\numberline {3.1.5}All nearest Neighbors}{7}{subsubsection.19}}
\@writefile{toc}{\contentsline {subsubsection}{\numberline {3.1.6}Fast kernel summation with trees}{8}{subsubsection.20}}
\newlabel{kernel_sum}{{3.2}{8}{Fast kernel summation with trees\relax }{equation.21}{}}
\@writefile{toc}{\contentsline {subsection}{\numberline {3.2}Kernel methods}{8}{subsection.25}}
\@writefile{toc}{\contentsline {subsubsection}{\numberline {3.2.1}Kernel principal component analysis}{9}{subsubsection.26}}
\@writefile{toc}{\contentsline {subsubsection}{\numberline {3.2.2}Customized Kernels}{10}{subsubsection.30}}
\@writefile{toc}{\contentsline {subsubsection}{\numberline {3.2.3}Stochastic optimization}{11}{subsubsection.35}}
\@writefile{toc}{\contentsline {subsubsection}{\numberline {3.2.4}Support Vector Machines}{11}{subsubsection.36}}
\@writefile{toc}{\setcounter {tocdepth}{2}}
\@writefile{toc}{\contentsline {section}{\numberline {\MakeUppercase {Section} 4}\MakeUppercase {Proposed Research}}{11}{section.37}}
\newlabel{proposed}{{4}{11}{Proposed Research\relax }{section.37}{}}
\@writefile{toc}{\contentsline {subsection}{\numberline {4.1}Large scale trees, for out of core memory}{11}{subsection.38}}
\@writefile{toc}{\contentsline {subsection}{\numberline {4.2}Memory layout of Kd-trees}{13}{subsection.39}}
\@writefile{toc}{\contentsline {subsection}{\numberline {4.3}System architecture}{13}{subsection.40}}
\@writefile{toc}{\contentsline {subsection}{\numberline {4.4}Memory Manager Architecture}{14}{subsection.41}}
\@writefile{toc}{\contentsline {subsection}{\numberline {4.5}Memory Mapped Files}{14}{subsection.42}}
\@writefile{toc}{\contentsline {subsection}{\numberline {4.6}User defined cache, with TPIE}{15}{subsection.43}}
\@writefile{toc}{\contentsline {subsection}{\numberline {3.2}Large scale trees, for out of core memory}{8}{subsection.20}}
\@writefile{toc}{\contentsline {subsection}{\numberline {3.3}Memory layout of Kd-trees}{9}{subsection.21}}
\@writefile{toc}{\contentsline {subsection}{\numberline {3.4}System architecture}{10}{subsection.22}}
\@writefile{toc}{\contentsline {subsection}{\numberline {3.5}Memory Manager Architecture}{10}{subsection.23}}
\@writefile{toc}{\contentsline {subsection}{\numberline {3.6}Memory Mapped Files}{11}{subsection.24}}
\@writefile{toc}{\contentsline {subsection}{\numberline {3.7}User defined cache, with TPIE}{11}{subsection.25}}
\@writefile{toc}{\contentsline {subsection}{\numberline {3.8}Kernel methods}{12}{subsection.26}}
\@writefile{toc}{\contentsline {subsubsection}{\numberline {3.8.1}Kernel principal component analysis}{13}{subsubsection.27}}
\@writefile{toc}{\contentsline {subsubsection}{\numberline {3.8.2}Customized Kernels}{14}{subsubsection.31}}
\citation{quatieri2002dts}
\citation{bentley1975bst}
\citation{moore2000ahu}
\citation{gray2000nbp}
\citation{ravindran:inr}
\@writefile{toc}{\contentsline {subsection}{\numberline {4.7}Customizing Kernels with fast gradient descent}{16}{subsection.44}}
\@writefile{toc}{\contentsline {subsection}{\numberline {4.8}Scalable Support Vector Machines}{16}{subsection.45}}
\@writefile{toc}{\contentsline {subsection}{\numberline {4.9}Computing Gaussian Kernel Matrix with trees}{16}{subsection.46}}
\@writefile{toc}{\contentsline {subsection}{\numberline {4.10}Intrinsic dimension of speech}{16}{subsection.47}}
\@writefile{toc}{\contentsline {subsection}{\numberline {4.11}Speech recognition with nearest neighbor search}{16}{subsection.48}}
\@writefile{toc}{\contentsline {subsection}{\numberline {3.9}Forming manifolds from data for clustering through non-parametric density estimation}{15}{subsection.36}}
\@writefile{toc}{\contentsline {subsection}{\numberline {3.10}Speech recognition with nearest neighbor search}{15}{subsection.37}}
\@writefile{toc}{\contentsline {subsection}{\numberline {3.11}Introduction}{15}{subsection.38}}
\citation{jansen2006ifa}
\citation{scholkopf:nca}
\citation{coifman2006dm}
\@writefile{toc}{\contentsline {subsection}{\numberline {4.12}Introduction}{17}{subsection.49}}
\citation{shamma_main}
\@writefile{toc}{\contentsline {subsection}{\numberline {3.12}Speech features}{16}{subsection.39}}
\citation{moore-tutorial}
\@writefile{toc}{\contentsline {subsection}{\numberline {4.13}Speech features}{18}{subsection.50}}
\@writefile{toc}{\contentsline {subsection}{\numberline {4.14}Fast N-body methods for feature comparison}{18}{subsection.51}}
\@writefile{lof}{\contentsline {figure}{\numberline {Figure\nobreakspace 4.1}{\ignorespaces A two dimensional kd-tree}}{19}{figure.53}}
\newlabel{kdtree}{{4.1}{19}{A two dimensional kd-tree\relax }{figure.53}{}}
\@writefile{lot}{\contentsline {table}{\numberline {Table\nobreakspace 4.1}{\ignorespaces Pseudo-code for the dual-tree all nearest neighbor algorithm}}{20}{table.52}}
\newlabel{dualtree_algorithm}{{4.1}{20}{Pseudo-code for the dual-tree all nearest neighbor algorithm\relax }{table.52}{}}
\@writefile{lof}{\contentsline {figure}{\numberline {Figure\nobreakspace 4.2}{\ignorespaces Simulation of the dual tree algorithm}}{20}{figure.54}}
\newlabel{dualkdtree}{{4.2}{20}{Simulation of the dual tree algorithm\relax }{figure.54}{}}
\@writefile{toc}{\contentsline {subsection}{\numberline {4.15}Kernel PCA on a manifold}{20}{subsection.55}}
\@writefile{toc}{\contentsline {subsection}{\numberline {4.16}Experimental Results}{21}{subsection.56}}
\@writefile{toc}{\contentsline {subsection}{\numberline {4.17}All nearest neighbor performance}{21}{subsection.57}}
\@writefile{toc}{\contentsline {subsection}{\numberline {4.18}Comparison of NRAF and MFCC}{22}{subsection.58}}
\@writefile{toc}{\contentsline {subsection}{\numberline {3.13}Fast N-body methods for feature comparison}{17}{subsection.40}}
\@writefile{lot}{\contentsline {table}{\numberline {Table\nobreakspace 3.1}{\ignorespaces Pseudo-code for the dual-tree all nearest neighbor algorithm}}{18}{table.41}}
\newlabel{dualtree_algorithm}{{3.1}{18}{Pseudo-code for the dual-tree all nearest neighbor algorithm\relax }{table.41}{}}
\@writefile{lof}{\contentsline {figure}{\numberline {Figure\nobreakspace 3.1}{\ignorespaces A two dimensional kd-tree}}{18}{figure.42}}
\newlabel{kdtree}{{3.1}{18}{A two dimensional kd-tree\relax }{figure.42}{}}
\@writefile{lof}{\contentsline {figure}{\numberline {Figure\nobreakspace 3.2}{\ignorespaces Simulation of the dual tree algorithm}}{18}{figure.43}}
\newlabel{dualkdtree}{{3.2}{18}{Simulation of the dual tree algorithm\relax }{figure.43}{}}
\@writefile{toc}{\contentsline {subsection}{\numberline {3.14}Kernel PCA on a manifold}{18}{subsection.44}}
\@writefile{toc}{\contentsline {subsection}{\numberline {3.15}Experimental Results}{19}{subsection.45}}
\@writefile{toc}{\contentsline {subsection}{\numberline {3.16}All nearest neighbor performance}{20}{subsection.46}}
\@writefile{toc}{\contentsline {subsection}{\numberline {3.17}Comparison of NRAF and MFCC}{20}{subsection.47}}
\citation{ravindran:inr}
\@writefile{toc}{\contentsline {subsection}{\numberline {4.19}Discussion}{23}{subsection.60}}
\@writefile{lot}{\contentsline {table}{\numberline {Table\nobreakspace 4.2}{\ignorespaces CPU time (in hours) for evaluating all 20 nearest neighbors on TIMIT database}}{23}{table.63}}
\newlabel{timing}{{4.19}{23}{Discussion\relax }{table.63}{}}
\@writefile{lof}{\contentsline {figure}{\numberline {Figure\nobreakspace 4.3}{\ignorespaces Principal Component Analysis of TIMIT (top) and NTIMIT (bottom) for MFCC (solid line) and for NRAF (dashed line).}}{24}{figure.61}}
\newlabel{pca}{{4.3}{24}{Principal Component Analysis of TIMIT (top) and NTIMIT (bottom) for MFCC (solid line) and for NRAF (dashed line)}{figure.61}{}}
\@writefile{lof}{\contentsline {figure}{\numberline {Figure\nobreakspace 4.4}{\ignorespaces Kernel PCA spectrum for TIMIT MFCC features}}{25}{figure.62}}
\newlabel{kernel_pca}{{4.4}{25}{Kernel PCA spectrum for TIMIT MFCC features\relax }{figure.62}{}}
\@writefile{lot}{\contentsline {table}{\numberline {Table\nobreakspace 4.3}{\ignorespaces Nearest k-neighbor classfier for TIMIT database}}{25}{table.64}}
\newlabel{kneighborTIMIT}{{4.19}{25}{Discussion\relax }{table.64}{}}
\@writefile{toc}{\contentsline {subsection}{\numberline {4.20}Acoustic Environment adaptation}{25}{subsection.68}}
\@writefile{toc}{\contentsline {subsection}{\numberline {3.18}Discussion}{21}{subsection.49}}
\@writefile{lof}{\contentsline {figure}{\numberline {Figure\nobreakspace 3.3}{\ignorespaces Principal Component Analysis of TIMIT (top) and NTIMIT (bottom) for MFCC (solid line) and for NRAF (dashed line).}}{22}{figure.50}}
\newlabel{pca}{{3.3}{22}{Principal Component Analysis of TIMIT (top) and NTIMIT (bottom) for MFCC (solid line) and for NRAF (dashed line)}{figure.50}{}}
\@writefile{lot}{\contentsline {table}{\numberline {Table\nobreakspace 3.2}{\ignorespaces CPU time (in hours) for evaluating all 20 nearest neighbors on TIMIT database}}{22}{table.52}}
\newlabel{timing}{{3.18}{22}{Discussion\relax }{table.52}{}}
\@writefile{toc}{\setcounter {tocdepth}{2}}
\@writefile{toc}{\contentsline {section}{\numberline {\MakeUppercase {Section} 5}\MakeUppercase {Work remaining to be done}}{25}{section.69}}
\newlabel{remains}{{5}{25}{Work remaining to be done\relax }{section.69}{}}
\@writefile{toc}{\contentsline {section}{\numberline {\MakeUppercase {Section} 4}\MakeUppercase {Proposed Research}}{22}{section.57}}
\newlabel{proposed}{{4}{22}{Proposed Research\relax }{section.57}{}}
\@writefile{lof}{\contentsline {figure}{\numberline {Figure\nobreakspace 3.4}{\ignorespaces Kernel PCA spectrum for TIMIT MFCC features}}{23}{figure.51}}
\newlabel{kernel_pca}{{3.4}{23}{Kernel PCA spectrum for TIMIT MFCC features\relax }{figure.51}{}}
\@writefile{lot}{\contentsline {table}{\numberline {Table\nobreakspace 3.3}{\ignorespaces Nearest k-neighbor classfier for TIMIT database}}{23}{table.53}}
\newlabel{kneighborTIMIT}{{3.18}{23}{Discussion\relax }{table.53}{}}
\@writefile{toc}{\contentsline {subsection}{\numberline {4.1}Customizing Kernels with fast gradient descent}{23}{subsection.58}}
\@writefile{lot}{\contentsline {table}{\numberline {Table\nobreakspace 3.4}{\ignorespaces Nearest k-neighbor classfier for NTIMIT database}}{23}{table.54}}
\newlabel{kneighborNTIMIT}{{3.18}{23}{Discussion\relax }{table.54}{}}
\@writefile{lot}{\contentsline {table}{\numberline {Table\nobreakspace 3.5}{\ignorespaces Nearest k-neighbor classfier for Leave One Speaker out Cross Validation in TIMIT database}}{24}{table.55}}
\newlabel{kneighborLOOCVTIMIT}{{3.18}{24}{Discussion\relax }{table.55}{}}
\@writefile{lot}{\contentsline {table}{\numberline {Table\nobreakspace 3.6}{\ignorespaces Nearest k-neighbor classfier for Leave One Speaker Out Cross validation in NTIMIT database}}{24}{table.56}}
\newlabel{kneighborLOOCVNTIMIT}{{3.18}{24}{Discussion\relax }{table.56}{}}
\@writefile{toc}{\contentsline {subsubsection}{\numberline {4.1.1}Fast kernel summation with trees}{24}{subsubsection.59}}
\newlabel{kernel_sum}{{4.1}{24}{Fast kernel summation with trees\relax }{equation.60}{}}
\@writefile{toc}{\contentsline {subsection}{\numberline {4.2}Computing Gaussian Kernel Matrix with trees}{25}{subsection.64}}
\@writefile{toc}{\contentsline {subsection}{\numberline {4.3}Intrinsic dimension of speech}{25}{subsection.65}}
\@writefile{toc}{\contentsline {subsection}{\numberline {4.4}Speech recognition with nearest neighbor search}{25}{subsection.66}}
\@writefile{toc}{\contentsline {subsection}{\numberline {4.5}Acoustic Environment adaptation}{25}{subsection.67}}
\@writefile{toc}{\setcounter {tocdepth}{2}}
\@writefile{toc}{\contentsline {section}{\numberline {\MakeUppercase {Section} 6}\MakeUppercase {Facilities and Equipment Needed}}{25}{section.70}}
\newlabel{facitlies}{{6}{25}{Facilities and Equipment Needed\relax }{section.70}{}}
\@writefile{lot}{\contentsline {table}{\numberline {Table\nobreakspace 4.4}{\ignorespaces Nearest k-neighbor classfier for NTIMIT database}}{25}{table.65}}
\newlabel{kneighborNTIMIT}{{4.19}{25}{Discussion\relax }{table.65}{}}
\@writefile{lot}{\contentsline {table}{\numberline {Table\nobreakspace 4.5}{\ignorespaces Nearest k-neighbor classfier for Leave One Speaker out Cross Validation in TIMIT database}}{26}{table.66}}
\newlabel{kneighborLOOCVTIMIT}{{4.19}{26}{Discussion\relax }{table.66}{}}
\@writefile{lot}{\contentsline {table}{\numberline {Table\nobreakspace 4.6}{\ignorespaces Nearest k-neighbor classfier for Leave One Speaker Out Cross validation in NTIMIT database}}{26}{table.67}}
\newlabel{kneighborLOOCVNTIMIT}{{4.19}{26}{Discussion\relax }{table.67}{}}
\@writefile{toc}{\setcounter {tocdepth}{2}}
\@writefile{toc}{\contentsline {section}{\numberline {\MakeUppercase {Appendix} A}\MakeUppercase {This is my first appendix}}{27}{section.71}}
\@writefile{toc}{\contentsline {section}{\numberline {\MakeUppercase {Section} 5}\MakeUppercase {Work remaining to be done}}{25}{section.68}}
\newlabel{remains}{{5}{25}{Work remaining to be done\relax }{section.68}{}}
\bibstyle{ieeetr}
\@writefile{toc}{\setcounter {tocdepth}{2}}
\@writefile{toc}{\contentsline {section}{\numberline {\MakeUppercase {Appendix} B}\MakeUppercase {This is my second appendix}}{28}{section.72}}
\@writefile{toc}{\contentsline {section}{\numberline {\MakeUppercase {Appendix} A}\MakeUppercase {This is my first appendix}}{26}{section.69}}
\bibdata{template}
\@writefile{toc}{\setcounter{tocdepth}{3}}
\@writefile{toc}{\contentsline {section}{\MakeUppercase {References}}{29}{section*.73}}
\@writefile{toc}{\contentsline {section}{\MakeUppercase {References}}{27}{section*.70}}
\@writefile{toc}{\setcounter {tocdepth}{2}}
+1 -1
View File
@@ -1,7 +1,7 @@
This is BibTeX, Version 0.99cThe top-level auxiliary file: proposal.aux
The style file: ieeetr.bst
I couldn't open database file template.bib
---line 109 of file proposal.aux
---line 103 of file proposal.aux
: \bibdata{template
: }
I'm skipping whatever remains of this command
Binary file not shown.
+4 -4
View File
@@ -1,4 +1,4 @@
\contentsline {figure}{\numberline {Figure\nobreakspace 4.1}{\ignorespaces A two dimensional kd-tree}}{19}{figure.53}
\contentsline {figure}{\numberline {Figure\nobreakspace 4.2}{\ignorespaces Simulation of the dual tree algorithm}}{20}{figure.54}
\contentsline {figure}{\numberline {Figure\nobreakspace 4.3}{\ignorespaces Principal Component Analysis of TIMIT (top) and NTIMIT (bottom) for MFCC (solid line) and for NRAF (dashed line).}}{24}{figure.61}
\contentsline {figure}{\numberline {Figure\nobreakspace 4.4}{\ignorespaces Kernel PCA spectrum for TIMIT MFCC features}}{25}{figure.62}
\contentsline {figure}{\numberline {Figure\nobreakspace 3.1}{\ignorespaces A two dimensional kd-tree}}{18}{figure.42}
\contentsline {figure}{\numberline {Figure\nobreakspace 3.2}{\ignorespaces Simulation of the dual tree algorithm}}{18}{figure.43}
\contentsline {figure}{\numberline {Figure\nobreakspace 3.3}{\ignorespaces Principal Component Analysis of TIMIT (top) and NTIMIT (bottom) for MFCC (solid line) and for NRAF (dashed line).}}{22}{figure.50}
\contentsline {figure}{\numberline {Figure\nobreakspace 3.4}{\ignorespaces Kernel PCA spectrum for TIMIT MFCC features}}{23}{figure.51}
+47 -49
View File
@@ -1,4 +1,4 @@
This is pdfeTeX, Version 3.141592-1.30.6-2.2 (MiKTeX 2.5) (preloaded format=latex 2007.2.26) 30 JUL 2007 18:04
This is pdfeTeX, Version 3.141592-1.30.6-2.2 (MiKTeX 2.5) (preloaded format=latex 2007.2.26) 30 JUL 2007 19:52
entering extended mode
**D:/nick/fastlib/u/nvasil/docs/proposal/proposal.tex
(D:\nick\fastlib\u\nvasil\docs\proposal\proposal.tex
@@ -547,61 +547,64 @@ rec-om-mended for speech MFCC,
[]
Section 3.
[3] [4] [5] [6] [7] [8] [9] [10]
Section 4.
[11] [12] [13] [14] [15]
[3] [4] [5] [6] [7] [8] [9] [10] [11] [12] [13]
Overfull \hbox (13.90266pt too wide) in paragraph at lines 625--625
[]\OT1/txr/bx/n/14.4 Forming man-i-folds from data for clus-ter-ing through non
-parametric
[]
LaTeX Warning: Citation `quatieri2002dts' on page 16 undefined on input line 70
0.
[14]
LaTeX Warning: Citation `quatieri2002dts' on page 15 undefined on input line 63
8.
LaTeX Warning: Citation `bentley1975bst' on page 16 undefined on input line 714
LaTeX Warning: Citation `bentley1975bst' on page 15 undefined on input line 652
.
LaTeX Warning: Citation `moore2000ahu' on page 16 undefined on input line 715.
LaTeX Warning: Citation `moore2000ahu' on page 15 undefined on input line 653.
LaTeX Warning: Citation `gray2000nbp' on page 16 undefined on input line 718.
LaTeX Warning: Citation `gray2000nbp' on page 15 undefined on input line 656.
LaTeX Warning: Citation `ravindran:inr' on page 16 undefined on input line 721.
LaTeX Warning: Citation `ravindran:inr' on page 15 undefined on input line 659.
[15]
LaTeX Warning: Citation `jansen2006ifa' on page 16 undefined on input line 669.
LaTeX Warning: Citation `scholkopf:nca' on page 16 undefined on input line 675.
LaTeX Warning: Citation `coifman2006dm' on page 16 undefined on input line 676.
LaTeX Warning: Citation `shamma_main' on page 16 undefined on input line 691.
[16]
LaTeX Warning: Citation `jansen2006ifa' on page 17 undefined on input line 731.
LaTeX Warning: Citation `scholkopf:nca' on page 17 undefined on input line 737.
LaTeX Warning: Citation `coifman2006dm' on page 17 undefined on input line 738.
[17]
LaTeX Warning: Citation `shamma_main' on page 18 undefined on input line 753.
LaTeX Warning: Citation `moore-tutorial' on page 18 undefined on input line 776
LaTeX Warning: Citation `moore-tutorial' on page 17 undefined on input line 714
.
[18]
LaTeX Font Info: Try loading font information for OT1+txtt on input line 803
LaTeX Font Info: Try loading font information for OT1+txtt on input line 741
.
("C:\Program Files\MiKTeX 2.5\tex\latex\txfonts\ot1txtt.fd"
("C:\Program Files\MiKTeX 2.5\tex\latex\txfonts\ot1txtt.fd"
File: ot1txtt.fd 2000/12/15 v3.1
)
File: kdtree.eps Graphic file (type eps)
<kdtree.eps>
File: dualtree_recursion.eps Graphic file (type eps)
<dualtree_recursion.eps> [19] [20] [21] [22]
<dualtree_recursion.eps> [17] [18] [19] [20]
LaTeX Warning: Citation `ravindran:inr' on page 23 undefined on input line 943.
LaTeX Warning: Citation `ravindran:inr' on page 21 undefined on input line 881.
File: timit_pca.eps Graphic file (type eps)
@@ -610,23 +613,18 @@ File: ntimit_pca.eps Graphic file (type eps)
<ntimit_pca.eps>
File: kernel_pca.eps Graphic file (type eps)
<kernel_pca.eps>
LaTeX Warning: `!h' float specifier changed to `!ht'.
Section 4.
[21] [22] [23] [24]
Section 5.
Section 6.
[23] [24]
Overfull \vbox (11.8856pt too high) has occurred while \output is active []
[25]
******** Body Ends ********
******** EndMatter Begins ********
[26]
Appendix A.
[27
]
Appendix B.
***************** References ******************
[28
[26
] (proposal.bbl
@@ -634,7 +632,7 @@ LaTeX Warning: Empty ``thebibliography'' environment on input line 3.
)
******** EndMatter Ends ********
[29
[27
] (proposal.aux)
@@ -643,12 +641,12 @@ LaTeX Warning: There were undefined references.
)
Here is how much of TeX's memory you used:
5349 strings out of 95420
68647 string characters out of 1187320
129426 words of memory out of 1125958
8213 multiletter control sequences out of 60000
5347 strings out of 95420
68627 string characters out of 1187320
129307 words of memory out of 1125839
8211 multiletter control sequences out of 60000
28619 words of font info for 89 fonts, out of 1000000 for 2000
14 hyphenation exceptions out of 8191
34i,8n,36p,341b,521s stack positions out of 5000i,500n,10000p,200000b,32768s
Output written on proposal.dvi (34 pages, 135100 bytes).
Output written on proposal.dvi (32 pages, 133300 bytes).
+6 -6
View File
@@ -1,6 +1,6 @@
\contentsline {table}{\numberline {Table\nobreakspace 4.1}{\ignorespaces Pseudo-code for the dual-tree all nearest neighbor algorithm}}{20}{table.52}
\contentsline {table}{\numberline {Table\nobreakspace 4.2}{\ignorespaces CPU time (in hours) for evaluating all 20 nearest neighbors on TIMIT database}}{23}{table.63}
\contentsline {table}{\numberline {Table\nobreakspace 4.3}{\ignorespaces Nearest k-neighbor classfier for TIMIT database}}{25}{table.64}
\contentsline {table}{\numberline {Table\nobreakspace 4.4}{\ignorespaces Nearest k-neighbor classfier for NTIMIT database}}{25}{table.65}
\contentsline {table}{\numberline {Table\nobreakspace 4.5}{\ignorespaces Nearest k-neighbor classfier for Leave One Speaker out Cross Validation in TIMIT database}}{26}{table.66}
\contentsline {table}{\numberline {Table\nobreakspace 4.6}{\ignorespaces Nearest k-neighbor classfier for Leave One Speaker Out Cross validation in NTIMIT database}}{26}{table.67}
\contentsline {table}{\numberline {Table\nobreakspace 3.1}{\ignorespaces Pseudo-code for the dual-tree all nearest neighbor algorithm}}{18}{table.41}
\contentsline {table}{\numberline {Table\nobreakspace 3.2}{\ignorespaces CPU time (in hours) for evaluating all 20 nearest neighbors on TIMIT database}}{22}{table.52}
\contentsline {table}{\numberline {Table\nobreakspace 3.3}{\ignorespaces Nearest k-neighbor classfier for TIMIT database}}{23}{table.53}
\contentsline {table}{\numberline {Table\nobreakspace 3.4}{\ignorespaces Nearest k-neighbor classfier for NTIMIT database}}{23}{table.54}
\contentsline {table}{\numberline {Table\nobreakspace 3.5}{\ignorespaces Nearest k-neighbor classfier for Leave One Speaker out Cross Validation in TIMIT database}}{24}{table.55}
\contentsline {table}{\numberline {Table\nobreakspace 3.6}{\ignorespaces Nearest k-neighbor classfier for Leave One Speaker Out Cross validation in NTIMIT database}}{24}{table.56}
+26 -27
View File
@@ -6,30 +6,29 @@
\BOOKMARK [1][-]{section.11}{Section 2 \204 Origin and history of the problem}{}
\BOOKMARK [1][-]{section.12}{Section 3 \204 Preliminary Research}{}
\BOOKMARK [2][-]{subsection.13}{The Nearest Neighbor Problem}{section.12}
\BOOKMARK [2][-]{subsection.25}{Kernel methods}{section.12}
\BOOKMARK [1][-]{section.37}{Section 4 \204 Proposed Research}{}
\BOOKMARK [2][-]{subsection.38}{Large scale trees, for out of core memory}{section.37}
\BOOKMARK [2][-]{subsection.39}{Memory layout of Kd-trees}{section.37}
\BOOKMARK [2][-]{subsection.40}{System architecture}{section.37}
\BOOKMARK [2][-]{subsection.41}{Memory Manager Architecture}{section.37}
\BOOKMARK [2][-]{subsection.42}{Memory Mapped Files}{section.37}
\BOOKMARK [2][-]{subsection.43}{User defined cache, with TPIE}{section.37}
\BOOKMARK [2][-]{subsection.44}{Customizing Kernels with fast gradient descent}{section.37}
\BOOKMARK [2][-]{subsection.45}{Scalable Support Vector Machines}{section.37}
\BOOKMARK [2][-]{subsection.46}{Computing Gaussian Kernel Matrix with trees}{section.37}
\BOOKMARK [2][-]{subsection.47}{Intrinsic dimension of speech}{section.37}
\BOOKMARK [2][-]{subsection.48}{Speech recognition with nearest neighbor search}{section.37}
\BOOKMARK [2][-]{subsection.49}{Introduction}{section.37}
\BOOKMARK [2][-]{subsection.50}{Speech features}{section.37}
\BOOKMARK [2][-]{subsection.51}{Fast N-body methods for feature comparison}{section.37}
\BOOKMARK [2][-]{subsection.55}{Kernel PCA on a manifold}{section.37}
\BOOKMARK [2][-]{subsection.56}{Experimental Results}{section.37}
\BOOKMARK [2][-]{subsection.57}{All nearest neighbor performance}{section.37}
\BOOKMARK [2][-]{subsection.58}{Comparison of NRAF and MFCC}{section.37}
\BOOKMARK [2][-]{subsection.60}{Discussion}{section.37}
\BOOKMARK [2][-]{subsection.68}{Acoustic Environment adaptation}{section.37}
\BOOKMARK [1][-]{section.69}{Section 5 \204 Work remaining to be done}{}
\BOOKMARK [1][-]{section.70}{Section 6 \204 Facilities and Equipment Needed}{}
\BOOKMARK [1][-]{section.71}{Appendix A \204 This is my first appendix}{}
\BOOKMARK [1][-]{section.72}{Appendix B \204 This is my second appendix}{}
\BOOKMARK [1][-]{section*.73}{References}{}
\BOOKMARK [2][-]{subsection.20}{Large scale trees, for out of core memory}{section.12}
\BOOKMARK [2][-]{subsection.21}{Memory layout of Kd-trees}{section.12}
\BOOKMARK [2][-]{subsection.22}{System architecture}{section.12}
\BOOKMARK [2][-]{subsection.23}{Memory Manager Architecture}{section.12}
\BOOKMARK [2][-]{subsection.24}{Memory Mapped Files}{section.12}
\BOOKMARK [2][-]{subsection.25}{User defined cache, with TPIE}{section.12}
\BOOKMARK [2][-]{subsection.26}{Kernel methods}{section.12}
\BOOKMARK [2][-]{subsection.36}{Forming manifolds from data for clustering through non-parametric density estimation}{section.12}
\BOOKMARK [2][-]{subsection.37}{Speech recognition with nearest neighbor search}{section.12}
\BOOKMARK [2][-]{subsection.38}{Introduction}{section.12}
\BOOKMARK [2][-]{subsection.39}{Speech features}{section.12}
\BOOKMARK [2][-]{subsection.40}{Fast N-body methods for feature comparison}{section.12}
\BOOKMARK [2][-]{subsection.44}{Kernel PCA on a manifold}{section.12}
\BOOKMARK [2][-]{subsection.45}{Experimental Results}{section.12}
\BOOKMARK [2][-]{subsection.46}{All nearest neighbor performance}{section.12}
\BOOKMARK [2][-]{subsection.47}{Comparison of NRAF and MFCC}{section.12}
\BOOKMARK [2][-]{subsection.49}{Discussion}{section.12}
\BOOKMARK [1][-]{section.57}{Section 4 \204 Proposed Research}{}
\BOOKMARK [2][-]{subsection.58}{Customizing Kernels with fast gradient descent}{section.57}
\BOOKMARK [2][-]{subsection.64}{Computing Gaussian Kernel Matrix with trees}{section.57}
\BOOKMARK [2][-]{subsection.65}{Intrinsic dimension of speech}{section.57}
\BOOKMARK [2][-]{subsection.66}{Speech recognition with nearest neighbor search}{section.57}
\BOOKMARK [2][-]{subsection.67}{Acoustic Environment adaptation}{section.57}
\BOOKMARK [1][-]{section.68}{Section 5 \204 Work remaining to be done}{}
\BOOKMARK [1][-]{section.69}{Appendix A \204 This is my first appendix}{}
\BOOKMARK [1][-]{section*.70}{References}{}
+529 -151
View File
@@ -358,143 +358,7 @@ In general the dual tree algorithm as we call it since the query and
the reference set lie on trees (they might share the same tree if
the query and the reference set are the same), gives linear
complexity over the number of data. This is empirical complexity and
\subsubsection{Fast kernel summation with trees}
Kernel summation is on of the most common tasks in kernel methods.
The general formula
\begin{equation}
\label{kernel_sum}
G(x)=\sum_{q=0}^{N} w_i\exp(\frac{(x-x_q)^2}{\sigma^2})
\end{equation}
appears in many cases, with the most interesting case the
eigenvalue/eigenvector computation. The kernel summation is a very
computationally intensive. For Gaussian kernels there are
approximations that can speed up computations orders of magnitude.
We will briefly discuss them here. Besides computation saving we can
save memory as well. The gaussian kernel has infinite support so the
kernel matrix $K, K(x,y)=\exp(\frac{(x-y)^2}{\sigma^2}$ is not
sparse. Numerical iterative methods for eigenvector computation are
based on the general iteration $x_n=x_{n-1}+\lambda Kx_{n-1}$, where
the multiplication of each row of $K$ with $x_{n-1}$ is a kernel
summation problem. So we don't really need to store the kernel
matrix, we just have to compute the kernel sum for every row using
$x_{n-1}$ as a the wight vector.
The gaussian can be approximated with the help of Hermite
polynomials. The Hermite polynomials $H_n(t)$ are defined by the
Rodrigues formula:
\begin{equation}
H_n(t)=(-1)^n\exp(t^2)D^n\exp(-t^2), t\in \Re
\end{equation}
where $D=d/dt$.
The gaussian can now be expressed as
\begin{equation}
\exp^(-(t-s)^2)=\sum_{n=0}^{\inf}\frac{s^n}{n!}h_n(t)
\end{equation}
where $h_n(t)$ is defined by
\begin{equation}
h_n(t)=\exp(-t^2)H_n(t).
\end{equation}
\subsection{Kernel methods}
The kernel or Gram matrix is an informational representation for a
set of data. There are a different interpretations for the kernel
matrix. The most mathematical interpretation is the dot product
matrix. Every element of the kernel matrix is the dot product
between the two data points. It is a positive semidefinite matrix
and it can be considered as the adjacency matrix of a graph where
the points are nodes. This graph can also be considered as an
approximation of the surface that data lie on. Positive
semidefinitness is the only property required for the gram matrix.
In most of the methods that we discuss in this section, we prefer
the interpretation of the dot product matrix, in the general Hilbert
sense. The dot products in most of the cases are non-linear.
\subsubsection{Kernel principal component analysis}
Kernel Principal Component Analysis (KPCA), is an extension of
principal component analysis where the data are mapped into a
different feature space through the kernel. In Principal Component
Analysis the goal is to find the principal components for the
covariance matrix $C$ for a given dataset $S={x_j, x_j\in \Re^d, j=1
\dots M}$, where $C=\frac{1}{M}\sum_{j=1}^{M}x_j x_j^T$. If we
define a mapping
\begin{equation}
\Phi : x \rightarrow X , x \in \Re^d, X\in \Re^f
\end{equation}
usually $f>>d$. So usually the nonlinear mapping, maps the data in a
higher dimensional space sometimes in an infinite dimensional space.
At this point it should be clarified that the mapping increases the
extrinsic dimension of the data not the intrinsic. In kernel PCA we
want to find the principal components of the covariance matrix
\begin{equation}
\tilde{C} = \frac{1}{M}\sum_{j=1}^{M}\phi(x_j)\phi(x_j)^T
\end{equation}
It turns out that the mapping $\Phi$ doesn't have to be known
explicitly, all we need is the kernel that represents the dot
product.
\begin{equation}
k(x, y) = \phi(x)\phi(y)^T
\end{equation}
For some kernels there is an analytic factorization while for others
there is not. Not every function can be a valid kernel (dot
product). There are conditions that must be fulfilled. Basically for
every dataset they must give a valid kernel matrix as described in
previous sections.
Kernel PCA has many cousins such as Laplacian eigenmaps, Diffusion
maps Isomap, Local Linear embedding
\subsubsection{Customized Kernels}
As we described in the previous section there are infinite functions
that can act as valid kernels. The most popular kernels are:
\begin{enumerate}
\item The gaussian kernel $k(x, y)=e^{-\frac{||x-y||^2}{\sigma^2)}}$
\item The polynomial kernel $k(x, y)=(xy)^n$
\item The epanechnikov kernel $k(x, y)=1-||x-y||^2/\sigma^2, ||x-y||<\sigma, else 0$
\item The k-nearest neighbor kernel
\end{enumerate}
Although these are valid kernels for any kernel method including
Kernel PCA they are not appropriate for every task. In some cases we
want to map data in higher dimensional spaces preserving distances
(isometry) while in other cases we want to make small distances
smaller and send high distances to infinity (clustering). The
gaussian kernel for example tends to behave better in clustering
since it fades out very quickly. On the other hand it fails to
unfold the data manifold so that the intrinsic dimension can be
estimated. Manifold unfolding is a very important procedure since it
can reveal the true dimension of the dataset.
In [?] the authors introduced the idea of building the kernel matrix
from scratch without using any fixed kernel. As we have already
mentioned the kernel matrix has to be a positive semidefinite
matrix. Initially we have to define a k neighborhood for every
point. So at first we have to connect points to their neighbors and
build a initialize the non zero values of the matrix to 1. Then we
maximize the trace of the kernel matrix under the constraint that
the the distances in the kernel matrix are preserved. This leads to
a semidefinite programming problem.
\subsubsection{Stochastic optimization}
Stochastic optimization or on line learning. On line svd is one key
aspect.
\subsubsection{Support Vector Machines}
\section{Proposed Research}
\label{proposed}
The goal of this thesis is to make kernel methods scalable.
Efficient computation of the kernel matrix for large amounts of
data. Special emphasis is given in high dimensional data $d>30$. One
of the goals is to use Kernel PCA for unfolding the speech manifold.
In other words estimate the intrinsic dimensionality of speech by
using big databases like timit, wall street journal and broadcast
news. Implicitly dimensionality reduction will lead to new feature
generation we will explain it in the following sections. The
ultimate goal is to do speech recognition by nearest neighbor
search.
\subsection{Large scale trees, for out of core memory}
Kd-trees have been proven to be be quite successful in
multidimensional indexing and turn out outperform several other
@@ -672,20 +536,461 @@ to address efficiently this. Pages are write protected in memory and
if a write occurs a signal is raised and the page is marked as
modified.
\subsection{Customizing Kernels with fast
gradient descent}
As we described in section ? it is possible to customize the kernel
for a given dataset, through semidefinite programming. Although this
method is very promising and gives very good results it is not
scalable. Semidefinite programming is solved with the intrinsic
point method which is of polynomial complexity.
\subsection{Scalable Support Vector Machines}
\subsection{Computing Gaussian Kernel Matrix with trees}
\subsection{Intrinsic dimension of speech}
\subsection{Kernel methods}
The kernel or Gram matrix is an informational representation for a
set of data. There are a different interpretations for the kernel
matrix. The most mathematical interpretation is the dot product
matrix. Every element of the kernel matrix is the dot product
between the two data points. It is a positive semidefinite matrix
and it can be considered as the adjacency matrix of a graph where
the points are nodes. This graph can also be considered as an
approximation of the surface that data lie on. Positive
semidefinitness is the only property required for the gram matrix.
In most of the methods that we discuss in this section, we prefer
the interpretation of the dot product matrix, in the general Hilbert
sense. The dot products in most of the cases are non-linear.
\subsubsection{Kernel principal component analysis}
Kernel Principal Component Analysis (KPCA), is an extension of
principal component analysis where the data are mapped into a
different feature space through the kernel. In Principal Component
Analysis the goal is to find the principal components for the
covariance matrix $C$ for a given dataset $S={x_j, x_j\in \Re^d, j=1
\dots M}$, where $C=\frac{1}{M}\sum_{j=1}^{M}x_j x_j^T$. If we
define a mapping
\begin{equation}
\Phi : x \rightarrow X , x \in \Re^d, X\in \Re^f
\end{equation}
usually $f>>d$. So usually the nonlinear mapping, maps the data in a
higher dimensional space sometimes in an infinite dimensional space.
At this point it should be clarified that the mapping increases the
extrinsic dimension of the data not the intrinsic. In kernel PCA we
want to find the principal components of the covariance matrix
\begin{equation}
\tilde{C} = \frac{1}{M}\sum_{j=1}^{M}\phi(x_j)\phi(x_j)^T
\end{equation}
It turns out that the mapping $\Phi$ doesn't have to be known
explicitly, all we need is the kernel that represents the dot
product.
\begin{equation}
k(x, y) = \phi(x)\phi(y)^T
\end{equation}
For some kernels there is an analytic factorization while for others
there is not. Not every function can be a valid kernel (dot
product). There are conditions that must be fulfilled. Basically for
every dataset they must give a valid kernel matrix as described in
previous sections.
Kernel PCA has many cousins such as Laplacian eigenmaps, Diffusion
maps Isomap, Local Linear embedding
\subsubsection{Customized Kernels}
As we described in the previous section there are infinite functions
that can act as valid kernels. The most popular kernels are:
\begin{enumerate}
\item The gaussian kernel $k(x, y)=e^{-\frac{||x-y||^2}{\sigma^2)}}$
\item The polynomial kernel $k(x, y)=(xy)^n$
\item The epanechnikov kernel $k(x, y)=1-||x-y||^2/\sigma^2, ||x-y||<\sigma, else 0$
\item The k-nearest neighbor kernel
\end{enumerate}
Although these are valid kernels for any kernel method including
Kernel PCA they are not appropriate for every task. In some cases we
want to map data in higher dimensional spaces preserving distances
(isometry) while in other cases we want to make small distances
smaller and send high distances to infinity (clustering). The
gaussian kernel for example tends to behave better in clustering
since it fades out very quickly. On the other hand it fails to
unfold the data manifold so that the intrinsic dimension can be
estimated. Manifold unfolding is a very important procedure since it
can reveal the true dimension of the dataset.
In [?] the authors introduced the idea of building the kernel matrix
from scratch without using any fixed kernel. As we have already
mentioned the kernel matrix has to be a positive semidefinite
matrix. Initially we have to define a k neighborhood for every
point. So at first we have to connect points to their neighbors and
build a initialize the non zero values of the matrix to 1. Then we
maximize the trace of the kernel matrix under the constraint that
the the distances in the kernel matrix are preserved. This leads to
a semidefinite programming problem.
\subsection{Parameter estimation for manifold learning, through density estimation}
Manifold learning turns out to be a very useful tool for many
applications of machine learning, such as classification.
Unfortunately the existing algorithms use ad hoc selection of the
parameters that define the geometry of the manifold. The parameter
choice affects significantly the performance of manifold learning
algorithms. Recent theoretical work has proven the equivalence
between the Mercer kernel learning methods and the kernel in kernel
density estimation. Based on this fact the problem of kernel
parameter estimation for manifold learning is addressed based on the
nonparametric statistical theory estimation. An automatic way of
determining the local bandwidths that define the geometry is
introduced. The results show that the automatic bandwidth selection
leads to improved clustering performance and reduces the
computational load versus ad hoc selection.
\subsection{INTRODUCTION}
Manifold learning has been active recently through the work of
\cite{Lafon} and \cite{Belkin}. Both works present algorithms for
computing the Laplace-Beltrami operator on a submanifold $\Gamma$
embedded in $\Re^d$. Lafon has also introduced a diffusion process
on the manifold based on a kernel $k(x,y)$ that defines the local
geometry on the points that belong to the manifold. The choice of
the kernel affects essentially the results of the diffusion process,
since it affects the notion of the neighborhood around the points
which is critical for the creation of the proximity graph of the
manifold. In their work Belkin and Lafon don't deal with the problem
of training the local bandwidth of the kernel, assuming that the
manifold has been sufficiently sampled. Jenssenl \cite{jenssen17lpd}
have proven that there is theoretical equivalence between the
Mercer kernel and the kernel in kernel density estimation. Based on
that fact we performed several experiments using the optimal
bandwidth for kernel density estimation. In this paper the Adaptive
Kernel Based Density Estimation Algorithm (AKDEA) \cite{Silverman}
was chosen. The proposed algorithm is a partial solution because it
doesn't deal with the final objective function directly, which in
some applications of manifold learning it could be clustering or
dimensionality reduction.
The connection between manifold
learning and Kernel Density Estimation (KDE) has also been addressed
in \cite{girolami2002osd} where the connection between manifold
learning and density estimation was also considered by using
orthogonal series density estimation rather than kernel density
estimation
In section 2 the diffusion operator is reviewed based on the work
of Lafon \cite{Lafon}. Section 3 outlines the adaptive kernel
estimation algorithm. The last section presents examples of the
diffusion operator when the kernel bandwidth is trained.
\subsection{Geometric diffusion on a manifold} Let $(\Gamma,\mu)$ be
a measure space, where $\Gamma$ is a finite set of $d$-dimensional
points and $\mu$ is a counting measure that represents the
distribution of the points on the data set. In other words $\Gamma$
is a submanifold of $\Re^{d}$. Assume that the geometry of $\Gamma$
is defined by a kernel $k(x,y)$. The kernel $k(x,y)$ measures the
degree of similarity between two points $x,y$. The kernel satisfies
the following conditions:
\begin{itemize}
\item $k$ is symmetric: $k(x,y)=k(y,x)$,
\item $k$ is positivity-preserving: for all $x,y$ in
$\Gamma, k(x,y)\geq 0$
\item $k$ is positive semi-definite: for all bounded
functions $f$ defined on $\Gamma$,
\[
\int_{\Gamma} \int_{\Gamma} k(x,y)f(x)f(y)d\mu(x)d\mu(y)\geq 0\]
\end{itemize}
In this paper it is assumed that $\Gamma$ is a subset of the
Euclidean space $\Re^{d}$. So for $x,y\in\Gamma$ the kernel
(similarity measure) is a function of the Euclidean distance
$\parallel x-y\parallel:$
\[
k(x,y)=\eta\left(\frac{\parallel x-y\parallel}{h}\right)
\]
In order to study the geometry of the submanifold $\Gamma$, an
oriented graph $G$ is formed. Every node corresponds to a data
point. Let $b(x,y)$ be the associated adjacency matrix, where
$b(x,y)=1$ if $x$ is in the neighborhood of $y$ and $b(x,y)=0$ if
$x$ is not in the neighborhood of $y$. The kernel $k$ defines the
notion of neighborhood between the points.
Let $u^2(x)=\int_{\Gamma}k(x,y)d\mu(y)$, then the normalized kernel
$a(x,y)=\frac{k(x,y)}{u(x)u(y)}$ is stochastic since:
\begin{equation}\label{kernel}
\int_{\Gamma}a(x,y)d\mu(y)=1.
\end{equation}
The kernel $a(x,y)$ can be considered as a transition matrix of a
Markov process on the submanifold $\Gamma$. So any admissible kernel
can be associated with a random walk on $\Gamma$. The operator
\[
Af(x)=\int_{\Gamma}a(x,y)f(y)d\mu(y)
\]
is called a diffusion operator. The eigenvalue analysis of the graph
with the adjacency matrix $a(x,y)$ is the basis for embedding the
data to a lower dimensional space.
\[
a(x,y)=\sum_{j=0}^{M} \lambda_{j}\phi_{j}(x)\phi_{j}(y)
\]
where $M$ is the number of data points, $x,y$ are nodes of the graph
and $\phi_j$ is the $j_{th}$ eigenvector of $a$ that corresponds to
the $\lambda_j$ eigenvalue.
\section{Non parametric density estimation}
The problem of density estimation of finite data set is very
significant for most of the well known machine learning algorithms.
Lafon has proven that the Laplace Beltrami operator on a manifold
can be computed without the knowledge of the density. Although the
density is not necessary for the computation of the operator it can
be useful in the estimation of the parameters of the kernel
[\ref{kernel}].
The density estimate $f$ of a finite set of data $X={x_i\ni
\Re^n,i=1,\dots,M}$ can be expressed with the following equation
\cite{Silverman}:
\begin{equation}\label{powerdensity}
f(x)=\frac{1}{M}\sum_{i=1}^{M}\frac{1}{h^d}k\left(\frac{\parallel x-x_i
\parallel}{h}\right)
\end{equation}
The shape of the density function varies significantly with
different values of $h$, which is the global bandwidth. Higher
values of $h$ lead to smoothed versions of the density. An
automatic way to determine the optimum value of $h$ is to maximize
the Leave One Out Cross Validation (LOOCV) criterion. In other words
the density is formed according to (\ref{powerdensity}) by leaving
one point out and then the probability of this point is evaluated
according to the density. The LOOCV can be mathematically expressed:
\begin{equation}
M(h)=\sum_{i=1}^{M}\log\sum_{i\neq j}f_j(x_i)
\end{equation}
where $f_j$ is the density computed according to
(\ref{powerdensity}) by leaving $x_j$ out.
Another restriction that can be imposed to the estimated density is
the minimization of the power of the second derivative of the
density $\int f''(x)^2$. This restriction tends to give smoother
densities and avoids the formation of isolated islands of points.
With the standard optimization techniques $h$ can be chosen.
The value of $h$ defines the neighborhood around a point through the
kernel $k$. For dense areas $h$ has to be small, while in sparse
areas $h$ must be high since the points are far away from each
other. This leads to the conclusion that $h$ should be different for
every point. The density can be written in the form:
\[
f(x)=\frac{1}{M}\sum_{i=1}^{M}\frac{1}{(\lambda_i h)^d}k(\frac{\parallel x-x_i
\parallel}{\lambda_i h})
\]
The $\lambda_i$ are called local bandwidths. The local bandwidths
can be estimated from the following algorithm:
\begin{itemize}
\item Find a pilot estimate $\tilde{f}(x)$ that satisfies
$\tilde{f(x_i)>0, i=1,\dots,M}$
\item The local bandwidth $\lambda_i$ is given by:
\[
\lambda_i=(\tilde{f}(x_i)/g)^{\frac{1}{d}}
\]
where g is the geometric mean of the $\tilde{f}(x_i)$:
\[
g=\sqrt[M]{\prod_{i=1}^{M}\tilde{f}(x_i)}
\]
\item The adaptive kernel estimate $f(x)$ is defined by:
\[
f(x)=\frac{1}{M}\sum_{i=1}^{M}\frac{1}{(\lambda_i
h)^d} k\left( \frac{\parallel x-x_i \parallel}{\lambda_i h} \right)
\]
\end{itemize}
The pilot estimate can be any ad hoc density estimate. A common
estimate is the \emph{kth nearest neighbor estimate} defined by:
\[
\tilde{f}(x)=\frac{1}{Md_k(x)}\sum_{i=1}^{M}k\left(\frac{
\parallel x-x_i \parallel}{d_k(x)}\right)
\]
where $d_k(x)$ is the distance between $x$ and the $kth$ nearest
neighbor. It turns out that the choice of $k$ doesn't affect the
values of $\lambda_i$. The adaptive kernel estimate can be improved
if we iterate the above procedure by replacing in every iteration
the pilot density estimate with the adaptive kernel estimate
computed in the last step.
\section{Eigenvalue analysis of the diffusion graph}
In the previous section we saw that the non-parametric estimation
can define the local bandwidth for every point. This information
can be used in the construction of the adjacency matrix $G$. The
local bandwidths as derived from the previous section will usually
give a graph that has connected components, or equivalently the
manifold is connected. A safe test to prove that is to look into
the eigenvalues also known as the spectrum of the graph. According
to \cite{Cvetkovic} the multiplicity of the first eigenvalue
reflects the number of disconnected components in the graph. There
are cases where the local bandwidths might not give a connected
manifold, then it is necessary to increase $h$ until the
multiplicity of the first eigenvalue becomes 1.
\subsection{Kernel choice}
There are many classes of kernels that can be used for the diffusion
operator on the manifold. The most popular is the Gaussian kernel:
\[
k(x)=\frac{1}{(2\pi)^{\frac{-d}{2}}}e^{(\frac{-1}{2}xx^T)}
\]
The problem with the gaussian kernel is that it is of infinite
support and it leads to adjacency matrices that are not sparse. It
would be preferable to use a kernel of finite support. The most
popular one is the Epanechnikov kernel \cite{Silverman}:
\begin{equation}
k(x)=\left\{
\begin{array}{ll}
\frac{1}{2}c_{d}^{-1}(d+2)(1-x^Tx) & x^Tx<1 \\
0 & x^Tx\geq 1 \\
\end{array}
\right.
\end{equation}
\subsection{Examples}
In this section two examples illustrated in the two dimensional
space are presented and a third one on real speech data. The first
example fig.~\ref{fig1} contains two clusters of data. Each cluster
contains 300 points uniformly distributed in a unity area square.
The Average Nearest Neighbor Distance (ANND) as a choice for the
bandwidth for the Epanechnikov kernel turns out to be inadequate
since it leads to a disconnected graph with 4 disjoint components.
This is probably because the two classes are sparsely sampled. The
AKDE based on maximization of LOOCV gives on the average local
bandwidths that are an order of magnitude higher than the ANND, and
lead to a smoother density. As it is shown in fig.~\ref{fig1}b the
density is not flat as it was expected to be and this is because of
the fact that the number of points that sample the manifold is
small. In fig.~\ref{fig1}c the first non trivial eigenvector is
depicted discriminating very sharply the two classes. In
fig.~\ref{fig1.1} the same data set is analyzed with a higher global
bandwidth. The results are very poor since there is high intra-class
deviation and the centroids of the two classes are closer.
The second example fig.~\ref{fig2} illustrates two classes of data
that lie on two halves of a circle, slightly separated. The clusters
are highly the sampled around the angle of $90^0$ degrees. This
fact makes very difficult to chose a universal value for the local
bandwidths since the notion of the locality is significantly
different in different areas of the manifold. In fig.~\ref{fig2}b
the density of the data points is shown, reflecting the high
concentration of points around the angle of $90^0$. In this case the
AKDE trains the local bandwidths $\lambda_i$ and the global $h$, so
that the LOOCV is maximized. The resulting values give a manifold
graph that is disconnected in two classes. A slight increase in the
global bandwidth makes the graph connected. The first non-trivial
eigenvector of the diffusion operator identifies very clearly the
two classes as shown in fig.~\ref{fig2}c. If the global bandwidth is
increased by an order of magnitude, then still the first non-trivial
eigenvector identifies the two classes, but the classes are not
compact and their centroids are closer fig.\ref{fig3}c. Another
disadvantage of keeping $h$ large is that the adjacency matrix
$a(x,y)$ fig.~\ref{fig5}b of the manifold graph becomes dense. This
increases tremendously the memory requirement and the computation of
the eigenvectors.
In this third example the effect of the local bandwidth on the
clustering problem of three phoneme classes is presented. The points
correspond to the Fast Fourier Transform of 20 msec speech frames
sampled at 8KHz, of the TIMIT database. 200 points from the three
phoneme classes /aa/,/ih/,/sh/ where chosen. In fig.~\ref{fig4}a the
data are plotted in 2 dimensions after the dimensionality reduction
process. It is clear that the diffusion process separates the
classes very well. In fig.~\ref{fig4}b where the bandwidth is an
order of magnitude larger than the optimal, the diffusion process
still separates the classes, but the clusters are not very compact.
When the local bandwidth is an order of magnitude smaller than the
optimal, the results are very poor. The corresponding graph
fig.~\ref{fig5}c shows that the points are isolated and not
connected.
\section{Discussion}
In this paper the problem of constructing the diffusion graph of a
manifold defined by data points in $\Re^d$ was addressed. The
parameters of the kernel that define the geometry over the manifold
were automatically trained through the estimation of the adaptive
kernel based density algorithm. Experiments on synthetic and real
data showed that if the global bandwidth is smaller than the one
that maximizes the LOOCV, then the graph that represents the
manifold is disconnected into many spurious clusters. On the other
hand if the bandwidth is higher, then the diffusion operator will
still identify the clusters, but their distance is smaller and they
are not as compact compared to the optimal bandwidth. Moreover high
bandwidth lead to dense graphs and as a consequence more
computations.
The approach followed is still heuristic, although it has an
underlying theoretical justification on the work described at
\cite{jenssen17lpd}. It is the first time that experimental results
are presented that show the effect of optimal kernel density
estimation on the manifold learning.
In terms of algorithmic
complexity, most of the operations performed are $O(n^2)$ with the
naive approach. Both kernel density estimation and graph formation
fall under the general case of N-body problem computations. The
speed can be significantly improved by using KD-trees as described
in \cite{gray2001nbp} and reduce the complexity up to $O(n)$.
\begin{figure}[h!]
\centerline{\includegraphics[height=4.5cm]{fig1.1b.eps}(a)}
\centerline{\includegraphics[height=4.5cm]{fig1.1c.eps}(b)}
\caption{The same data with fig.~\ref{fig1}a analyzed with a
larger global bandwidth. (a) The estimated density, (b) The first
non-trivial eigenvector }
\label{fig1.1}
\end{figure}
\bibliographystyle{IEEEbib}
\bibliography{main}
\begin{figure}[tb]
\centerline{\includegraphics[height=4.4cm]{fig1a.eps}(a)}
\centerline{\includegraphics[height=4.4cm]{fig1b.eps}(b)}
\centerline{\includegraphics[height=4.4cm]{fig1c.eps}(c)}
\caption{(a) Two clusters of uniformly distributed 2-D points,
(b) The estimated density, (c) The first non-trivial eigenvector }
\label{fig1}
\end{figure}
\begin{figure}[tb]
\centerline{\includegraphics[height=6.0cm]{fig2a.eps}(a)}
\centerline{\includegraphics[height=6.0cm]{fig2b.eps}(b)}
\centerline{\includegraphics[height=6.0cm]{fig2c.eps}(c)}
\caption{(a) Two clusters of non-uniformly distributed 2-D points,
(b)The estimated density, (c) The first non-trivial eigenvector
\vspace{1.2cm}}
\label{fig2}
\end{figure}
\begin{figure}[tb]
\centerline{\includegraphics[height=6cm]{fig3a.eps}(a)}
\centerline{\includegraphics[height=6cm]{fig3b.eps}(b)}
\centerline{\includegraphics[height=6cm]{fig3c.eps}(c)}
\caption{The same data set with fig.~\ref{fig2} processed with a
bandwidth one order of magnitude larger. Although the first
eigenvector still separates the two classes, the classes are closer
and not compact (a) Two clusters of non-uniformly distributes 2-D
points, (b)The estimated density, (c) The first non trivial
eigenvector }
\label{fig3}
\end{figure}
\begin{figure}[tb]
\centerline{\includegraphics[height=6cm]{fig4a.eps}(a)}
\centerline{\includegraphics[height=6cm]{fig4c.eps}(b)}
\centerline{\includegraphics[height=6cm]{fig4e.eps}(c)}
\caption{(a) The 3 phoneme classes after the dimensionality reduction
with the optimal bandwidth, plotted in two dimensions, (b) The 3
phoneme classes with a larger bandwidth,(c) The 3 phoneme classes
with a smaller bandwidth}
\label{fig4}
\end{figure}
\begin{figure}[tb]
\centerline{\includegraphics[height=6cm]{fig4b.eps}(a)}
\centerline{\includegraphics[height=6cm]{fig4d.eps}(b)}
\centerline{\includegraphics[height=6cm]{fig4f.eps}(c)}
\caption{The diffusion graph for three different local bandwidths
(a)Optimal, (b)An order of magnitude larger, (c)An order of
magnitude smaller}
\label{fig5}
\end{figure}
\subsection{Speech recognition with nearest neighbor search}
In this paper we present a computationally efficient means of
@@ -1043,18 +1348,91 @@ MFCC features}{8cm}
} \label{kneighborLOOCVNTIMIT}
\end{table}
\section{Proposed Research}
\label{proposed}
The goal of this thesis is to make kernel methods scalable.
Efficient computation of the kernel matrix for large amounts of
data. Special emphasis is given in high dimensional data $d>30$. One
of the goals is to use Kernel PCA for unfolding the speech manifold.
In other words estimate the intrinsic dimensionality of speech by
using big databases like timit, wall street journal and broadcast
news. Implicitly dimensionality reduction will lead to new feature
generation we will explain it in the following sections. The
ultimate goal is to do speech recognition by nearest neighbor
search.
\subsection{Customizing Kernels with fast
gradient descent}
As we described in section ? it is possible to customize the kernel
for a given dataset, through semidefinite programming. Although this
method is very promising and gives very good results it is not
scalable. Semidefinite programming is solved with the intrinsic
point method which is of polynomial complexity.
\subsubsection{Fast kernel summation with trees}
Kernel summation is on of the most common tasks in kernel methods.
The general formula
\begin{equation}
\label{kernel_sum}
G(x)=\sum_{q=0}^{N} w_i\exp(\frac{(x-x_q)^2}{\sigma^2})
\end{equation}
appears in many cases, with the most interesting case the
eigenvalue/eigenvector computation. The kernel summation is a very
computationally intensive. For Gaussian kernels there are
approximations that can speed up computations orders of magnitude.
We will briefly discuss them here. Besides computation saving we can
save memory as well. The gaussian kernel has infinite support so the
kernel matrix $K, K(x,y)=\exp(\frac{(x-y)^2}{\sigma^2}$ is not
sparse. Numerical iterative methods for eigenvector computation are
based on the general iteration $x_n=x_{n-1}+\lambda Kx_{n-1}$, where
the multiplication of each row of $K$ with $x_{n-1}$ is a kernel
summation problem. So we don't really need to store the kernel
matrix, we just have to compute the kernel sum for every row using
$x_{n-1}$ as a the wight vector.
The gaussian can be approximated with the help of Hermite
polynomials. The Hermite polynomials $H_n(t)$ are defined by the
Rodrigues formula:
\begin{equation}
H_n(t)=(-1)^n\exp(t^2)D^n\exp(-t^2), t\in \Re
\end{equation}
where $D=d/dt$.
The gaussian can now be expressed as
\begin{equation}
\exp^(-(t-s)^2)=\sum_{n=0}^{\inf}\frac{s^n}{n!}h_n(t)
\end{equation}
where $h_n(t)$ is defined by
\begin{equation}
h_n(t)=\exp(-t^2)H_n(t).
\end{equation}
\subsection{Computing Gaussian Kernel Matrix with trees}
Direct computation of the kernel(Gaussian) kernel leads to a
non-sparse matrix. There is a way though to compute the
inverse/eigenvalues without explicitly computing the kernel matrix.
Several optimizations can be done for high dimensions
\subsection{Intrinsic dimension of speech}
\subsection{Speech recognition with nearest neighbor search}
Implement a merge reduce scheme to create different descriptions of
speech. Use nearest neighbor algorithms to replace gaussian mixtures
in hidden markov models.
\subsection{Acoustic Environment adaptation}
\section{Work remaining to be done}
\label{remains}
\section{Facilities and Equipment Needed}
\label{facitlies}
\end{Body}
\begin{EndMatter}
\appendix{This is my first appendix}
\appendix{This is my second appendix}
%% Below are some optional environments that can be used here
\bibliographystyle{ieeetr}
\references % Generates the bibliography page
+191 -150
View File
@@ -156,7 +156,37 @@ literature dimensionality is used without the prefix intrinsic or
extrinsic, causing confusion. The extrinsic dimensionality is just
the number of variables needed to describe a point, i.e. the
extrinsic dimensionality of $x\in \Re^d$ is just $d$. Each variable
of the d-dimensional vector will be called dimension $1\dots d$
of the d-dimensional vector will be called dimension $1\dots d$. In
reality some of the dimensions might be linearly or non-linearly
dependent on other dimensions. So less than $d$ dimensions are
sufficient to describe a point. We can define the intrinsic
dimensionality of a dataset as the minimum number of "independent"
variables needed to describe a point in space. Still the term is not
very clear, for the scope of this thesis we mean that there is no
algorithm that can predict a dimension from the others. Principal
Component Analysis for example finds the most important dimensions
on a dataset that are linearly independent. In practice because of
noise all dimensions appear to be strictly independent. Modern
kernel methods can capture non-lineal dependencies between
dimensions. Kernel PCA is the father of non-linear dimensionality
reduction methods.
Although Kernel PCA seems to be a very powerful tool it suffers from
polynomial complexity. It requires $O(N^2)$ kernel computations,
which make it unsuitable for large scale problems. In this thesis we
will show how we are going to use multidimensional trees to reduce
the complexity, along with other approximations to bring it down to
linear.
The driving force for investigating dimensionality reduction was the
estimation of the true dimension of speech. Many representations
have been recommended for speech MFCC, NRAFF, LPC, etc. The true
dimension is yet unknown. Meanwhile there is theoretical evidence
that speech signals are embedded in high dimensional manifolds. In
this thesis we will attempt to estimate the true dimensionality of
speech based on several representations. We will see in the future
sections how this can help in developing new directions in speech
recognition.
@@ -328,143 +358,7 @@ In general the dual tree algorithm as we call it since the query and
the reference set lie on trees (they might share the same tree if
the query and the reference set are the same), gives linear
complexity over the number of data. This is empirical complexity and
\subsubsection{Fast kernel summation with trees}
Kernel summation is on of the most common tasks in kernel methods.
The general formula
\begin{equation}
\label{kernel_sum}
G(x)=\sum_{q=0}^{N} w_i\exp(\frac{(x-x_q)^2}{\sigma^2})
\end{equation}
appears in many cases, with the most interesting case the
eigenvalue/eigenvector computation. The kernel summation is a very
computationally intensive. For Gaussian kernels there are
approximations that can speed up computations orders of magnitude.
We will briefly discuss them here. Besides computation saving we can
save memory as well. The gaussian kernel has infinite support so the
kernel matrix $K, K(x,y)=\exp(\frac{(x-y)^2}{\sigma^2}$ is not
sparse. Numerical iterative methods for eigenvector computation are
based on the general iteration $x_n=x_{n-1}+\lambda Kx_{n-1}$, where
the multiplication of each row of $K$ with $x_{n-1}$ is a kernel
summation problem. So we don't really need to store the kernel
matrix, we just have to compute the kernel sum for every row using
$x_{n-1}$ as a the wight vector.
The gaussian can be approximated with the help of Hermite
polynomials. The Hermite polynomials $H_n(t)$ are defined by the
Rodrigues formula:
\begin{equation}
H_n(t)=(-1)^n\exp(t^2)D^n\exp(-t^2), t\in \Re
\end{equation}
where $D=d/dt$.
The gaussian can now be expressed as
\begin{equation}
\exp^(-(t-s)^2)=\sum_{n=0}^{\inf}\frac{s^n}{n!}h_n(t)
\end{equation}
where $h_n(t)$ is defined by
\begin{equation}
h_n(t)=\exp(-t^2)H_n(t).
\end{equation}
\subsection{Kernel methods}
The kernel or Gram matrix is an informational representation for a
set of data. There are a different interpretations for the kernel
matrix. The most mathematical interpretation is the dot product
matrix. Every element of the kernel matrix is the dot product
between the two data points. It is a positive semidefinite matrix
and it can be considered as the adjacency matrix of a graph where
the points are nodes. This graph can also be considered as an
approximation of the surface that data lie on. Positive
semidefinitness is the only property required for the gram matrix.
In most of the methods that we discuss in this section, we prefer
the interpretation of the dot product matrix, in the general Hilbert
sense. The dot products in most of the cases are non-linear.
\subsubsection{Kernel principal component analysis}
Kernel Principal Component Analysis (KPCA), is an extension of
principal component analysis where the data are mapped into a
different feature space through the kernel. In Principal Component
Analysis the goal is to find the principal components for the
covariance matrix $C$ for a given dataset $S={x_j, x_j\in \Re^d, j=1
\dots M}$, where $C=\frac{1}{M}\sum_{j=1}^{M}x_j x_j^T$. If we
define a mapping
\begin{equation}
\Phi : x \rightarrow X , x \in \Re^d, X\in \Re^f
\end{equation}
usually $f>>d$. So usually the nonlinear mapping, maps the data in a
higher dimensional space sometimes in an infinite dimensional space.
At this point it should be clarified that the mapping increases the
extrinsic dimension of the data not the intrinsic. In kernel PCA we
want to find the principal components of the covariance matrix
\begin{equation}
\tilde{C} = \frac{1}{M}\sum_{j=1}^{M}\phi(x_j)\phi(x_j)^T
\end{equation}
It turns out that the mapping $\Phi$ doesn't have to be known
explicitly, all we need is the kernel that represents the dot
product.
\begin{equation}
k(x, y) = \phi(x)\phi(y)^T
\end{equation}
For some kernels there is an analytic factorization while for others
there is not. Not every function can be a valid kernel (dot
product). There are conditions that must be fulfilled. Basically for
every dataset they must give a valid kernel matrix as described in
previous sections.
Kernel PCA has many cousins such as Laplacian eigenmaps, Diffusion
maps Isomap, Local Linear embedding
\subsubsection{Customized Kernels}
As we described in the previous section there are infinite functions
that can act as valid kernels. The most popular kernels are:
\begin{enumerate}
\item The gaussian kernel $k(x, y)=e^{-\frac{||x-y||^2}{\sigma^2)}}$
\item The polynomial kernel $k(x, y)=(xy)^n$
\item The epanechnikov kernel $k(x, y)=1-||x-y||^2/\sigma^2, ||x-y||<\sigma, else 0$
\item The k-nearest neighbor kernel
\end{enumerate}
Although these are valid kernels for any kernel method including
Kernel PCA they are not appropriate for every task. In some cases we
want to map data in higher dimensional spaces preserving distances
(isometry) while in other cases we want to make small distances
smaller and send high distances to infinity (clustering). The
gaussian kernel for example tends to behave better in clustering
since it fades out very quickly. On the other hand it fails to
unfold the data manifold so that the intrinsic dimension can be
estimated. Manifold unfolding is a very important procedure since it
can reveal the true dimension of the dataset.
In [?] the authors introduced the idea of building the kernel matrix
from scratch without using any fixed kernel. As we have already
mentioned the kernel matrix has to be a positive semidefinite
matrix. Initially we have to define a k neighborhood for every
point. So at first we have to connect points to their neighbors and
build a initialize the non zero values of the matrix to 1. Then we
maximize the trace of the kernel matrix under the constraint that
the the distances in the kernel matrix are preserved. This leads to
a semidefinite programming problem.
\subsubsection{Stochastic optimization}
Stochastic optimization or on line learning. On line svd is one key
aspect.
\subsubsection{Support Vector Machines}
\section{Proposed Research}
\label{proposed}
The goal of this thesis is to make kernel methods scalable.
Efficient computation of the kernel matrix for large amounts of
data. Special emphasis is given in high dimensional data $d>30$. One
of the goals is to use Kernel PCA for unfolding the speech manifold.
In other words estimate the intrinsic dimensionality of speech by
using big databases like timit, wall street journal and broadcast
news. Implicitly dimensionality reduction will lead to new feature
generation we will explain it in the following sections. The
ultimate goal is to do speech recognition by nearest neighbor
search.
\subsection{Large scale trees, for out of core memory}
Kd-trees have been proven to be be quite successful in
multidimensional indexing and turn out outperform several other
@@ -642,19 +536,93 @@ to address efficiently this. Pages are write protected in memory and
if a write occurs a signal is raised and the page is marked as
modified.
\subsection{Customizing Kernels with fast
gradient descent}
As we described in section ? it is possible to customize the kernel
for a given dataset, through semidefinite programming. Although this
method is very promising and gives very good results it is not
scalable. Semidefinite programming is solved with the intrinsic
point method which is of polynomial complexity.
\subsection{Scalable Support Vector Machines}
\subsection{Computing Gaussian Kernel Matrix with trees}
\subsection{Intrinsic dimension of speech}
\subsection{Kernel methods}
The kernel or Gram matrix is an informational representation for a
set of data. There are a different interpretations for the kernel
matrix. The most mathematical interpretation is the dot product
matrix. Every element of the kernel matrix is the dot product
between the two data points. It is a positive semidefinite matrix
and it can be considered as the adjacency matrix of a graph where
the points are nodes. This graph can also be considered as an
approximation of the surface that data lie on. Positive
semidefinitness is the only property required for the gram matrix.
In most of the methods that we discuss in this section, we prefer
the interpretation of the dot product matrix, in the general Hilbert
sense. The dot products in most of the cases are non-linear.
\subsubsection{Kernel principal component analysis}
Kernel Principal Component Analysis (KPCA), is an extension of
principal component analysis where the data are mapped into a
different feature space through the kernel. In Principal Component
Analysis the goal is to find the principal components for the
covariance matrix $C$ for a given dataset $S={x_j, x_j\in \Re^d, j=1
\dots M}$, where $C=\frac{1}{M}\sum_{j=1}^{M}x_j x_j^T$. If we
define a mapping
\begin{equation}
\Phi : x \rightarrow X , x \in \Re^d, X\in \Re^f
\end{equation}
usually $f>>d$. So usually the nonlinear mapping, maps the data in a
higher dimensional space sometimes in an infinite dimensional space.
At this point it should be clarified that the mapping increases the
extrinsic dimension of the data not the intrinsic. In kernel PCA we
want to find the principal components of the covariance matrix
\begin{equation}
\tilde{C} = \frac{1}{M}\sum_{j=1}^{M}\phi(x_j)\phi(x_j)^T
\end{equation}
It turns out that the mapping $\Phi$ doesn't have to be known
explicitly, all we need is the kernel that represents the dot
product.
\begin{equation}
k(x, y) = \phi(x)\phi(y)^T
\end{equation}
For some kernels there is an analytic factorization while for others
there is not. Not every function can be a valid kernel (dot
product). There are conditions that must be fulfilled. Basically for
every dataset they must give a valid kernel matrix as described in
previous sections.
Kernel PCA has many cousins such as Laplacian eigenmaps, Diffusion
maps Isomap, Local Linear embedding
\subsubsection{Customized Kernels}
As we described in the previous section there are infinite functions
that can act as valid kernels. The most popular kernels are:
\begin{enumerate}
\item The gaussian kernel $k(x, y)=e^{-\frac{||x-y||^2}{\sigma^2)}}$
\item The polynomial kernel $k(x, y)=(xy)^n$
\item The epanechnikov kernel $k(x, y)=1-||x-y||^2/\sigma^2, ||x-y||<\sigma, else 0$
\item The k-nearest neighbor kernel
\end{enumerate}
Although these are valid kernels for any kernel method including
Kernel PCA they are not appropriate for every task. In some cases we
want to map data in higher dimensional spaces preserving distances
(isometry) while in other cases we want to make small distances
smaller and send high distances to infinity (clustering). The
gaussian kernel for example tends to behave better in clustering
since it fades out very quickly. On the other hand it fails to
unfold the data manifold so that the intrinsic dimension can be
estimated. Manifold unfolding is a very important procedure since it
can reveal the true dimension of the dataset.
In [?] the authors introduced the idea of building the kernel matrix
from scratch without using any fixed kernel. As we have already
mentioned the kernel matrix has to be a positive semidefinite
matrix. Initially we have to define a k neighborhood for every
point. So at first we have to connect points to their neighbors and
build a initialize the non zero values of the matrix to 1. Then we
maximize the trace of the kernel matrix under the constraint that
the the distances in the kernel matrix are preserved. This leads to
a semidefinite programming problem.
\subsection{Forming manifolds from data for clustering through
non-parametric density estimation}
\subsection{Speech recognition with nearest neighbor search}
@@ -1013,18 +981,91 @@ MFCC features}{8cm}
} \label{kneighborLOOCVNTIMIT}
\end{table}
\section{Proposed Research}
\label{proposed}
The goal of this thesis is to make kernel methods scalable.
Efficient computation of the kernel matrix for large amounts of
data. Special emphasis is given in high dimensional data $d>30$. One
of the goals is to use Kernel PCA for unfolding the speech manifold.
In other words estimate the intrinsic dimensionality of speech by
using big databases like timit, wall street journal and broadcast
news. Implicitly dimensionality reduction will lead to new feature
generation we will explain it in the following sections. The
ultimate goal is to do speech recognition by nearest neighbor
search.
\subsection{Customizing Kernels with fast
gradient descent}
As we described in section ? it is possible to customize the kernel
for a given dataset, through semidefinite programming. Although this
method is very promising and gives very good results it is not
scalable. Semidefinite programming is solved with the intrinsic
point method which is of polynomial complexity.
\subsubsection{Fast kernel summation with trees}
Kernel summation is on of the most common tasks in kernel methods.
The general formula
\begin{equation}
\label{kernel_sum}
G(x)=\sum_{q=0}^{N} w_i\exp(\frac{(x-x_q)^2}{\sigma^2})
\end{equation}
appears in many cases, with the most interesting case the
eigenvalue/eigenvector computation. The kernel summation is a very
computationally intensive. For Gaussian kernels there are
approximations that can speed up computations orders of magnitude.
We will briefly discuss them here. Besides computation saving we can
save memory as well. The gaussian kernel has infinite support so the
kernel matrix $K, K(x,y)=\exp(\frac{(x-y)^2}{\sigma^2}$ is not
sparse. Numerical iterative methods for eigenvector computation are
based on the general iteration $x_n=x_{n-1}+\lambda Kx_{n-1}$, where
the multiplication of each row of $K$ with $x_{n-1}$ is a kernel
summation problem. So we don't really need to store the kernel
matrix, we just have to compute the kernel sum for every row using
$x_{n-1}$ as a the wight vector.
The gaussian can be approximated with the help of Hermite
polynomials. The Hermite polynomials $H_n(t)$ are defined by the
Rodrigues formula:
\begin{equation}
H_n(t)=(-1)^n\exp(t^2)D^n\exp(-t^2), t\in \Re
\end{equation}
where $D=d/dt$.
The gaussian can now be expressed as
\begin{equation}
\exp^(-(t-s)^2)=\sum_{n=0}^{\inf}\frac{s^n}{n!}h_n(t)
\end{equation}
where $h_n(t)$ is defined by
\begin{equation}
h_n(t)=\exp(-t^2)H_n(t).
\end{equation}
\subsection{Computing Gaussian Kernel Matrix with trees}
Direct computation of the kernel(Gaussian) kernel leads to a
non-sparse matrix. There is a way though to compute the
inverse/eigenvalues without explicitly computing the kernel matrix.
Several optimizations can be done for high dimensions
\subsection{Intrinsic dimension of speech}
\subsection{Speech recognition with nearest neighbor search}
Implement a merge reduce scheme to create different descriptions of
speech. Use nearest neighbor algorithms to replace gaussian mixtures
in hidden markov models.
\subsection{Acoustic Environment adaptation}
\section{Work remaining to be done}
\label{remains}
\section{Facilities and Equipment Needed}
\label{facitlies}
\end{Body}
\begin{EndMatter}
\appendix{This is my first appendix}
\appendix{This is my second appendix}
%% Below are some optional environments that can be used here
\bibliographystyle{ieeetr}
\references % Generates the bibliography page
+29 -34
View File
@@ -15,42 +15,37 @@
\contentsline {subsubsection}{\numberline {3.1.3}Ball trees}{5}{subsubsection.16}
\contentsline {subsubsection}{\numberline {3.1.4}Nearest Neighbor Algorithm}{6}{subsubsection.18}
\contentsline {subsubsection}{\numberline {3.1.5}All nearest Neighbors}{7}{subsubsection.19}
\contentsline {subsubsection}{\numberline {3.1.6}Fast kernel summation with trees}{8}{subsubsection.20}
\contentsline {subsection}{\numberline {3.2}Kernel methods}{8}{subsection.25}
\contentsline {subsubsection}{\numberline {3.2.1}Kernel principal component analysis}{9}{subsubsection.26}
\contentsline {subsubsection}{\numberline {3.2.2}Customized Kernels}{10}{subsubsection.30}
\contentsline {subsubsection}{\numberline {3.2.3}Stochastic optimization}{11}{subsubsection.35}
\contentsline {subsubsection}{\numberline {3.2.4}Support Vector Machines}{11}{subsubsection.36}
\contentsline {subsection}{\numberline {3.2}Large scale trees, for out of core memory}{8}{subsection.20}
\contentsline {subsection}{\numberline {3.3}Memory layout of Kd-trees}{9}{subsection.21}
\contentsline {subsection}{\numberline {3.4}System architecture}{10}{subsection.22}
\contentsline {subsection}{\numberline {3.5}Memory Manager Architecture}{10}{subsection.23}
\contentsline {subsection}{\numberline {3.6}Memory Mapped Files}{11}{subsection.24}
\contentsline {subsection}{\numberline {3.7}User defined cache, with TPIE}{11}{subsection.25}
\contentsline {subsection}{\numberline {3.8}Kernel methods}{12}{subsection.26}
\contentsline {subsubsection}{\numberline {3.8.1}Kernel principal component analysis}{13}{subsubsection.27}
\contentsline {subsubsection}{\numberline {3.8.2}Customized Kernels}{14}{subsubsection.31}
\contentsline {subsection}{\numberline {3.9}Forming manifolds from data for clustering through non-parametric density estimation}{15}{subsection.36}
\contentsline {subsection}{\numberline {3.10}Speech recognition with nearest neighbor search}{15}{subsection.37}
\contentsline {subsection}{\numberline {3.11}Introduction}{15}{subsection.38}
\contentsline {subsection}{\numberline {3.12}Speech features}{16}{subsection.39}
\contentsline {subsection}{\numberline {3.13}Fast N-body methods for feature comparison}{17}{subsection.40}
\contentsline {subsection}{\numberline {3.14}Kernel PCA on a manifold}{18}{subsection.44}
\contentsline {subsection}{\numberline {3.15}Experimental Results}{19}{subsection.45}
\contentsline {subsection}{\numberline {3.16}All nearest neighbor performance}{20}{subsection.46}
\contentsline {subsection}{\numberline {3.17}Comparison of NRAF and MFCC}{20}{subsection.47}
\contentsline {subsection}{\numberline {3.18}Discussion}{21}{subsection.49}
\setcounter {tocdepth}{2}
\contentsline {section}{\numberline {\MakeUppercase {Section} 4}\MakeUppercase {Proposed Research}}{11}{section.37}
\contentsline {subsection}{\numberline {4.1}Large scale trees, for out of core memory}{11}{subsection.38}
\contentsline {subsection}{\numberline {4.2}Memory layout of Kd-trees}{13}{subsection.39}
\contentsline {subsection}{\numberline {4.3}System architecture}{13}{subsection.40}
\contentsline {subsection}{\numberline {4.4}Memory Manager Architecture}{14}{subsection.41}
\contentsline {subsection}{\numberline {4.5}Memory Mapped Files}{14}{subsection.42}
\contentsline {subsection}{\numberline {4.6}User defined cache, with TPIE}{15}{subsection.43}
\contentsline {subsection}{\numberline {4.7}Customizing Kernels with fast gradient descent}{16}{subsection.44}
\contentsline {subsection}{\numberline {4.8}Scalable Support Vector Machines}{16}{subsection.45}
\contentsline {subsection}{\numberline {4.9}Computing Gaussian Kernel Matrix with trees}{16}{subsection.46}
\contentsline {subsection}{\numberline {4.10}Intrinsic dimension of speech}{16}{subsection.47}
\contentsline {subsection}{\numberline {4.11}Speech recognition with nearest neighbor search}{16}{subsection.48}
\contentsline {subsection}{\numberline {4.12}Introduction}{17}{subsection.49}
\contentsline {subsection}{\numberline {4.13}Speech features}{18}{subsection.50}
\contentsline {subsection}{\numberline {4.14}Fast N-body methods for feature comparison}{18}{subsection.51}
\contentsline {subsection}{\numberline {4.15}Kernel PCA on a manifold}{20}{subsection.55}
\contentsline {subsection}{\numberline {4.16}Experimental Results}{21}{subsection.56}
\contentsline {subsection}{\numberline {4.17}All nearest neighbor performance}{21}{subsection.57}
\contentsline {subsection}{\numberline {4.18}Comparison of NRAF and MFCC}{22}{subsection.58}
\contentsline {subsection}{\numberline {4.19}Discussion}{23}{subsection.60}
\contentsline {subsection}{\numberline {4.20}Acoustic Environment adaptation}{25}{subsection.68}
\contentsline {section}{\numberline {\MakeUppercase {Section} 4}\MakeUppercase {Proposed Research}}{22}{section.57}
\contentsline {subsection}{\numberline {4.1}Customizing Kernels with fast gradient descent}{23}{subsection.58}
\contentsline {subsubsection}{\numberline {4.1.1}Fast kernel summation with trees}{24}{subsubsection.59}
\contentsline {subsection}{\numberline {4.2}Computing Gaussian Kernel Matrix with trees}{25}{subsection.64}
\contentsline {subsection}{\numberline {4.3}Intrinsic dimension of speech}{25}{subsection.65}
\contentsline {subsection}{\numberline {4.4}Speech recognition with nearest neighbor search}{25}{subsection.66}
\contentsline {subsection}{\numberline {4.5}Acoustic Environment adaptation}{25}{subsection.67}
\setcounter {tocdepth}{2}
\contentsline {section}{\numberline {\MakeUppercase {Section} 5}\MakeUppercase {Work remaining to be done}}{25}{section.69}
\contentsline {section}{\numberline {\MakeUppercase {Section} 5}\MakeUppercase {Work remaining to be done}}{25}{section.68}
\setcounter {tocdepth}{2}
\contentsline {section}{\numberline {\MakeUppercase {Section} 6}\MakeUppercase {Facilities and Equipment Needed}}{25}{section.70}
\setcounter {tocdepth}{2}
\contentsline {section}{\numberline {\MakeUppercase {Appendix} A}\MakeUppercase {This is my first appendix}}{27}{section.71}
\setcounter {tocdepth}{2}
\contentsline {section}{\numberline {\MakeUppercase {Appendix} B}\MakeUppercase {This is my second appendix}}{28}{section.72}
\contentsline {section}{\numberline {\MakeUppercase {Appendix} A}\MakeUppercase {This is my first appendix}}{26}{section.69}
\setcounter {tocdepth}{3}
\contentsline {section}{\MakeUppercase {References}}{29}{section*.73}
\contentsline {section}{\MakeUppercase {References}}{27}{section*.70}
\setcounter {tocdepth}{2}