X=[227 237 234 235 235 235 235 235 235 235 235 235 235 235 235 235 235 235 234 237 227 237 255 255 255 255 255 255 255 255 255 255 255 255 255 255 255 255 255 255 255 237 237 255 254 254 254 254 254 254 254 254 254 254 255 255 254 254 254 254 254 255 237 237 255 254 255 255 255 255 255 254 254 255 255 254 244 255 255 255 255 254 255 237 237 255 254 255 255 255 255 255 255 255 228 200 167 112 209 255 254 255 254 255 237 237 255 254 255 255 255 255 254 227 167 92 0 37 0 106 240 255 254 254 255 237 237 255 254 255 255 255 254 192 67 0 0 120 218 137 64 224 255 254 254 255 237 237 255 254 255 254 255 234 67 0 0 111 230 243 124 70 224 255 254 254 255 237 237 255 254 255 254 255 247 167 161 208 240 255 224 77 78 221 255 254 254 255 237 237 255 254 255 255 254 255 252 255 255 231 235 191 0 85 231 255 254 254 255 237 237 255 254 255 255 254 255 255 234 157 41 82 71 0 182 255 254 254 254 255 237 237 255 254 255 254 255 240 166 68 61 121 62 0 3 169 252 255 255 255 255 237 237 255 254 254 255 224 124 0 154 237 221 78 0 24 0 164 228 238 246 255 237 237 255 254 255 228 96 21 166 241 214 82 31 183 209 145 0 82 107 202 255 235 237 255 255 242 120 0 155 197 133 0 106 215 254 255 251 215 137 147 234 255 236 237 255 255 228 0 0 14 0 56 161 234 255 254 254 255 255 255 255 255 255 237 237 255 255 248 187 113 141 204 228 254 255 254 255 255 254 254 254 254 254 255 237 237 255 254 255 255 245 248 255 255 255 254 255 255 255 255 255 255 255 254 255 237 237 255 254 254 254 255 255 254 254 254 254 254 254 254 254 254 254 254 254 255 237 237 255 255 255 255 255 255 255 255 255 255 255 255 255 255 255 255 255 255 255 237 227 237 234 235 235 235 235 234 235 234 234 234 234 234 235 235 235 235 234 237 227]; imagesc(X), colormap(gray), axis('image') [U,S,V] = svds(double(X),5); XX = abs(U*S*V'); imagesc(XX); axis('image')