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

wuming1429

一、使用曲线拟合模拟病人等待时间的模拟matlab程序实现：
2 x$ q, C+ f7 T[A]=xlsread('C:\MATLAB7\work\问题3数据5');
len=length(A);
i=1:len;
' P5 A- l3 _) u4 gj=A(:,5);
figure(1);
plot(i,j,'.-');
& b- N8 k" A% m2 Phold on;
i=i';
  U( X4 T0 e' Z' z' K. O9 Mp1=polyfit(i,j,2)
f1=polyval(p1,i);
plot(i,j,'.',i,f1,'-');
8 T$ p2 B3 t& t* T8 S- [# q: Uhold on;
title('病人平均等待时间与病人到达时间的关系');
- ]2 q/ R+ P6 Lxlabel('病人门诊时间距离7 月12 日的时间差');ylabel('病人平均等待时间
');
axis equal;
) S: w" {1 U# ~figure(2);
, e3 U6 J1 k% d+ B) R* j+ q# C/ M5 GA1=[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
5 X% M: ?7 o3 _, `0 W32 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
# L. ^/ E; E9 b$ ]0 F11.5 13.25 14 13 12 11 13 13 12 12.25 13 15.5 15 14 13];
" L% N% Q; ]+ m' X, w' Y# |subplot(1,2,1);
( C: E) F* S9 U# v5 |8 ]& ^plot(A1,B1,'.-');
! |9 c! h' U5 M5 W( x. l" F& nhold on;
" J- Y9 x8 p) X! t9 [2 mp2=polyfit(A1,B1,2)
. ~- b2 I! p- a2 C/ l' df2=polyval(p2,A1);
9 D0 k  M% t2 ]6 w% Xplot(A1,B1,'.',A1,f2,'-');
hold on;
title('白内障病人情况');
xlabel('白内障病人门诊时间距7月12日的时间差');ylabel('白内障病人平均
等待时间');
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 ];
" R0 p: |5 P. U0 S; r' D+ bB2=[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);
plot(A2,B2,'.-');
hold on;
p3=polyfit(A2,B2,2)
f3=polyval(p3,A2);
2 a" s3 @' {# {plot(A2,B2,'.',A2,f3,'-');
hold on;
; M8 |" Y- H" s$ u$ ptitle('白内障双眼病人情况');
. t6 h0 N0 V4 k* F' ~xlabel('白内障双眼病人门诊时间距7月12日的时间差');ylabel('白内障双眼
& r' t- S' M, u% y' ]9 w$ t3 D病人平均等待时间');
axis equal;
figure(3);
- X6 ], P4 ^+ B) [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
& t1 T  \; [) S' K4 u. }28 29 30 31 32 33 34 35 36 37 38 ];
8 D( h/ y0 a( X7 F% eB3=[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
- Q# X, ^  {/ ?& {7 N# 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);
: G3 ~, B1 H1 W& O4 Nplot(A3,B3,'.-');
hold on;
p4=polyfit(A3,B3,2)
f4=polyval(p4,A3);
plot(A3,B3,'.',A3,f4,'-');
/ Y& m# v; m; ^$ shold on;
title('视网膜病人情况');
7 Z8 N8 Z# Y* k) y+ \xlabel('门诊时间距7 月12 日的时间差');ylabel('视网膜病人平均等待时间
) k4 U, L; t* A( ~6 l');
9 Q0 _2 C- j& c" M& f3 {6 g/ haxis 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
- M! h) P( D) W) D6 @* `13 12 12 13 13 12];
  R. C* J8 k, o# ^: ?  y3 Tsubplot(2,2,2);
plot(A4,B4,'.-');
hold on;
1 R" C: f% k6 A1 J6 yp5=polyfit(A4,B4,2)
f5=polyval(p5,A4);
plot(A4,B4,'.',A4,f5,'-');
hold on;
* Z4 D- E  \" F0 J& Q. Stitle('青光眼病人情况');
. @! w+ V0 i$ ~! `+ uxlabel('门诊时间距7 月12日的时间差');ylabel('视网膜病人平均等待时间
');
axis equal;
0 ?' F3 O- H1 ^) s0 x. |' e2 mA4=[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
0 w( u6 v" F$ Z% X7 L9 i- V37];
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];
, ?3 H+ v* j! ?0 ]% |subplot(2,2,2);
/ X, G. z2 @, D. c6 m/ F/ g6 ]1 z2 ?plot(A4,B4,'.-');
hold on;
p5=polyfit(A4,B4,2)
4 _$ }: J' w' e1 B4 [f5=polyval(p5,A4);
plot(A4,B4,'.',A4,f5,'-');
, U" m- ]% ~( E5 w5 P5 G1 K- thold on;
% g7 ^3 j7 o& Vtitle('青光眼病人情况');
( `5 O- A$ Q* e9 X6 w: Zxlabel('门诊时间距7 月12 日的时间差');ylabel('青光眼病人平均等待时间
');
7 r$ v& E+ Z# n( O2 E* Q7 faxis 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,'.-');
; R: R7 I5 ]9 E& x' f7 J3 Ghold on;
& \3 a) w5 w) x  d& bp6=polyfit(A5,B5,2)
4 b8 `) x6 {4 U; j; c) mf6=polyval(p6,A5);
  [5 N; m8 H5 L6 T( w: F+ o  Nplot(A5,B5,'.',A5,f6,'-');
hold on;
' A3 w3 {: R. H& a9 o/ r1 qtitle('外伤病人情况');
7 d/ S- d, f5 T* M1 V' wxlabel('外伤病人门诊时间距离7月12日的时间差');ylabel('外伤病人平均等
( b4 _  _' E' ?$ G: c" z, i待时间');
axis equal;
' l* _3 h: _1 P二、层次分析法的一致性指标判定的matlab 实现:
: g  e' n6 D; B2 T/ n: c5 tfunction H=ff(A)
t=eig(A)
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(
3 m$ p( C+ Z  z" t4 U. d" T, T! b- C9)+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))
- w! d/ l5 I$ s  f' i6 z(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))
7 a' ~. x7 G, f# C' z2 `5 ](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))]
6 M% C/ \& v, e6 Y0 S* ~3 t

linmatsas

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

迷途的羊

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

linmatsas

回复 迷途的羊 的帖子

是自己编的还用别人运行么……不过要是楼主自己编的……小弟就佩服死啦
$ X# w& s4 E3 n# g   

紫辰

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

chenrye121

佩服两个字…………

水木年华zzu

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

wuming1429

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

hhao

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

abc321

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