数学建模社区-数学中国
标题:
求熵值法matlab程序！！！
[打印本页]
作者:
w13king
时间:
2009-9-5 11:28
标题:
求熵值法matlab程序！！！
熵值法有固定的程序可以套用的么？
0 W) j2 [% @& E2 T3 n* W
, h0 |6 y o7 E4 O. F2 a. `. ~
% R5 ]9 t; ]& w% V4 C7 v
各位好心的童鞋。。。
1 X$ y2 R4 k7 B; r0 Z0 s- {1 S0 [
第一次接触这个。。。
7 {7 p) {0 V' K0 ]9 E' v( I
不会。。。
% _- T/ n* F* v1 O/ P8 z2 R0 o4 O
希望大家帮帮忙！
7 i) B& l$ p* R2 B3 c
谢谢~~
作者:
aileenrain
时间:
2009-9-5 11:54
我也想知道！
作者:
东方明珠-WDZYQ
时间:
2009-9-5 12:21
希望能够有啊
作者:
w13king
时间:
2009-9-5 14:19
大家都不知道么？
! y7 M) ]: I' ? w
不会的吧...
! @: i3 I% m% W S! p' S
在线等高人解答啊~~~
作者:
熙123
时间:
2009-9-5 15:01
我也不知道呀！！无力！！
作者:
ykl126
时间:
2009-9-5 15:36
大哥呀，下次说清楚一些好不
作者:
w13king
时间:
2009-9-5 15:49
就是那个 纯净水安全监控问题
, N, R- Q# i$ N3 a
这是我们的大概方法。。。
$ [# K+ l7 A' e, F! A) E7 v% U
1.对于问题1，通过分析首先确定各危害指标的权值，这里我们采取两种确定方法，即主成分分析和专家评定，考虑到因专家评定带有主观因素与事实可能不符，进一步用主成分分析法，依照客观数据给出结果，对两种方法进行综合考虑，取其平均值作为最终的权重。在此基础上建立了纯净水安全风险检测的科学评价模型，并对该城区所有批次的纯净水进行评判排序。
; u0 h' `1 C9 p1 z/ N& m$ v
2.对于问题2，在问题一模型的基础上，按照每个公司各项危害指标的平均值最为该公司的危险指标，利用问题1的模型求出各公司的风险度，并进行排序。依据各公司的各项危害指标及风险度，分析公司共产品的主要可能的危害因素。
+ G$ [0 v7 ~. `2 q+ M. f
3.首先筛选出不合格的样品，根据问题一中建立的科学评价模型对其进行成因分析。利用熵值法求出各危害因素在生产流通环节上所占的权重。评价各个危害指标的分布规律。最后，通过四类危害指标的分析，讨论出A、B、D公司的管理情况。
' [' s3 a% V9 k& c# ?% }
0 g0 B2 a3 e6 q9 T; W3 h% b- ^ H
就是不知道那个熵值法怎么编程啊。。。
作者:
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
* k" S2 R% h" w4 q
89.98 7 5 0 7.1875
U, e9 K3 [, |; n9 T- t, p1 `
82.73 10 5 2 8.75
% o( c( H7 O/ s" a9 R+ d2 T: g$ `
76.82 5 8 0 6.25
+ w: ^. p. s: A$ ~
76.91 7.6 10 0 9.375
$ J! m5 T: L' `
85.22 10 0 6 7.5
+ r: n, Y i! ~" Y6 V% ~2 E, a
80.68 5.5 8 0 8.125
+ W' N; y6 e( S; R
84.59 3 0 0 4.6875
8 ? @8 c4 S$ m# `
86.05 9.5 5 3 4.375
2 ] g9 T: U. T7 r0 S: J# @$ _
83.89 10 0 0 5.3125
: W" j8 A, V" Z; T2 t! L. t9 f
81.43 1 8 0 7.5
4 G. k% C0 g- ^9 E9 @) C% O' L
87.15 0.5 0 0 5.625
" P- B1 t# M# }+ U
73.22 1 0 0 3.75
0 P6 k; J1 \, g6 {3 Q' v% B6 \
92.82 10 5 4 9.0625];
- b5 V! X O \9 l, L$ _
[ma,na]=size(A);
- j( a+ J1 y# P& K/ T( W
B=zeros(1,na);
; q1 _2 @: b2 }$ Y5 ^
P=zeros(1,na);
1 Y2 y4 f) n, A1 d4 ?5 {3 i
G=zeros(1,na);
+ `( |9 k5 `2 x$ B4 A" y6 z8 T
S=zeros(1,ma);
/ G1 t8 K% C1 [8 I( y
g=0;
8 ?2 R! {# C" b. U: f- ^9 A% c' O
a=0.5;
4 j! U+ i2 v- u; o( {' S
MAX=max(A);
3 j& I( z; x2 C& { S( ?
MIN=min(A);
' X" y9 i5 l$ D8 t0 z# t5 D) I
for i=1:ma
& e6 ^4 F( ^: x( i4 l
for j=1:na
8 F( t% E$ g& X4 @! P% s- w
A(i,j)=((A(i,j)-MIN(j))/(MAX(j)-MIN(j)))*a+(1-a); %标准化矩阵A
6 ^4 e8 r. t2 f8 [$ F6 `; q- T+ ~
end
j- b# p+ _5 W: s* c9 |+ C
end
# ?- F7 ?, ]; X6 `% [
for j=1:na
. c/ N# ?" {! @6 a9 u3 m# ^
for i=1:ma
, `3 T, Z. Y5 h O) |' B
B(j)=B(j)+A(i,j);
# m$ Q, y1 p6 F: j
end
( r+ j0 D7 ]: E& E ?* b
end
3 h6 Q8 o9 e& f# O; `% }
for i=1:ma
) p, y0 `1 k2 h5 s) n
for j=1:na
* t+ V, v w+ S o3 ~4 G* j; B" ]# U
A(i,j)=A(i,j)/B(j)
: V3 w* Z! e* \4 w3 c# s
end
3 D( P6 t8 I/ U7 G9 h
end
& r8 A6 K* a6 \4 V! T! ]' {. q
for j=1:na
$ q5 x9 k- \% x. J& w5 t
for i=1:ma
6 W! T4 [7 v) T
P(j)=P(j)+A(i,j)*log(A(i,j));
. V3 {5 l! N4 U0 L; \' `$ ^! E, \
end
/ u" w3 I" R/ a2 I* R, H# R
end
# w' Q# `2 ?2 |4 b$ a( N5 c
k=1/log(na);
$ t* W& L6 h% |3 e4 Q& }0 @
P=P*(-k)
3 i6 `9 ^/ P a: I5 j: ]
for j=1:na
A* D& `8 w9 E( G* @2 i5 _
G(j)=1-P(j)
7 Z) j7 l+ V/ K( O0 x9 p
g=g+G(j);
' f7 [% w3 y- c( x0 _
end
# O$ i( \: r6 s, M4 e5 K
for j=1:na
9 F4 D8 [. l- ?9 x9 ]
b(j)=G(j)/g
7 u7 a4 |; ^2 v, o: f6 h
end
, x6 o% `5 v3 Y( O0 p) B
for i=1:ma
( Q5 R8 Z- L9 C8 M6 {6 E2 \+ u4 c
for j=1:na
4 r# t5 P0 q5 P7 x$ V: w
S(i)=b(j)*A(i,j)
7 K7 k/ w) L1 P) Z- H3 K9 O& O2 n# g
end
& b5 y/ ]7 B5 m$ Y2 y* W3 k8 V1 ~ H
end
5 t2 K& b5 ?' G3 l: `" [( I
8 @/ K1 S' t! j: i5 z1 b
q d( E& P+ P- ^0 k
. G! l% F i5 {3 K- }
) y2 Q0 o/ Q) x# R8 ^1 |$ @+ `
% {6 E! U1 m% P b0 d; F [, n/ h
作者:
Reciprocal
时间:
2014-8-24 08:51
谢谢了楼主
作者:
逻辑学
时间:
2014-9-3 10:52
商权法是用来做排序的
欢迎光临 数学建模社区-数学中国 (http://www.madio.net/)
Powered by Discuz! X2.5
大家都不知道么？! y7 M) ]: I' ? w