QQ登录

只需要一步,快速开始

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

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

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

5

主题

4

听众

23

积分

升级  18.95%

该用户从未签到

群组华南理工大学

跳转到指定楼层
1#
发表于 2010-8-1 20:45 |只看该作者 |倒序浏览
|招呼Ta 关注Ta
一、使用曲线拟合模拟病人等待时间的模拟matlab程序实现:
$ b4 R& j2 Q- b- q[A]=xlsread('C:\MATLAB7\work\问题3数据5');7 q' s! }! K4 p' E8 Y: V3 q; q# u
len=length(A);) K" U) o3 L! l. v
i=1:len;
- |; D0 W8 x2 Y; zj=A(:,5);
6 A3 c  ?: p) y! P1 g" ffigure(1);5 l9 R. K3 b8 X
plot(i,j,'.-');
# S8 Q$ y- \" \, z6 ^) jhold on;
2 L" Q. M* ?1 ]% m0 h% o' o& G; ni=i';: I' Y+ W$ O/ l! P! j
p1=polyfit(i,j,2)2 d$ B" m' T8 B( G" Y
f1=polyval(p1,i);( a! W1 ^5 b5 m3 E  }9 S: b
plot(i,j,'.',i,f1,'-');, Y7 f5 S6 \3 C; L) N$ F
hold on;
( N; k! {* U' R# n/ R' A7 otitle('病人平均等待时间与病人到达时间的关系');/ w/ l& y  ^. _) Z
xlabel('病人门诊时间距离7 月12 日的时间差');ylabel('病人平均等待时间0 ^! r8 T; [, [8 J" C+ \/ A
');3 `/ X$ p% J4 J' z9 S& r7 l
axis equal;
3 r1 ~6 [- a, y. L4 nfigure(2);  v  ]+ T* y! ?* J; y$ S  B; c5 N- p' K
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
8 p% {9 a  u6 h( ?8 U$ C32 34 35 36 37 38 39 40 41 42];
" ?8 c( _: M3 E, V  GB1=[12 12 12 12 12 12.5 13 13 13 13 13.4 14 14 13 12.5 11 11 12 12 12 11( r, Q+ J0 E6 y
11.5 13.25 14 13 12 11 13 13 12 12.25 13 15.5 15 14 13];: M, f3 ~# @' O% |
subplot(1,2,1);" J7 Y) R! S4 o( y6 P$ \
plot(A1,B1,'.-');# R( Q! o! P. e5 b
hold on;
' N: M4 _5 _9 B. mp2=polyfit(A1,B1,2)0 E. R! B1 e; K* \9 z9 ]& K; K
f2=polyval(p2,A1);% z' U8 _: @- q$ S
plot(A1,B1,'.',A1,f2,'-');
5 |& L% p3 R" T* ]% j7 d3 {0 Rhold on;
6 _. G# F8 P2 f0 j7 c+ rtitle('白内障病人情况');5 A2 w) L4 e7 V9 W
xlabel('白内障病人门诊时间距7月12日的时间差');ylabel('白内障病人平均. S0 r6 Y# W! @3 N8 Q
等待时间');( b; g0 B0 W2 B# V
axis equal;
9 A. L( ~" ?2 h7 WA2=[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
. s# w6 c/ c) S+ E* u6 h% y0 ?4 Y0 }31 32 33 34 35 36 37 ];
! I6 {; o. M; y$ W# L/ r6 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- n. B. w! x' |  @" d0 I  J
12 12 13.25 14 14 13 12 12 12.4 13.33 13 13 13 12];
  @9 I$ f6 M1 B5 H5 Xsubplot(1,2,2);
5 @, J) R9 _+ {! S+ u9 B* bplot(A2,B2,'.-');
: i9 s; e4 @+ M' ]# r0 Dhold on;
( {. }: V/ v4 }, V4 K' hp3=polyfit(A2,B2,2)* q5 I5 g$ s7 ~( M* S: Q
f3=polyval(p3,A2);+ b% C9 M1 _' `. K& K
plot(A2,B2,'.',A2,f3,'-');# A' X7 \  |1 w8 W9 E. O0 G( z
hold on;
1 J# r) h1 w; x7 t/ l: [  wtitle('白内障双眼病人情况');
; N( E3 a; z/ |% T0 J& {xlabel('白内障双眼病人门诊时间距7月12日的时间差');ylabel('白内障双眼
- ^4 [# ?2 f/ b7 u0 M病人平均等待时间');
3 |! p; @/ Z9 @! B9 Yaxis equal;
! I9 n0 X+ w5 _figure(3);9 O* ]7 y) E; }: t
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
! `" ~3 C; Z/ K+ p( G+ F1 r8 q2 a* J5 v28 29 30 31 32 33 34 35 36 37 38 ];
* Z0 X6 ^; Y( v; RB3=[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 H, s$ ?( v2 B5 d+ T- {11.67 12 12.25 12 11 11.75 13 14 14 13 12 12 12.33 13 13 13 12.5 12 12];
8 {0 U7 N! q$ i6 p' f8 H, Ssubplot(2,2,1);
) I: s* r+ ~  r- ?) jplot(A3,B3,'.-');
/ @0 J/ l* k. c1 X) X- N- rhold on;! J( W+ [3 l- D1 j1 h! I
p4=polyfit(A3,B3,2)
/ X5 `: l5 I9 of4=polyval(p4,A3);
* ]$ p5 y+ m; N3 W# O" t7 }plot(A3,B3,'.',A3,f4,'-');
4 D/ P$ y' F$ X9 thold on;
: a% c& m- U& K! O7 A, stitle('视网膜病人情况');
5 s" k3 l* c3 }' Q% b- |xlabel('门诊时间距7 月12 日的时间差');ylabel('视网膜病人平均等待时间
* S3 |+ T# i& h' E, D');
- Q3 l( z" T# r3 }axis equal;! j; a1 |* f$ U* S" n
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$ F$ V4 A! U* R2 P! I  E
37];( b* P$ ~+ L( _" v
B4=[12 13 12 12 12 12 13.5 13 13.5 14 13 12 12 11 10.25 12 11 11 15 15$ }5 v, ~; M" L
13 12 12 13 13 12];: N8 i# [2 A* x3 g# k4 F3 A. @
subplot(2,2,2);
2 }0 N2 P* n6 q' f; Oplot(A4,B4,'.-');5 Y, l5 X- S) r4 V
hold on;
1 F- I; T3 K: g2 o, k7 m. d$ @! {* ap5=polyfit(A4,B4,2)5 H  M4 P  `: @; a! K& c
f5=polyval(p5,A4);! C. z0 V/ s- }. P* `: f, h5 _
plot(A4,B4,'.',A4,f5,'-');
- @! g. \, T4 ^hold on;: w3 Q- {% Z2 V/ O
title('青光眼病人情况');$ [; Z% r) d! n) g9 T2 S
xlabel('门诊时间距7 月12日的时间差');ylabel('视网膜病人平均等待时间# I3 S) l; r+ Y- z. C2 A% w- F7 A; U
');: {( j% v2 d# g3 V" P& _
axis equal;
9 U5 f& j9 U: }( ~# X  J+ b! yA4=[1 2 3 4 5 6 7 10 11 12 14 15 16 17 19 20 22 24 26 27 29 30 32 33 344 u! U+ n$ `) }. l! ~$ O1 h
37];
8 Y- `# i3 Z5 S/ b  m8 X( n6 `1 iB4=[12 13 12 12 12 12 13.5 13 13.5 14 13 12 12 11 10.25 12 11 11 15 15
/ Z9 {3 P6 S9 L0 ?" U13 12 12 13 13 12];
1 g4 z% J" d- Q. {" l0 X2 osubplot(2,2,2);! ~6 t$ ~* {. G0 G5 U; o0 p+ F6 a2 b
plot(A4,B4,'.-');& h& C5 G  i/ h
hold on;7 w) A3 ]3 a( W3 v0 ^' Z! F) d& H
p5=polyfit(A4,B4,2)
3 ~  ]9 ]3 W' W$ T4 Ef5=polyval(p5,A4);) K0 a! i% u6 ^  U) x. I" V$ R5 @4 M
plot(A4,B4,'.',A4,f5,'-');
: _1 y1 N2 o7 ~* ~hold on;3 a) ^  v- z8 g1 J  F- D* A
title('青光眼病人情况');' L$ i4 m, |% Y! [! i( |
xlabel('门诊时间距7 月12 日的时间差');ylabel('青光眼病人平均等待时间1 Y8 g% C" I2 S. a' W5 k
');
9 p% `9 C' a  D4 E9 R( saxis equal;
: I4 Z& i& Q: X( U! k1 E6 [1 W, EA5=[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];. b7 i& l' T. w
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];
% T2 J# ]( e. ~subplot(2,2,3);
% L- g! T$ S. O& r, Lplot(A5,B5,'.-');* S7 U/ |: C- v$ Y, f/ q; E
hold on;" i3 A, V, D3 S3 \
p6=polyfit(A5,B5,2)
! ?' C6 n- @4 g- Z' If6=polyval(p6,A5);8 ]8 O% M0 G/ \  ~
plot(A5,B5,'.',A5,f6,'-');
! V6 {# u( f, F8 K& E3 _! Fhold on;
" m6 X1 B1 V7 t' `( A) S& e! Ztitle('外伤病人情况');
! H7 c  L" k2 N3 f% n/ P! C9 Axlabel('外伤病人门诊时间距离7月12日的时间差');ylabel('外伤病人平均等
! k6 K: F) J+ }1 ]待时间');
  f+ V7 {$ r) p; Y5 d  |( D; Baxis equal;0 I" G/ D  _5 E. J: z9 }
二、层次分析法的一致性指标判定的matlab 实现:
6 Q9 ~4 G( ^% ~( Q/ w" u0 y, qfunction H=ff(A)
% ]. b: M: f' m4 M# I3 w# O6 Y# }t=eig(A)
" b, r$ ]) D7 K3 i8 yw=max(t), O% v, Z& }) |3 ]
q=(w-4)/3/.93 d/ z& d6 t( S8 z8 w/ 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(
; S4 O) Z9 Q0 a$ y9)+A(10)+A(11)+A(12)+A(13)+A(14)+A(15)+A(16))
" Z$ p1 n4 o: |2 x$ k(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)
9 [/ H! S- T2 Y4 J8 h+A(10)+A(11)+A(12)+A(13)+A(14)+A(15)+A(16))
7 U  F) W8 s0 A(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)
! J2 n6 t# ]) q0 X3 |- a: h7 N, i/ o0 c+A(10)+A(11)+A(12)+A(13)+A(14)+A(15)+A(16))8 {! `( s) N  d0 m& r
(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)8 ^# c" g1 R# W2 E( m: X" M1 Q
+A(10)+A(11)+A(12)+A(13)+A(14)+A(15)+A(16))]2 e! C6 d- E: Y" O

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

    回复 迷途的羊 的帖子, e5 b! @& ~2 P, T3 I

    , n; s! K4 i( X. L8 A5 B4 I$ J是自己编的还用别人运行么……不过要是楼主自己编的……小弟就佩服死啦
    ) S# `' M4 T7 a- a7 z! ^9 e3 ~   
    回复

    使用道具 举报

    紫辰 实名认证       

    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');4 t+ H' |/ r) S. X; T
    在不同的地方运行上面一句要改成本地数据保存的路径的
    回复

    使用道具 举报

    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 23:45 , Processed in 0.597747 second(s), 107 queries .

    回顶部