recognizedigits.m

this is digit recognition using som

% Y is label.
% I is input image, which can be assumed to be 28-by-28.
function Y = recognizeDigits(I)
  if size(I,1) ~= 28 | size(I,2) ~= 28
    error('I is of incorrect size');
  end
W1 =[  2.8040
    2.6566
    2.7139
    8.5952
    5.6010
    4.7382
    6.1281
   -1.3221
    3.7724
    2.4032
    7.7519
    4.7580
    7.3120
    9.9566
    5.3142
    0.7607
    7.8809
    0.3704
    5.8282
    9.2741
    2.1743
    7.9325
    0.5968
    3.3670
    5.7413
    0.3866
    4.8349
    5.4046
    0.5891
    6.5681
    7.8684
   -1.9800
    9.0608
    7.9031
    7.4736
    0.5161
    9.2113
    4.3898
    0.6197
    1.2997
    2.2024
    0.3511
    7.2531
    6.6083
    8.6932
    2.7309
    1.7473
    5.3583
    2.5410
    2.3759
    7.7399
    1.4362
    2.0397
    8.1061
    9.0581
    6.1178
    0.3456
    2.9989
    0.5835
   -1.8519
    4.8170
    4.9191
    2.8994
    3.6224
    9.5508
    7.7459
    8.2537
    2.8278
    0.2948
    6.2317
    3.5219
    4.4311
    6.6821
    1.4218
    8.1630
   -2.4035
    5.3564
    7.6380
    3.7858
   -2.5787
    3.2510
    7.3132
    7.4341
    0.8294
    3.7846
    1.2955
    3.1717
    5.2990
    4.0388
    4.4778
    0.2406
    0.0046
    9.1209
    1.0865
    8.6997
    1.5515
    7.4561
    7.2111
    1.3292
    5.8722
    2.0152
    8.1214
    8.8167
    3.1274
    5.5393
    8.8772
    5.2448
    7.3746
    0.4788
    6.1460
    7.6640
    3.6871
    0.6148
    9.4888
    2.1584
   -1.7857
    3.6759
    8.3890
    3.0795
   -2.8696
    7.2370
    0.0326
    5.1113
    4.5074
    4.8077
    4.0878
    3.1849
    6.9187
    3.2753
    7.2380
    5.4587
    1.3213
    1.8130
    2.6363
    7.8070
    4.3484
    7.3669
    6.0515
    0.9956
    5.5910
    7.3358
    4.0452
    9.8558
    8.8634
    7.4874
    0.5339
    7.8754
    5.1017
    5.4655
    6.5213
    9.1547
    2.7080
    3.7510
    0.8377
    7.1037
    5.6947
    9.3262
    0.7672
    6.4229
    4.8681
    4.4695
    9.3312
    1.6478
    7.5518
    9.5297
    7.4850
    0.8724
    8.0529
    3.0268
    4.9181
    5.9312
    6.7746
    2.8738
    2.1121
    9.1585
    9.4835
    1.0872
    4.7754
    6.7885
    0.5096
    4.5195
    4.5780
    1.2654
    8.5321
    2.1869
    7.5688
    2.7092
    8.7914
    1.7413
    1.3301
    7.6766
    3.3675
    3.9665
    2.8311
    8.1355
    4.0019
    5.5294
    5.5602
    5.3735
    3.6149
    8.6910
    6.2489
    7.9836
    1.1813
    7.6809
    7.6885
    2.5807
   -2.5568
    4.1407
    5.1327
    7.9450
    2.5517
    4.6633
    0.6981
    3.2787
   -2.9131
    7.6549
    9.3102
    0.5640
    8.9736
    1.5147
    2.2004
    5.7275
    2.3502
    3.6827
    4.1417
    9.2076
   -2.7622
    2.6180
    6.7448
    6.1195
    7.4744
    6.1842
    8.3918
    2.8309
    7.2824
    5.7600
    2.6458
    5.0156
    9.7651
    1.9559
    4.5367
    9.8126
   -2.6299
    7.7741
    6.5701
    7.5122
    5.9429
    4.6950
    9.1604
    2.2781
    4.9607
    6.6533
    2.9218
    8.3962
    7.4018
    2.5729
    9.2214
    8.1462
    0.1895
    5.7669
    9.5434
    2.2051
    7.4364
    3.1649
    9.2263
    0.8407
    7.0257
    6.3160
    2.0035
    8.2353
    1.4342
    9.1140
    8.0970
    4.9014
   -2.7873
    2.6888
    6.8203
    9.0810
    9.4816
    2.5290
    4.5290
    6.1578
    7.9649
    7.5992
    2.0759
    2.8043
    3.8006
    7.1074
    4.0708
    4.4269
    2.7734
    7.9001
    1.3785
    2.2109
   -2.8907
    9.9530
    8.4636
    0.7216
    9.7057
    0.6135
    7.5155
    3.0163
    6.9369
    7.8188
    3.4249
    3.3352
    5.7531
    2.4149
    6.5201
    1.1953
    7.0278
    8.8797
    6.7918
    6.9946
    5.6583
    9.5115
    0.7284
    7.3941
    7.6052
    8.6266
    7.3421
    0.6375
    3.9351
    2.8415
    5.4150
    0.7945
    6.3904
    6.4866
    7.0537
    6.6116
    7.4288
    8.5443
    2.5728
    4.8740
    2.1790
    2.7121
    3.9065
    0.9025
   -2.4068
    8.2011
    6.0879
    5.4753
    4.4647
    8.5715
    3.9084
    9.0727
   -0.3732
    5.2606
    3.8824
    8.8064
    7.4626
    3.1164
    8.4757
    0.0699
    7.6104
    7.0075
    8.6252
    5.6116
    7.6767
    8.1887
    5.5216
    7.5991
    4.9743
    7.8502
    8.2380
    5.9525
    6.6448
    5.4524
    6.6868
    1.9627
   -2.6885
    0.9408
    3.5278
    3.9121
    7.8704
    6.6184
    1.4677
    0.7773
    5.2342
    6.2582
    9.7777
    3.3669
    3.7351
    9.2690
    8.1782
    7.4292
    6.2356
    6.4147
    7.0688
    9.5069
    7.0785
    0.5188
    5.5426
    9.2001
    6.1129
    6.8934
    9.9963
    8.4999
    8.0470
    1.7170
    1.7909
    1.6562
   -2.4081
    7.2562
    9.5498
    0.0409
    6.0772
    3.3471
    1.9286
    3.6066
    4.2867
    3.5496
    5.0254
    1.4661
    5.2607
    4.0361
    8.5121
    3.8364
    8.5506
    6.7239
    6.9328
    5.2540
    2.8613
    6.1398
    7.0012
    3.4389
    3.8402
    1.5981
    9.1369
    2.9893
    5.3582
    1.9234
    7.3581
    8.4960
   -3.0151
    8.1249
    3.7375
    9.9472
    8.5822
    8.4120
    6.0622
    5.3003
    3.4192
    1.4921
    3.3161
    0.7360
    4.2889
    8.1517
    5.5308
    4.8973
    6.8405
    0.2849
    2.2997
    6.8023
    1.0799
    5.0190
    7.2313
    3.1990
   -1.1584
    3.7065
    8.9272
    3.8686
    0.2378
    8.5206
    3.2246
    5.5471
    4.7602
    0.6018
    2.4167
    5.1626
    5.5569
    0.0987
    6.3788
    2.0614
    4.5659
    5.3107
    6.1359
    6.2078
    3.2503
    9.3570
    5.7908
    8.5863
    4.4440
    0.6905
    4.5291
    0.6314
    5.0796
    2.2373
    4.4434
    8.6602
    5.2264
    3.7797
    8.1673
    2.0926
    3.2757
    5.1776
    7.9030
    2.3204
    6.2203
    4.5596
    1.6316
    3.2115
    9.2321
    7.5080
    1.5928
    7.5491
    5.9220
    1.9073
   10.0161
    8.0561
    4.9655
    0.7453
    0.7891
    3.6366
    5.4438
    1.9384
    7.2296
    9.6883
   -0.2304
    4.5396
    8.9071
    3.1777
   -1.9435
    9.0150
    9.8542
    3.1173
    3.5301
    1.0401
    5.0230
    0.7481
    4.1588
    0.8848
    0.3373
    9.1148
    6.7415
    9.7029
    1.9692
    4.9383
   -1.5063
    9.3991
    0.5778
    4.2698
    4.4221
    3.5652
    0.8261
    5.1061
    7.5152
    7.1089
    9.6494
    2.9480
    6.5857
    0.0480
    9.5329
    2.3775
    0.4386
    2.3354
    5.2981
    2.9034
    4.4667
    6.9597
    9.5275
    6.6186
    9.0017
    9.3006
    1.9844
    0.9865
    4.4427
    0.6802
    6.4286
    1.4199
    9.6409
    8.4047
    8.7210
    0.9437
    1.5583
    1.4123
    5.7456
    6.4343
    3.8937
    3.3874
    2.7119
    5.9949
    5.7125
    9.8437
    3.0334
    5.2864
    3.2420
    9.0600
    1.7374
    9.7946
    2.5319
    2.2006
    1.6930
    7.0677
    6.6321
    1.5808
    2.1063
    7.6931
    7.9623
    6.9409
    1.3190
    3.8095
    1.7311
    3.9355
    6.9160
    4.7437
    3.1904
    5.4081
    2.2636
    8.7574
    8.6969
    4.1640
    4.2419
    1.1649
    6.7001
    1.3483
    1.4678
    0.5063
   -2.2009
    9.5121
    8.8173
    1.2814
    6.5906
    1.7528
    5.1620
    9.6742
    9.7526
    1.6391
    9.8657
    6.9958
    5.6914
    5.4579
    1.0082
    7.6631
    8.5494
    6.1822
    8.9411
    8.0330
    2.4733
    5.9898
    5.6725
    4.2672
    4.9758
    8.0923
    0.5392
    7.0270
    3.7645
    3.2372
    8.0061
    0.5777
    1.9688
    0.4001
    1.6258
    3.4551
    2.7296
    5.7799
    2.1760
    1.1023
    2.1295
    1.7016
    8.9850
    2.3014
    1.4612
    4.9366
    3.5598
    1.0509
    5.3080
    3.2898
    5.8461
    9.6487
    9.3614
    1.3887
    5.5414
    6.6050
    7.0912
    6.2541
    3.9225
    8.8026
    5.6713
    1.8300
    2.8234
    2.7696
    8.8176
    9.1495
    7.1371
    6.0253
    6.5439
    1.3527
    1.4379
    8.6613
    2.2776
    5.4568
    2.7282
    5.4373
    3.5973
    5.3587
    2.0725
    8.4522
    2.7841
    3.7953
    2.4727
    6.2877
    4.0947
    4.3181
    2.9354
    2.2235
    5.0313
    3.7904
    8.6597
    4.9790
   10.0852
    9.9243
    9.7664
    4.5073
    4.0209
    1.3306
    6.3018
    9.4757
    9.7412
    1.6937
    2.1855
    7.6311
    7.5820
    4.4569
    9.9285
    5.1867
    8.1228
    5.7914
    7.7922
    3.5019
    2.8933
    8.1102
    0.2534
    9.7954
   -0.4233
    6.8846
    3.9226
    3.4472
    6.2646
    3.5964
    9.2709
    2.0618
    1.4918