数学建模社区-数学中国
标题:
求熵值法matlab程序！！！
[打印本页]
作者:
w13king
时间:
2009-9-5 11:28
标题:
求熵值法matlab程序！！！
熵值法有固定的程序可以套用的么？
. n, v" z: J8 I8 [/ M/ N. d
' g j6 k$ I! `, u! Q# `0 C
1 i) k7 t5 m8 K
各位好心的童鞋。。。
5 r1 ^' e/ m% \/ X" A1 X
第一次接触这个。。。
% v V; B+ N- F$ n
不会。。。
% o, r7 v- u1 d" c
希望大家帮帮忙！
7 F, p% @' L/ c/ w, n' e+ H
谢谢~~
作者:
aileenrain
时间:
2009-9-5 11:54
我也想知道！
作者:
东方明珠-WDZYQ
时间:
2009-9-5 12:21
希望能够有啊
作者:
w13king
时间:
2009-9-5 14:19
大家都不知道么？
x& j: Q& d3 P( u
不会的吧...
/ P6 U& d& P0 A. G( j8 A9 `: F1 G
在线等高人解答啊~~~
作者:
熙123
时间:
2009-9-5 15:01
我也不知道呀！！无力！！
作者:
ykl126
时间:
2009-9-5 15:36
大哥呀，下次说清楚一些好不
作者:
w13king
时间:
2009-9-5 15:49
就是那个 纯净水安全监控问题
2 c% f+ o- y/ a; q7 L
这是我们的大概方法。。。
4 L3 }* v4 U0 h; |
1.对于问题1，通过分析首先确定各危害指标的权值，这里我们采取两种确定方法，即主成分分析和专家评定，考虑到因专家评定带有主观因素与事实可能不符，进一步用主成分分析法，依照客观数据给出结果，对两种方法进行综合考虑，取其平均值作为最终的权重。在此基础上建立了纯净水安全风险检测的科学评价模型，并对该城区所有批次的纯净水进行评判排序。
" t) C3 o; z) O9 ]* B, B# }
2.对于问题2，在问题一模型的基础上，按照每个公司各项危害指标的平均值最为该公司的危险指标，利用问题1的模型求出各公司的风险度，并进行排序。依据各公司的各项危害指标及风险度，分析公司共产品的主要可能的危害因素。
+ v0 w1 U% i# t# f6 s$ t1 w* e' \
3.首先筛选出不合格的样品，根据问题一中建立的科学评价模型对其进行成因分析。利用熵值法求出各危害因素在生产流通环节上所占的权重。评价各个危害指标的分布规律。最后，通过四类危害指标的分析，讨论出A、B、D公司的管理情况。
" C- E3 o0 h4 Z# Q$ u m
/ b7 L$ o" A6 [2 a( I
就是不知道那个熵值法怎么编程啊。。。
作者:
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
) A* {# Y( T. ^# s7 v! s
89.98 7 5 0 7.1875
% w# w+ s1 n8 f) s. }$ C
82.73 10 5 2 8.75
$ v5 C1 Z* ~& z. S: O0 ?4 W5 t5 k
76.82 5 8 0 6.25
% O- T) Y9 [& v. h2 B) k, F5 Q
76.91 7.6 10 0 9.375
4 `$ ? c+ w. g6 O+ M
85.22 10 0 6 7.5
8 f) m7 u) i3 C) G4 r
80.68 5.5 8 0 8.125
5 W" M6 y: f4 |) E& w# e0 t
84.59 3 0 0 4.6875
' s! v" ^; J; y6 [# Q5 U
86.05 9.5 5 3 4.375
: n# C7 `' E% j( g, c: S) T4 T
83.89 10 0 0 5.3125
' W" b8 z0 s4 b/ g$ G# I2 I
81.43 1 8 0 7.5
2 X1 D6 ]9 C' d# P- K5 e
87.15 0.5 0 0 5.625
/ q$ F2 U" e5 Y' Z4 v
73.22 1 0 0 3.75
7 r: u5 ?' _$ l/ x4 Z1 n5 k7 E
92.82 10 5 4 9.0625];
; j9 l6 ^% t: D
[ma,na]=size(A);
# S% b5 `. W) ?8 V4 H6 k3 {; c5 ]
B=zeros(1,na);
4 f/ X) L0 z3 v! E
P=zeros(1,na);
/ Q: r& t) |4 z0 Q' I8 `
G=zeros(1,na);
& r. H x+ ^( [
S=zeros(1,ma);
- V. f& H1 Q% A" m4 x# _
g=0;
. V- U: o6 `0 l, N; D! y
a=0.5;
! n6 O- _+ M5 C( j2 q9 N% E
MAX=max(A);
! F' Y R* S% d. p% U
MIN=min(A);
, o5 Y0 R' D4 \9 `9 J8 d5 i! O
for i=1:ma
6 V! H# D, { l
for j=1:na
+ Q( A9 G/ @5 {0 a4 D
A(i,j)=((A(i,j)-MIN(j))/(MAX(j)-MIN(j)))*a+(1-a); %标准化矩阵A
1 U Y$ t: m1 Q
end
4 ^2 ]1 A( c5 e% C7 _+ P4 F4 |2 w
end
! U5 b6 L/ {, c5 ^
for j=1:na
/ k+ A5 D0 G3 L. z/ S6 A2 R" Z5 R
for i=1:ma
- F5 u: K' h" R7 t! q5 ]& S( X
B(j)=B(j)+A(i,j);
" b# m4 Z, I) W# m* x1 o
end
+ @: H- D4 l% D. @! u
end
+ } t, b3 v& w* Y5 E
for i=1:ma
* Y9 |2 \9 @+ F2 J8 F" v* A
for j=1:na
4 F7 P# H6 g& b- @- i. R
A(i,j)=A(i,j)/B(j)
& O* E8 d3 ~. h" u* g! @
end
; f r0 J5 q8 p9 p& N1 m2 s
end
4 ]4 j& t& {+ u* g
for j=1:na
! y& M' K, U/ ?
for i=1:ma
1 R* v! b2 C/ ~
P(j)=P(j)+A(i,j)*log(A(i,j));
& n$ A- ~. a8 a/ g9 q0 z3 x
end
6 X# Y- t3 y0 Z
end
0 ^$ M4 @' K: I$ ]
k=1/log(na);
- ?) m( m) j, t* i; ?5 e
P=P*(-k)
$ V& v, m f, _3 c: L
for j=1:na
( ~/ J3 D1 _5 \- a m7 g1 F7 d5 M) J& k
G(j)=1-P(j)
2 e7 E- } |/ V
g=g+G(j);
$ E2 H% D' d4 h+ B
end
. ~9 G. ?" r4 J
for j=1:na
; f" s# x$ L4 \+ V6 n6 Q4 z6 A0 K$ C
b(j)=G(j)/g
6 Y0 a% V. h+ [9 R
end
& a7 w* ]6 X& R5 g+ e
for i=1:ma
8 y! |5 \% _3 H( }5 K4 L$ z6 ]; J
for j=1:na
1 N2 Z9 M% F9 |$ G
S(i)=b(j)*A(i,j)
$ y9 p5 i6 t9 E d3 @* s5 P# v) V
end
. @( A' r; a9 n' t* @
end
7 ~- C& g5 P3 M
: {- C0 a; J: r
+ t; ?) W' k8 L7 h, r( ]; b N1 Q$ q
3 k" v% Z. a+ u3 ]/ s+ A7 M& e Y
$ D3 F/ L; O M" S
' P, C' `; K; y. t- B/ V: A
作者:
Reciprocal
时间:
2014-8-24 08:51
谢谢了楼主
作者:
逻辑学
时间:
2014-9-3 10:52
商权法是用来做排序的
欢迎光临 数学建模社区-数学中国 (http://www.madio.net/)
Powered by Discuz! X2.5
大家都不知道么？ x& j: Q& d3 P( u