数学建模社区-数学中国
标题:
求熵值法matlab程序！！！
[打印本页]
作者:
w13king
时间:
2009-9-5 11:28
标题:
求熵值法matlab程序！！！
熵值法有固定的程序可以套用的么？
0 q. i- @* X% C6 u0 v
7 T' o+ | D: \5 b9 a+ [
+ |. @3 y! f" o/ o9 c
各位好心的童鞋。。。
3 I: O5 V p# z) l |, L
第一次接触这个。。。
, }- K7 z6 d& R6 j3 l$ U! Y9 q6 N, p
不会。。。
3 w1 h, Y1 i/ w% G4 i
希望大家帮帮忙！
, T, K C6 I; d& w7 y
谢谢~~
作者:
aileenrain
时间:
2009-9-5 11:54
我也想知道！
作者:
东方明珠-WDZYQ
时间:
2009-9-5 12:21
希望能够有啊
作者:
w13king
时间:
2009-9-5 14:19
大家都不知道么？
8 k6 d3 n& P" {" Q& F3 b1 N
不会的吧...
! k U( F* ~4 `2 Q) A
在线等高人解答啊~~~
作者:
熙123
时间:
2009-9-5 15:01
我也不知道呀！！无力！！
作者:
ykl126
时间:
2009-9-5 15:36
大哥呀，下次说清楚一些好不
作者:
w13king
时间:
2009-9-5 15:49
就是那个 纯净水安全监控问题
3 g, `& O" P& B, T: F9 o
这是我们的大概方法。。。
" B. A& n s, q! a
1.对于问题1，通过分析首先确定各危害指标的权值，这里我们采取两种确定方法，即主成分分析和专家评定，考虑到因专家评定带有主观因素与事实可能不符，进一步用主成分分析法，依照客观数据给出结果，对两种方法进行综合考虑，取其平均值作为最终的权重。在此基础上建立了纯净水安全风险检测的科学评价模型，并对该城区所有批次的纯净水进行评判排序。
0 {; Q( b' V. b2 d6 K; c
2.对于问题2，在问题一模型的基础上，按照每个公司各项危害指标的平均值最为该公司的危险指标，利用问题1的模型求出各公司的风险度，并进行排序。依据各公司的各项危害指标及风险度，分析公司共产品的主要可能的危害因素。
' G# u) j4 ~' N- \5 @5 M
3.首先筛选出不合格的样品，根据问题一中建立的科学评价模型对其进行成因分析。利用熵值法求出各危害因素在生产流通环节上所占的权重。评价各个危害指标的分布规律。最后，通过四类危害指标的分析，讨论出A、B、D公司的管理情况。
3 K: j% _* @3 p! i% N# J
/ ?1 b! h7 l3 h) g5 B
就是不知道那个熵值法怎么编程啊。。。
作者:
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
, ]( l! Y/ T+ l6 }9 N! @5 P9 ~
89.98 7 5 0 7.1875
1 x5 V0 ~( t, ?% {5 P
82.73 10 5 2 8.75
; p. n* k& ^- w5 {
76.82 5 8 0 6.25
# N' Y! S0 K2 n$ v# G5 @" K
76.91 7.6 10 0 9.375
8 C2 ?& l8 Y: d- Y1 L5 z. M
85.22 10 0 6 7.5
8 i) ^! V: J2 z+ x! Z8 g1 q
80.68 5.5 8 0 8.125
9 h3 z6 [( [9 C1 r
84.59 3 0 0 4.6875
; l$ o/ g8 m3 X6 x. ?& Z' D8 V/ Z
86.05 9.5 5 3 4.375
: J- }3 E9 R# y. m |! C Y: L
83.89 10 0 0 5.3125
0 d* q3 s4 h$ i2 F0 z3 t% a
81.43 1 8 0 7.5
% C( ?0 U% G; P; d4 b
87.15 0.5 0 0 5.625
3 Q8 ^( H9 ^" K% X# x: A
73.22 1 0 0 3.75
2 n/ Z4 Q( m: B+ V; n7 j
92.82 10 5 4 9.0625];
0 m' H9 X% g) I1 n( ^8 }
[ma,na]=size(A);
& H; ~! K* G' p+ v! r
B=zeros(1,na);
" O8 p$ C* O5 Y, s! X
P=zeros(1,na);
% G2 u1 `- C5 m5 T4 {
G=zeros(1,na);
+ k6 b% P6 ?) j" E
S=zeros(1,ma);
; b: z, G7 U+ s4 [+ Q* G. i" J( j9 f% a
g=0;
% A3 O, z5 T% [) q: k/ u* B) d
a=0.5;
$ J1 ?) @/ B9 O% Z4 C: J# B# p% S/ J. U
MAX=max(A);
+ {1 y1 S( q3 W& Z; w. M
MIN=min(A);
9 o \2 ~& i3 [2 p* Z* E C
for i=1:ma
& `" f* [9 y9 @4 a/ W2 W
for j=1:na
0 P- |$ _# |7 h% @) A& B
A(i,j)=((A(i,j)-MIN(j))/(MAX(j)-MIN(j)))*a+(1-a); %标准化矩阵A
5 Y$ @" L& W' F* M! C/ a* y
end
* e0 E1 [- t5 ]' ]0 l
end
; H( i! u6 v0 d& N4 ?
for j=1:na
* ]9 q6 J5 M, P' h8 J& b) X7 T
for i=1:ma
& T! r. u2 I% `
B(j)=B(j)+A(i,j);
: _2 t1 j- ]5 U" k" ~ l
end
& i8 w, b/ j# o0 A2 l1 Q+ \3 I
end
' Y' i0 R& G5 U% k
for i=1:ma
9 H3 E6 D% O- B0 U7 _6 w4 G
for j=1:na
' i7 ]: i+ I) G' J# I2 a; ^- F: k/ F; l D; T
A(i,j)=A(i,j)/B(j)
k$ K* j" F6 U) y, j3 T6 C
end
3 t- M5 O8 H3 o2 Y
end
, Q* i3 }3 f+ j$ V) t; s
for j=1:na
: v, j' [. q/ v3 w9 B- P, R
for i=1:ma
2 ]7 A7 {( Z" f; }6 l6 ?
P(j)=P(j)+A(i,j)*log(A(i,j));
! W1 W4 a* j h t
end
, ~) W. s2 y7 v
end
# z% a( @3 W/ n1 U' V% k
k=1/log(na);
; @- a! m" ]2 D5 W$ X
P=P*(-k)
& Y8 L# `$ \2 L! w) Y
for j=1:na
1 h5 @, S4 c: [2 L! b
G(j)=1-P(j)
3 }& Y& i% l! O! @) }* u1 Z
g=g+G(j);
8 _/ t' V& M) y( G8 X4 S+ f
end
9 e! I7 A& g& J
for j=1:na
* [; y$ `( \0 C; I3 ]2 B- `9 M
b(j)=G(j)/g
; a3 z, Q0 V" v1 B, K* s
end
4 o8 _5 y) Q8 \6 ?$ o
for i=1:ma
4 I8 z% w+ Z: }8 i9 B, l5 }2 G
for j=1:na
5 C: O% k- l2 V; v
S(i)=b(j)*A(i,j)
' @! h5 b- @% \9 o/ B
end
7 j1 d0 y8 P% s/ U* d5 `& q
end
2 D1 K8 c/ L; f
- D3 |! v7 z/ m+ f- T k$ Y0 S
. F+ A. \0 F& |+ j) l" }
* ]. K, k# Z) y9 I/ E* D+ t0 W
L. d* j& W9 t& o
1 W+ g5 t' `8 p. g) u
作者:
Reciprocal
时间:
2014-8-24 08:51
谢谢了楼主
作者:
逻辑学
时间:
2014-9-3 10:52
商权法是用来做排序的
欢迎光临 数学建模社区-数学中国 (http://www.madio.net/)
Powered by Discuz! X2.5
大家都不知道么？