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

wuming1429

一、使用曲线拟合模拟病人等待时间的模拟matlab程序实现：
[A]=xlsread('C:\MATLAB7\work\问题3数据5');
len=length(A);
i=1:len;
/ I2 D. a; ?$ c; b; A" q6 l1 t1 Sj=A(:,5);
1 j3 q" d& Z0 f& q7 A$ ?" B$ gfigure(1);
( g( Y7 J) n) W% @plot(i,j,'.-');
hold on;
i=i';
( Z  m2 n4 I6 y2 q; lp1=polyfit(i,j,2)
f1=polyval(p1,i);
* a5 i* W6 r& _% c- Pplot(i,j,'.',i,f1,'-');
hold on;
- e3 D- M$ w$ w/ Gtitle('病人平均等待时间与病人到达时间的关系');
xlabel('病人门诊时间距离7 月12 日的时间差');ylabel('病人平均等待时间
7 q! i& j0 D* \" [$ V! s1 |- [; c7 k. }');
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
( y# |8 M2 j" z6 e1 ]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];
8 D2 x5 _* b0 z, E6 w7 b& Jsubplot(1,2,1);
plot(A1,B1,'.-');
hold on;
p2=polyfit(A1,B1,2)
  d8 F/ u$ I+ d, f8 C- X6 r2 \f2=polyval(p2,A1);
( s. E* i( g: Y1 }! Uplot(A1,B1,'.',A1,f2,'-');
( N7 x. P; i( m5 s, Z5 R* M- G  Q9 nhold on;
7 a0 K7 l$ Y9 ^# `$ z  G2 ztitle('白内障病人情况');
xlabel('白内障病人门诊时间距7月12日的时间差');ylabel('白内障病人平均
等待时间');
axis equal;
! z: i! s# `' D0 V' D8 ?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
: U* L1 C& P& N% [+ Q31 32 33 34 35 36 37 ];
# N& ~! o' L! H7 o" z) \2 YB2=[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
& e% H* a* Q0 F/ K/ u1 f2 H12 12 13.25 14 14 13 12 12 12.4 13.33 13 13 13 12];
subplot(1,2,2);
, t: I! V- |3 Gplot(A2,B2,'.-');
4 C4 Q$ F2 z  p1 O5 _hold on;
p3=polyfit(A2,B2,2)
8 ~  T3 ?/ `( j! h/ h+ H6 vf3=polyval(p3,A2);
# a2 i2 l. [9 y6 [3 J, `plot(A2,B2,'.',A2,f3,'-');
. _2 X  a$ R5 p$ T; ^hold on;
title('白内障双眼病人情况');
1 r8 J9 F& {* u  {9 Oxlabel('白内障双眼病人门诊时间距7月12日的时间差');ylabel('白内障双眼
0 M0 z8 w2 m# G0 z病人平均等待时间');
- O6 d6 \, R! t( Caxis 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
' K  H( f$ u$ Y1 B: o; d; V6 M11.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);
' z0 @! x" b8 ^: dplot(A3,B3,'.-');
hold on;
p4=polyfit(A3,B3,2)
f4=polyval(p4,A3);
4 I& Z. w- h- F4 ?9 x1 e& Pplot(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];
% u& l0 c4 k+ lB4=[12 13 12 12 12 12 13.5 13 13.5 14 13 12 12 11 10.25 12 11 11 15 15
5 w% a- c* w1 f" X! ^13 12 12 13 13 12];
8 b( a6 U; b5 t1 D6 ^, csubplot(2,2,2);
plot(A4,B4,'.-');
" Y( V" K( m2 s3 M: R0 Nhold on;
p5=polyfit(A4,B4,2)
; q4 z7 i. z( {3 Q" |# Ff5=polyval(p5,A4);
$ u. M) s6 v, N  X4 mplot(A4,B4,'.',A4,f5,'-');
' k) i% i) c  A& j7 y$ g3 ^3 E  uhold on;
title('青光眼病人情况');
- k6 l0 _' {, E! k" wxlabel('门诊时间距7 月12日的时间差');ylabel('视网膜病人平均等待时间
');
. K6 |6 @7 \# W0 R' L7 _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
4 h: s/ v+ n/ M1 A37];
( v4 c; E* a) N( h- a2 q2 Z8 U9 OB4=[12 13 12 12 12 12 13.5 13 13.5 14 13 12 12 11 10.25 12 11 11 15 15
9 }: Z# o) U1 p6 [5 R' k13 12 12 13 13 12];
subplot(2,2,2);
% `# m0 t/ A4 K5 C* L9 G2 xplot(A4,B4,'.-');
% _1 [: n0 @' W" jhold on;
; Y! a8 n! ]1 g! ^7 tp5=polyfit(A4,B4,2)
) Q8 X  g; ]: A# ?f5=polyval(p5,A4);
plot(A4,B4,'.',A4,f5,'-');
. d8 D& t. r8 s! f$ khold on;
8 a) o$ E9 g' z6 F" W  r" c8 `title('青光眼病人情况');
- z9 c/ L: k: Oxlabel('门诊时间距7 月12 日的时间差');ylabel('青光眼病人平均等待时间
');
% g5 z) K. u! {+ ~1 P3 ]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];
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,'.-');
' w+ K  \6 W: x) mhold on;
; H; ]' d3 a' A, {. @8 Cp6=polyfit(A5,B5,2)
f6=polyval(p6,A5);
4 ]+ w1 x3 G2 oplot(A5,B5,'.',A5,f6,'-');
hold on;
) K% ]  Z9 N2 [, E5 R7 Rtitle('外伤病人情况');
$ b) o; b; X% T: F% dxlabel('外伤病人门诊时间距离7月12日的时间差');ylabel('外伤病人平均等
待时间');
axis equal;
$ r0 l( W" S/ N0 h# [+ M& X9 `二、层次分析法的一致性指标判定的matlab 实现:
9 r2 Q% Z* ~, r/ Q6 E; c8 wfunction H=ff(A)
+ _. G9 W" b/ d" o: rt=eig(A)
- J5 Y# Z8 f3 M& ]; y( a6 @w=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(
0 z8 P1 l. s' d+ p9)+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)
; {# N$ }4 r3 O0 {# ~+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))]
  K9 L5 D) F5 o; M: A

linmatsas

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

迷途的羊

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

linmatsas

回复 迷途的羊 的帖子
: Y6 _6 g! u& v0 }' p$ ~
是自己编的还用别人运行么……不过要是楼主自己编的……小弟就佩服死啦
4 ?* z. W' P, Q' L) @( g   

紫辰

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

chenrye121

佩服两个字…………

水木年华zzu

[A]=xlsread('C:\MATLAB7\work\问题3数据5');
( h; y" H- W, O0 H. f6 e! O在不同的地方运行上面一句要改成本地数据保存的路径的

wuming1429

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

hhao

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

abc321

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