- 在线时间
- 0 小时
- 最后登录
- 2010-11-19
- 注册时间
- 2010-7-31
- 听众数
- 4
- 收听数
- 0
- 能力
- 0 分
- 体力
- 10 点
- 威望
- 0 点
- 阅读权限
- 20
- 积分
- 23
- 相册
- 0
- 日志
- 0
- 记录
- 0
- 帖子
- 42
- 主题
- 5
- 精华
- 0
- 分享
- 0
- 好友
- 0
升级   18.95% 该用户从未签到
 群组: 华南理工大学 |
一、使用曲线拟合模拟病人等待时间的模拟matlab程序实现:
2 x$ q, C+ f7 T[A]=xlsread('C:\MATLAB7\work\问题3数据5');/ V6 q8 E* N; i& M9 i+ D
len=length(A);# J5 v ]5 h6 X- |& C9 }1 h9 }
i=1:len;
' P5 A- l3 _) u4 gj=A(:,5);+ _ {/ x1 Q9 ]$ M% z% R
figure(1);: C9 t+ x& a0 @
plot(i,j,'.-');
& b- N8 k" A% m2 Phold on;$ W7 D2 x) c7 Z( |
i=i';
U( X4 T0 e' Z' z' K. O9 Mp1=polyfit(i,j,2)2 E$ ^ i, @! W7 k D
f1=polyval(p1,i);5 L* ?7 A6 B- s4 h! l
plot(i,j,'.',i,f1,'-');
8 T$ p2 B3 t& t* T8 S- [# q: Uhold on;3 t [2 G8 w! ] @' i# T
title('病人平均等待时间与病人到达时间的关系');
- ]2 q/ R+ P6 Lxlabel('病人门诊时间距离7 月12 日的时间差');ylabel('病人平均等待时间$ q. O( q, u0 f- C% ^
');* ?) t& L1 u, e" |" o
axis equal;
) S: w" {1 U# ~figure(2);
, e3 U6 J1 k% d+ B) R* j+ q# C/ M5 GA1=[1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 17 18 21 22 23 24 25 26 28 29 31
5 X% M: ?7 o3 _, `0 W32 34 35 36 37 38 39 40 41 42];3 b: v& y6 }9 E" W
B1=[12 12 12 12 12 12.5 13 13 13 13 13.4 14 14 13 12.5 11 11 12 12 12 11
# L. ^/ E; E9 b$ ]0 F11.5 13.25 14 13 12 11 13 13 12 12.25 13 15.5 15 14 13];
" L% N% Q; ]+ m' X, w' Y# |subplot(1,2,1);
( C: E) F* S9 U# v5 |8 ]& ^plot(A1,B1,'.-');
! |9 c! h' U5 M5 W( x. l" F& nhold on;
" J- Y9 x8 p) X! t9 [2 mp2=polyfit(A1,B1,2)
. ~- b2 I! p- a2 C/ l' df2=polyval(p2,A1);
9 D0 k M% t2 ]6 w% Xplot(A1,B1,'.',A1,f2,'-');! B4 s& _9 X, r, \4 D; i' h
hold on;: D' m( b, L$ ~3 F6 Z
title('白内障病人情况');4 }: k9 f" `$ R8 T H
xlabel('白内障病人门诊时间距7月12日的时间差');ylabel('白内障病人平均) m, b, F; B) ~
等待时间');) Z- @ w7 V" I1 t5 ^2 v% S
axis equal;7 b: W2 p, f# l( f$ J
A2=[2 3 4 5 6 7 8 9 11 13 15 16 17 18 19 20 21 22 23 25 26 27 28 29 30 q: _* L& u6 X t
31 32 33 34 35 36 37 ];
" R0 p: |5 P. U0 S; r' D+ bB2=[12 12.5 12 12.5 12.5 13 13 13.25 14 14 12.5 12 11 11 10 11.75 12 11.67( w |; X, q- H7 B# a3 d' r2 m" X
12 12 13.25 14 14 13 12 12 12.4 13.33 13 13 13 12];* w- i6 e1 v9 e/ s( e% ^- R
subplot(1,2,2);9 [ e# z, b5 J5 |1 c. I
plot(A2,B2,'.-');% |8 J5 o) C' J6 [
hold on;1 ?) ^8 P- n- B
p3=polyfit(A2,B2,2) Z% m4 q2 u' m- E: \- K6 T- L
f3=polyval(p3,A2);
2 a" s3 @' {# {plot(A2,B2,'.',A2,f3,'-');& E7 z+ u( N. A/ I, n" L) y
hold on;
; M8 |" Y- H" s$ u$ ptitle('白内障双眼病人情况');
. t6 h0 N0 V4 k* F' ~xlabel('白内障双眼病人门诊时间距7月12日的时间差');ylabel('白内障双眼
& r' t- S' M, u% y' ]9 w$ t3 D病人平均等待时间');! H$ l" q0 w N6 J
axis equal;' a6 m+ O q' c* w$ T: ]8 }& U8 C) i
figure(3);
- X6 ], P4 ^+ B) [A3=[1 2 3 4 5 6 7 8 10 11 12 13 15 16 17 18 19 20 21 22 23 24 25 26 27
& t1 T \; [) S' K4 u. }28 29 30 31 32 33 34 35 36 37 38 ];
8 D( h/ y0 a( X7 F% eB3=[12.33 12.5 12.3 12 11.67 12.33 12 13 13 14.2 14.25 14 13 12 11 11 10.75
- Q# X, ^ {/ ?& {7 N# B11.67 12 12.25 12 11 11.75 13 14 14 13 12 12 12.33 13 13 13 12.5 12 12];& Q# N+ M; R% ?; p, U$ C
subplot(2,2,1);
: G3 ~, B1 H1 W& O4 Nplot(A3,B3,'.-');8 N& q! |) F/ V6 u2 J4 Q* j7 h% Z
hold on;& v+ U* r. z/ P7 v. P
p4=polyfit(A3,B3,2)- U! j. ]1 b$ ^% R/ x* {/ w
f4=polyval(p4,A3);" T! k% K& E/ _$ |
plot(A3,B3,'.',A3,f4,'-');
/ Y& m# v; m; ^$ shold on;& a. i l# j% `: e
title('视网膜病人情况');
7 Z8 N8 Z# Y* k) y+ \xlabel('门诊时间距7 月12 日的时间差');ylabel('视网膜病人平均等待时间
) k4 U, L; t* A( ~6 l');
9 Q0 _2 C- j& c" M& f3 {6 g/ haxis equal;" { c0 U" [9 g1 z$ R% D( H
A4=[1 2 3 4 5 6 7 10 11 12 14 15 16 17 19 20 22 24 26 27 29 30 32 33 34, q! g5 u- R5 O: X* @
37];" ?3 ]# V- j1 T/ r/ m
B4=[12 13 12 12 12 12 13.5 13 13.5 14 13 12 12 11 10.25 12 11 11 15 15
- M! h) P( D) W) D6 @* `13 12 12 13 13 12];
R. C* J8 k, o# ^: ? y3 Tsubplot(2,2,2);1 L" @& P8 Q. C" F' I1 t
plot(A4,B4,'.-');" ~* E7 z: L) ?2 J! Q E7 l
hold on;
1 R" C: f% k6 A1 J6 yp5=polyfit(A4,B4,2)6 _% e: e) ? r/ t
f5=polyval(p5,A4);) l/ k4 k3 i% L
plot(A4,B4,'.',A4,f5,'-');3 }1 @$ g) V) _& {1 B5 ^
hold on;
* Z4 D- E \" F0 J& Q. Stitle('青光眼病人情况');
. @! w+ V0 i$ ~! `+ uxlabel('门诊时间距7 月12日的时间差');ylabel('视网膜病人平均等待时间6 Z# q" T h9 g6 v% u
');$ k3 q) N8 O2 L- J- I
axis equal;
0 ?' F3 O- H1 ^) s0 x. |' e2 mA4=[1 2 3 4 5 6 7 10 11 12 14 15 16 17 19 20 22 24 26 27 29 30 32 33 34
0 w( u6 v" F$ Z% X7 L9 i- V37]; w9 ?% C. i% u
B4=[12 13 12 12 12 12 13.5 13 13.5 14 13 12 12 11 10.25 12 11 11 15 159 H- ]8 J. v5 c. J) [, h
13 12 12 13 13 12];
, ?3 H+ v* j! ?0 ]% |subplot(2,2,2);
/ X, G. z2 @, D. c6 m/ F/ g6 ]1 z2 ?plot(A4,B4,'.-');! c, U6 w0 i; R) M; a. Z
hold on;# j$ N) _% L! q- Y! `
p5=polyfit(A4,B4,2)
4 _$ }: J' w' e1 B4 [f5=polyval(p5,A4);& k) l3 D6 L+ e! F- K/ P
plot(A4,B4,'.',A4,f5,'-');
, U" m- ]% ~( E5 w5 P5 G1 K- thold on;
% g7 ^3 j7 o& Vtitle('青光眼病人情况');
( `5 O- A$ Q* e9 X6 w: Zxlabel('门诊时间距7 月12 日的时间差');ylabel('青光眼病人平均等待时间! @0 l$ b. S+ t+ |- U- o* T7 I
');
7 r$ v& E+ Z# n( O2 E* Q7 faxis equal;' C% x' r# ]$ \" }2 h" w2 a
A5=[1 43 5 6 7 8 9 10 11 15 16 18 19 21 22 23 26 28 29 30 32 36 38 39 40];/ q- L* c# J7 ?" ]; p* t8 C
B5=[1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1];. u1 Y- {1 U) W
subplot(2,2,3);( `3 x i Z1 j f( A& s
plot(A5,B5,'.-');
; R: R7 I5 ]9 E& x' f7 J3 Ghold on;
& \3 a) w5 w) x d& bp6=polyfit(A5,B5,2)
4 b8 `) x6 {4 U; j; c) mf6=polyval(p6,A5);
[5 N; m8 H5 L6 T( w: F+ o Nplot(A5,B5,'.',A5,f6,'-');/ Z2 c) T: ]4 O0 F; _" Q- D# w* k
hold on;
' A3 w3 {: R. H& a9 o/ r1 qtitle('外伤病人情况');
7 d/ S- d, f5 T* M1 V' wxlabel('外伤病人门诊时间距离7月12日的时间差');ylabel('外伤病人平均等
( b4 _ _' E' ?$ G: c" z, i待时间');0 D4 Q4 s% c- R9 D6 Q3 r5 z8 l+ ?
axis equal;
' l* _3 h: _1 P二、层次分析法的一致性指标判定的matlab 实现:
: g e' n6 D; B2 T/ n: c5 tfunction H=ff(A)/ c" w* C4 q1 u' U+ L4 T
t=eig(A). J, J5 n) p) A+ Y" i
w=max(t)5 T* _7 j6 e) W w9 r
q=(w-4)/3/.9. B% S) V4 [7 @$ V3 B# F
s=[(A(1)+A(5)+A(9)+A(13))/(A(1)+A(2)+A(3)+A(4)+A(5)+A(6)+A(7)+A(8)+A(
3 m$ p( C+ Z z" t4 U. d" T, T! b- C9)+A(10)+A(11)+A(12)+A(13)+A(14)+A(15)+A(16))8 f! c# w% L+ e7 U
(A(2)+A(6)+A(10)+A(14))/(A(1)+A(2)+A(3)+A(4)+A(5)+A(6)+A(7)+A(8)+A(9)/ e- T0 [" K7 }4 Z/ O* z
+A(10)+A(11)+A(12)+A(13)+A(14)+A(15)+A(16))
- w! d/ l5 I$ s f' i6 z(A(3)+A(7)+A(11)+A(15))/(A(1)+A(2)+A(3)+A(4)+A(5)+A(6)+A(7)+A(8)+A(9)1 @: \+ g/ k1 Z, F' E! ?# V- [
+A(10)+A(11)+A(12)+A(13)+A(14)+A(15)+A(16))
7 a' ~. x7 G, f# C' z2 `5 ](A(4)+A(8)+A(12)+A(16))/(A(1)+A(2)+A(3)+A(4)+A(5)+A(6)+A(7)+A(8)+A(9)- J% D2 p5 s* ?9 s5 D
+A(10)+A(11)+A(12)+A(13)+A(14)+A(15)+A(16))]
6 M% C/ \& v, e6 Y0 S* ~3 t |
zan
|