请问MATLAB编程高手啊谁能将这个程序运行啊

wuming1429

一、使用曲线拟合模拟病人等待时间的模拟matlab程序实现：
+ b* a& n- \* r& u1 O[A]=xlsread('C:\MATLAB7\work\问题3数据5');
len=length(A);
- V" J7 ^6 ?8 K+ Di=1:len;
j=A(:,5);
figure(1);
8 e3 ~* g$ ^0 t7 \% vplot(i,j,'.-');
1 U; R6 ]4 C' ?5 G' n( P2 xhold on;
. a& k5 d2 D+ \' si=i';
p1=polyfit(i,j,2)
f1=polyval(p1,i);
( M/ G+ U6 S1 L. n8 l! Yplot(i,j,'.',i,f1,'-');
/ y2 ]; n7 }  l9 G' Chold on;
title('病人平均等待时间与病人到达时间的关系');
  `( H1 A! p6 r5 n! S. nxlabel('病人门诊时间距离7 月12 日的时间差');ylabel('病人平均等待时间
');
, N9 w6 P- F2 ^3 C& o8 y6 C! Faxis 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
& S. d. `2 R' Z2 T! h/ y/ P+ g32 34 35 36 37 38 39 40 41 42];
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
, y  J5 ?* \1 F# Z11.5 13.25 14 13 12 11 13 13 12 12.25 13 15.5 15 14 13];
subplot(1,2,1);
) Y* y% U' Y- }  o/ iplot(A1,B1,'.-');
0 s9 X& C: ]" Bhold on;
9 @& U+ e4 O" I" s/ U. a! u, Op2=polyfit(A1,B1,2)
f2=polyval(p2,A1);
+ I9 B3 @8 v! H4 M0 [$ W$ iplot(A1,B1,'.',A1,f2,'-');
hold on;
title('白内障病人情况');
7 S. F1 ]" p' D5 Oxlabel('白内障病人门诊时间距7月12日的时间差');ylabel('白内障病人平均
) P% ]) ~7 C& B/ f5 _. x- o: T等待时间');
axis equal;
: N" `9 g; P# m5 CA2=[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
# M* G& c# w0 T' U  u4 `( v6 r5 k31 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
1 w! U8 _9 \& O12 12 13.25 14 14 13 12 12 12.4 13.33 13 13 13 12];
subplot(1,2,2);
3 {  p, k$ _7 T0 z+ u( w& F/ |8 Rplot(A2,B2,'.-');
: J1 [7 `% o9 ]; hhold on;
p3=polyfit(A2,B2,2)
, l2 B# E; }5 gf3=polyval(p3,A2);
plot(A2,B2,'.',A2,f3,'-');
hold on;
$ F! H1 l  M0 N8 ttitle('白内障双眼病人情况');
8 x1 P5 T/ e8 `) G4 Yxlabel('白内障双眼病人门诊时间距7月12日的时间差');ylabel('白内障双眼
病人平均等待时间');
axis equal;
figure(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
28 29 30 31 32 33 34 35 36 37 38 ];
0 G9 a; L* t7 t! D  A4 `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];
4 ~8 a: I2 M# U  q/ E4 _5 C4 F) F# C( gsubplot(2,2,1);
" f1 |$ o/ N; a& |& H) [: ]7 s/ b0 Nplot(A3,B3,'.-');
9 R$ Z4 J% |: ]: ]1 d$ N* L) F2 Y2 N* @! Khold on;
p4=polyfit(A3,B3,2)
6 v4 M! }* R$ Y# A2 wf4=polyval(p4,A3);
6 j* q3 U8 X! S9 s9 Yplot(A3,B3,'.',A3,f4,'-');
hold on;
title('视网膜病人情况');
xlabel('门诊时间距7 月12 日的时间差');ylabel('视网膜病人平均等待时间
5 ]2 X4 ]- F1 |3 Y+ ^');
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
37];
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];
9 ]1 G# P0 }1 h% j; z. ^/ psubplot(2,2,2);
) }" b0 A% d4 y( eplot(A4,B4,'.-');
hold on;
' Q* B+ l# ^& E  f3 E! ?; Sp5=polyfit(A4,B4,2)
f5=polyval(p5,A4);
plot(A4,B4,'.',A4,f5,'-');
  \# W' e8 x% B4 ~+ mhold on;
title('青光眼病人情况');
0 m' w& Y4 P7 _8 mxlabel('门诊时间距7 月12日的时间差');ylabel('视网膜病人平均等待时间
');
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
37];
' y0 W& W) K4 i/ Z8 c* U  qB4=[12 13 12 12 12 12 13.5 13 13.5 14 13 12 12 11 10.25 12 11 11 15 15
5 h; G  `' H- i3 [13 12 12 13 13 12];
% [% N" B+ O) {( B8 w' Bsubplot(2,2,2);
plot(A4,B4,'.-');
hold on;
p5=polyfit(A4,B4,2)
$ ~3 \% `5 ~8 r  Cf5=polyval(p5,A4);
plot(A4,B4,'.',A4,f5,'-');
' o4 X% _- v1 P0 Rhold on;
title('青光眼病人情况');
2 I9 }' U. L7 u+ |. ]xlabel('门诊时间距7 月12 日的时间差');ylabel('青光眼病人平均等待时间
');
# z: z" ?& w; {axis equal;
1 z/ ?) f# \- B$ u8 qA5=[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];
" _8 ~9 V' g7 v4 IB5=[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];
subplot(2,2,3);
$ J$ Y4 ]" \6 F# c3 r9 Z- H- Iplot(A5,B5,'.-');
hold on;
0 F- H9 o, \# }$ a) ^, J. Pp6=polyfit(A5,B5,2)
7 w4 \( M  e: S3 }) R6 a7 of6=polyval(p6,A5);
plot(A5,B5,'.',A5,f6,'-');
hold on;
7 B( c7 D$ m6 Q$ H# I3 Ktitle('外伤病人情况');
; n8 X+ `# l7 \( Exlabel('外伤病人门诊时间距离7月12日的时间差');ylabel('外伤病人平均等
4 n2 }5 F' r' X8 V; t2 }* k  l待时间');
axis equal;
二、层次分析法的一致性指标判定的matlab 实现:
function H=ff(A)
t=eig(A)
" q6 ~3 M" B7 d5 p* \) V4 _w=max(t)
& O' U  X/ B, }8 Wq=(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(
7 C; l" V: X: y( E) m4 S+ v9)+A(10)+A(11)+A(12)+A(13)+A(14)+A(15)+A(16))
(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)
0 M  H4 z3 j3 T+A(10)+A(11)+A(12)+A(13)+A(14)+A(15)+A(16))
# F( d+ z, s* u1 P' t& 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)
1 u6 {% [4 v' \3 I: \8 n1 I) b+A(10)+A(11)+A(12)+A(13)+A(14)+A(15)+A(16))
( D, L/ {" h! }* z3 `. ^9 ^) ](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)
# ^1 ?. H& m* K0 G$ Q- g) v+A(10)+A(11)+A(12)+A(13)+A(14)+A(15)+A(16))]
$ Y9 s- X: N' {0 c

dongdong111

其实代码质量不高。。。。明显可视化设置没到位。然后后面的AHP写的也很乱还有bug

ilava

试试运气啦~~~~~~~~~~~

qiuhanty

声明一下：本人看贴和回贴的规则，好贴必看，精华贴必回。

hhao

我来了~~~~~~~~~ 闪人~~~~~~~~~~~~~~~~  

Hyacinth

我来了~~~~~~~~~ 闪人~~~~~~~~~~~~~~~~  

liguoli

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

hab0458

楼主那种裂纸欲出的大手笔，竟使我忍不住一次次的翻开楼主的帖子……   

星际丑男

强烈支持。楼主万岁

casper

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