{VERSION 3 0 "IBM INTEL NT" "3.0" }
{USTYLETAB {CSTYLE "Maple Input" -1 0 "Courier" 0 1 255 0 0 1 0 1 0 0 
1 0 0 0 0 }{CSTYLE "2D Math" -1 2 "Times" 0 1 0 0 0 0 0 0 2 0 0 0 0 0 
0 }{CSTYLE "2D Comment" 2 18 "" 0 1 0 0 0 0 0 0 0 0 0 0 0 0 0 }
{CSTYLE "2D Output" 2 20 "" 0 1 0 0 255 1 0 0 0 0 0 0 0 0 0 }{PSTYLE "
Normal" -1 0 1 {CSTYLE "" -1 -1 "" 0 1 0 0 0 0 0 0 0 0 0 0 0 0 0 }0 0 
0 -1 -1 -1 0 0 0 0 0 0 -1 0 }{PSTYLE "Maple Output" 0 11 1 {CSTYLE "" 
-1 -1 "" 0 1 0 0 0 0 0 0 0 0 0 0 0 0 0 }3 3 0 -1 -1 -1 0 0 0 0 0 0 -1 
0 }{PSTYLE "" 11 12 1 {CSTYLE "" -1 -1 "" 0 1 0 0 0 0 0 0 0 0 0 0 0 0 
0 }1 0 0 -1 -1 -1 0 0 0 0 0 0 -1 0 }{PSTYLE "Maple Plot" 0 13 1 
{CSTYLE "" -1 -1 "" 0 1 0 0 0 0 0 0 0 0 0 0 0 0 0 }3 0 0 -1 -1 -1 0 0 
0 0 0 0 -1 0 }{PSTYLE "" 0 256 1 {CSTYLE "" -1 -1 "" 0 1 0 0 0 0 0 0 
0 0 0 0 0 0 0 }3 0 0 -1 -1 -1 0 0 0 0 0 0 -1 0 }}
{SECT 0 {EXCHG {PARA 256 "" 0 "" {TEXT -1 33 "Demonstration of image p
rocessing" }}{PARA 0 "" 0 "" {TEXT -1 0 "" }}{PARA 0 "" 0 "" {TEXT -1 
0 "" }}{PARA 0 "" 0 "" {TEXT -1 0 "" }}{PARA 0 "" 0 "" {TEXT -1 21 "A \+
simple image matrix" }}{PARA 0 "" 0 "" {TEXT -1 0 "" }}{PARA 0 "> " 0 
"" {MPLTEXT 1 0 11 "Digits:=10;" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%'
DigitsG\"#5" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 17 "A:=matrix(20
,20);" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%\"AG-%&arrayG6%;\"\"\"\"#?F
(7\"" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 58 "for i to 20 do\n   \+
for j to 20 do\n      A[i,j]:=0\n   od\nod;" }}}{EXCHG {PARA 0 "> " 0 
"" {MPLTEXT 1 0 123 "for j from 5 to 15 do\n   for i to 5 do\n      A[
i,j]:=1\n   od;\n   for i from 15 to 20 do\n      A[i,j]:=1   \n   od
\nod;\n      " }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 121 "for i fro
m 6 to 14 do\n   for j from 5 to 8 do\n      A[i,j]:=1\n   od;\n   for
 j from 12 to 15 do\n      A[i,j]:=1\n   od;\nod;" }}}{EXCHG {PARA 0 "
> " 0 "" {MPLTEXT 1 0 9 "evalm(A);" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#
-%'MATRIXG6#7676\"\"!F(F(F(\"\"\"F)F)F)F)F)F)F)F)F)F)F(F(F(F(F(F'F'F'F
'76F(F(F(F(F)F)F)F)F(F(F(F)F)F)F)F(F(F(F(F(F*F*F*F*F*F*F*F*F'F'F'F'F'F
'" }}}{EXCHG {PARA 0 "" 0 "" {TEXT -1 0 "" }}{PARA 0 "" 0 "" {TEXT -1 
0 "" }}{PARA 0 "" 0 "" {TEXT -1 28 "Singular value decomposition" }}
{PARA 0 "" 0 "" {TEXT -1 0 "" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 
0 13 "with(linalg):" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 21 "s:=e
valf(Svd(A,U,V));" }}{PARA 12 "" 1 "" {XPPMATH 20 "6#>%\"sG-%'VECTORG6
#76$\"0Gsx7O3M\"!#8$\"04B(e2ONO!#9$\"0h$QNcbE9!#G$\"0&z)Gz$3YE!#H$\"05
2NHT-]%!#V$\"0SaW%)Q9E#!#c$\"0M(zn_8DW!#e$\"0!o5TgQIO!#s$\"0hf\"4!*G:;
!#&)$\"0%42*y)G#>#!#()$\"0(>k*)[dx7!$,\"\"\"!FJFJFJFJFJFJFJFJ" }}}
{EXCHG {PARA 0 "" 0 "" {TEXT -1 0 "" }}{PARA 0 "" 0 "" {TEXT -1 0 "" }
}{PARA 0 "" 0 "" {TEXT -1 48 "Only the first 2 singular values are non
zero, so" }}{PARA 0 "" 0 "" {TEXT -1 13 "         A = " }{XPPEDIT 18 
0 "u[1]" "6#&%\"uG6#\"\"\"" }{TEXT -1 1 "*" }{XPPEDIT 18 0 "s[1]" "6#&
%\"sG6#\"\"\"" }{TEXT -1 1 "*" }{XPPEDIT 18 0 "transpose(v[1])" "6#-%*
transposeG6#&%\"vG6#\"\"\"" }{TEXT -1 0 "" }{TEXT -1 3 " + " }
{XPPEDIT 18 0 "u[2]" "6#&%\"uG6#\"\"#" }{TEXT -1 1 "*" }{XPPEDIT 18 0 
"s[2]" "6#&%\"sG6#\"\"#" }{TEXT -1 1 "*" }{XPPEDIT 18 0 "transpose(v[2
])" "6#-%*transposeG6#&%\"vG6#\"\"#" }{TEXT -1 0 "" }}{PARA 0 "" 0 "" 
{TEXT -1 0 "" }}{PARA 0 "" 0 "" {TEXT -1 0 "" }}{PARA 0 "" 0 "" {TEXT 
-1 22 "Get the column vectors" }}{PARA 0 "" 0 "" {TEXT -1 0 "" }}}
{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 13 "u1:=col(U,1):" }}}{EXCHG 
{PARA 0 "> " 0 "" {MPLTEXT 1 0 13 "u2:=col(U,2):" }}}{EXCHG {PARA 0 ">
 " 0 "" {MPLTEXT 1 0 13 "v1:=col(V,1):" }}}{EXCHG {PARA 0 "> " 0 "" 
{MPLTEXT 1 0 13 "v2:=col(V,2):" }}}{EXCHG {PARA 0 "" 0 "" {TEXT -1 0 "
" }}{PARA 0 "" 0 "" {TEXT -1 22 "Rank one approximation" }}{PARA 0 "" 
0 "" {TEXT -1 0 "" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 57 "B:=eva
lm(s[1]*u1&*transpose(v1)):   # s1*u1*transpose(v1)" }}}{EXCHG {PARA 
0 "" 0 "" {TEXT -1 0 "" }}{PARA 0 "" 0 "" {TEXT -1 22 "Rank two approx
imation" }}{PARA 0 "" 0 "" {TEXT -1 0 "" }}}{EXCHG {PARA 0 "> " 0 "" 
{MPLTEXT 1 0 38 "C:=evalm( B+ s[2]*u2&*transpose(v2) ):" }}}{EXCHG 
{PARA 0 "" 0 "" {TEXT -1 0 "" }}{PARA 0 "" 0 "" {TEXT -1 19 "Observe t
he results" }}{PARA 0 "" 0 "" {TEXT -1 0 "" }}}{EXCHG {PARA 0 "> " 0 "
" {MPLTEXT 1 0 10 "Digits:=2;" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%'Di
gitsG\"\"#" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 9 "evalm(C);" }}
{PARA 12 "" 1 "" {XPPMATH 20 "6#-%'MATRIXG6#7676\"\"!$!#7!#:$\"#B!#;F(
$\"#5!\"\"F/F/F/F/F/F/F/F/F/F/F(F(F(F(F(F'F'F'F'76F($!#SF.$!#QF.F(F/F/
F/F/$!\"&F+F7F7F/F/F/F/F(F(F(F(F(F2F2F2F2F2F2F2F2F'F'F'F'F'F'" }}}
{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 9 "evalm(B);" }}{PARA 12 "" 1 "
" {XPPMATH 20 "6#-%'MATRIXG6#7676\"\"!$!#$*!#;$!#=!#<F($\"#6!\"\"F/F/F
/$\"#l!\"#F2F2F/F/F/F/F(F(F(F(F(F'F'F'F'76F($!#wF+$!#:F.F($\"#))F4F:F:
F:$\"#`F4F<F<F:F:F:F:F(F(F(F(F(F5F5F5F5F5F5F5F5F'F'F'F'F'F'" }}}
{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 33 "plots[matrixplot](A,style=pa
tch);" }}{PARA 13 "" 1 "" {GLPLOT3D 380 380 380 {PLOTDATA 3 "6%-%%GRID
G6%;\"\"\"\"#?F&7676\"\"!F+F+F+$F'F+F,F,F,F,F,F,F,F,F,F,F+F+F+F+F+F*F*
F*F*76F+F+F+F+F,F,F,F,F+F+F+F,F,F,F,F+F+F+F+F+F-F-F-F-F-F-F-F-F*F*F*F*
F*F*-%&STYLEG6#%&PATCHG-%+AXESLABELSG6%%$rowG%'columnG%\"AG" 1 2 0 1 
0 2 1 1 1 2 1.000000 1.000000 1.000000 0 }}}}{EXCHG {PARA 0 "> " 0 "" 
{MPLTEXT 1 0 33 "plots[matrixplot](B,style=patch);" }}{PARA 13 "" 1 "
" {GLPLOT3D 380 380 380 {PLOTDATA 3 "6%-%%GRIDG6%;\"\"\"\"#?F&7676\"\"
!$!0lQ&3_5!H*!#H$!0q$)>3*4-=!#IF+$\"0XAE\">Mz5!#9F2F2F2$\"0-;$o!o3Z'!#
:F5F5F2F2F2F2F+F+F+F+F+76F+$!0xQ&3_5!H*F.$!0r$)>3*4-=F1F+$\"0YAE\">Mz5
F4F=F=F=$\"05;$o!o3Z'F7F?F?F=F=F=F=F+F+F+F+F+76F+$!0iQ&3_5!H*F.$!0p$)>
3*4-=F1F+$\"0WAE\">Mz5F4FFFFFF$\"0*fJo!o3Z'F7FHFHFFFFFFFFF+F+F+F+F+FAF
A76F+$!0QA\"=W([e(F.$!0K3@f<8Z\"F1F+$\"0y\\8*)\\A\"))F7FOFOFO$\"0xl'fz
6$G&F7FQFQFOFOFOFOF+F+F+F+F+FJFJFJFJFJFJFJFJFAFAFAFAFAFA-%&STYLEG6#%&P
ATCHG-%+AXESLABELSG6%%$rowG%'columnG%\"AG" 1 2 0 1 0 2 1 1 1 2 
1.000000 0.000000 1.000000 0 }}}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 
0 33 "plots[matrixplot](C,style=patch);" }}{PARA 13 "" 1 "" {GLPLOT3D 
380 380 380 {PLOTDATA 3 "6%-%%GRIDG6%;\"\"\"\"#?F&7676\"\"!$!0)*=IQFm;
\"!#G$\"05K-%yIrA!#HF+$\"0,++++++\"!#9F2F2F2$\"0+++++++\"F4F5F5F2F2F2F
2F+F+F+F+F+76F+$!0**=IQFm;\"F.F/F+$\"0-++++++\"F4F:F:F:F2F2F2F:F:F:F:F
+F+F+F+F+76F+F,F/F+$\"0(**************!#:F=F=F=$\"0)**************F?F@
F@F=F=F=F=F+F+F+F+F+F<F<76F+$!0Vu+(\\uFSF1$!0%*y#4!fq\"QF1F+F5F5F5F5$!
\"&F?FGFGF5F5F5F5F+F+F+F+F+FBFBFBFBFBFBFBFBF<F<F<F<F<F<-%&STYLEG6#%&PA
TCHG-%+AXESLABELSG6%%$rowG%'columnG%\"AG" 1 2 0 1 0 2 1 1 1 2 
1.000000 0.000000 1.000000 0 }}}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 
0 0 "" }}}}{MARK "3 0 0" 62 }{VIEWOPTS 1 1 0 1 1 1803 }