QQ登录

只需要一步,快速开始

 注册地址  找回密码
查看: 14105|回复: 35
打印 上一主题 下一主题

[代码资源] 请问MATLAB编程高手啊谁能将这个程序运行啊

[复制链接]
字体大小: 正常 放大

5

主题

4

听众

23

积分

升级  18.95%

该用户从未签到

群组华南理工大学

跳转到指定楼层
1#
发表于 2010-8-1 20:45 |只看该作者 |正序浏览
|招呼Ta 关注Ta
一、使用曲线拟合模拟病人等待时间的模拟matlab程序实现:
+ b* a& n- \* r& u1 O[A]=xlsread('C:\MATLAB7\work\问题3数据5');) D* V5 ?: v% `' p9 o  U: i  [  ?
len=length(A);
- V" J7 ^6 ?8 K+ Di=1:len;2 _8 ]  B! g" s; j3 n& b
j=A(:,5);# L  r1 Q4 M, L+ j* K) x
figure(1);
8 e3 ~* g$ ^0 t7 \% vplot(i,j,'.-');
1 U; R6 ]4 C' ?5 G' n( P2 xhold on;
. a& k5 d2 D+ \' si=i';  G* f- ?( `' r1 w! S
p1=polyfit(i,j,2)! m4 x; j+ a8 f: ^4 G) `
f1=polyval(p1,i);
( M/ G+ U6 S1 L. n8 l! Yplot(i,j,'.',i,f1,'-');
/ y2 ]; n7 }  l9 G' Chold on;: g8 u+ ]% A  r" T9 }
title('病人平均等待时间与病人到达时间的关系');
  `( H1 A! p6 r5 n! S. nxlabel('病人门诊时间距离7 月12 日的时间差');ylabel('病人平均等待时间1 a/ l2 _/ u0 E& v6 h4 ?! y# ]5 H
');
, N9 w6 P- F2 ^3 C& o8 y6 C! Faxis equal;2 z0 p* j( b- {
figure(2);1 `% L6 J8 e6 |" I) P5 i
A1=[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
& S. d. `2 R' Z2 T! h/ y/ P+ g32 34 35 36 37 38 39 40 41 42];- b9 _3 v7 V) Q6 u8 f; i1 `5 ?
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
, y  J5 ?* \1 F# Z11.5 13.25 14 13 12 11 13 13 12 12.25 13 15.5 15 14 13];2 U. f3 T' O  W; ]& G
subplot(1,2,1);
) Y* y% U' Y- }  o/ iplot(A1,B1,'.-');
0 s9 X& C: ]" Bhold on;
9 @& U+ e4 O" I" s/ U. a! u, Op2=polyfit(A1,B1,2)$ O0 V1 J" ]& G# ?; U( j
f2=polyval(p2,A1);
+ I9 B3 @8 v! H4 M0 [$ W$ iplot(A1,B1,'.',A1,f2,'-');# R' D" `4 o( I  A& J3 S9 _
hold on;2 I( {; o* e: X) |
title('白内障病人情况');
7 S. F1 ]" p' D5 Oxlabel('白内障病人门诊时间距7月12日的时间差');ylabel('白内障病人平均
) P% ]) ~7 C& B/ f5 _. x- o: T等待时间');7 N1 A! ~: }/ Y& R- g
axis equal;
: N" `9 g; P# m5 CA2=[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
# M* G& c# w0 T' U  u4 `( v6 r5 k31 32 33 34 35 36 37 ];4 i, C' o0 w/ J; I- K/ B
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
1 w! U8 _9 \& O12 12 13.25 14 14 13 12 12 12.4 13.33 13 13 13 12];2 D8 e; q* W; c  b0 p# h
subplot(1,2,2);
3 {  p, k$ _7 T0 z+ u( w& F/ |8 Rplot(A2,B2,'.-');
: J1 [7 `% o9 ]; hhold on;" T& t) @5 O6 ^+ P7 H7 {  ~
p3=polyfit(A2,B2,2)
, l2 B# E; }5 gf3=polyval(p3,A2);6 X9 {: _) s5 D/ ~/ T7 O7 a
plot(A2,B2,'.',A2,f3,'-');+ x% P2 J% e  B  |' Q# c6 G
hold on;
$ F! H1 l  M0 N8 ttitle('白内障双眼病人情况');
8 x1 P5 T/ e8 `) G4 Yxlabel('白内障双眼病人门诊时间距7月12日的时间差');ylabel('白内障双眼: E6 G* k2 C9 K; p
病人平均等待时间');0 E5 u4 R9 G! y
axis equal;% T) x) M* f  H' f  ]/ \1 Q7 G; l
figure(3);; H/ n3 |* W- b0 P
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 276 r) N9 x) S' n' d( [$ o0 I+ o4 c9 }
28 29 30 31 32 33 34 35 36 37 38 ];
0 G9 a; L* t7 t! D  A4 `B3=[12.33 12.5 12.3 12 11.67 12.33 12 13 13 14.2 14.25 14 13 12 11 11 10.753 U& V2 G1 C# G6 e( I+ Q2 O
11.67 12 12.25 12 11 11.75 13 14 14 13 12 12 12.33 13 13 13 12.5 12 12];
4 ~8 a: I2 M# U  q/ E4 _5 C4 F) F# C( gsubplot(2,2,1);
" f1 |$ o/ N; a& |& H) [: ]7 s/ b0 Nplot(A3,B3,'.-');
9 R$ Z4 J% |: ]: ]1 d$ N* L) F2 Y2 N* @! Khold on;7 y# u8 J* ^- @, R2 z# q3 h$ n2 A
p4=polyfit(A3,B3,2)
6 v4 M! }* R$ Y# A2 wf4=polyval(p4,A3);
6 j* q3 U8 X! S9 s9 Yplot(A3,B3,'.',A3,f4,'-');4 K/ }4 ~9 q( Z( o7 x5 k
hold on;3 F. g* g3 N( u' F7 w/ m
title('视网膜病人情况');" C+ B' B, S3 t1 Y5 ^( e# y
xlabel('门诊时间距7 月12 日的时间差');ylabel('视网膜病人平均等待时间
5 ]2 X4 ]- F1 |3 Y+ ^');$ [1 \! H! l" ~2 G
axis equal;: H+ T- P$ F) _5 e
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& c# B1 A4 ]$ s+ ?6 @! a8 @
37];( U9 C1 i: v& o1 R
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  R3 G' r) Q8 h
13 12 12 13 13 12];
9 ]1 G# P0 }1 h% j; z. ^/ psubplot(2,2,2);
) }" b0 A% d4 y( eplot(A4,B4,'.-');0 ~. K+ S/ \/ E( P" L& [: Z
hold on;
' Q* B+ l# ^& E  f3 E! ?; Sp5=polyfit(A4,B4,2)/ n  C+ C) d: r, L
f5=polyval(p5,A4);3 }% R9 o6 r* C& B+ b( S0 P0 }
plot(A4,B4,'.',A4,f5,'-');
  \# W' e8 x% B4 ~+ mhold on;4 c# |( b( x) J7 Q* G- \7 y
title('青光眼病人情况');
0 m' w& Y4 P7 _8 mxlabel('门诊时间距7 月12日的时间差');ylabel('视网膜病人平均等待时间6 v) J7 g0 Y1 k
');  q5 W2 Y% [/ z7 p9 j
axis equal;6 F9 K9 L0 @4 _% C) j) r
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 341 h% Z2 n( g: N
37];
' y0 W& W) K4 i/ Z8 c* U  qB4=[12 13 12 12 12 12 13.5 13 13.5 14 13 12 12 11 10.25 12 11 11 15 15
5 h; G  `' H- i3 [13 12 12 13 13 12];
% [% N" B+ O) {( B8 w' Bsubplot(2,2,2);. E6 H( b! w* i* ^: y
plot(A4,B4,'.-');, j% T+ G4 n) i* y& j& ?5 E. X
hold on;  _- {$ E7 m1 |1 O. O8 A
p5=polyfit(A4,B4,2)
$ ~3 \% `5 ~8 r  Cf5=polyval(p5,A4);1 W/ h: {( v: M. H/ C
plot(A4,B4,'.',A4,f5,'-');
' o4 X% _- v1 P0 Rhold on;/ F' K) ^  G5 X9 C6 y  Q
title('青光眼病人情况');
2 I9 }' U. L7 u+ |. ]xlabel('门诊时间距7 月12 日的时间差');ylabel('青光眼病人平均等待时间. p7 d3 N; x0 O" X* O' k/ L
');
# z: z" ?& w; {axis equal;
1 z/ ?) f# \- B$ u8 qA5=[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];
" _8 ~9 V' g7 v4 IB5=[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];4 c" a- z# B# a2 e5 B7 _
subplot(2,2,3);
$ J$ Y4 ]" \6 F# c3 r9 Z- H- Iplot(A5,B5,'.-');! I2 {+ E  }/ v$ }, O: ~5 F
hold on;
0 F- H9 o, \# }$ a) ^, J. Pp6=polyfit(A5,B5,2)
7 w4 \( M  e: S3 }) R6 a7 of6=polyval(p6,A5);4 O4 S# z, B# u; c9 B; [
plot(A5,B5,'.',A5,f6,'-');* h5 S1 g3 R+ w, g$ \
hold on;
7 B( c7 D$ m6 Q$ H# I3 Ktitle('外伤病人情况');
; n8 X+ `# l7 \( Exlabel('外伤病人门诊时间距离7月12日的时间差');ylabel('外伤病人平均等
4 n2 }5 F' r' X8 V; t2 }* k  l待时间');/ O- `; ]7 ]' d0 [
axis equal;6 R/ f; r* d) Z2 }
二、层次分析法的一致性指标判定的matlab 实现:8 ]3 D  ~2 [' L4 O$ n, F- M
function H=ff(A)! S3 h% f4 @" v, |4 Y; G
t=eig(A)
" q6 ~3 M" B7 d5 p* \) V4 _w=max(t)
& O' U  X/ B, }8 Wq=(w-4)/3/.9& N! E. |# C( G: f8 h
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(
7 C; l" V: X: y( E) m4 S+ v9)+A(10)+A(11)+A(12)+A(13)+A(14)+A(15)+A(16))% z: t( f( ^6 j) a: n
(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)
0 M  H4 z3 j3 T+A(10)+A(11)+A(12)+A(13)+A(14)+A(15)+A(16))
# F( d+ z, s* u1 P' t& E(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 u6 {% [4 v' \3 I: \8 n1 I) b+A(10)+A(11)+A(12)+A(13)+A(14)+A(15)+A(16))
( D, L/ {" h! }* z3 `. ^9 ^) ](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)
# ^1 ?. H& m* K0 G$ Q- g) v+A(10)+A(11)+A(12)+A(13)+A(14)+A(15)+A(16))]
$ Y9 s- X: N' {0 c

B2009.doc

688.5 KB, 下载次数: 1, 下载积分: 体力 -2 点

09年b题

2009B题0315(1).pdf

352.05 KB, 下载次数: 3, 下载积分: 体力 -2 点

问题三数据五.xls

68.5 KB, 下载次数: 3, 下载积分: 体力 -2 点

zan
转播转播0 分享淘帖0 分享分享0 收藏收藏0 支持支持0 反对反对0 微信微信

22

主题

5

听众

854

积分

  • TA的每日心情

    2013-6-5 16:56
  • 签到天数: 202 天

    [LV.7]常住居民III

    群组学术交流A

    其实代码质量不高。。。。明显可视化设置没到位。然后后面的AHP写的也很乱还有bug
    回复

    使用道具 举报

    ilava        

    0

    主题

    2

    听众

    74

    积分

    升级  72.63%

    该用户从未签到

    新人进步奖

    回复

    使用道具 举报

    qiuhanty        

    0

    主题

    2

    听众

    32

    积分

    升级  28.42%

    该用户从未签到

    新人进步奖

    回复

    使用道具 举报

    hhao        

    0

    主题

    2

    听众

    46

    积分

    升级  43.16%

    该用户从未签到

    新人进步奖

    回复

    使用道具 举报

    Hyacinth        

    0

    主题

    2

    听众

    73

    积分

    升级  71.58%

    该用户从未签到

    新人进步奖

    回复

    使用道具 举报

    liguoli        

    0

    主题

    2

    听众

    71

    积分

    升级  69.47%

    该用户从未签到

    新人进步奖

    楼主的帖子实在是写得太好了。可是我立刻想到,这么好的帖子,倘若别人看不到,那么不是浪费楼主的心血吗?经过痛苦的思想斗争,我终于下定决心,牺牲小我,奉献大我。我要拿出这帖子奉献给世人赏阅,我要把这个帖子一直往上顶,往上顶!顶到所有人都看到为止!  
    回复

    使用道具 举报

    hab0458        

    0

    主题

    2

    听众

    45

    积分

    升级  42.11%

    该用户从未签到

    网络挑战赛参赛者

    新人进步奖

    回复

    使用道具 举报

    0

    主题

    2

    听众

    72

    积分

    升级  70.53%

    该用户从未签到

    新人进步奖

    回复

    使用道具 举报

    casper        

    0

    主题

    2

    听众

    72

    积分

    升级  70.53%

    该用户从未签到

    新人进步奖

    我基本上是采用看英语文章的办法,先泛读,再精读,再一句一句看,最后再提纲挈领,总算是明白一点了,当然,也可能还是领悟错了。最后要说的一句话是:楼主,你很牛叉,希望你不是真的有病。   
    回复

    使用道具 举报

    您需要登录后才可以回帖 登录 | 注册地址

    qq
    收缩
    • 电话咨询

    • 04714969085
    fastpost

    关于我们| 联系我们| 诚征英才| 对外合作| 产品服务| QQ

    手机版|Archiver| |繁體中文 手机客户端  

    蒙公网安备 15010502000194号

    Powered by Discuz! X2.5   © 2001-2013 数学建模网-数学中国 ( 蒙ICP备14002410号-3 蒙BBS备-0002号 )     论坛法律顾问:王兆丰

    GMT+8, 2026-4-12 22:38 , Processed in 0.545502 second(s), 109 queries .

    回顶部