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

wuming1429

一、使用曲线拟合模拟病人等待时间的模拟matlab程序实现：
[A]=xlsread('C:\MATLAB7\work\问题3数据5');
len=length(A);
. x) Z; O5 k6 i" q+ @5 v; {) ~i=1:len;
( o  J& `; f6 N9 Y& Nj=A(:,5);
figure(1);
plot(i,j,'.-');
hold on;
i=i';
: }. v3 S' i1 Rp1=polyfit(i,j,2)
5 \  V+ s! D* F; S/ i( e5 c: rf1=polyval(p1,i);
: w( U2 {9 a& E9 I) O) J- f, tplot(i,j,'.',i,f1,'-');
hold on;
title('病人平均等待时间与病人到达时间的关系');
6 y# y5 V# I+ `xlabel('病人门诊时间距离7 月12 日的时间差');ylabel('病人平均等待时间
');
axis equal;
8 j, u; E( t# ^' {& Yfigure(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
32 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
- F$ k- O  L8 P8 I, h" W1 f11.5 13.25 14 13 12 11 13 13 12 12.25 13 15.5 15 14 13];
subplot(1,2,1);
1 j% S5 l6 R" Iplot(A1,B1,'.-');
; P# X" E! Y$ U( [hold on;
p2=polyfit(A1,B1,2)
0 F% m4 c6 t6 q0 G* d" t, Ff2=polyval(p2,A1);
plot(A1,B1,'.',A1,f2,'-');
hold on;
1 w5 t; f! C$ `8 g: ^6 q7 d+ z2 m5 Wtitle('白内障病人情况');
7 Q5 o* b0 B7 Fxlabel('白内障病人门诊时间距7月12日的时间差');ylabel('白内障病人平均
等待时间');
- i' g9 m6 {$ H6 J: }( \' J) taxis equal;
8 F5 K$ G/ X! p3 d; L& C/ TA2=[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
# @; O6 z0 e/ i$ ^- f" I31 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
4 A* C1 g6 z& D9 _# d3 ?12 12 13.25 14 14 13 12 12 12.4 13.33 13 13 13 12];
subplot(1,2,2);
plot(A2,B2,'.-');
: v: g6 ~" e8 a/ S3 [hold on;
p3=polyfit(A2,B2,2)
f3=polyval(p3,A2);
plot(A2,B2,'.',A2,f3,'-');
- S1 Q: U7 Q% [) {hold on;
9 p( ], b* ]4 B. |; Z: P) ?title('白内障双眼病人情况');
xlabel('白内障双眼病人门诊时间距7月12日的时间差');ylabel('白内障双眼
9 r" {; n0 f3 l: X( a病人平均等待时间');
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 ];
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];
0 [4 p& j+ W0 ^) F: }subplot(2,2,1);
. Z& X* \+ h! }- Z3 M8 |5 k  gplot(A3,B3,'.-');
hold on;
p4=polyfit(A3,B3,2)
; v$ F/ x" Y. M4 \" P: ef4=polyval(p4,A3);
: y3 @. R$ f$ o& Jplot(A3,B3,'.',A3,f4,'-');
hold on;
title('视网膜病人情况');
xlabel('门诊时间距7 月12 日的时间差');ylabel('视网膜病人平均等待时间
');
axis equal;
6 T  E! X& I3 H) Q9 a" ?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
1 o5 p, w5 _' @- x  G' H8 A13 12 12 13 13 12];
7 y5 O* J: @) K) Xsubplot(2,2,2);
* C  ~( K5 ?9 Mplot(A4,B4,'.-');
hold on;
/ d, Z4 _. t; d. Mp5=polyfit(A4,B4,2)
f5=polyval(p5,A4);
8 `1 a8 u4 ~6 V4 Oplot(A4,B4,'.',A4,f5,'-');
hold on;
title('青光眼病人情况');
xlabel('门诊时间距7 月12日的时间差');ylabel('视网膜病人平均等待时间
');
axis equal;
3 N/ ]3 L" d2 @8 ]% a( P0 NA4=[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
" C9 ?1 X/ Q' G% i: [1 J37];
/ @8 M. y; j' e, I% xB4=[12 13 12 12 12 12 13.5 13 13.5 14 13 12 12 11 10.25 12 11 11 15 15
2 }# }* S- |  G! o7 S6 r6 f+ x13 12 12 13 13 12];
: ?( z2 U, E# g& F; @/ usubplot(2,2,2);
3 ], x% x4 q( q1 z2 oplot(A4,B4,'.-');
hold on;
2 B5 @# M' f/ L+ c2 rp5=polyfit(A4,B4,2)
f5=polyval(p5,A4);
2 u3 |" n( W4 X1 M: zplot(A4,B4,'.',A4,f5,'-');
7 d5 R* c6 t2 G# I* {hold on;
8 p# o& f0 _% z1 Ctitle('青光眼病人情况');
xlabel('门诊时间距7 月12 日的时间差');ylabel('青光眼病人平均等待时间
4 p- l# h$ M' S2 X: w- W');
axis 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];
8 R( Q" H8 N' r( U" yB5=[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);
plot(A5,B5,'.-');
hold on;
p6=polyfit(A5,B5,2)
f6=polyval(p6,A5);
plot(A5,B5,'.',A5,f6,'-');
, m2 G& m0 t/ ]2 x+ shold on;
title('外伤病人情况');
0 A4 J6 h6 T8 a- E. m& d$ Y0 R* uxlabel('外伤病人门诊时间距离7月12日的时间差');ylabel('外伤病人平均等
待时间');
, [$ a- o+ [- H+ Gaxis equal;
二、层次分析法的一致性指标判定的matlab 实现:
& H0 u  w' h* W, Mfunction H=ff(A)
1 j1 h) U' G0 v6 s' w+ `t=eig(A)
  V9 T3 H. h5 r* Cw=max(t)
q=(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(
: Q7 G6 |" O4 a* e, A& y. f9)+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)
) _4 h) d" e! H* a' W+A(10)+A(11)+A(12)+A(13)+A(14)+A(15)+A(16))
/ S& x1 Y% M  g. Y0 V(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)
' M) X6 Y$ p- E2 K+A(10)+A(11)+A(12)+A(13)+A(14)+A(15)+A(16))
(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)
" h8 F/ j6 @5 q* ~6 \' W: P+A(10)+A(11)+A(12)+A(13)+A(14)+A(15)+A(16))]
5 V& }% ^4 N' S0 h6 s8 P9 p

linmatsas

运行不了嘛？那就是错了呗……

迷途的羊

你很强大哦，编这么 长的程序

linmatsas

回复 迷途的羊 的帖子
4 H  i" X- e, ]! v
8 E  @2 f5 I& Y' q3 d; u是自己编的还用别人运行么……不过要是楼主自己编的……小弟就佩服死啦
   

紫辰

你可以在不同 的MATLAB的版本试试看吧

chenrye121

佩服两个字…………

水木年华zzu

[A]=xlsread('C:\MATLAB7\work\问题3数据5');
; I8 q+ t4 G% a+ o在不同的地方运行上面一句要改成本地数据保存的路径的

wuming1429

数据要转换成矩阵，不知怎么弄啊？汉字数据读取不了啊

hhao

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

abc321

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