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

wuming1429

一、使用曲线拟合模拟病人等待时间的模拟matlab程序实现：
[A]=xlsread('C:\MATLAB7\work\问题3数据5');
4 h: C$ B% g9 }, n  u1 s3 H; _1 ylen=length(A);
i=1:len;
. u; q: `; d2 D) m) J4 ej=A(:,5);
8 v5 V/ t( A0 R9 H* G2 s- h# Z  Tfigure(1);
) L2 V: T/ H  }8 ]3 Uplot(i,j,'.-');
hold on;
i=i';
p1=polyfit(i,j,2)
2 g# u5 q0 G. o  Ef1=polyval(p1,i);
plot(i,j,'.',i,f1,'-');
( L# j  o. c* W9 l0 Rhold on;
title('病人平均等待时间与病人到达时间的关系');
( ?% @: l4 ]8 x' X+ l- J7 Exlabel('病人门诊时间距离7 月12 日的时间差');ylabel('病人平均等待时间
  ~* X1 _- f- p; ^');
axis equal;
" X2 f, u$ b0 [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
3 J, B. g3 `  K: H+ @" S32 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
9 @: e( M0 \6 I9 ?+ ]3 I11.5 13.25 14 13 12 11 13 13 12 12.25 13 15.5 15 14 13];
subplot(1,2,1);
plot(A1,B1,'.-');
hold on;
. q: F' [/ w% _p2=polyfit(A1,B1,2)
: C) e* o# I% g/ hf2=polyval(p2,A1);
2 I4 o3 C! E5 q1 X- n7 Q# mplot(A1,B1,'.',A1,f2,'-');
hold on;
5 N7 o1 S# S- F) e. W/ H1 Etitle('白内障病人情况');
' j# y; P. Z: |1 ^; `" Wxlabel('白内障病人门诊时间距7月12日的时间差');ylabel('白内障病人平均
等待时间');
0 k. d4 c7 q( H" }axis equal;
, J& P$ C$ z0 KA2=[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 ];
' |2 z0 R: Z2 R! pB2=[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
* U# f" h) S( t1 a12 12 13.25 14 14 13 12 12 12.4 13.33 13 13 13 12];
subplot(1,2,2);
8 p4 R3 F2 z# }- jplot(A2,B2,'.-');
0 C4 D: y2 Q, L) J3 [) Jhold on;
p3=polyfit(A2,B2,2)
* `/ J) M+ L  Y( Lf3=polyval(p3,A2);
plot(A2,B2,'.',A2,f3,'-');
  r& D' e8 n4 ]4 T* ]' m3 w% B- Chold on;
title('白内障双眼病人情况');
0 T0 x5 w$ Y7 |( Gxlabel('白内障双眼病人门诊时间距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 ];
: A) Q9 n* v1 _7 y6 DB3=[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
0 j- ~7 q$ v3 a1 u( A; a, B11.67 12 12.25 12 11 11.75 13 14 14 13 12 12 12.33 13 13 13 12.5 12 12];
subplot(2,2,1);
plot(A3,B3,'.-');
hold on;
  O: c" e0 J* z, b& z- tp4=polyfit(A3,B3,2)
2 O; o# t# R3 p: H4 b2 ]f4=polyval(p4,A3);
plot(A3,B3,'.',A3,f4,'-');
hold on;
title('视网膜病人情况');
xlabel('门诊时间距7 月12 日的时间差');ylabel('视网膜病人平均等待时间
');
- y/ q3 M* i% {2 `0 l  E9 f3 Laxis 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
4 |0 x( p7 ?3 E; ?6 i5 U13 12 12 13 13 12];
subplot(2,2,2);
plot(A4,B4,'.-');
8 J3 h/ V, Z' L" D: ^hold on;
5 q' X. ]& I, X2 z$ B# ^' B% Up5=polyfit(A4,B4,2)
f5=polyval(p5,A4);
plot(A4,B4,'.',A4,f5,'-');
! Q9 q8 Y' x* z' b# vhold on;
title('青光眼病人情况');
2 O( B' \7 k! H* P; t  d5 D, \8 Cxlabel('门诊时间距7 月12日的时间差');ylabel('视网膜病人平均等待时间
$ l1 V# D( R: i9 w7 Q');
axis equal;
% K1 }2 T$ E" c3 y' r, wA4=[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
% x; [- r5 B/ O" f7 b" ^8 A37];
% c6 Q  Y, T' M( f2 _* WB4=[12 13 12 12 12 12 13.5 13 13.5 14 13 12 12 11 10.25 12 11 11 15 15
& t4 F& B6 M; R13 12 12 13 13 12];
subplot(2,2,2);
plot(A4,B4,'.-');
hold on;
3 F# E, l: w% C: fp5=polyfit(A4,B4,2)
f5=polyval(p5,A4);
plot(A4,B4,'.',A4,f5,'-');
hold on;
! I+ h8 H% J# r0 v, b3 U8 F$ Htitle('青光眼病人情况');
" \0 d* Z9 H- W6 Y' x. Gxlabel('门诊时间距7 月12 日的时间差');ylabel('青光眼病人平均等待时间
');
$ g) b* P+ Y6 C% C3 g6 A' t  n: W2 maxis 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];
7 i2 e0 x+ l$ L/ H0 \$ _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];
% G4 P+ s  y6 C) j$ |: m, O2 Isubplot(2,2,3);
plot(A5,B5,'.-');
hold on;
( x- s+ I; Z3 p6 w; Tp6=polyfit(A5,B5,2)
f6=polyval(p6,A5);
plot(A5,B5,'.',A5,f6,'-');
hold on;
% I! T0 C( s3 Q: V2 T3 @# Y. ztitle('外伤病人情况');
xlabel('外伤病人门诊时间距离7月12日的时间差');ylabel('外伤病人平均等
待时间');
axis equal;
二、层次分析法的一致性指标判定的matlab 实现:
% m$ a% `% A. _; u7 c- m+ s, Nfunction H=ff(A)
t=eig(A)
w=max(t)
7 U& f) D0 E8 n, d. c0 G" S2 y" @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(
9)+A(10)+A(11)+A(12)+A(13)+A(14)+A(15)+A(16))
: u2 i- M- `3 ]* i4 C! L& I(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(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)
+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)
! x8 S4 R# e/ t2 W! T! F+A(10)+A(11)+A(12)+A(13)+A(14)+A(15)+A(16))]
7 {6 o; f* b( K" l5 N  r

linmatsas

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

迷途的羊

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

linmatsas

回复 迷途的羊 的帖子
' |  B6 z; b. ?5 d
是自己编的还用别人运行么……不过要是楼主自己编的……小弟就佩服死啦
   

紫辰

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

chenrye121

佩服两个字…………

水木年华zzu

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

wuming1429

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

hhao

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

abc321

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