数学建模社区-数学中国
标题:
求熵值法matlab程序!!!
[打印本页]
作者:
w13king
时间:
2009-9-5 11:28
标题:
求熵值法matlab程序!!!
熵值法有固定的程序可以套用的么?
& J& l* }- N1 M+ c- ~# q
' t. @6 U1 J: @
7 ?% j/ u5 H, W$ N. u8 i) \6 X
各位好心的童鞋。。。
7 s: {! k% \# L) _3 _
第一次接触这个。。。
# K" ?9 _1 u5 }8 g5 v$ `9 _; D
不会。。。
& ~; c, b1 @, p0 g* \
希望大家帮帮忙!
/ r1 @* B# Q1 ]& u
谢谢~~
作者:
aileenrain
时间:
2009-9-5 11:54
我也想知道!
作者:
东方明珠-WDZYQ
时间:
2009-9-5 12:21
希望能够有啊
作者:
w13king
时间:
2009-9-5 14:19
大家都不知道么?
5 h- I8 k+ O0 ?& F
不会的吧...
3 y8 d2 m% ]. a4 Q0 l0 I3 I: }) ~
在线等高人解答啊~~~
作者:
熙123
时间:
2009-9-5 15:01
我也不知道呀!!无力!!
作者:
ykl126
时间:
2009-9-5 15:36
大哥呀,下次说清楚一些好不
作者:
w13king
时间:
2009-9-5 15:49
就是那个 纯净水安全监控问题
* A; d9 G$ ]8 {, V
这是我们的大概方法。。。
, ~9 S5 L4 X1 Z) X
1.对于问题1,通过分析首先确定各危害指标的权值,这里我们采取两种确定方法,即主成分分析和专家评定,考虑到因专家评定带有主观因素与事实可能不符,进一步用主成分分析法,依照客观数据给出结果,对两种方法进行综合考虑,取其平均值作为最终的权重。在此基础上建立了纯净水安全风险检测的科学评价模型,并对该城区所有批次的纯净水进行评判排序。
1 \( I$ X: ^6 w& f5 t2 R
2.对于问题2,在问题一模型的基础上,按照每个公司各项危害指标的平均值最为该公司的危险指标,利用问题1的模型求出各公司的风险度,并进行排序。依据各公司的各项危害指标及风险度,分析公司共产品的主要可能的危害因素。
& M+ Q% L) E' K0 P l
3.首先筛选出不合格的样品,根据问题一中建立的科学评价模型对其进行成因分析。利用熵值法求出各危害因素在生产流通环节上所占的权重。评价各个危害指标的分布规律。最后,通过四类危害指标的分析,讨论出A、B、D公司的管理情况。
, a* Y" b) a7 c$ A3 E$ @: u' y
$ p2 q4 r Y; L
就是不知道那个熵值法怎么编程啊。。。
作者:
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
! A0 E# W/ ^8 k6 Z
89.98 7 5 0 7.1875
) T" h( U3 I y
82.73 10 5 2 8.75
! L" i1 i0 W" l
76.82 5 8 0 6.25
/ R& E; v; m3 H
76.91 7.6 10 0 9.375
1 t6 B$ e9 Y# r4 `" v
85.22 10 0 6 7.5
6 W5 I' v$ j, P7 M3 I+ a2 E( D' [
80.68 5.5 8 0 8.125
; Y% A3 O* N, h- \" \
84.59 3 0 0 4.6875
% }& |. v( ^- B; ~8 z
86.05 9.5 5 3 4.375
7 P$ h$ ~7 y8 b
83.89 10 0 0 5.3125
$ x1 @7 A( l9 D' g' x
81.43 1 8 0 7.5
% ^+ c, w8 H. U! `8 Y* F
87.15 0.5 0 0 5.625
; K0 S+ D+ c8 u% X5 T8 a
73.22 1 0 0 3.75
' I5 k; J% l* V/ W" d! Q( b
92.82 10 5 4 9.0625];
+ `9 h1 ?+ n# p( j8 r1 E6 E6 B
[ma,na]=size(A);
% ^: x7 ]5 ^3 d* s. D; d5 }
B=zeros(1,na);
9 x% ]- R! o4 t# k
P=zeros(1,na);
! r! s8 c* i2 h$ H- A6 d
G=zeros(1,na);
( @9 b) s% x" |$ Y2 ?
S=zeros(1,ma);
4 u( ~2 @+ Y$ {6 b- _$ n
g=0;
* e7 l( K" r- N( H
a=0.5;
. [2 I$ E" a+ s
MAX=max(A);
9 A3 U5 F& L$ q; }3 R, H0 @; B
MIN=min(A);
1 w3 I6 G7 w y4 X0 A9 i
for i=1:ma
; o" {+ Q8 [* U" C& g
for j=1:na
- I D2 N s3 ^
A(i,j)=((A(i,j)-MIN(j))/(MAX(j)-MIN(j)))*a+(1-a); %标准化矩阵A
- G$ ?; D8 a" I {1 c
end
! A0 ~. m6 ~4 q. M& C
end
4 `: p9 e `' X
for j=1:na
) Z" Q3 Y* Q4 M1 n
for i=1:ma
1 ]/ M8 p( n% D% {' }8 n8 }
B(j)=B(j)+A(i,j);
( T1 l$ {8 @% A n
end
2 m1 c0 s, e4 u9 Z: F% h
end
! Y% z5 f1 h8 [
for i=1:ma
i" @9 J7 ^( T: `
for j=1:na
0 \$ i& _5 u: X4 m
A(i,j)=A(i,j)/B(j)
( z' l" ]6 _. s5 L
end
9 k4 `. C1 m' t. F) ~. g( d
end
7 J- T' ^% q$ z0 a. }& m
for j=1:na
" a5 U# O+ v4 N# B: ^$ y) i7 m
for i=1:ma
$ }5 D6 R4 ~+ a* _; H, f
P(j)=P(j)+A(i,j)*log(A(i,j));
; H- ?. | \3 }6 A
end
6 u. D: ~$ [2 p- T$ B
end
: L f$ e! a8 I& g8 G
k=1/log(na);
0 q4 u+ u# x6 p0 }% U
P=P*(-k)
1 k0 t( v' I) q% ]: o# h( ~6 M
for j=1:na
2 g: u& y6 D4 r- E+ T6 B
G(j)=1-P(j)
8 k$ X5 z( d. j: R$ V1 ]
g=g+G(j);
' G6 x8 M/ m5 t$ b/ A
end
+ P" G. f3 |0 a7 W. _
for j=1:na
- @* j! r/ N) C1 F
b(j)=G(j)/g
6 I- ~0 G0 A, m7 y; o' \7 f$ E
end
: T) k; K5 I6 v3 K9 E$ Z$ z0 c$ p d( g
for i=1:ma
s+ s) K( k$ f% o8 p9 i( y
for j=1:na
/ V- _% u0 k; |# V6 |$ [
S(i)=b(j)*A(i,j)
0 X, q. I0 o" _& H( W/ R( j
end
u0 l7 }& f. f3 z+ Z
end
1 B% z8 _; z# O/ J% x: s# P* }1 ?
5 [/ ~7 P# L- L
) @" v( C( G5 Y0 O$ Q
4 v# e& C2 ^! x# M8 Q2 b
8 v& x+ v9 T8 [/ N1 ~0 T
[8 N& I# Z# k
作者:
Reciprocal
时间:
2014-8-24 08:51
谢谢了楼主
作者:
逻辑学
时间:
2014-9-3 10:52
商权法是用来做排序的
欢迎光临 数学建模社区-数学中国 (http://www.madio.net/)
Powered by Discuz! X2.5