- 在线时间
- 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程序实现:
8 W* m: p6 ~3 U7 x6 g3 w" u[A]=xlsread('C:\MATLAB7\work\问题3数据5');
# p0 O# {+ P# A" w1 g8 i4 slen=length(A);& K+ A) u+ ]% }" `/ @! m( T
i=1:len; u0 C+ _" N* X! Y: o& i4 k
j=A(:,5);" m, ^2 A3 v' U. m
figure(1);: W# q" F) U7 o K5 v1 t. F
plot(i,j,'.-');% e: p; s7 n8 \5 @: U( H
hold on;8 t, |+ K, O- x6 Y* w9 K7 p
i=i';
6 g: \2 ^4 R% y* e c" i: Q9 `/ }p1=polyfit(i,j,2)
3 U. i) o6 e+ y( A' j; `( f& Qf1=polyval(p1,i);
& S4 E P. j4 i4 j5 Qplot(i,j,'.',i,f1,'-');
" `5 \+ s$ w8 b. Z- Whold on;
* g& @- E8 y$ ktitle('病人平均等待时间与病人到达时间的关系');1 g$ E7 |8 q5 m" P% { ~! B8 f! I
xlabel('病人门诊时间距离7 月12 日的时间差');ylabel('病人平均等待时间; Q$ w) p0 s2 b
');- Q" m1 G* P, v$ e. R4 p# Y) z
axis equal;
( ^% a9 z0 ]' q8 a$ |- R- v+ j" r8 Kfigure(2);
7 L4 i$ y- e4 L; v) p9 KA1=[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. O) W) o% M: ~( M) H2 `7 a; x+ T
32 34 35 36 37 38 39 40 41 42];
3 w5 P, ]5 g9 ^ s2 E, W6 wB1=[12 12 12 12 12 12.5 13 13 13 13 13.4 14 14 13 12.5 11 11 12 12 12 11
" Y+ h3 i6 Z3 }! J4 X5 p11.5 13.25 14 13 12 11 13 13 12 12.25 13 15.5 15 14 13];/ V1 H. e9 M0 I/ C$ ^
subplot(1,2,1);# X& {! `: }) N/ F0 m( z2 z
plot(A1,B1,'.-');- O* ~6 l V8 ^7 _. z
hold on;
$ B1 F+ |# v7 _; T' ^p2=polyfit(A1,B1,2)4 c5 C( h+ L! C; f c5 K$ k0 T
f2=polyval(p2,A1);9 e5 S" T4 Z* l. w% d$ U
plot(A1,B1,'.',A1,f2,'-');% p- w+ {& H& B/ \. S+ x& K
hold on;
1 ^7 r6 _" {1 @title('白内障病人情况');4 C( k" m9 S6 K
xlabel('白内障病人门诊时间距7月12日的时间差');ylabel('白内障病人平均) O4 f4 G+ e8 F2 ^4 A% I& p
等待时间');
/ z! ~' u" A- C$ m; Daxis equal;
) w4 ?' E) \4 i2 H. Q( L/ F3 gA2=[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
# v: h. F; ?/ Y4 ]31 32 33 34 35 36 37 ];! Q& T0 f1 k8 I
B2=[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
) x* j# H3 F; g" [12 12 13.25 14 14 13 12 12 12.4 13.33 13 13 13 12];
8 W L. D' h- H* x; r( C# I+ }subplot(1,2,2);
7 \$ h$ d" t! z' s1 R9 l- M% h2 uplot(A2,B2,'.-');
' I& N" m$ H- V A2 |! Jhold on;
0 h+ B5 ?) g2 x& G- Z6 y& dp3=polyfit(A2,B2,2)7 ^0 ^; S/ Y1 Q5 x( Z5 z# g3 ^3 @
f3=polyval(p3,A2);. Q: X* o8 z6 D, u: f6 h
plot(A2,B2,'.',A2,f3,'-');
: i# R6 x* z$ j* K% i: Uhold on;
0 F& p- ?- k; l( A. v2 B% ~title('白内障双眼病人情况');4 O. `* ~7 B: _; P1 ~8 P; J
xlabel('白内障双眼病人门诊时间距7月12日的时间差');ylabel('白内障双眼9 Z; N6 C' b2 m8 }- ~' J* l; ~
病人平均等待时间');
( y/ t. K' w. D8 l& A( yaxis equal; D) q- A4 r0 o5 P) n& Q
figure(3);
) V% O; c6 a+ `3 s( tA3=[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* j5 m$ f6 C( f5 t. c1 Y5 _
28 29 30 31 32 33 34 35 36 37 38 ];
& v3 _% K8 c2 a. x6 C+ 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
% ^6 \! J8 z2 Q: A11.67 12 12.25 12 11 11.75 13 14 14 13 12 12 12.33 13 13 13 12.5 12 12];
" ~3 U9 L/ A( B$ [( Asubplot(2,2,1);
& a* E* ~& u/ n& Eplot(A3,B3,'.-');
, y0 s# }0 u4 K F/ E' v, ~7 Bhold on;, u- I9 c+ X5 r& E/ i1 g# h
p4=polyfit(A3,B3,2)
9 X4 {& p; c( _9 @& D! ~f4=polyval(p4,A3);
% Z+ M/ u8 v O- M' ~- b# i0 Q, Nplot(A3,B3,'.',A3,f4,'-');
+ J' N8 o9 m! P3 j& V# ?9 qhold on;
' a0 k( p; G/ Jtitle('视网膜病人情况');) B- X, h; K, t" h" z
xlabel('门诊时间距7 月12 日的时间差');ylabel('视网膜病人平均等待时间
4 g1 r9 l/ p# t' V6 d');3 A( c/ [" x$ {! b) ~+ T! S1 t0 z
axis equal;
3 }% `8 ?9 x8 z0 v% ?% a. z1 A4 ^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 347 n( a, x. F: V# V5 ]! w( g
37];- ]7 i7 R: C0 i& ?
B4=[12 13 12 12 12 12 13.5 13 13.5 14 13 12 12 11 10.25 12 11 11 15 15
% S: z, C' z/ }- _2 {13 12 12 13 13 12];
4 P* ~" R9 X# S) S2 V$ C3 q2 ^subplot(2,2,2);+ q* m4 K4 {1 ?/ |6 j0 K
plot(A4,B4,'.-');
0 \/ K" y/ x8 {6 Ohold on;
7 p9 t9 M4 F' k( Sp5=polyfit(A4,B4,2)
9 }( A# u9 u8 n0 p) Rf5=polyval(p5,A4);+ [5 j7 t: q2 f1 h/ Q
plot(A4,B4,'.',A4,f5,'-');
" L' a# |1 r8 f+ I$ O( z% nhold on;. I4 O* `: {' @# M4 Q% P
title('青光眼病人情况');
R! @( x+ |; [1 Fxlabel('门诊时间距7 月12日的时间差');ylabel('视网膜病人平均等待时间0 f g! f, c% S/ x, i
');9 H7 N7 P& ]( H% @$ |" l5 ]
axis equal;
Q. K+ [8 l3 a7 BA4=[1 2 3 4 5 6 7 10 11 12 14 15 16 17 19 20 22 24 26 27 29 30 32 33 347 P1 F# l& I2 r' D0 q
37];
0 }) D/ @. P8 M, J9 b' [* x O9 [B4=[12 13 12 12 12 12 13.5 13 13.5 14 13 12 12 11 10.25 12 11 11 15 15- f) u" e% b1 }& P+ g
13 12 12 13 13 12];
8 g5 B7 Q* y* Qsubplot(2,2,2);' {2 O. x G+ t2 J' J. C
plot(A4,B4,'.-');5 H: J l/ P9 m
hold on;# \+ L$ x! Z6 D5 p3 Q
p5=polyfit(A4,B4,2)
. y$ I' L% `: u0 i' df5=polyval(p5,A4);. v$ S2 B @4 J1 A; Q1 L/ V
plot(A4,B4,'.',A4,f5,'-');) W; w) `" W1 z- | P
hold on;- k+ ?, ^9 g7 q
title('青光眼病人情况');1 m" O: n' O6 p2 v8 ]
xlabel('门诊时间距7 月12 日的时间差');ylabel('青光眼病人平均等待时间3 {. {( ~0 c& M# t) V9 M
');8 @/ Y0 S0 B0 y( @7 w' D, h3 ~
axis equal;& u1 p, F# V: V, J
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];
+ e3 Q- ?8 `4 b+ BB5=[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];: M/ y# I8 H1 y8 G3 f. d
subplot(2,2,3);2 N0 w* L" J3 D0 K- W
plot(A5,B5,'.-');! |' K- h$ J) d/ c
hold on;. E/ ]3 U: W. | H( U
p6=polyfit(A5,B5,2)
9 \4 P7 U4 m& ]5 L1 |" K- M4 cf6=polyval(p6,A5);
9 `9 ?4 e' e5 d$ ]6 g% K( Nplot(A5,B5,'.',A5,f6,'-');/ V6 |6 ]7 H* Q0 L! f' `1 v
hold on;1 k0 s! P: P3 X M, u- C, Y
title('外伤病人情况');, m% [- ~( I' b2 Y) [0 Z& j4 A: W
xlabel('外伤病人门诊时间距离7月12日的时间差');ylabel('外伤病人平均等; Y* c' e' }5 \ Z
待时间');3 P, g E" e/ P) q% P
axis equal;, o; d' r6 ~7 Z B$ `
二、层次分析法的一致性指标判定的matlab 实现:% ~6 s7 V% b1 ]# v) x' D. ~
function H=ff(A)
9 ]4 R- U5 ? p8 ^ N" Qt=eig(A)& d4 W% w- f, L$ s* n- ^ }* \
w=max(t)* Q" c4 I4 P* s6 t$ g+ B+ s
q=(w-4)/3/.9& b$ }% T9 N+ M, R6 A' @
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(
+ W4 O! O+ Z {# R/ b9)+A(10)+A(11)+A(12)+A(13)+A(14)+A(15)+A(16))
! J0 @. x8 H- m1 i7 A* ~(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)
! {/ X6 e9 K: i2 b- A/ m+A(10)+A(11)+A(12)+A(13)+A(14)+A(15)+A(16))3 K0 u. G! Z4 y3 ~0 u! V
(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)# l2 D- Q' ]( @: @- x
+A(10)+A(11)+A(12)+A(13)+A(14)+A(15)+A(16))/ W a1 B0 ~$ I
(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)
! W% H' u7 p& r8 ?. W+A(10)+A(11)+A(12)+A(13)+A(14)+A(15)+A(16))]2 R9 l {. g8 E
|
zan
|