数学建模社区-数学中国

标题: 请问MATLAB编程高手啊谁能将这个程序运行啊 [打印本页]

---

作者: wuming1429    时间: 2010-8-1 20:45
标题: 请问MATLAB编程高手啊谁能将这个程序运行啊
一、使用曲线拟合模拟病人等待时间的模拟matlab程序实现：
[A]=xlsread('C:\MATLAB7\work\问题3数据5');
len=length(A);
i=1:len;
+ B; Y  s! ^. ]  U, y4 v6 o# qj=A(:,5);
" D0 N( X, Z; Lfigure(1);
plot(i,j,'.-');
( x: @& s; k. h. chold on;
i=i';
0 M2 e% I+ P* F$ t: U" u# E" Fp1=polyfit(i,j,2)
) H2 S( x8 J9 U- E3 u) P, ?f1=polyval(p1,i);
8 C/ ], P6 J5 d) X" V* J# U4 Rplot(i,j,'.',i,f1,'-');
hold on;
title('病人平均等待时间与病人到达时间的关系');
2 q" e1 Q( p# _' q5 D, z1 y4 Oxlabel('病人门诊时间距离7 月12 日的时间差');ylabel('病人平均等待时间
5 a" Y/ X. ^, g) B: U');
axis equal;
figure(2);
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
- \/ V8 n* m+ e% a6 E32 34 35 36 37 38 39 40 41 42];
3 {# s# Y  s8 g" @+ L: lB1=[12 12 12 12 12 12.5 13 13 13 13 13.4 14 14 13 12.5 11 11 12 12 12 11
11.5 13.25 14 13 12 11 13 13 12 12.25 13 15.5 15 14 13];
) h" A; z5 U1 E# j" r, ?" vsubplot(1,2,1);
; v- C& I% Z$ d' ~/ C: dplot(A1,B1,'.-');
' ]/ B* P1 M: h4 T" G8 jhold on;
) G* J5 b( r9 n0 @1 Ep2=polyfit(A1,B1,2)
f2=polyval(p2,A1);
3 W- {6 M  P2 z6 l0 V* v9 {8 cplot(A1,B1,'.',A1,f2,'-');
/ F9 ^( W4 x' c3 l$ f1 Zhold on;
title('白内障病人情况');
8 ?1 _. N1 K( N4 h( k2 dxlabel('白内障病人门诊时间距7月12日的时间差');ylabel('白内障病人平均
) `# x/ E9 f$ R等待时间');
1 p& I  G5 {" S; p; A; j9 L% Aaxis equal;
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
31 32 33 34 35 36 37 ];
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
12 12 13.25 14 14 13 12 12 12.4 13.33 13 13 13 12];
" l! g* F% f* t7 D4 D3 ssubplot(1,2,2);
- d/ i# X) ?. X/ s, wplot(A2,B2,'.-');
hold on;
p3=polyfit(A2,B2,2)
  P( g4 ]: F! u4 ?& ?f3=polyval(p3,A2);
0 w. c8 `# Y+ v( v; r7 R0 L" a6 w4 Mplot(A2,B2,'.',A2,f3,'-');
hold on;
title('白内障双眼病人情况');
3 t3 R9 c8 F  R+ v. i7 H& S* C6 sxlabel('白内障双眼病人门诊时间距7月12日的时间差');ylabel('白内障双眼
, d" S4 `) _8 a. [' n# X4 w病人平均等待时间');
axis equal;
  Z; c( b6 ~/ T  r2 c) P- v/ P/ @+ _: v! Zfigure(3);
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
- ?% Y" D4 ~/ p28 29 30 31 32 33 34 35 36 37 38 ];
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.75
11.67 12 12.25 12 11 11.75 13 14 14 13 12 12 12.33 13 13 13 12.5 12 12];
3 m' F1 n8 F4 @3 P6 ssubplot(2,2,1);
0 T5 t9 y4 Z  \: ]  M- g6 rplot(A3,B3,'.-');
hold on;
p4=polyfit(A3,B3,2)
f4=polyval(p4,A3);
5 _- s1 p2 }+ M, Gplot(A3,B3,'.',A3,f4,'-');
$ A' A; _4 ~. H! hhold on;
+ {1 x0 B3 b  i- ]7 B( [) F( C# \title('视网膜病人情况');
! X5 K4 `5 X6 N* l2 \1 d) oxlabel('门诊时间距7 月12 日的时间差');ylabel('视网膜病人平均等待时间
" g3 |1 u  k! E  u+ W' ]# `& d# b+ s');
axis equal;
1 q4 F: x# i7 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 34
37];
# s2 j# D2 N" K8 j" h8 {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) b) x+ H9 n9 B13 12 12 13 13 12];
/ s( R% ^8 P+ s% m, \! I- o$ esubplot(2,2,2);
plot(A4,B4,'.-');
hold on;
p5=polyfit(A4,B4,2)
f5=polyval(p5,A4);
, B) S$ w" y- _" I: }plot(A4,B4,'.',A4,f5,'-');
hold on;
9 P7 d  @1 x. ~2 j9 |% |title('青光眼病人情况');
xlabel('门诊时间距7 月12日的时间差');ylabel('视网膜病人平均等待时间
: n" X! z7 P6 X! a' R& L1 p');
axis equal;
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
3 `  x* V2 P9 {3 }$ K3 Q) ^9 K37];
B4=[12 13 12 12 12 12 13.5 13 13.5 14 13 12 12 11 10.25 12 11 11 15 15
13 12 12 13 13 12];
4 C/ r# p7 l: s3 M3 D& g8 \3 gsubplot(2,2,2);
plot(A4,B4,'.-');
- f. M5 E4 ?4 z0 xhold on;
p5=polyfit(A4,B4,2)
4 a0 [- o7 V+ D. `. Y0 a/ r. yf5=polyval(p5,A4);
plot(A4,B4,'.',A4,f5,'-');
( R0 L* `) r$ F: D$ x) D) Dhold on;
3 @. O5 R& D5 v6 b% L: u+ vtitle('青光眼病人情况');
" _" k1 }: [1 e0 q+ lxlabel('门诊时间距7 月12 日的时间差');ylabel('青光眼病人平均等待时间
');
  J7 ]; y1 X" l+ O, Vaxis equal;
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];
2 ]( Z' @8 s6 e% b. _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];
, L9 V4 p0 j6 ]# ]8 \) H# Osubplot(2,2,3);
plot(A5,B5,'.-');
- a# ?: L4 e7 H2 Zhold on;
p6=polyfit(A5,B5,2)
f6=polyval(p6,A5);
plot(A5,B5,'.',A5,f6,'-');
hold on;
- ^% r+ j8 k8 T8 ptitle('外伤病人情况');
% ~: J8 B0 s9 Z% A; }( Q. cxlabel('外伤病人门诊时间距离7月12日的时间差');ylabel('外伤病人平均等
待时间');
axis equal;
2 v) Q5 b- [; ]- w8 c9 k' [二、层次分析法的一致性指标判定的matlab 实现:
$ y6 Q0 V; o( f% f/ zfunction H=ff(A)
t=eig(A)
w=max(t)
: k* K6 q3 A; mq=(w-4)/3/.9
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(
9)+A(10)+A(11)+A(12)+A(13)+A(14)+A(15)+A(16))
9 N; [5 G) M# P' E0 t(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)
+A(10)+A(11)+A(12)+A(13)+A(14)+A(15)+A(16))
' A$ T# C* P: r7 O% c9 y(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)
; @' O% O  m% S3 s* N- k7 r$ }5 b+A(10)+A(11)+A(12)+A(13)+A(14)+A(15)+A(16))
% A9 a1 F$ b0 `( v(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)
+A(10)+A(11)+A(12)+A(13)+A(14)+A(15)+A(16))]
8 a, `* M# a: }& }8 \( K

B2009.doc

2010-8-1 20:44 上传

点击文件名下载附件

下载积分: 体力 -2 点

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

09年b题

2009B题0315(1).pdf

2010-8-1 20:44 上传

点击文件名下载附件

下载积分: 体力 -2 点

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

问题三数据五.xls

2010-8-1 20:45 上传

点击文件名下载附件

下载积分: 体力 -2 点

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

---

作者: linmatsas    时间: 2010-8-1 21:39
运行不了嘛？那就是错了呗……

---

作者: 迷途的羊    时间: 2010-8-1 21:53
你很强大哦，编这么 长的程序

---

作者: linmatsas    时间: 2010-8-1 22:38
回复 迷途的羊 的帖子

是自己编的还用别人运行么……不过要是楼主自己编的……小弟就佩服死啦
/ M/ Q- x, v2 w. v# e+ q/ t2 K0 {   

---

作者: 紫辰    时间: 2010-8-2 00:06
你可以在不同 的MATLAB的版本试试看吧

---

作者: chenrye121    时间: 2010-8-2 14:06
佩服两个字…………

---

作者: 水木年华zzu    时间: 2010-8-2 16:46
[A]=xlsread('C:\MATLAB7\work\问题3数据5');
在不同的地方运行上面一句要改成本地数据保存的路径的

---

作者: wuming1429    时间: 2010-8-2 19:06
数据要转换成矩阵，不知怎么弄啊？汉字数据读取不了啊

---

作者: hhao    时间: 2010-8-26 11:33
我基本上是采用看英语文章的办法，先泛读，再精读，再一句一句看，最后再提纲挈领，总算是明白一点了，当然，也可能还是领悟错了。最后要说的一句话是：楼主，你很牛叉，希望你不是真的有病。   

---

作者: abc321    时间: 2010-8-26 11:34
试试运气啦~~~~~~~~~~~

---

作者: snrl    时间: 2010-8-26 11:35
楼主，你写得实在是太好了。我惟一能做的，就只有把这个帖子顶上去这件事了

---

作者: chshfxfx    时间: 2010-8-26 11:36
顶顶更健康，越顶吃的越香。

---

作者: icm    时间: 2010-8-26 11:38
顶顶更健康，越顶吃的越香。

---

作者: spss    时间: 2010-8-26 11:46
楼主，你写得实在是太好了。我惟一能做的，就只有把这个帖子顶上去这件事了

---

作者: william745    时间: 2010-8-26 13:52
来报道!!!!!!!!!!!

---

作者: huxiao9026    时间: 2010-8-26 18:59
哦～～

---

作者: racheltong    时间: 2010-8-26 19:17
顶顶更健康，越顶吃的越香。

---

作者: xphoenix    时间: 2010-8-27 00:00
我来了~~~~~~~~~ 闪人~~~~~~~~~~~~~~~~  

---

作者: lxmh815    时间: 2010-8-27 12:00
留个脚印```````

---

作者: nn58123    时间: 2010-8-27 15:00
声明一下：本人看贴和回贴的规则，好贴必看，精华贴必回。

---

作者: yunwuya    时间: 2010-8-27 20:00
强烈支持。楼主万岁

---

作者: huhu0079    时间: 2010-8-28 12:00
声明一下：本人看贴和回贴的规则，好贴必看，精华贴必回。

---

作者: lwx193520    时间: 2010-8-28 15:00
强烈支持。楼主万岁

---

作者: zjf822    时间: 2010-8-28 20:00
顶顶更健康，越顶吃的越香。

---

作者: diaohaiq    时间: 2010-8-28 23:59
我要把这个帖子一直往上顶，往上顶！

---

作者: qiuqirong    时间: 2010-8-29 08:00
试试运气啦~~~~~~~~~~~

---

作者: ruiqiu    时间: 2010-8-29 12:00
我基本上是采用看英语文章的办法，先泛读，再精读，再一句一句看，最后再提纲挈领，总算是明白一点了，当然，也可能还是领悟错了。最后要说的一句话是：楼主，你很牛叉，希望你不是真的有病。   

---

作者: casper    时间: 2010-8-29 15:00
我基本上是采用看英语文章的办法，先泛读，再精读，再一句一句看，最后再提纲挈领，总算是明白一点了，当然，也可能还是领悟错了。最后要说的一句话是：楼主，你很牛叉，希望你不是真的有病。   

---

作者: 星际丑男    时间: 2010-8-29 20:00
强烈支持。楼主万岁

---

作者: hab0458    时间: 2010-8-30 12:00
楼主那种裂纸欲出的大手笔，竟使我忍不住一次次的翻开楼主的帖子……   

---

作者: liguoli    时间: 2010-8-30 15:00
楼主的帖子实在是写得太好了。可是我立刻想到，这么好的帖子，倘若别人看不到，那么不是浪费楼主的心血吗？经过痛苦的思想斗争，我终于下定决心，牺牲小我，奉献大我。我要拿出这帖子奉献给世人赏阅，我要把这个帖子一直往上顶，往上顶！顶到所有人都看到为止！  

---

作者: Hyacinth    时间: 2010-8-30 20:00
我来了~~~~~~~~~ 闪人~~~~~~~~~~~~~~~~  

---

作者: hhao    时间: 2010-8-31 08:00
我来了~~~~~~~~~ 闪人~~~~~~~~~~~~~~~~  

---

作者: qiuhanty    时间: 2010-8-31 12:00
声明一下：本人看贴和回贴的规则，好贴必看，精华贴必回。

---

作者: ilava    时间: 2010-8-31 15:00
试试运气啦~~~~~~~~~~~

---

作者: dongdong111    时间: 2012-6-19 15:31
其实代码质量不高。。。。明显可视化设置没到位。然后后面的AHP写的也很乱还有bug

---

欢迎光临 数学建模社区-数学中国 (http://www.madio.net/)

Powered by Discuz! X2.5