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

wuming1429

一、使用曲线拟合模拟病人等待时间的模拟matlab程序实现：
[A]=xlsread('C:\MATLAB7\work\问题3数据5');
" q* X4 u, B1 o; M* ~len=length(A);
i=1:len;
# @( G2 O( w7 I5 i- r5 g) s* d& G" ej=A(:,5);
3 o# b' A" S1 P- A+ X+ ]; K! _$ a/ `figure(1);
plot(i,j,'.-');
hold on;
3 k( E! p& }" n7 v! Z3 pi=i';
p1=polyfit(i,j,2)
, b6 d) Y: K) w$ Y  d" h* X5 K( jf1=polyval(p1,i);
plot(i,j,'.',i,f1,'-');
hold on;
title('病人平均等待时间与病人到达时间的关系');
xlabel('病人门诊时间距离7 月12 日的时间差');ylabel('病人平均等待时间
" r! G% M( Q# N');
) n) P$ i$ u- z% G/ M5 b3 V8 kaxis equal;
& m" d/ }8 g. a, d: H0 _9 pfigure(2);
& B* k4 S" K8 ^& d) F6 T. @) N& jA1=[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
11.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;
p2=polyfit(A1,B1,2)
0 D3 D# v0 g6 h5 \3 f* cf2=polyval(p2,A1);
plot(A1,B1,'.',A1,f2,'-');
6 Y  z8 ], \: h) b2 E4 B, Zhold on;
title('白内障病人情况');
& z+ L8 a4 m) {' t0 axlabel('白内障病人门诊时间距7月12日的时间差');ylabel('白内障病人平均
" N5 u1 Z, G0 e# \) q& z- K等待时间');
) p' d/ R6 }9 Laxis 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
5 x" ]$ |% e+ Z- f. b. j9 R31 32 33 34 35 36 37 ];
7 b6 c1 F! h7 _" h, K. q' qB2=[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
* o" n/ X5 h& o( m12 12 13.25 14 14 13 12 12 12.4 13.33 13 13 13 12];
subplot(1,2,2);
* d( r4 k7 N3 H: ~+ Lplot(A2,B2,'.-');
# N) Y, r' J* o% w) o" e, shold on;
p3=polyfit(A2,B2,2)
f3=polyval(p3,A2);
# r! Y! P9 n7 A4 qplot(A2,B2,'.',A2,f3,'-');
" w1 j& @7 V$ a% k, v, zhold on;
title('白内障双眼病人情况');
+ }5 d! }8 A5 `$ I% m# x: rxlabel('白内障双眼病人门诊时间距7月12日的时间差');ylabel('白内障双眼
病人平均等待时间');
( v& A+ p; I" D) X; G) s; {axis equal;
$ P  E$ g7 B8 Q' [! a& D9 ]# }figure(3);
# o- z7 U/ k, T6 WA3=[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
1 z9 R  b+ H. h, C9 H9 O; l, K28 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
% Y. N0 \' ]8 |3 U: l1 @11.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,'.-');
5 N6 O5 Q! t! ahold on;
p4=polyfit(A3,B3,2)
f4=polyval(p4,A3);
% L7 x, ~; M- `4 c) Q' c- Hplot(A3,B3,'.',A3,f4,'-');
hold on;
title('视网膜病人情况');
xlabel('门诊时间距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];
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];
/ J  e' C# z. Z. e  Ksubplot(2,2,2);
& \: ^+ T0 E& U: x+ Wplot(A4,B4,'.-');
! Z# G4 X3 _3 _/ ^$ dhold on;
) y: R: z4 y" L9 xp5=polyfit(A4,B4,2)
f5=polyval(p5,A4);
& U- ?9 M! a/ J' n! E5 Cplot(A4,B4,'.',A4,f5,'-');
6 P6 F: ]* K1 F* C+ C" {  ]' }hold on;
# o) d* ^* {4 P) O7 vtitle('青光眼病人情况');
xlabel('门诊时间距7 月12日的时间差');ylabel('视网膜病人平均等待时间
');
/ a- g9 W! k, \  s# y. saxis 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];
subplot(2,2,2);
$ o% y, N) R; `1 yplot(A4,B4,'.-');
hold on;
( Z' c7 y/ o1 G5 mp5=polyfit(A4,B4,2)
# J" M. Q- p1 O9 r$ \f5=polyval(p5,A4);
plot(A4,B4,'.',A4,f5,'-');
1 a4 |- Y6 w- khold on;
title('青光眼病人情况');
xlabel('门诊时间距7 月12 日的时间差');ylabel('青光眼病人平均等待时间
' n# _2 P) P- G" Q. ]');
axis equal;
  S) f# a7 d/ w  A& f) r& D  j" oA5=[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];
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];
# G( r8 N& b! _( Y1 U( Lsubplot(2,2,3);
plot(A5,B5,'.-');
hold on;
p6=polyfit(A5,B5,2)
# I+ d; f8 i3 f0 f( Lf6=polyval(p6,A5);
7 W6 g. p9 P" T( G. p' {plot(A5,B5,'.',A5,f6,'-');
hold on;
title('外伤病人情况');
xlabel('外伤病人门诊时间距离7月12日的时间差');ylabel('外伤病人平均等
待时间');
axis equal;
4 A2 |: `' \" P8 R1 P# i二、层次分析法的一致性指标判定的matlab 实现:
& Q: P9 v, [: q; O! t% ~+ ffunction H=ff(A)
t=eig(A)
w=max(t)
* o0 S2 P' C7 W( Z0 w& 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(
( D! r) g8 S0 t4 y9)+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)
" U1 [' ]; M) F+ _+A(10)+A(11)+A(12)+A(13)+A(14)+A(15)+A(16))
& b1 m, k9 Z: X5 n- A(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)
+A(10)+A(11)+A(12)+A(13)+A(14)+A(15)+A(16))]

linmatsas

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

迷途的羊

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

linmatsas

回复 迷途的羊 的帖子

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

紫辰

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

chenrye121

佩服两个字…………

水木年华zzu

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

wuming1429

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

hhao

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

abc321

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