数学建模社区-数学中国
标题:
求熵值法matlab程序！！！
[打印本页]
作者:
w13king
时间:
2009-9-5 11:28
标题:
求熵值法matlab程序！！！
熵值法有固定的程序可以套用的么？
5 o; x' B$ ~ F2 t' y
2 I/ Y) c) H5 u/ a# s. q
7 T; }+ D; G0 V/ n" m
各位好心的童鞋。。。
# L* @+ o3 z5 P" b& ^& r
第一次接触这个。。。
- e9 d3 U1 `5 w! R8 x
不会。。。
/ N0 q1 ?* U, @
希望大家帮帮忙！
' Y- N1 L3 h3 J P2 ?$ I" L5 ^
谢谢~~
作者:
aileenrain
时间:
2009-9-5 11:54
我也想知道！
作者:
东方明珠-WDZYQ
时间:
2009-9-5 12:21
希望能够有啊
作者:
w13king
时间:
2009-9-5 14:19
大家都不知道么？
9 [% ^% {( O; u7 ], P
不会的吧...
; t8 s# u- B4 P/ ? P& ], Z: ~
在线等高人解答啊~~~
作者:
熙123
时间:
2009-9-5 15:01
我也不知道呀！！无力！！
作者:
ykl126
时间:
2009-9-5 15:36
大哥呀，下次说清楚一些好不
作者:
w13king
时间:
2009-9-5 15:49
就是那个 纯净水安全监控问题
" Y. ~2 }# `$ P6 T+ a
这是我们的大概方法。。。
& r* N# w8 e% U# }
1.对于问题1，通过分析首先确定各危害指标的权值，这里我们采取两种确定方法，即主成分分析和专家评定，考虑到因专家评定带有主观因素与事实可能不符，进一步用主成分分析法，依照客观数据给出结果，对两种方法进行综合考虑，取其平均值作为最终的权重。在此基础上建立了纯净水安全风险检测的科学评价模型，并对该城区所有批次的纯净水进行评判排序。
5 t+ c) m" T1 g1 Y5 i( Y6 q
2.对于问题2，在问题一模型的基础上，按照每个公司各项危害指标的平均值最为该公司的危险指标，利用问题1的模型求出各公司的风险度，并进行排序。依据各公司的各项危害指标及风险度，分析公司共产品的主要可能的危害因素。
2 U: J4 {0 Z" S7 ^7 j
3.首先筛选出不合格的样品，根据问题一中建立的科学评价模型对其进行成因分析。利用熵值法求出各危害因素在生产流通环节上所占的权重。评价各个危害指标的分布规律。最后，通过四类危害指标的分析，讨论出A、B、D公司的管理情况。
n3 K" \; C- E
( M# r( M) R1 F' w9 K
就是不知道那个熵值法怎么编程啊。。。
作者:
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
! K1 A( u: x, V9 e* p0 O
89.98 7 5 0 7.1875
' q0 s3 j7 d; b. Y! F
82.73 10 5 2 8.75
X" h3 d" ?$ o" \+ J& z$ i
76.82 5 8 0 6.25
! q% ]3 k" h0 r1 j* j6 L& j
76.91 7.6 10 0 9.375
8 `/ q u$ H) ^8 E; j5 n( I' B
85.22 10 0 6 7.5
/ A1 S6 k: |# X! g7 n2 S0 R( g
80.68 5.5 8 0 8.125
' }! P( o U0 w A7 w; |
84.59 3 0 0 4.6875
8 M+ u0 s) z. x2 J( W( I
86.05 9.5 5 3 4.375
, F+ f9 w- X* _8 D( C1 V- |% S9 M7 [
83.89 10 0 0 5.3125
8 W5 y; s6 G' e# p- }5 N
81.43 1 8 0 7.5
; m- G l; L/ K+ }
87.15 0.5 0 0 5.625
' H7 a4 j: y5 a( }
73.22 1 0 0 3.75
( d* S' p/ s5 A0 m
92.82 10 5 4 9.0625];
, O- c4 @& e5 ~' J- {* s
[ma,na]=size(A);
7 F4 R! n2 q% j7 _. d+ s: E
B=zeros(1,na);
3 _/ I7 E. N. T: Q" Z. t8 |
P=zeros(1,na);
+ t" n, O9 B9 a5 X1 K. {
G=zeros(1,na);
9 W) R4 }. u) F+ }+ m6 t! v% ]' ^+ ^
S=zeros(1,ma);
& V* ~% a9 U6 C6 g8 T5 P
g=0;
2 S5 d9 k! S* ^4 S7 _5 q
a=0.5;
# h" c2 e. j4 V
MAX=max(A);
+ x v! A0 y, i+ z
MIN=min(A);
- |) }: q6 v S. x
for i=1:ma
# q/ M; ^- ?$ n/ l
for j=1:na
# y9 s6 Y+ I. R) V# E
A(i,j)=((A(i,j)-MIN(j))/(MAX(j)-MIN(j)))*a+(1-a); %标准化矩阵A
% G8 x+ V! a. G; J# ]( H" b" S
end
! B. L5 b5 M4 u9 r" r
end
2 }! q/ C C. w/ q! K& L( i8 c) w
for j=1:na
, @" ]9 |6 `% t9 C& {/ g9 T7 ~; E! ]
for i=1:ma
# b$ Q+ N+ C# \. e! J
B(j)=B(j)+A(i,j);
: {5 ^6 ~. U: ^! x0 A
end
/ [- n G5 s% q# w/ [
end
$ Y, @- {9 {$ A
for i=1:ma
! w5 N, r' v# t# B
for j=1:na
3 U' H. I- H3 J8 @: d
A(i,j)=A(i,j)/B(j)
0 ^0 {3 i( H9 F4 Z( }
end
; E6 l% C- K6 W% N- C* T
end
0 Z% D ]3 S: X# l
for j=1:na
/ \9 P" E7 @% H9 h& G' S9 V" M6 y$ B
for i=1:ma
8 d. N h& |1 t" a; B2 F
P(j)=P(j)+A(i,j)*log(A(i,j));
S" r. ?% R1 B I+ g
end
. j: m+ O2 f) T
end
: X; M/ x3 u* [- u$ f" H, I" ^, v
k=1/log(na);
! l$ p5 w" {3 q
P=P*(-k)
7 B+ a- w* _6 q8 Z5 v! j
for j=1:na
: Y; x4 s3 i, d$ Z3 L) @/ `
G(j)=1-P(j)
+ y3 O$ }5 d4 Q, h* L' M2 X
g=g+G(j);
, } o' X( G) c2 ?" |: {+ |4 `- X
end
4 C. d4 T- i; y) H
for j=1:na
4 s. p+ T1 W( @# {/ p6 n5 }8 [
b(j)=G(j)/g
% x( q7 q0 c: Z4 U- C
end
0 F: B9 ^+ F& Q2 l+ |8 k; ?, g
for i=1:ma
+ P v. e' h7 f. h% `9 b) F" Q4 u
for j=1:na
8 Z3 W8 e9 T7 ]; u0 X1 l
S(i)=b(j)*A(i,j)
- \$ u* W# K% L( M9 h2 D
end
; V, v. `& |! B& S
end
% n2 M6 r( o3 t5 l' O
8 T8 M4 ~; d) i; p
3 ]* U; s; d/ F9 R4 X! ]
6 Z& S8 f8 ^2 @" l& h" ~* o" U8 }
- e% W* d5 R! K
0 E7 G% K/ P- r
作者:
Reciprocal
时间:
2014-8-24 08:51
谢谢了楼主
作者:
逻辑学
时间:
2014-9-3 10:52
商权法是用来做排序的
欢迎光临 数学建模社区-数学中国 (http://www.madio.net/)
Powered by Discuz! X2.5
大家都不知道么？