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

wuming1429

一、使用曲线拟合模拟病人等待时间的模拟matlab程序实现：
1 @7 b- ?/ g) {6 i2 f9 d[A]=xlsread('C:\MATLAB7\work\问题3数据5');
: S' _0 E& f4 N8 G& Llen=length(A);
, j  \, J. H0 v% x  }' ii=1:len;
j=A(:,5);
figure(1);
8 D# n4 }& A; {4 u, C% rplot(i,j,'.-');
hold on;
i=i';
p1=polyfit(i,j,2)
0 Y" j; l! a* t& Wf1=polyval(p1,i);
plot(i,j,'.',i,f1,'-');
& O& l) ?$ X- E; @& [hold on;
4 j0 D7 |1 q& d: J5 i' t+ k$ R# utitle('病人平均等待时间与病人到达时间的关系');
4 s7 P6 L! w$ B# ]; Xxlabel('病人门诊时间距离7 月12 日的时间差');ylabel('病人平均等待时间
  Q3 ^: {0 @8 Q4 a4 E');
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
# f0 r9 K5 R7 f2 r4 N! y- h+ E) j32 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
11.5 13.25 14 13 12 11 13 13 12 12.25 13 15.5 15 14 13];
- e" z  b1 e7 S# z/ i8 Zsubplot(1,2,1);
. S' O- g* O; j0 `plot(A1,B1,'.-');
hold on;
! a' v9 K6 k) J/ C0 [, U! I7 np2=polyfit(A1,B1,2)
f2=polyval(p2,A1);
plot(A1,B1,'.',A1,f2,'-');
) i0 Q% Y' [1 a  C& R  Dhold on;
title('白内障病人情况');
3 d2 X0 h4 `% R" z, V; oxlabel('白内障病人门诊时间距7月12日的时间差');ylabel('白内障病人平均
等待时间');
axis equal;
5 g; w  r0 _3 _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 ];
4 `7 i! E& u  U6 H# Y. A0 G% U* `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
. \. C3 @9 ~$ x12 12 13.25 14 14 13 12 12 12.4 13.33 13 13 13 12];
subplot(1,2,2);
( o: C+ o4 C/ Q3 f- C2 Lplot(A2,B2,'.-');
# F2 G  r( j; s6 hhold on;
p3=polyfit(A2,B2,2)
# \$ b9 V9 m& m; W; C3 k& rf3=polyval(p3,A2);
plot(A2,B2,'.',A2,f3,'-');
hold on;
title('白内障双眼病人情况');
1 n; R# C* j4 o' r$ M# gxlabel('白内障双眼病人门诊时间距7月12日的时间差');ylabel('白内障双眼
9 b( Z, G3 d' ~# z病人平均等待时间');
axis equal;
figure(3);
+ V; T  z! T$ E- G* x1 OA3=[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 ];
& z5 s% `! s1 Z' 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
1 w& @) A$ q8 W; a11.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);
5 Z0 j* Z2 u+ l7 r8 T7 C! X. e. rplot(A3,B3,'.-');
; ?) t$ z# q7 I2 f& J5 shold on;
p4=polyfit(A3,B3,2)
9 D9 z  r8 g' o6 S' f9 ff4=polyval(p4,A3);
  h4 R  L7 P( vplot(A3,B3,'.',A3,f4,'-');
/ Y' {* f# M) H" N* Q' Phold on;
, |# v* e: N4 r: b+ {5 [title('视网膜病人情况');
xlabel('门诊时间距7 月12 日的时间差');ylabel('视网膜病人平均等待时间
) P$ q& ]0 _" c1 o  ^& ~' E' ~: E, J');
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
/ T( j3 r! |2 }: i6 x* H. M37];
' U) C! T: R1 b' w) ]  B5 h6 KB4=[12 13 12 12 12 12 13.5 13 13.5 14 13 12 12 11 10.25 12 11 11 15 15
  s; {% R, c6 o6 U( T# _8 p13 12 12 13 13 12];
: Z" J8 H6 F+ |: z8 m4 Hsubplot(2,2,2);
plot(A4,B4,'.-');
+ r, M  D( K, z# _3 ~& u( ^' ^hold on;
. v" Z6 O* P% V5 M$ B) [4 `p5=polyfit(A4,B4,2)
f5=polyval(p5,A4);
  H/ E' |! r% H: N3 T, {: h: }plot(A4,B4,'.',A4,f5,'-');
' E. M2 y9 o5 \/ p+ q6 z% [/ e/ Dhold on;
title('青光眼病人情况');
' b5 L3 F2 G- N/ Ixlabel('门诊时间距7 月12日的时间差');ylabel('视网膜病人平均等待时间
/ u  a; j; r8 g- H5 k');
0 N3 f2 o$ V3 M5 o% Y, \axis equal;
- Y7 `! [* Z' N/ x% l7 P  dA4=[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
" T$ ^) @% L% x9 Q37];
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];
0 O' b$ x' C; Q" y. D9 jsubplot(2,2,2);
( e* y' [1 L; R  n2 `. u) mplot(A4,B4,'.-');
hold on;
  x/ c6 s! ^3 Vp5=polyfit(A4,B4,2)
0 n$ O# A1 U5 f6 x- d& lf5=polyval(p5,A4);
, t( U, a1 G& E2 qplot(A4,B4,'.',A4,f5,'-');
hold on;
title('青光眼病人情况');
+ ^" ?$ s, o: k3 B1 t8 Bxlabel('门诊时间距7 月12 日的时间差');ylabel('青光眼病人平均等待时间
');
axis equal;
4 e3 B- V# ]6 B7 c4 ~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];
" ~! _' J. P# K+ _- s1 ?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];
subplot(2,2,3);
# y1 y$ v3 h1 l" f8 [plot(A5,B5,'.-');
2 h/ H% e" C7 s6 t4 |8 x% l4 n  Zhold on;
" M* i* u. z; w4 I$ F: Cp6=polyfit(A5,B5,2)
& M$ P8 k7 p; {6 C6 f2 Mf6=polyval(p6,A5);
2 G% U8 p4 G$ g; f: Y% o  ^* pplot(A5,B5,'.',A5,f6,'-');
1 C5 U& M1 M4 h7 r" B* Ohold on;
title('外伤病人情况');
xlabel('外伤病人门诊时间距离7月12日的时间差');ylabel('外伤病人平均等
待时间');
axis equal;
二、层次分析法的一致性指标判定的matlab 实现:
function H=ff(A)
, b4 O" u  T5 o& Ct=eig(A)
' i( s) @7 a" J# ~+ tw=max(t)
q=(w-4)/3/.9
4 j) g. a3 T% T% Rs=[(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(
8 s6 X) U2 A6 F& 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)
+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)
9 q# J0 N- I3 ?! y4 ~( y+A(10)+A(11)+A(12)+A(13)+A(14)+A(15)+A(16))]
3 z) F0 @3 P& U6 q

linmatsas

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

迷途的羊

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

linmatsas

回复 迷途的羊 的帖子

是自己编的还用别人运行么……不过要是楼主自己编的……小弟就佩服死啦
   

紫辰

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

chenrye121

佩服两个字…………

水木年华zzu

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

wuming1429

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

hhao

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

abc321

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