21.06.2014 Views

Numerical Methods Contents - SAM

Numerical Methods Contents - SAM

Numerical Methods Contents - SAM

SHOW MORE
SHOW LESS

You also want an ePaper? Increase the reach of your titles

YUMPU automatically turns print PDFs into web optimized ePapers that Google loves.

✬<br />

✩<br />

✬<br />

✩<br />

Lemma 5.3.4 (Ncut and Rayleigh quotient). (→ [36, Sect. 2])<br />

With z ∈ {−1, 1} N according to (5.3.9) there holds<br />

∑<br />

d i<br />

Ncut(X) = yT Ay<br />

y T Dy , y := (1 + z) − β(1 − z) , β := z i >0<br />

∑ .<br />

d i<br />

✫<br />

generalized Rayleigh quotient ρ A,D (y)<br />

z i 0d ∑ d i = β<br />

1+β<br />

. Also observe<br />

i<br />

y T Dy = (1 + z) T D(1 + z) + β 2 (1 − z) T D(1 − z) =<br />

4( ∑ d i + β 2∑ d i ) = 4β ∑<br />

z i >0 z i N<br />

1 x ∗ i } . (5.3.15)<br />

i=1<br />

mean value of entries of x ∗<br />

Ôº¿ º¿<br />

Ôº¿ º¿<br />

3 xs = reshape (V ( : , 2 ) ,m, n ) ; Xidx = find ( xs >(sum(sum( xs ) ) / ( n∗m) ) ) ;<br />

Code 5.3.13: 1st stage of segmentation of greyscale image<br />

1 P = imread ( ’ image .pbm ’ ) ; [m, n ] = size (P) ; [ A,D] = imgsegmat (P) ;<br />

2 [ V, E ] = eig ( f u l l (A+D) , f u l l (D) ) ; % grossly inefficient !

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!