QQ登录

只需要一步,快速开始

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

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

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

5

主题

4

听众

23

积分

升级  18.95%

该用户从未签到

群组华南理工大学

跳转到指定楼层
1#
发表于 2010-8-1 20:45 |只看该作者 |倒序浏览
|招呼Ta 关注Ta
一、使用曲线拟合模拟病人等待时间的模拟matlab程序实现:
9 a7 r9 K$ {8 i2 U[A]=xlsread('C:\MATLAB7\work\问题3数据5');
& H/ o  b9 R( N) A; |0 llen=length(A);3 t4 r3 v! v+ e
i=1:len;
* \$ o. P1 y  a6 \j=A(:,5);  Q6 ^% W% e9 j+ t9 f; h
figure(1);8 [- Z) T8 A2 L& I- e
plot(i,j,'.-');0 ?+ f0 n8 }3 ^. l
hold on;
2 e4 V9 U4 b7 W" E4 u% P7 g1 Oi=i';% C6 _: X0 k9 K- _0 m! @- l
p1=polyfit(i,j,2)
6 W) N! Z3 _! U8 H" ^f1=polyval(p1,i);7 K6 e3 `! @! N! _; r% \
plot(i,j,'.',i,f1,'-');
3 j! ^% N7 u8 f$ l; Nhold on;
: d$ i3 j& a( Z5 `. Y; {0 Btitle('病人平均等待时间与病人到达时间的关系');* o& P$ y/ k. l% D; W+ B! ~4 S( R
xlabel('病人门诊时间距离7 月12 日的时间差');ylabel('病人平均等待时间
) U7 C3 d1 _# ^: L7 w1 D');7 I. C) K: \5 |" o+ _) j  G8 ~
axis equal;
+ V4 `8 j1 G5 c; U6 ofigure(2);5 w% b) ~* r' @% z
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 319 C" h7 ]& z& F' s, Q' E
32 34 35 36 37 38 39 40 41 42];
! z9 r& z8 a1 iB1=[12 12 12 12 12 12.5 13 13 13 13 13.4 14 14 13 12.5 11 11 12 12 12 11
, `; X, o0 T- E) X11.5 13.25 14 13 12 11 13 13 12 12.25 13 15.5 15 14 13];
' a) R4 \9 i; W) F" p) v3 jsubplot(1,2,1);
3 d3 J: p  s" Iplot(A1,B1,'.-');
  i5 Y5 y$ B- Y, U7 P" q4 ]3 {3 k0 Whold on;, }) c8 y$ l9 y3 F" G  Q4 ~, v
p2=polyfit(A1,B1,2)3 `& ?1 Z" `; f, `) U
f2=polyval(p2,A1);3 g/ L4 i5 O# I) @  e
plot(A1,B1,'.',A1,f2,'-');
3 N& n, s' A' n5 nhold on;) q- w+ X6 j$ C% R: v3 a% F& e
title('白内障病人情况');
. n; ]* n: X' Y% k( cxlabel('白内障病人门诊时间距7月12日的时间差');ylabel('白内障病人平均
+ @0 [' \$ o( V) P$ \2 `0 m等待时间');' E3 v9 A/ d9 J( w" g# Z
axis equal;
& @4 g0 R( ~' ?6 s  C- Q" d( ?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
3 C, x* E1 z. n( i31 32 33 34 35 36 37 ];6 O8 y9 c2 p; A. g2 C& v. F% w
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
- p, h, z! M* I12 12 13.25 14 14 13 12 12 12.4 13.33 13 13 13 12];6 ?# H$ r& f1 c3 }! j4 d+ y, e, ?
subplot(1,2,2);+ H. R; ^! r8 w5 F( P
plot(A2,B2,'.-');( i" x' B. K! X+ p0 q
hold on;
6 j1 M/ k5 A; I3 v. xp3=polyfit(A2,B2,2)
) p5 j) X- D3 U& b; af3=polyval(p3,A2);; b' H" D% {4 E! D7 m  q
plot(A2,B2,'.',A2,f3,'-');! V. Z% D# V# Z) q7 ?2 |' Y8 P
hold on;
5 _3 Y1 E* V0 S1 J- |title('白内障双眼病人情况');0 l) q/ |' N5 e3 x: r; U2 F
xlabel('白内障双眼病人门诊时间距7月12日的时间差');ylabel('白内障双眼
7 i9 [( {9 |" K- }3 |9 f病人平均等待时间');+ `" a- M) k7 w
axis equal;
. n3 I4 z: |2 N- j4 A6 y% U* vfigure(3);
1 I% d* n2 U1 b5 A2 v, c- lA3=[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. j2 k. i8 R% `. q6 m, B' r3 [' x
28 29 30 31 32 33 34 35 36 37 38 ];
  C8 J* ?- Z; T  oB3=[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" F4 z  t% Q3 j' L- v+ L" Q
11.67 12 12.25 12 11 11.75 13 14 14 13 12 12 12.33 13 13 13 12.5 12 12];+ k( W( V1 X: I+ j
subplot(2,2,1);6 k3 O( A. R" o" j3 B
plot(A3,B3,'.-');# Q3 W% S+ c# C5 s* i
hold on;/ H5 R$ C! f5 e3 B3 ]
p4=polyfit(A3,B3,2)' Y% L; S6 d5 D
f4=polyval(p4,A3);
  J2 B1 g9 K2 `4 G4 Y+ `1 _9 Q7 Splot(A3,B3,'.',A3,f4,'-');
  B8 x7 P/ }/ I0 L7 dhold on;
! X6 ~0 M3 [0 etitle('视网膜病人情况');1 v6 I" i% J1 a2 J
xlabel('门诊时间距7 月12 日的时间差');ylabel('视网膜病人平均等待时间9 V* o4 z- B" J( p
');
( U1 k+ V: ~& T6 y; Q; Xaxis equal;3 }, h6 L1 v) U. H" Q1 \
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+ g9 P" I8 U( L) t' a
37];* U  I) T% W4 d( b, f4 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
/ B4 E9 B. ?! X+ w! Z6 G13 12 12 13 13 12];
( H" X1 p" |& O5 U* Z5 [subplot(2,2,2);
5 v% s5 r8 a' E5 H0 m% }6 y+ y3 `plot(A4,B4,'.-');$ _& y8 W; K- l* z
hold on;9 O1 d, |6 a: j
p5=polyfit(A4,B4,2)
- c4 `0 a; `, d! X6 L) Z. y& p  |f5=polyval(p5,A4);
9 A. @$ w. p8 W/ [) Bplot(A4,B4,'.',A4,f5,'-');
& @& ]4 D9 J& Vhold on;
) v% D" K0 v' I1 i2 L) [7 z% Utitle('青光眼病人情况');9 ~4 o& f" R# d
xlabel('门诊时间距7 月12日的时间差');ylabel('视网膜病人平均等待时间
- \! ~( |* e: r5 T( T6 N');3 G' f4 B- P  |% u" C0 G5 @. U- ~
axis equal;
; ~0 r# c9 `! G& z, rA4=[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
/ w$ W' d; d- ^+ _5 r% d; d+ h37];( h) V( Z- t: {- X6 Z$ A3 s
B4=[12 13 12 12 12 12 13.5 13 13.5 14 13 12 12 11 10.25 12 11 11 15 15
4 [2 W  g) p8 {13 12 12 13 13 12];& N  K& n  U+ F4 T8 b3 R2 {% u
subplot(2,2,2);" l0 S7 w. a6 ?+ W* r8 p( w
plot(A4,B4,'.-');/ E+ O3 G+ M& r
hold on;2 n% _% h. B  C$ T; f7 n) a
p5=polyfit(A4,B4,2)9 y0 A0 c0 `! b# |; {( I9 x
f5=polyval(p5,A4);/ b1 h0 j4 o3 |; Y/ T
plot(A4,B4,'.',A4,f5,'-');
1 `4 X! |) d; |3 Khold on;1 \( i- U; m) O
title('青光眼病人情况');
7 N& |: p; c) u; nxlabel('门诊时间距7 月12 日的时间差');ylabel('青光眼病人平均等待时间
3 N% L* G& i/ z# \1 S" w');" j6 m5 q0 s- d: J, ?% {7 I" W
axis equal;' m! K5 R( }6 x) _/ N5 |- @7 {
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];
3 J7 w/ j, b' i! ^5 x! H8 lB5=[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];5 P6 G$ ]6 x+ q6 v" s* z6 [
subplot(2,2,3);
/ \3 s) j, Z0 Nplot(A5,B5,'.-');
4 I' Q/ S1 m" S  L4 x  E! T' q% I4 U. zhold on;7 L% k% h5 ?9 a! l* P. }% h( f
p6=polyfit(A5,B5,2)" H8 s$ S( v" u6 {8 U  E6 |0 D) `
f6=polyval(p6,A5);* T  U. ~4 F; d( z( F
plot(A5,B5,'.',A5,f6,'-');& u  P, N$ q5 n. [9 M
hold on;
8 K% t: F; u7 c7 N( Rtitle('外伤病人情况');3 L: a+ R# `/ z. V* P1 D
xlabel('外伤病人门诊时间距离7月12日的时间差');ylabel('外伤病人平均等' i" M( Q; w# d9 Z5 N9 ]
待时间');" V9 H2 ~3 R& l2 ~, O* g1 V# q2 W
axis equal;, D4 E2 k" j3 z/ n# G
二、层次分析法的一致性指标判定的matlab 实现:
, f2 {  y4 o* `+ H, r0 h$ |function H=ff(A); C, ]% S- ]" H' f% ?
t=eig(A)& L- G$ h' e3 Q2 Q) M: X( _+ Y
w=max(t)
( n- c/ E, [- S4 d: fq=(w-4)/3/.92 C, L" e0 M  l" y, _& n
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(! m. q) m) a+ W; o. }  w9 \( U1 C$ e
9)+A(10)+A(11)+A(12)+A(13)+A(14)+A(15)+A(16))' C' ~# A: [+ Z6 j
(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)
3 S6 x; y. N* G- t% B7 K! H+A(10)+A(11)+A(12)+A(13)+A(14)+A(15)+A(16))
- ^: V' K- J# s+ 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)
/ ^# y( b; B3 d+A(10)+A(11)+A(12)+A(13)+A(14)+A(15)+A(16))) a7 F( R- r- p' e4 p, s0 ]8 Z1 b, n
(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)
+ l  a" V. t- U  P+A(10)+A(11)+A(12)+A(13)+A(14)+A(15)+A(16))]
7 f- H( r, L  X" x& b

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 微信微信
linmatsas 实名认证       

53

主题

13

听众

3592

积分

逍遥游

  • TA的每日心情
    奋斗
    2014-12-2 09:53
  • 签到天数: 54 天

    [LV.5]常住居民I

    自我介绍
    额。。。。世界上最讨厌的事情就是自我介绍。。。

    邮箱绑定达人 新人进步奖 发帖功臣 最具活力勋章

    群组Matlab讨论组

    群组数学建模

    群组小草的客厅

    群组2012数学一考研交流

    群组C 语言讨论组

    回复

    使用道具 举报

    3

    主题

    6

    听众

    264

    积分

    升级  82%

  • TA的每日心情
    开心
    2023-4-16 10:24
  • 签到天数: 6 天

    [LV.2]偶尔看看I

    自我介绍
    路漫漫其修远兮,吾将上下而求索!

    社区QQ达人 新人进步奖

    群组Matlab讨论组

    回复

    使用道具 举报

    linmatsas 实名认证       

    53

    主题

    13

    听众

    3592

    积分

    逍遥游

  • TA的每日心情
    奋斗
    2014-12-2 09:53
  • 签到天数: 54 天

    [LV.5]常住居民I

    自我介绍
    额。。。。世界上最讨厌的事情就是自我介绍。。。

    邮箱绑定达人 新人进步奖 发帖功臣 最具活力勋章

    群组Matlab讨论组

    群组数学建模

    群组小草的客厅

    群组2012数学一考研交流

    群组C 语言讨论组

    回复 迷途的羊 的帖子
    ' o! x  y/ n( O6 q) }
    ; Q' _+ J$ R; N4 E7 [  t% f是自己编的还用别人运行么……不过要是楼主自己编的……小弟就佩服死啦) Z- ]- J' ^: p# [8 s, \+ a
       
    回复

    使用道具 举报

    紫辰 实名认证       

    12

    主题

    16

    听众

    1304

    积分

    升级  30.4%

  • TA的每日心情
    擦汗
    2013-2-5 09:29
  • 签到天数: 35 天

    [LV.5]常住居民I

    自我介绍
    200 字节以内

    不支持自定义 Discuz! 代码

    群组学术交流A

    群组数学建模保研联盟

    群组Matlab讨论组

    群组湖南大学数学建模

    群组学术交流B

    回复

    使用道具 举报

    1

    主题

    2

    听众

    128

    积分

    升级  14%

  • TA的每日心情

    2014-9-9 21:30
  • 签到天数: 3 天

    [LV.2]偶尔看看I

    自我介绍
    平凡的我!

    平凡的快乐!
    回复

    使用道具 举报

    21

    主题

    7

    听众

    3435

    积分

    升级  47.83%

  • TA的每日心情

    2014-5-25 20:58
  • 签到天数: 20 天

    [LV.4]偶尔看看III

    新人进步奖 优秀斑竹奖

    群组Matlab讨论组

    群组小草的客厅

    群组数学趣味、游戏、IQ等

    群组C 语言讨论组

    群组我行我数

    [A]=xlsread('C:\MATLAB7\work\问题3数据5');
    * W0 ^( S0 t: F在不同的地方运行上面一句要改成本地数据保存的路径的
    回复

    使用道具 举报

    5

    主题

    4

    听众

    23

    积分

    升级  18.95%

    该用户从未签到

    群组华南理工大学

    回复

    使用道具 举报

    hhao        

    0

    主题

    2

    听众

    46

    积分

    升级  43.16%

    该用户从未签到

    新人进步奖

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

    使用道具 举报

    abc321        

    0

    主题

    2

    听众

    30

    积分

    升级  26.32%

    该用户从未签到

    新人进步奖

    回复

    使用道具 举报

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

    qq
    收缩
    • 电话咨询

    • 04714969085
    fastpost

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

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

    蒙公网安备 15010502000194号

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

    GMT+8, 2026-7-31 19:37 , Processed in 0.647164 second(s), 107 queries .

    回顶部