数学建模社区-数学中国
标题:
求熵值法matlab程序！！！
[打印本页]
作者:
w13king
时间:
2009-9-5 11:28
标题:
求熵值法matlab程序！！！
熵值法有固定的程序可以套用的么？
. r# x: F! T% c# k' t
# G5 ~& v- a$ n. V' C: g
8 s V4 m9 E# f5 t, m
各位好心的童鞋。。。
, n- _+ W4 h7 |$ r2 ^8 u* t1 f2 V
第一次接触这个。。。
i9 c, N: V" i7 q* I' k4 O5 |: ]
不会。。。
6 r! m: o' N! y& \6 x
希望大家帮帮忙！
0 W5 R3 f) `% _5 r. o* f" q
谢谢~~
作者:
aileenrain
时间:
2009-9-5 11:54
我也想知道！
作者:
东方明珠-WDZYQ
时间:
2009-9-5 12:21
希望能够有啊
作者:
w13king
时间:
2009-9-5 14:19
大家都不知道么？
, U3 T& d6 O3 H! P
不会的吧...
- p4 s7 ? ~. y% n) Z5 ~
在线等高人解答啊~~~
作者:
熙123
时间:
2009-9-5 15:01
我也不知道呀！！无力！！
作者:
ykl126
时间:
2009-9-5 15:36
大哥呀，下次说清楚一些好不
作者:
w13king
时间:
2009-9-5 15:49
就是那个 纯净水安全监控问题
( Z! P4 ]4 X4 O& i4 R
这是我们的大概方法。。。
8 V/ x( M- X3 x! a4 w/ v
1.对于问题1，通过分析首先确定各危害指标的权值，这里我们采取两种确定方法，即主成分分析和专家评定，考虑到因专家评定带有主观因素与事实可能不符，进一步用主成分分析法，依照客观数据给出结果，对两种方法进行综合考虑，取其平均值作为最终的权重。在此基础上建立了纯净水安全风险检测的科学评价模型，并对该城区所有批次的纯净水进行评判排序。
* q/ B; y+ R) [& B5 o
2.对于问题2，在问题一模型的基础上，按照每个公司各项危害指标的平均值最为该公司的危险指标，利用问题1的模型求出各公司的风险度，并进行排序。依据各公司的各项危害指标及风险度，分析公司共产品的主要可能的危害因素。
$ n/ ]6 o7 @+ v, k6 b# h( y
3.首先筛选出不合格的样品，根据问题一中建立的科学评价模型对其进行成因分析。利用熵值法求出各危害因素在生产流通环节上所占的权重。评价各个危害指标的分布规律。最后，通过四类危害指标的分析，讨论出A、B、D公司的管理情况。
4 q6 p) I$ B' x- Q& S |' e
9 [" M- F5 ^* P% V; d1 _
就是不知道那个熵值法怎么编程啊。。。
作者:
crystal9
时间:
2011-8-21 18:00
希望可以有啊。。。我也要。。。
作者:
小p才
时间:
2011-8-22 16:39
真是好东西啊。。。。。
作者:
谢凌风
时间:
2012-8-12 18:55
求高人指点。
作者:
benfifi
时间:
2013-8-12 16:51
A=[91.45 9.8 0 0 8.125
: Q" y4 K" L, s# a7 {
89.98 7 5 0 7.1875
8 g) A/ }8 c( O# `% g
82.73 10 5 2 8.75
3 P1 x9 k( x) f- x$ H
76.82 5 8 0 6.25
; T% q- m3 r4 z9 x
76.91 7.6 10 0 9.375
- @8 i' t( s# a0 E
85.22 10 0 6 7.5
( R7 `9 x5 `: K; R# u0 d6 C: P
80.68 5.5 8 0 8.125
) y b4 j% ~& ?4 _- x( k
84.59 3 0 0 4.6875
- T9 @$ G9 V8 \7 `5 \
86.05 9.5 5 3 4.375
$ V* Y- m) {* G$ ~3 M- x$ ~
83.89 10 0 0 5.3125
# q% b! |; ?- M- {& P5 c
81.43 1 8 0 7.5
, H8 P" k. f }% ?7 N# y& O R
87.15 0.5 0 0 5.625
1 r8 D* \# c" Y6 r8 G" D
73.22 1 0 0 3.75
4 W; `! _4 S+ x* {/ m$ s5 O3 N8 h
92.82 10 5 4 9.0625];
' j6 N- A5 O& d& j0 y2 Z" e
[ma,na]=size(A);
- g, n" A8 v8 ~* ?& h3 d3 l
B=zeros(1,na);
! l6 c ]& z0 p& d6 @9 v. F
P=zeros(1,na);
1 ] s9 U; I! K k
G=zeros(1,na);
8 U! n: |& P& d$ c+ V) M; b% i' i* P
S=zeros(1,ma);
. u2 L- j) [1 @+ Q
g=0;
/ _, i4 Q- f: C" ^+ L/ K9 i
a=0.5;
: f4 k6 l: k* \2 k' t! w. U- O6 H
MAX=max(A);
' j" o5 \* Z. G
MIN=min(A);
0 P2 c2 o: X& G5 D h2 H# q
for i=1:ma
S& V0 p t/ t& m/ H' w( F4 w% |( O/ _
for j=1:na
. C) n+ v; K w" m, C- c" @
A(i,j)=((A(i,j)-MIN(j))/(MAX(j)-MIN(j)))*a+(1-a); %标准化矩阵A
+ T' A5 C7 J7 o/ R6 s- \0 L3 Z
end
- Y1 U* K# S9 E% j
end
# C7 F$ F0 X9 \+ |! M8 ^5 J
for j=1:na
M" _' ?& N2 u! r7 n0 f
for i=1:ma
& ]1 k& `3 I/ [! z+ d) k
B(j)=B(j)+A(i,j);
1 O0 b* u! v) `. o# v! e4 D1 r
end
; A) o! l* `/ ^; y! o" @' @( W
end
) ]! K, N% p- v# }
for i=1:ma
* P/ c8 J) k# y5 E+ A
for j=1:na
+ D% J& x* j( @! x6 |& e; h. p e
A(i,j)=A(i,j)/B(j)
6 K0 T8 H' j# N4 n8 |1 w7 ?
end
+ U; ?9 W( k3 b7 K9 B q- O
end
) L" R% Q6 {4 ]
for j=1:na
& g+ e" X! s; y# I }4 R( n1 Z
for i=1:ma
& ^$ N+ K6 x4 V, o" H2 G5 J" C% e& {
P(j)=P(j)+A(i,j)*log(A(i,j));
3 _7 {6 i8 g$ H/ a4 v
end
! K8 m3 p4 z4 L: M8 D8 J+ g/ k
end
$ { s& `# H% T; r( J t2 u# m1 r
k=1/log(na);
3 D: @1 m9 s/ J9 e1 N
P=P*(-k)
' n4 W* A5 q& G4 t. U
for j=1:na
) O% j+ T% |) C$ [+ J( D5 s
G(j)=1-P(j)
2 n+ j9 ]6 G8 J% Z
g=g+G(j);
w! ?5 Z( J [( O7 q0 V- w5 D5 l
end
" m% Y4 c5 A1 C
for j=1:na
- w `9 m" v* m0 C' V& w+ s ?
b(j)=G(j)/g
8 I6 |+ e3 i2 E9 l# t) j& `8 v0 P5 x
end
( I* {$ _; s% Y) p
for i=1:ma
9 ?$ X- B9 L. p* R+ G8 G" G) ~
for j=1:na
" ]' A2 ?+ L' C; c
S(i)=b(j)*A(i,j)
% f( v( \) t6 i; O: `% _ d
end
5 G' b y& ]+ S# t; k
end
) T6 m* F8 b4 E% q
3 B8 Y0 @, Y# R, S
$ g3 B; r3 I7 }# _8 I
, N5 n: ^$ s( C. d5 Y' ~5 K' k5 Z
* ?7 P& q3 {, G2 w5 s) V$ E, D
: z& \' u; h' l# Q$ T
作者:
Reciprocal
时间:
2014-8-24 08:51
谢谢了楼主
作者:
逻辑学
时间:
2014-9-3 10:52
商权法是用来做排序的
欢迎光临 数学建模社区-数学中国 (http://www.madio.net/)
Powered by Discuz! X2.5
大家都不知道么？