QQ登录

只需要一步,快速开始

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

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

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

5

主题

4

听众

23

积分

升级  18.95%

该用户从未签到

群组华南理工大学

跳转到指定楼层
1#
发表于 2010-8-1 20:45 |只看该作者 |倒序浏览
|招呼Ta 关注Ta
一、使用曲线拟合模拟病人等待时间的模拟matlab程序实现:
# g. \- K0 {" y[A]=xlsread('C:\MATLAB7\work\问题3数据5');3 p% Y4 ]. I. }
len=length(A);) i9 l- @: q  Y' M. G
i=1:len;
" F7 X3 F; O3 h/ n1 zj=A(:,5);! M; u  m3 ]8 B3 s$ q) ~4 Y8 f) L$ y; W! [
figure(1);
& T0 G  S. k+ d6 ?plot(i,j,'.-');
# r+ C) _, ^1 i* W- j9 Ihold on;8 ^% p) h) [- \, F# G
i=i';0 Y( M: _0 L. ~. X
p1=polyfit(i,j,2)- K+ o! I3 j5 A3 D7 v) P& g  e, c
f1=polyval(p1,i);
+ z( Z/ t( B+ O1 z  Eplot(i,j,'.',i,f1,'-');
2 W/ ^$ T# @, l' g$ khold on;' d4 ~$ ]/ M0 S6 Y$ L
title('病人平均等待时间与病人到达时间的关系');, z9 N5 B) p3 ]/ Y" @
xlabel('病人门诊时间距离7 月12 日的时间差');ylabel('病人平均等待时间
0 t- Y/ @: J/ w! z; m+ A3 e& n');
+ m8 ?5 y% L4 ~0 Q9 O# {axis equal;
( }% x! k  C$ }! g2 x6 ufigure(2);$ \+ Y. w) c! \+ E/ Q  M! ]' A  t
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
1 C' \4 @5 n7 ~* r" W3 D$ F: L32 34 35 36 37 38 39 40 41 42];
4 n. i* x8 ?) G, |1 A; F0 R' \) |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- ^( N2 A( p" Y7 M
11.5 13.25 14 13 12 11 13 13 12 12.25 13 15.5 15 14 13];- v/ O/ }/ n( {) o
subplot(1,2,1);
5 ]4 X+ y4 R+ x% T# P+ hplot(A1,B1,'.-');+ D& j7 G. G. k+ @& c/ D
hold on;8 W/ e, B' I9 E4 S: n% w
p2=polyfit(A1,B1,2)! u& P. o" D3 x) ~$ `
f2=polyval(p2,A1);
) b7 n" x- `9 x" I& Zplot(A1,B1,'.',A1,f2,'-');
9 c- R) G, V: ]/ @0 ~# U% [hold on;, r, l: t4 z+ D* T! h! @
title('白内障病人情况');
; E9 {5 K; f; ?- L% P+ G5 Dxlabel('白内障病人门诊时间距7月12日的时间差');ylabel('白内障病人平均
; L! l# @  ^# |4 N4 j等待时间');$ X, {4 p8 |, n- @; a  A  O
axis equal;
6 N; \% o! Q2 O' IA2=[2 3 4 5 6 7 8 9 11 13 15 16 17 18 19 20 21 22 23 25 26 27 28 29 302 p  @$ N! K9 l/ C5 S
31 32 33 34 35 36 37 ];( v' ^9 D0 J' G- O5 g* v! h
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
5 F% e3 E* `% D1 c; E, F12 12 13.25 14 14 13 12 12 12.4 13.33 13 13 13 12];: Q- Z9 T0 N" A  ]" D
subplot(1,2,2);$ J+ F! D  o" g1 L; b$ U
plot(A2,B2,'.-');
+ T0 D' V3 ?7 F" e# b- K6 \hold on;' F4 j) i' ~) ^* y: ~
p3=polyfit(A2,B2,2)+ }9 C- X0 W0 t+ V7 j5 `9 P
f3=polyval(p3,A2);( ^# s! @$ F  j( m
plot(A2,B2,'.',A2,f3,'-');
! n- R* }8 `3 {5 v' \hold on;
6 L3 J/ V7 b4 \# z/ @2 ytitle('白内障双眼病人情况');- s+ U5 ^' v1 j7 `* T
xlabel('白内障双眼病人门诊时间距7月12日的时间差');ylabel('白内障双眼
2 J2 u" u; Z6 k病人平均等待时间');
& @! |5 X$ D. W/ V* haxis equal;
( {9 u4 L1 W; o" ?% A3 Ifigure(3);
6 Z8 Q- a; e4 K7 c. qA3=[1 2 3 4 5 6 7 8 10 11 12 13 15 16 17 18 19 20 21 22 23 24 25 26 275 J( B/ r) y: `7 G- J+ h
28 29 30 31 32 33 34 35 36 37 38 ];
6 G5 m7 J! I; W5 ]9 TB3=[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
) A7 m5 E6 ]0 G5 q: Y/ M7 o11.67 12 12.25 12 11 11.75 13 14 14 13 12 12 12.33 13 13 13 12.5 12 12];( P$ V4 s, l+ G
subplot(2,2,1);
4 U. R9 F6 c8 R$ D" y  h, W3 b; splot(A3,B3,'.-');1 H. z. Q' S& i* T4 X+ b# m2 ~0 H
hold on;
5 f. H: R( x' g" K% jp4=polyfit(A3,B3,2)/ e9 j( X# G$ O7 N; c/ T
f4=polyval(p4,A3);& F$ A8 V/ Z4 l, s5 W% U
plot(A3,B3,'.',A3,f4,'-');
3 ~% A, \. s. A' F# L* shold on;
- s# k  e6 @$ u3 Otitle('视网膜病人情况');
% m. G4 {3 {" Y) ^4 g$ oxlabel('门诊时间距7 月12 日的时间差');ylabel('视网膜病人平均等待时间
- e' c6 m6 \* Z3 p' n) w3 B');
5 C4 }( h' \, o7 [% h' ~/ ^, jaxis equal;2 I1 V) H1 t- l" n8 M
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 _( s' W$ S, P( C0 {$ p
37];
7 d" a5 w8 m! R* }  E& s+ `1 gB4=[12 13 12 12 12 12 13.5 13 13.5 14 13 12 12 11 10.25 12 11 11 15 15. F, Y. B+ Z' P: z
13 12 12 13 13 12];0 f% K4 F: O* m8 w
subplot(2,2,2);, E+ D9 Z8 u" D* d- V3 Q
plot(A4,B4,'.-');- F- D. X" A! n% a+ b2 Z2 o  v
hold on;3 s* y9 M( U3 x8 }
p5=polyfit(A4,B4,2)
) x: C" `& F; c: N- P" I9 Sf5=polyval(p5,A4);! R# y, M% n% G) n1 A/ k! K
plot(A4,B4,'.',A4,f5,'-');3 C8 C8 @3 ^, ^% W0 g8 F
hold on;; i- t6 q/ l4 p/ t9 j" t; q* A
title('青光眼病人情况');
5 a4 D/ {) [* ~! }xlabel('门诊时间距7 月12日的时间差');ylabel('视网膜病人平均等待时间
0 K1 `: P; ]6 m" [: Q');8 |) l8 E" a8 r
axis equal;" D- a0 q% ]' B4 @6 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 34
2 T! H# |# u& E4 l37];
8 |+ f' `6 y- {% [# dB4=[12 13 12 12 12 12 13.5 13 13.5 14 13 12 12 11 10.25 12 11 11 15 15
. F& y0 F1 N- v' U2 i* P3 P13 12 12 13 13 12];
  Z  N! ?- `1 W7 K$ x* u8 Bsubplot(2,2,2);
* ^; {4 z/ m& G) _# v( {* Tplot(A4,B4,'.-');: B7 G) e/ `( o
hold on;
+ }1 z. p0 Q# l7 V; rp5=polyfit(A4,B4,2)0 D+ R" v* R5 e/ C9 H
f5=polyval(p5,A4);4 F0 c! |4 t$ C7 z9 Z) a8 A3 y
plot(A4,B4,'.',A4,f5,'-');
2 K$ w$ i" T5 G. D8 y8 b. Qhold on;
' d, a4 T+ T" ?4 F: `title('青光眼病人情况');. H* e: z6 _+ s' K7 P3 Z
xlabel('门诊时间距7 月12 日的时间差');ylabel('青光眼病人平均等待时间
; B1 [% z; f' b: R');
1 i0 {7 M1 A2 a+ T+ [* z  p9 D9 ^3 ^axis equal;6 X  V& \5 W. t+ |& |1 d) l
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];
! {/ {( ^4 A4 \7 ~/ q) V( Y+ }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];) c& x7 B" S! E$ `- n
subplot(2,2,3);. B" w4 S% ^9 Y1 m" J: d
plot(A5,B5,'.-');1 P: E5 c; r& _8 T
hold on;
) e8 H: [* z% o3 V0 T2 i7 Ep6=polyfit(A5,B5,2)7 s8 _. }' R5 Y
f6=polyval(p6,A5);! R0 }# q( L  l6 U6 a" _
plot(A5,B5,'.',A5,f6,'-');% f1 t' Z' y% P# O/ C5 n
hold on;
" y0 }+ n! U* p7 i, otitle('外伤病人情况');0 j% t) X! X' C- T: R- g/ F
xlabel('外伤病人门诊时间距离7月12日的时间差');ylabel('外伤病人平均等$ Y+ W- p1 P0 K) z
待时间');
+ B) Y7 Q% x% A$ v4 G" t3 Yaxis equal;
' Z; R6 P* v5 H  r, _% `5 L二、层次分析法的一致性指标判定的matlab 实现:# ]$ ~# u; b( O% W  \! S
function H=ff(A)4 W7 z6 K) h, C9 ]& y
t=eig(A)
1 d: F$ y6 H4 Sw=max(t), O' ?1 W$ T4 e+ \) R" n
q=(w-4)/3/.90 w  Q, N# K% O' |! C3 t
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(
2 i  x" @8 ^1 ^9)+A(10)+A(11)+A(12)+A(13)+A(14)+A(15)+A(16))6 p+ u( N* i' t+ P: E
(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)
* Z) T' `1 L; O! g* A6 }) {0 x% A- W+A(10)+A(11)+A(12)+A(13)+A(14)+A(15)+A(16))  {9 Z- c1 t0 N, |
(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)
4 @! k* E" u# ?+A(10)+A(11)+A(12)+A(13)+A(14)+A(15)+A(16))  y! K$ n8 y" O0 W
(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)
/ P4 o2 }) Y; n7 r1 I. P+A(10)+A(11)+A(12)+A(13)+A(14)+A(15)+A(16))]
  D& g" p. N8 K0 [, v6 z5 P6 ]+ M

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

听众

3591

积分

逍遥游

  • 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

    听众

    3591

    积分

    逍遥游

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

    [LV.5]常住居民I

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

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

    群组Matlab讨论组

    群组数学建模

    群组小草的客厅

    群组2012数学一考研交流

    群组C 语言讨论组

    回复 迷途的羊 的帖子
    * w/ h: S3 Y' b* W
    ( Q& o4 u2 j/ K8 F% A# Y是自己编的还用别人运行么……不过要是楼主自己编的……小弟就佩服死啦8 a! p9 r- ]- O. }
       
    回复

    使用道具 举报

    紫辰 实名认证       

    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');
    ! Y2 O5 d, c' Q9 [$ I8 q- r5 @在不同的地方运行上面一句要改成本地数据保存的路径的
    回复

    使用道具 举报

    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-5-26 04:42 , Processed in 3.601714 second(s), 108 queries .

    回顶部