归一化matlab代码

>> A=[1 1/2 4 3 3;2 1 7 5 5;1/4 1/7 1 1/2 1/3;1/3 1/5 2 1 1;1/3 1/5 3 1 1];

>> B1=[1 2 5;1/2 1 2;1/5 1/2 1];

>> B2=[1 1/3 1/8;3 1 1/3;8 3 1];

>> B3=[1 1 3;1 1 3;1/3 1/3 1];

>> B4=[1 3 4;1/3 1 1;1/4 1 1];

>> B5=[1 1 1/4;1 1 1/4;4 4 1];

>> [x,y]=eig(A)

x =

  -0.4658             0.4419 + 0.2711i   0.4419 - 0.2711i  -0.3672 + 0.2415i  -0.3672 - 0.2415i

  -0.8409             0.7773             0.7773             0.8575             0.8575          

  -0.0951            -0.0200 - 0.1557i  -0.0200 + 0.1557i  -0.0190 + 0.00i  -0.0190 - 0.00i

  -0.1733            -0.0261 + 0.0783i  -0.0261 - 0.0783i  -0.0737 - 0.2146i  -0.0737 + 0.2146i

  -0.1920            -0.2829 + 0.1247i  -0.2829 - 0.1247i   0.0748 + 0.1185i   0.0748 - 0.1185i

y =

   5.0721                  0                  0                  0                  0          

        0            -0.0307 + 0.6008i        0                  0                  0          

        0                  0            -0.0307 - 0.6008i        0                  0          

        0                  0                  0            -0.0053 + 0.0548i        0          

        0                  0                  0                  0            -0.0053 - 0.0548i

>> lamda=y(1,1)

lamda =

    5.0721

>> ci1=(lamda-5)/4;cr1=ci1/1.12

cr1 =

    0.0161

>> w1=x(:,1)/sum(x(:,1))

w1 =

    0.2636

    0.4758

    0.0538

    0.0981

    0.1087

>> [x,y]=eig(B1)

x =

  -0.02            -0.02            -0.02          

  -0.4132             0.2066 + 0.3578i   0.2066 - 0.3578i

  -0.1918             0.0959 - 0.1661i   0.0959 + 0.1661i

y =

   3.0055                  0                  0          

        0            -0.0028 + 0.1290i        0          

        0                  0            -0.0028 - 0.1290i

>> lamda=y(1,1)

lamda =

    3.0055

>> ci21=(lamda-3)/2;cr21=ci21/0.58

cr21 =

    0.0048

>> w21=x(:,1)/sum(x(:,1))

w21 =

    0.5954

    0.27

    0.1283

>> [x,y]=eig(B2)

x =

  -0.1128            -0.05 - 0.0977i  -0.05 + 0.0977i

  -0.3255            -0.1627 + 0.2819i  -0.1627 - 0.2819i

  -0.9388             0.9388             0.9388          

y =

   3.0015                  0                  0          

        0            -0.0008 + 0.0680i        0          

        0                  0            -0.0008 - 0.0680i

>> lamda=y(1,1)

lamda =

    3.0015

>> ci22=(lamda-3)/2;cr22=ci22/0.58

cr22 =

    0.0013

>> w22=x(:,1)/sum(x(:,1))

w22 =

    0.0819

    0.2363

    0.6817

>> [x,y]=eig(B3)

x =

   -0.8847    0.6882   -0.2058

    0.4423    0.6882   -0.9057

    0.1474    0.2294    0.3705

y =

   -0.0000         0         0

         0    3.0000         0

         0         0    0.0000

>> lamda=y(2,2)

lamda =

    3.0000

>> ci23=(lamda-3)/2;cr23=ci23/0.58

cr23 =

 -7.6567e-016

>> w23=x(:,2)/sum(x(:,2))

w23 =

    0.4286

    0.4286

    0.1429

>> [x,y]=eig(B4)

x =

  -0.9255             0.9255             0.9255          

  -0.2803            -0.1401 - 0.2427i  -0.1401 + 0.2427i

  -0.2547            -0.1273 + 0.2205i  -0.1273 - 0.2205i

y =

   3.0092                  0                  0          

        0            -0.0046 + 0.1663i        0          

        0                  0            -0.0046 - 0.1663i

>> lamda=y(1,1)

lamda =

    3.0092

>> ci24=(lamda-3)/2;cr24=ci24/0.58

cr24 =

    0.0079

>> w24=x(:,1)/sum(x(:,1))

w24 =

    0.6337

    0.1919

    0.1744

>> [x,y]=eig(B5)

x =

   -0.43    0.2357         0

    0.2182    0.2357   -0.2425

    0.8729    0.9428    0.9701

y =

   -0.0000         0         0

         0    3.0000         0

         0         0         0

>> lamda=y(2,2)

lamda =

     3

>> ci25=(lamda-3)/2;cr25=ci25/0.58

cr25 =

     0

>> w25=x(:,2)/sum(x(:,2))

w25 =

    0.1667

    0.1667

    0.6667

>> ci=[ci21 ci22 ci23 ci24 ci25];

>> cr=ci*w1/sum(0.58*w1)

cr =

    0.0027

>> w=[w21 w22 w23 w24 w25]*w1

w =

    0.2993

    0.2453

    0.4554