culture.m

基于自组织特征映射（SOM）的系统辨识算法实现

P=[
    82.95 88.20 75.90 71.65 40.61 81.31 80.47;
    77.15 72.40 75.84 60.24 58.04 72.65 66.65;
    83.52 47.72 45.01 66.47 42.90 84.32 80.55;
    72.70 69.65 62.82 53.70 58.99 69.69 58.80;
    54.88 61.15 52.81 66.52 53.49 76.49 70.09;
    48.37 51.40 47.96 58.47 62.22 58.90 53.84;
    65.20 50.75 44.01 50.96 52.36 71.16 46.54;
    47.41 54.96 63.60 54.27 41.63 35.83 54.76;
    45.14 42.04 39.61 67.85 29.13 53.27 57.74;
    41.52 46.47 59.52 47.81 50.97 39.27 47.35;
    47.01 51.35 49.78 33.20 66.42 36.99 44.57;
    33.94 31.14 36.88 73.07 44.17 47.69 61.74;
    50.75 60.52 52.66 43.28 44.89 29.01 38.41;
    51.36 46.53 55.46 49.48 37.88 37.94 37.73;
    42.80 44.06 52.18 45.07 46.81 39.43 44.99;
    36.83 48.59 43.13 34.75 61.39 45.23 44.31;
    43.94 46.58 52.39 43.44 55.50 34.00 36.74;
    38.83 42.43 52.48 33.43 52.76 42.86 49.03;
    44.30 44.14 44.02 39.61 60.28 40.91 36.61;
    37.24 39.36 48.25 48.34 59.60 34.97 39.46;
    41.90 37.69 38.48 48.46 49.99 48.42 40.71;
    35.24 33.56 33.40 49.08 50.51 46.95 45.15;
    36.89 35.76 33.47 71.86 24.58 42.06 45.36;
    35.84 44.04 56.53 39.74 49.94 31.96 33.90;
    49.63 43.27 31.17 37.19 50.77 38.03 34.48;
    41.23 37.49 35.16 40.64 39.37 51.31 34.78;
    41.06 39.44 32.72 36.53 54.96 39.00 36.12;
    38.61 35.99 36.53 45.18 37.01 28.11 44.76;
    38.17 28.08 43.87 38.57 40.62 37.89 38.79;
    53.88 47.11 40.26 28.68 19.29 35.09 36.46;
    39.02 34.47 35.52 36.22 40.75 33.10 37.79;
    34.95 27.81 32.47 38.37 44.50 40.65 34.38;
    36.39 42.31 40.24 31.02 39.14 27.78 33.53;
    34.61 29.52 31.10 20.67 45.17 45.51 43.56;
    29.93 37.28 39.78 27.42 46.84 32.93 31.32;
    30.28 34.74 27.17 22.51 23.18 29.91 27.58]';
net=newsom(minmax(P),[3 2]);
a=[100 150 200];
for i=1:3
    net.trainParam.epochs=a(i);
    net=train(net,P);
    y=sim(net,P);
    y=vec2ind(y)
end