QQ登录

只需要一步,快速开始

 注册地址  找回密码
12
返回列表 发新帖
楼主: w13king
打印 上一主题 下一主题

求熵值法matlab程序!!!

[复制链接]
字体大小: 正常 放大
benfifi        

0

主题

4

听众

76

积分

升级  74.74%

  • TA的每日心情
    开心
    2013-11-10 09:36
  • 签到天数: 38 天

    [LV.5]常住居民I

    自我介绍
    一名普通大学生,偏爱数学!
    11#
    发表于 2013-8-12 16:51 |只看该作者
    |招呼Ta 关注Ta
    A=[91.45        9.8        0        0        8.125
    9 z4 g* `- Y% z% M( v/ k89.98        7        5        0        7.1875
    2 C. K  c! `& O* w! {2 m1 W0 M' V( Z$ n82.73        10        5        2        8.75
    9 X# M# L, x5 \& W2 k76.82        5        8        0        6.25
    & Q) P3 a0 K3 Q$ V7 h76.91        7.6        10        0        9.375
    ) D0 {' m4 h& G4 D3 j6 E85.22        10        0        6        7.5: ~& y/ m7 G9 M& s5 S
    80.68        5.5        8        0        8.1258 v& q8 c; }% g; T7 W7 s% _
    84.59        3        0        0        4.6875
    4 H  l5 A) R, {" Q8 f7 l86.05        9.5        5        3        4.375: t' B, @; C" I: J8 d
    83.89        10        0        0        5.3125
    % T2 a8 D/ y! Q& H) [81.43        1        8        0        7.5( |7 t1 S& d, f1 h; V  B+ B$ x
    87.15        0.5        0        0        5.625
    + I: L& I' k4 H5 k3 Z73.22        1        0        0        3.754 x/ A# H( M. f* D$ u- l
    92.82        10        5        4        9.0625];6 s/ b" V1 m0 V! W- u. S, C$ K
    [ma,na]=size(A);
    ; U* S$ _; x2 @* z, w6 N' gB=zeros(1,na);
    5 V( m1 e) _4 g0 AP=zeros(1,na);  D7 d! _* |2 j% Y
    G=zeros(1,na);/ T- W$ O4 q6 {) H" r0 K
    S=zeros(1,ma);! S  Z3 n2 N' w
    g=0;
    5 K3 p" d5 P# ja=0.5; ; P3 m/ t0 T9 M3 e" A
    MAX=max(A);& E# |) s$ |& I# M" p. x$ m6 f
    MIN=min(A);
    8 N* K4 ^9 F# s. A' N6 n4 Wfor i=1:ma     
    2 g- Z: I2 A8 t5 X1 `$ a; b4 J& n    for j=1:na      & N! w: h+ W$ C: W+ @
            A(i,j)=((A(i,j)-MIN(j))/(MAX(j)-MIN(j)))*a+(1-a);   %标准化矩阵A
    " a$ a& o% T' l7 N' [    end
    0 i/ V& F; t9 O" mend
    : K) u  C/ I  t3 p. s3 L" ]for j=1:na     
    ; |/ C4 j3 C1 \    for i=1:ma     6 q( V1 m' K* p3 k3 Q1 c9 {+ _
            B(j)=B(j)+A(i,j);( X7 c3 b8 ~: ~) ?$ J; P
        end& ?; A9 Z! F0 d4 f* ]  ~2 e
    end " B, Z' \7 z& y, k2 W: u
    for i=1:ma  / k7 s( p' e% M9 ^
        for j=1:na    % \, s0 |: {* K/ @5 y
            A(i,j)=A(i,j)/B(j)
    ( f7 \  F$ _* }    end
    9 n: O  `' R) e- q' Mend
    ( K# K3 c! b& L" {for j=1:na
      ]! Q: j2 V0 B% g6 t; A! `+ r    for i=1:ma
    ) d6 W3 \% J% }# f' J        P(j)=P(j)+A(i,j)*log(A(i,j));% e& z2 s+ I! I0 f
        end 5 W6 R, J8 ~3 D1 o+ X9 l
    end
    . F, D6 H6 o: o: t  lk=1/log(na);
      g- u/ I7 p" l9 L4 ]P=P*(-k)
    0 `9 h- }( a1 X: o/ ^+ |for j=1:na   $ F$ ]/ n$ r0 L/ Z+ U# C) g, v
        G(j)=1-P(j): g3 }& J6 k$ }0 E. W
        g=g+G(j); , f# W% b5 J: O1 {" u5 e( v
    end
    0 X& o  Q# |: K/ V9 q) yfor j=1:na  ( P/ a( h+ s6 P" {& g4 I) W0 M' O
        b(j)=G(j)/g8 q8 e7 S8 T  I$ p3 A0 g3 E$ S
    end 7 s! Q- A$ _+ L* Y
    for i=1:ma
    / [: |$ ~; V0 Y- |    for j=1:na, [, j) W$ ~( K' {
            S(i)=b(j)*A(i,j)) b. \* i, i% ~, W+ H& j
        end
    + Y2 r, A. c% m6 tend* a7 H, v( E- d
            
    3 ]1 i/ ?0 }; ^9 x* H0 q
      [2 |) x' q9 Y1 W* s5 n5 g! U" J) l) R
    & b. {( F( P& S' A4 ^# r  e

    ) _. Y) K- V1 p5 l
    回复

    使用道具 举报

    38

    主题

    13

    听众

    444

    积分

    升级  48%

  • TA的每日心情
    开心
    2015-3-30 23:54
  • 签到天数: 71 天

    [LV.6]常住居民II

    自我介绍
    仰望星空,脚踏实地

    社区QQ达人

    群组第六届国赛赛前冲刺培

    群组全国大学生数学建模竞

    回复

    使用道具 举报

    逻辑学        

    0

    主题

    1

    听众

    2

    积分

    升级  40%

    该用户从未签到

    自我介绍
    开始建模,↖(^ω^)↗

    群组2014数学建模国赛备战

    回复

    使用道具 举报

    您需要登录后才可以回帖 登录 | 注册地址

    qq
    收缩
    • 电话咨询

    • 04714969085
    fastpost

    关于我们| 联系我们| 诚征英才| 对外合作| 产品服务| QQ

    手机版|Archiver| |繁體中文 手机客户端  

    蒙公网安备 15010502000194号

    Powered by Discuz! X2.5   © 2001-2013 数学建模网-数学中国 ( 蒙ICP备14002410号-3 蒙BBS备-0002号 )     论坛法律顾问:王兆丰

    GMT+8, 2026-7-10 10:09 , Processed in 1.176976 second(s), 68 queries .

    回顶部