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.1250 [' t1 P" k+ |- [7 H) {7 N) ?
    89.98        7        5        0        7.18756 Z9 G/ z( `" r3 a  T) y
    82.73        10        5        2        8.75, ]& B; z) x7 C6 Y5 I, o. n5 e; C
    76.82        5        8        0        6.25
    2 p/ j/ s* R' o+ k76.91        7.6        10        0        9.375
    4 H3 u! m2 `1 r85.22        10        0        6        7.5
    ) n  O( `$ \% L* z; n' l. S80.68        5.5        8        0        8.125
    : j" o+ Q) N; O; q. u0 ~8 j84.59        3        0        0        4.6875
    4 D0 b) d. F6 n. P( }86.05        9.5        5        3        4.375/ x6 A9 i* @4 r& f8 E: C
    83.89        10        0        0        5.3125
    4 \6 P9 b2 s/ e! [81.43        1        8        0        7.5
    & e: E7 M# O) j  ^8 ~87.15        0.5        0        0        5.625  ^9 H4 r7 b* ^  v1 J
    73.22        1        0        0        3.75
    2 z7 A& y' P4 `92.82        10        5        4        9.0625];
    3 T0 S1 f& G2 `, v[ma,na]=size(A);; j  t8 _# e* z2 M% g
    B=zeros(1,na);# m, G9 Z+ g8 ?1 |: ]) M6 {
    P=zeros(1,na);3 `' z6 ~% k3 G; q, b9 n/ Y
    G=zeros(1,na);
    # V: X3 o/ d( s9 v* aS=zeros(1,ma);
    " ~! v4 n; `' O3 p0 d! G; tg=0;
    8 T; c5 X$ ~& p9 u3 K' Ra=0.5;
    0 v  s$ H3 A8 Q3 @" i$ GMAX=max(A);
    : ?; k2 O3 \$ F: d1 H/ ^MIN=min(A); 7 U/ p) m* P- E6 G5 r7 H' a
    for i=1:ma       _3 C1 L5 \2 W2 L
        for j=1:na      
    : b- a, V/ F# K3 A$ q1 M) O        A(i,j)=((A(i,j)-MIN(j))/(MAX(j)-MIN(j)))*a+(1-a);   %标准化矩阵A& E" N5 D$ f( l8 L' o5 \8 d
        end
    5 L4 ?6 p  z! Lend
    - Q* P! ]- M9 B* q& efor j=1:na     5 m) q9 }3 H: I
        for i=1:ma     
    8 t3 n" \" x( D; A! c: ^! a        B(j)=B(j)+A(i,j);
    & C7 b9 k, y8 M! v    end5 K( L/ S# p; @9 ^8 P; g" f& @. I
    end % q# |; \) w/ Q) h1 L
    for i=1:ma  
    4 B0 G4 s4 l1 w7 `    for j=1:na    + w& n( i1 ~( ]. w8 x6 o" f8 I. b. c
            A(i,j)=A(i,j)/B(j)
    9 ?4 z% `9 z0 M# ?6 x2 {    end
    0 N0 u4 J) l* C1 E1 N+ Uend ' Y) g' x: z# t7 a* d0 J2 b
    for j=1:na / T7 w5 a5 O. Z4 v4 c: |; I
        for i=1:ma
    3 j# d. q+ k  I7 Q: N7 j        P(j)=P(j)+A(i,j)*log(A(i,j));
    . }) p  R; |" ?5 N    end 9 v4 ?& Q6 n6 Y0 j9 n
    end8 G" |( Y) D; T7 S; i2 U
    k=1/log(na); 8 y2 E- e2 n! n
    P=P*(-k)
    : J& S+ Y' a$ f+ {) j2 pfor j=1:na   
    1 u% x% s( \" F( d    G(j)=1-P(j)# K2 }  \2 ?8 N) n/ L
        g=g+G(j); 2 _9 q# ]4 }0 }0 n3 y/ H5 J
    end ) {0 ?; n! a7 U- Y" u2 [" n# [, U
    for j=1:na  8 s" z" ~2 A  s( G2 V/ `2 o+ }
        b(j)=G(j)/g4 B3 \: z# Y! [8 q
    end , b: f1 [6 ^2 Y( A/ K
    for i=1:ma8 y: x; E+ J% \; q) b6 K
        for j=1:na2 F7 W# N1 g- x5 i$ p4 A& v7 d
            S(i)=b(j)*A(i,j)
    # Y. U9 H5 I3 j- u" c4 v- T    end; C: Q2 s; i! w4 L( C4 H7 W
    end+ {5 \) D8 ?5 J. U- n1 k$ ~  n  P3 q
            
    4 ?. |* T3 _, s* S5 X; I9 V: L  W5 E6 {7 d3 J# T; F
    9 N/ V& A% D$ J
    / Q. {3 L( |- I. W+ J

    % E' G  u0 C/ h3 @2 `4 l% h
    回复

    使用道具 举报

    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-10-10 03:40 , Processed in 0.439648 second(s), 67 queries .

    回顶部