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

wuming1429

一、使用曲线拟合模拟病人等待时间的模拟matlab程序实现：
, {: @5 i6 D, [- t[A]=xlsread('C:\MATLAB7\work\问题3数据5');
len=length(A);
i=1:len;
1 m6 w# J" l: Q/ Ej=A(:,5);
: j5 E# d0 P8 A* h+ Z6 \figure(1);
$ C8 i; B+ K& r  X; \  gplot(i,j,'.-');
hold on;
i=i';
p1=polyfit(i,j,2)
% y$ D6 |( W! I! gf1=polyval(p1,i);
$ y& ~- R1 Z6 t4 Y/ n2 D& Gplot(i,j,'.',i,f1,'-');
% O1 j7 R: ~& l. z0 [hold on;
! Z2 D# }8 H) o' Q. ?+ ptitle('病人平均等待时间与病人到达时间的关系');
xlabel('病人门诊时间距离7 月12 日的时间差');ylabel('病人平均等待时间
');
7 l1 q! ^6 p5 Z, E, `. Eaxis equal;
figure(2);
: I. j8 j( J3 M1 BA1=[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
/ w0 n6 Z* t" m( ~& A32 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);
6 m* f; K/ k8 Y2 W$ [/ T9 Bplot(A1,B1,'.-');
hold on;
p2=polyfit(A1,B1,2)
) g: S5 {. `/ ]f2=polyval(p2,A1);
7 L( `7 t- q, U' iplot(A1,B1,'.',A1,f2,'-');
/ A% x0 l: g; b+ Jhold on;
title('白内障病人情况');
xlabel('白内障病人门诊时间距7月12日的时间差');ylabel('白内障病人平均
等待时间');
. W+ y# k4 i, z1 J8 Maxis 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
) z. e* [% t* i8 X5 n! j( y31 32 33 34 35 36 37 ];
, U4 c6 x' j0 }" aB2=[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];
! I; D- Q+ m! A$ ssubplot(1,2,2);
plot(A2,B2,'.-');
hold on;
p3=polyfit(A2,B2,2)
; g2 r' n2 K$ o7 z! w& lf3=polyval(p3,A2);
plot(A2,B2,'.',A2,f3,'-');
hold on;
title('白内障双眼病人情况');
1 l3 V6 G" p( w- m, nxlabel('白内障双眼病人门诊时间距7月12日的时间差');ylabel('白内障双眼
病人平均等待时间');
' z; B$ Y7 m! T  g# H; M* Q6 p$ Z( Baxis 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 ];
$ p' s4 a6 v/ C  jB3=[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];
subplot(2,2,1);
plot(A3,B3,'.-');
- F  O) p; Q# ]8 c& }hold on;
; \( J! @, k0 y' j9 ~! vp4=polyfit(A3,B3,2)
' ^, ]2 ~- r! M; I0 bf4=polyval(p4,A3);
) Y$ _. V! O3 f! c! [4 Zplot(A3,B3,'.',A3,f4,'-');
hold on;
title('视网膜病人情况');
. e: u' w9 K3 z7 ^xlabel('门诊时间距7 月12 日的时间差');ylabel('视网膜病人平均等待时间
' G& j- U- P2 n% s2 Q' t');
axis equal;
' Y3 g+ i1 j" a9 q! Y/ S* tA4=[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) F3 y1 `# V# ]* I4 W! OB4=[12 13 12 12 12 12 13.5 13 13.5 14 13 12 12 11 10.25 12 11 11 15 15
  [+ R$ _3 n1 \+ S( Y13 12 12 13 13 12];
7 _6 ~1 d$ s9 b) P2 ssubplot(2,2,2);
plot(A4,B4,'.-');
hold on;
; o! d3 o2 M  c1 P; ~# G3 Bp5=polyfit(A4,B4,2)
8 Z; @) W/ {9 ]6 k1 \5 qf5=polyval(p5,A4);
plot(A4,B4,'.',A4,f5,'-');
hold on;
+ L+ C) H) {4 h+ ?' \title('青光眼病人情况');
xlabel('门诊时间距7 月12日的时间差');ylabel('视网膜病人平均等待时间
');
axis equal;
; D7 H5 [4 Q3 p' a# VA4=[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
) e( i6 F  m0 K) X1 D9 p; q% U5 [1 N' ~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
1 J+ W; Y; {/ R1 k- j! C13 12 12 13 13 12];
subplot(2,2,2);
plot(A4,B4,'.-');
8 ]+ W4 O, S! }! g9 }: ?hold on;
" W. i! u1 Z! {+ {p5=polyfit(A4,B4,2)
f5=polyval(p5,A4);
5 W0 @# e# J. z9 k/ {$ Lplot(A4,B4,'.',A4,f5,'-');
# O1 ^: R% \9 hhold on;
1 s& Q+ F4 K1 \+ Q2 ^3 d% Etitle('青光眼病人情况');
: G/ k1 d4 [4 A; N& Kxlabel('门诊时间距7 月12 日的时间差');ylabel('青光眼病人平均等待时间
5 n. i0 t" V$ W7 _; N. [2 ~0 J! q');
4 ~/ p  f, q# N) W8 b/ f: p; _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];
+ V9 B% Q+ A3 {0 L6 _: ~6 @" d" ~subplot(2,2,3);
plot(A5,B5,'.-');
hold on;
p6=polyfit(A5,B5,2)
f6=polyval(p6,A5);
( i8 i# q$ ~5 x: r4 s3 t: {plot(A5,B5,'.',A5,f6,'-');
hold on;
title('外伤病人情况');
( d! }5 F8 `; a) k  t& T! T4 N; Axlabel('外伤病人门诊时间距离7月12日的时间差');ylabel('外伤病人平均等
4 n8 Q, ?4 H! }: j7 ]待时间');
* s; e0 Y. N. `! Paxis equal;
二、层次分析法的一致性指标判定的matlab 实现:
% w2 Z! N% X. ]! V2 @% Gfunction H=ff(A)
6 v( }3 V5 _! k! n! Yt=eig(A)
/ O6 [" K: V, O& i/ r- [' c+ Vw=max(t)
4 g4 O4 `$ ?: q+ w& pq=(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(
- ]/ q- J( }- p0 c1 U' d- y4 f! N9)+A(10)+A(11)+A(12)+A(13)+A(14)+A(15)+A(16))
: A3 Q$ }+ n  o( d# u(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)
& d2 m1 `. H! L7 t- ?% q' }, x+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)
  T( y5 C. B4 l( e  j  h! p7 g+A(10)+A(11)+A(12)+A(13)+A(14)+A(15)+A(16))]

linmatsas

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

迷途的羊

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

linmatsas

回复 迷途的羊 的帖子

) V1 `0 K5 B2 s* }4 {/ @+ F是自己编的还用别人运行么……不过要是楼主自己编的……小弟就佩服死啦
6 M8 m" _% g$ _7 r" P0 o+ p8 `   

紫辰

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

chenrye121

佩服两个字…………

水木年华zzu

[A]=xlsread('C:\MATLAB7\work\问题3数据5');
  z' \& w$ i  T. I" G& L* b在不同的地方运行上面一句要改成本地数据保存的路径的

wuming1429

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

hhao

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

abc321

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