QQ登录

只需要一步,快速开始

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

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

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

5

主题

4

听众

23

积分

升级  18.95%

该用户从未签到

群组华南理工大学

跳转到指定楼层
1#
发表于 2010-8-1 20:45 |只看该作者 |倒序浏览
|招呼Ta 关注Ta
一、使用曲线拟合模拟病人等待时间的模拟matlab程序实现:
, {: @5 i6 D, [- t[A]=xlsread('C:\MATLAB7\work\问题3数据5');# i4 g( i( N; E6 `3 Q
len=length(A);7 ~, h9 L2 Y, G; q& s8 N7 t* ^
i=1:len;
1 m6 w# J" l: Q/ Ej=A(:,5);
: j5 E# d0 P8 A* h+ Z6 \figure(1);
$ C8 i; B+ K& r  X; \  gplot(i,j,'.-');8 B7 q$ h; o, s* M
hold on;* K/ K( u  q  p+ f  H6 \+ Y) J1 u
i=i';, B  A0 H4 T0 Q& y# O
p1=polyfit(i,j,2)
% y$ D6 |( W! I! gf1=polyval(p1,i);
$ y& ~- R1 Z6 t4 Y/ n2 D& Gplot(i,j,'.',i,f1,'-');
% O1 j7 R: ~& l. z0 [hold on;
! Z2 D# }8 H) o' Q. ?+ ptitle('病人平均等待时间与病人到达时间的关系');1 W+ S1 J. P9 p
xlabel('病人门诊时间距离7 月12 日的时间差');ylabel('病人平均等待时间$ [! J0 _- u5 `& H5 p
');
7 l1 q! ^6 p5 Z, E, `. Eaxis equal;9 Q1 w% P. F/ H2 ?; `
figure(2);
: I. j8 j( J3 M1 BA1=[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
/ w0 n6 Z* t" m( ~& A32 34 35 36 37 38 39 40 41 42];0 M. Y" @" H$ P
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, e. B' r* p! }. Q1 p9 X6 R
11.5 13.25 14 13 12 11 13 13 12 12.25 13 15.5 15 14 13];' t, E- l' b# h% @- ~) W
subplot(1,2,1);
6 m* f; K/ k8 Y2 W$ [/ T9 Bplot(A1,B1,'.-');( T5 n& h% s$ t& ]
hold on;6 v9 x2 b. ?# K& c2 i
p2=polyfit(A1,B1,2)
) g: S5 {. `/ ]f2=polyval(p2,A1);
7 L( `7 t- q, U' iplot(A1,B1,'.',A1,f2,'-');
/ A% x0 l: g; b+ Jhold on;( ^& f+ q( a- l, ^: _* e& A
title('白内障病人情况');- F" r" d% X/ X
xlabel('白内障病人门诊时间距7月12日的时间差');ylabel('白内障病人平均5 c) {$ y( x: D( ~
等待时间');
. W+ y# k4 i, z1 J8 Maxis equal;6 [+ Y. F# i0 ^) H  z8 q* ~
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
) z. e* [% t* i8 X5 n! j( y31 32 33 34 35 36 37 ];
, U4 c6 x' j0 }" aB2=[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, c2 S" v, I# v  `! S. Q9 }; S3 p" g
12 12 13.25 14 14 13 12 12 12.4 13.33 13 13 13 12];
! I; D- Q+ m! A$ ssubplot(1,2,2);8 P5 V. k  X  y% y1 b* Z
plot(A2,B2,'.-');. T) ]6 Q" Z6 [; C# X
hold on;& H: P. Y: S% L% f/ m+ l
p3=polyfit(A2,B2,2)
; g2 r' n2 K$ o7 z! w& lf3=polyval(p3,A2);+ H% w" t1 |; ~$ M5 P9 ~  w+ b
plot(A2,B2,'.',A2,f3,'-');# E! l7 ?- A& u- p" Q8 M4 j
hold on;5 C5 n& K- G7 x2 ~) C
title('白内障双眼病人情况');
1 l3 V6 G" p( w- m, nxlabel('白内障双眼病人门诊时间距7月12日的时间差');ylabel('白内障双眼1 |8 k5 Y5 U5 D+ {
病人平均等待时间');
' z; B$ Y7 m! T  g# H; M* Q6 p$ Z( Baxis equal;3 P; Y: L" R( o) I
figure(3);, `) {3 Z8 H# X$ ]
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 277 [* S+ D$ R1 n
28 29 30 31 32 33 34 35 36 37 38 ];
$ p' s4 a6 v/ C  jB3=[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, K8 W/ @7 ]* J. `
11.67 12 12.25 12 11 11.75 13 14 14 13 12 12 12.33 13 13 13 12.5 12 12];9 H. M3 C! `) c: D. q+ U6 R- D1 n
subplot(2,2,1);  N$ N; z' m6 t- J$ R: P1 [' n$ J
plot(A3,B3,'.-');
- F  O) p; Q# ]8 c& }hold on;
; \( J! @, k0 y' j9 ~! vp4=polyfit(A3,B3,2)
' ^, ]2 ~- r! M; I0 bf4=polyval(p4,A3);
) Y$ _. V! O3 f! c! [4 Zplot(A3,B3,'.',A3,f4,'-');0 Q0 ^6 P" I( Z8 v% f
hold on;! Q5 I3 `2 Z  Q' k3 n" n
title('视网膜病人情况');
. e: u' w9 K3 z7 ^xlabel('门诊时间距7 月12 日的时间差');ylabel('视网膜病人平均等待时间
' G& j- U- P2 n% s2 Q' t');& w: Y7 s% K5 J) X5 e
axis equal;
' Y3 g+ i1 j" a9 q! Y/ S* tA4=[1 2 3 4 5 6 7 10 11 12 14 15 16 17 19 20 22 24 26 27 29 30 32 33 343 i5 M; n1 L/ o9 r$ x# w7 m
37];
" m) F3 y1 `# V# ]* I4 W! OB4=[12 13 12 12 12 12 13.5 13 13.5 14 13 12 12 11 10.25 12 11 11 15 15
  [+ R$ _3 n1 \+ S( Y13 12 12 13 13 12];
7 _6 ~1 d$ s9 b) P2 ssubplot(2,2,2);) r. n& S0 A; v/ P3 H9 E5 D
plot(A4,B4,'.-');5 }& T& {0 R  _! e3 k% P/ d
hold on;
; o! d3 o2 M  c1 P; ~# G3 Bp5=polyfit(A4,B4,2)
8 Z; @) W/ {9 ]6 k1 \5 qf5=polyval(p5,A4);& {4 e& t6 B" b1 y
plot(A4,B4,'.',A4,f5,'-');. Z* ]% w! y; A, q( ]) c8 n  |0 d
hold on;
+ L+ C) H) {4 h+ ?' \title('青光眼病人情况');# C8 O4 m! _! Z  N4 G
xlabel('门诊时间距7 月12日的时间差');ylabel('视网膜病人平均等待时间- N. Y6 n2 x' M0 Q; b8 o
');" c2 W, T/ B4 W1 A0 T
axis equal;
; D7 H5 [4 Q3 p' a# VA4=[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
) e( i6 F  m0 K) X1 D9 p; q% U5 [1 N' ~37];/ W! C7 o! X1 `! l! k6 L
B4=[12 13 12 12 12 12 13.5 13 13.5 14 13 12 12 11 10.25 12 11 11 15 15
1 J+ W; Y; {/ R1 k- j! C13 12 12 13 13 12];2 o. c: M$ V* X+ ^& X
subplot(2,2,2);& i9 p2 ~# N5 N* Y, j8 \
plot(A4,B4,'.-');
8 ]+ W4 O, S! }! g9 }: ?hold on;
" W. i! u1 Z! {+ {p5=polyfit(A4,B4,2)3 Q4 [) n: o* Z& d  \
f5=polyval(p5,A4);
5 W0 @# e# J. z9 k/ {$ Lplot(A4,B4,'.',A4,f5,'-');
# O1 ^: R% \9 hhold on;
1 s& Q+ F4 K1 \+ Q2 ^3 d% Etitle('青光眼病人情况');
: G/ k1 d4 [4 A; N& Kxlabel('门诊时间距7 月12 日的时间差');ylabel('青光眼病人平均等待时间
5 n. i0 t" V$ W7 _; N. [2 ~0 J! q');
4 ~/ p  f, q# N) W8 b/ f: p; _axis equal;+ w* Y/ y' o( \4 p. _$ H8 L0 s
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];5 Q, j  Y8 U9 _
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];
+ V9 B% Q+ A3 {0 L6 _: ~6 @" d" ~subplot(2,2,3);- h) T: G, R, N+ o: ]0 W& ^
plot(A5,B5,'.-');4 E9 Q# k1 a. r& B. S# r
hold on;3 n0 w) v+ u, t7 C  Y1 Y# U0 F
p6=polyfit(A5,B5,2)# H4 a8 m4 w- k* w& W( V
f6=polyval(p6,A5);
( i8 i# q$ ~5 x: r4 s3 t: {plot(A5,B5,'.',A5,f6,'-');  j& U" R# O1 h5 r7 E8 R5 _
hold on;9 M8 b8 s9 n3 f$ p9 ]" g# ]3 r4 T2 u
title('外伤病人情况');
( d! }5 F8 `; a) k  t& T! T4 N; Axlabel('外伤病人门诊时间距离7月12日的时间差');ylabel('外伤病人平均等
4 n8 Q, ?4 H! }: j7 ]待时间');
* s; e0 Y. N. `! Paxis equal;- O/ X) {) C; Q5 L6 v; p
二、层次分析法的一致性指标判定的matlab 实现:
% w2 Z! N% X. ]! V2 @% Gfunction H=ff(A)
6 v( }3 V5 _! k! n! Yt=eig(A)
/ O6 [" K: V, O& i/ r- [' c+ Vw=max(t)
4 g4 O4 `$ ?: q+ w& pq=(w-4)/3/.9" i- A) [0 b7 R4 W$ S  ^3 y
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(
- ]/ q- J( }- p0 c1 U' d- y4 f! N9)+A(10)+A(11)+A(12)+A(13)+A(14)+A(15)+A(16))
: A3 Q$ }+ n  o( d# 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)' r' H* V5 @& {+ i
+A(10)+A(11)+A(12)+A(13)+A(14)+A(15)+A(16)): i& K! v. n; |; c
(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)
& d2 m1 `. H! L7 t- ?% q' }, x+A(10)+A(11)+A(12)+A(13)+A(14)+A(15)+A(16))+ N9 Y$ X( q) _0 s: q/ q
(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)
  T( y5 C. B4 l( e  j  h! p7 g+A(10)+A(11)+A(12)+A(13)+A(14)+A(15)+A(16))]4 u# d; ~" ?/ F. H7 s/ K; n

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 语言讨论组

    回复 迷途的羊 的帖子& o3 P9 n8 ^% c+ m  G2 L# x% \

    ) V1 `0 K5 B2 s* }4 {/ @+ F是自己编的还用别人运行么……不过要是楼主自己编的……小弟就佩服死啦
    6 M8 m" _% g$ _7 r" P0 o+ p8 `   
    回复

    使用道具 举报

    紫辰 实名认证       

    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');
      z' \& w$ i  T. I" G& L* b在不同的地方运行上面一句要改成本地数据保存的路径的
    回复

    使用道具 举报

    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-11 15:14 , Processed in 0.494475 second(s), 107 queries .

    回顶部