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

wuming1429

一、使用曲线拟合模拟病人等待时间的模拟matlab程序实现：
- O0 e$ z2 F  g; A/ Z3 L[A]=xlsread('C:\MATLAB7\work\问题3数据5');
7 N' i6 ]1 s' n5 U8 Llen=length(A);
5 L* S: P( V5 M6 W- q) `i=1:len;
9 e8 Z$ \/ x0 p4 d5 a1 V% Lj=A(:,5);
& _' e  B9 N, ?+ \7 Dfigure(1);
" ]& D7 n, A0 [; }: ~plot(i,j,'.-');
) t5 X: a- b" ], Bhold on;
) L, X1 C* E$ Y0 V/ Mi=i';
p1=polyfit(i,j,2)
9 b3 S  i/ S4 [- Lf1=polyval(p1,i);
plot(i,j,'.',i,f1,'-');
, g" Z- T* u3 B+ [/ M+ |. Whold on;
title('病人平均等待时间与病人到达时间的关系');
xlabel('病人门诊时间距离7 月12 日的时间差');ylabel('病人平均等待时间
');
$ ~( ?( V1 |7 b. x& _5 V* ]axis equal;
/ S* U+ ^' r4 m$ @figure(2);
3 Y* m4 N! ?  r' d; m1 G& D1 SA1=[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
9 P* t6 j  J0 h1 g9 S32 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];
+ ~  _/ L' I# P2 `+ \, isubplot(1,2,1);
plot(A1,B1,'.-');
hold on;
5 x: }+ d! @" I  Qp2=polyfit(A1,B1,2)
f2=polyval(p2,A1);
% z& D8 T, l) ^plot(A1,B1,'.',A1,f2,'-');
hold on;
title('白内障病人情况');
xlabel('白内障病人门诊时间距7月12日的时间差');ylabel('白内障病人平均
等待时间');
" b! F% O. {) j8 ~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
6 ^. N' B( X. y! T2 i31 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);
plot(A2,B2,'.-');
1 P2 Y0 M0 Q" v/ ^0 @5 @+ phold on;
p3=polyfit(A2,B2,2)
3 Q4 g" X4 y# \f3=polyval(p3,A2);
' i! `6 Q( b- l5 J8 cplot(A2,B2,'.',A2,f3,'-');
; V5 \. ?3 Y4 M- A' X: m' qhold on;
! ^# n1 J' r" M& Atitle('白内障双眼病人情况');
7 e: t2 F# b. P5 u' h# M% pxlabel('白内障双眼病人门诊时间距7月12日的时间差');ylabel('白内障双眼
$ `8 C$ n4 A. H2 o6 r, `, }5 W病人平均等待时间');
, F) \6 d: {# W# }) T9 o% |axis equal;
figure(3);
% y$ n* Y% }3 O# H2 U' l- IA3=[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
+ q2 G! F/ A6 L; H4 C2 N. J# z2 L28 29 30 31 32 33 34 35 36 37 38 ];
# p, A: S+ A/ ^/ MB3=[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
6 \' m4 [% T. g. k0 ^11.67 12 12.25 12 11 11.75 13 14 14 13 12 12 12.33 13 13 13 12.5 12 12];
; [* m/ p, w$ H4 X( |4 X7 t, E3 gsubplot(2,2,1);
plot(A3,B3,'.-');
hold on;
p4=polyfit(A3,B3,2)
& J- l( V) x" J* [& g2 Kf4=polyval(p4,A3);
8 b$ `6 W- ?4 k/ P5 ^( ]# L6 G# ^plot(A3,B3,'.',A3,f4,'-');
' ]2 a9 F: r# J! Y+ T# X- }4 Vhold on;
title('视网膜病人情况');
xlabel('门诊时间距7 月12 日的时间差');ylabel('视网膜病人平均等待时间
% t- _7 c; k, d6 z" b! u');
axis equal;
! B: y' [7 r4 Y5 q+ j' S/ v/ m' xA4=[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];
# {9 x  f" p  i6 N$ l7 gB4=[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;
p5=polyfit(A4,B4,2)
5 D$ v$ W1 a9 Ff5=polyval(p5,A4);
. f6 S: B" q" a0 w! G1 {2 D" iplot(A4,B4,'.',A4,f5,'-');
; @' U1 d3 s" U( A$ @1 K; o! Whold on;
0 d2 r, C; ^3 e3 K' H* @) G9 Htitle('青光眼病人情况');
6 y# ^* t0 T; |6 v( z* jxlabel('门诊时间距7 月12日的时间差');ylabel('视网膜病人平均等待时间
');
axis equal;
* h( j3 {2 Q% r# N, f1 [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];
' d7 I; t1 @1 \: q+ RB4=[12 13 12 12 12 12 13.5 13 13.5 14 13 12 12 11 10.25 12 11 11 15 15
& _1 `7 u) n1 r" j9 F13 12 12 13 13 12];
* P( J& L5 Y$ {& j: L- |1 r/ a" Fsubplot(2,2,2);
5 S6 n/ U9 Z: T: |% tplot(A4,B4,'.-');
hold on;
) V8 k5 f. @( g$ rp5=polyfit(A4,B4,2)
, S. {* H' k9 Z3 P, V, xf5=polyval(p5,A4);
plot(A4,B4,'.',A4,f5,'-');
hold on;
title('青光眼病人情况');
xlabel('门诊时间距7 月12 日的时间差');ylabel('青光眼病人平均等待时间
* J* u% b; [- j0 t7 E) `( @');
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);
& I* j9 x! T  l; c+ wplot(A5,B5,'.-');
hold on;
  P( [5 C! q: Z2 O* d2 s" b& C! ]p6=polyfit(A5,B5,2)
8 v, _  P9 g. s6 |/ l8 Y# Sf6=polyval(p6,A5);
plot(A5,B5,'.',A5,f6,'-');
) _! }$ b) [5 e- f8 p2 C( B5 W: {, uhold on;
title('外伤病人情况');
xlabel('外伤病人门诊时间距离7月12日的时间差');ylabel('外伤病人平均等
待时间');
axis equal;
" Q2 \; W! }8 |3 J# {% y二、层次分析法的一致性指标判定的matlab 实现:
" n( d: K8 l* i! y' Bfunction H=ff(A)
/ _( _2 y: w  @3 G# Kt=eig(A)
4 o0 w2 N/ T  T6 yw=max(t)
# W2 t1 y! B% U* Hq=(w-4)/3/.9
+ G% E6 e$ }# _+ w' F3 _3 u2 z" ]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(
4 N6 b; R! ^0 K1 f' v# x" C! Z/ J9)+A(10)+A(11)+A(12)+A(13)+A(14)+A(15)+A(16))
8 M0 X' R+ s# {7 u6 X(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)
! F# b$ E6 G# L8 n, Y$ q, @+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)
+A(10)+A(11)+A(12)+A(13)+A(14)+A(15)+A(16))]
$ ?' R6 t5 ^. [

linmatsas

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

迷途的羊

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

linmatsas

回复 迷途的羊 的帖子
2 s5 q. F) |; S" Z' q
8 ^1 }- b4 `/ t是自己编的还用别人运行么……不过要是楼主自己编的……小弟就佩服死啦
* {. Q4 l6 x& w! N. G3 }0 ]1 m   

紫辰

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

chenrye121

佩服两个字…………

水木年华zzu

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

wuming1429

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

hhao

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

abc321

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