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

wuming1429

一、使用曲线拟合模拟病人等待时间的模拟matlab程序实现：
[A]=xlsread('C:\MATLAB7\work\问题3数据5');
len=length(A);
' g: M" F0 S9 D( B. C: Fi=1:len;
* t" }3 p  b" C6 Xj=A(:,5);
figure(1);
plot(i,j,'.-');
hold on;
' \, }) C$ Z  j9 `2 ]) f3 @i=i';
p1=polyfit(i,j,2)
f1=polyval(p1,i);
plot(i,j,'.',i,f1,'-');
hold on;
4 j! Z% i3 k8 }& Y, Ztitle('病人平均等待时间与病人到达时间的关系');
9 M+ c! z* ~- L* b3 r5 @! _- Q. uxlabel('病人门诊时间距离7 月12 日的时间差');ylabel('病人平均等待时间
');
axis equal;
figure(2);
  h# T' |" z( \! ~. H7 t6 z" cA1=[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
, J8 K, \7 ^( k; f4 X32 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];
- K. A8 i# N( A, ?2 b* j7 Xsubplot(1,2,1);
; j5 u  O- A5 D0 e8 e4 ^plot(A1,B1,'.-');
1 s0 F9 o* ^1 s& @  phold on;
p2=polyfit(A1,B1,2)
f2=polyval(p2,A1);
1 U) |) Y2 `: U. @* o9 m. ^/ z6 \plot(A1,B1,'.',A1,f2,'-');
$ }/ G$ i/ L0 D0 x; Zhold on;
title('白内障病人情况');
xlabel('白内障病人门诊时间距7月12日的时间差');ylabel('白内障病人平均
  m5 d/ S% f6 ?& p, G" T4 [9 H等待时间');
, N* A! w& X5 P7 A5 p# J8 j. @axis 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
31 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
12 12 13.25 14 14 13 12 12 12.4 13.33 13 13 13 12];
subplot(1,2,2);
3 E0 F+ N) J4 Q7 ?- D% {plot(A2,B2,'.-');
hold on;
p3=polyfit(A2,B2,2)
8 ~( @! J/ [/ e$ V+ x& T1 Q8 yf3=polyval(p3,A2);
plot(A2,B2,'.',A2,f3,'-');
hold on;
title('白内障双眼病人情况');
xlabel('白内障双眼病人门诊时间距7月12日的时间差');ylabel('白内障双眼
病人平均等待时间');
6 o7 m  m- T/ n( l9 d% {7 U: faxis 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
$ v& G, ?% ~, B28 29 30 31 32 33 34 35 36 37 38 ];
6 P% D# A. ]4 U2 LB3=[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];
- }3 E8 D( [  l9 Zsubplot(2,2,1);
$ G$ c# M- D. R, F" W) Gplot(A3,B3,'.-');
hold on;
p4=polyfit(A3,B3,2)
f4=polyval(p4,A3);
plot(A3,B3,'.',A3,f4,'-');
. \4 H6 {- l0 G4 V7 X6 Ihold on;
title('视网膜病人情况');
xlabel('门诊时间距7 月12 日的时间差');ylabel('视网膜病人平均等待时间
');
/ I, N8 b6 B6 H3 Paxis 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];
% ^- M/ y- V" m+ G3 Q7 {, YB4=[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);
plot(A4,B4,'.-');
hold on;
! N0 ^- _- M6 y- ^  Zp5=polyfit(A4,B4,2)
* H9 _% [8 @! u' X( if5=polyval(p5,A4);
plot(A4,B4,'.',A4,f5,'-');
hold on;
6 F2 Y5 Z' }# U; V( `3 ptitle('青光眼病人情况');
xlabel('门诊时间距7 月12日的时间差');ylabel('视网膜病人平均等待时间
1 g; z" r8 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
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];
* ~1 _3 p3 c9 z+ Hsubplot(2,2,2);
3 n/ N$ r6 A/ N* j3 gplot(A4,B4,'.-');
hold on;
p5=polyfit(A4,B4,2)
  H4 D. F9 G5 F- ~2 @3 Z" Zf5=polyval(p5,A4);
plot(A4,B4,'.',A4,f5,'-');
: F3 \1 |# m( u. i6 uhold on;
title('青光眼病人情况');
xlabel('门诊时间距7 月12 日的时间差');ylabel('青光眼病人平均等待时间
');
7 y/ b% i6 F* a8 yaxis equal;
( O. D- g) T- l" K2 AA5=[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];
  G# i* S8 F  J9 _: ?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);
plot(A5,B5,'.-');
, F6 x2 g7 Y- M; `# m4 ohold on;
p6=polyfit(A5,B5,2)
f6=polyval(p6,A5);
7 J) i- r6 z  hplot(A5,B5,'.',A5,f6,'-');
" e9 f0 d0 ]8 L: Q/ Z1 A3 Nhold on;
title('外伤病人情况');
xlabel('外伤病人门诊时间距离7月12日的时间差');ylabel('外伤病人平均等
0 d$ C0 j+ {, p* c5 P待时间');
axis equal;
! j( \8 ]8 b( B( ^: j/ d二、层次分析法的一致性指标判定的matlab 实现:
/ Y2 ?! ]7 R4 Z. xfunction H=ff(A)
t=eig(A)
7 S$ i: z9 a( c- F* y1 xw=max(t)
, k- N* ]2 _. |5 j5 v4 b- l! D: T$ Cq=(w-4)/3/.9
6 i; h- {7 Y% g) H. u9 A+ A( q$ js=[(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))
(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)
0 g3 A* d& u1 B' }+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

回复 迷途的羊 的帖子
5 G! `- i' ^5 g7 L8 x' B2 f
是自己编的还用别人运行么……不过要是楼主自己编的……小弟就佩服死啦
! i: y/ g, z3 Z9 C( b8 f   

紫辰

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

chenrye121

佩服两个字…………

水木年华zzu

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

wuming1429

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

hhao

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

abc321

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