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

wuming1429

一、使用曲线拟合模拟病人等待时间的模拟matlab程序实现：
  N0 p7 @( N& J: X[A]=xlsread('C:\MATLAB7\work\问题3数据5');
7 Z0 `( u( |: H: f6 nlen=length(A);
7 s/ [% k2 z: U6 K. @/ e2 qi=1:len;
$ ]. g! m9 L" J7 u6 k! n  aj=A(:,5);
  A5 {- s" {, @figure(1);
plot(i,j,'.-');
, L7 k, d; y0 P6 l4 m* |' Bhold on;
" w5 p8 w1 K* z1 fi=i';
p1=polyfit(i,j,2)
7 j5 b. W5 Q* W. q0 `8 i: x% ff1=polyval(p1,i);
- ]  r$ ~) m! [- R0 Y6 u4 Uplot(i,j,'.',i,f1,'-');
hold on;
title('病人平均等待时间与病人到达时间的关系');
xlabel('病人门诊时间距离7 月12 日的时间差');ylabel('病人平均等待时间
" R! \& c  ?# X4 z1 ]' F');
* T& {. E' w6 Taxis equal;
+ Q9 v3 \$ V! Y1 v2 A. E& I% {/ Z$ kfigure(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
32 34 35 36 37 38 39 40 41 42];
7 A( b6 }$ D2 I* f+ O; PB1=[12 12 12 12 12 12.5 13 13 13 13 13.4 14 14 13 12.5 11 11 12 12 12 11
1 g& s' n$ h& w/ A( Q. E; E. I& L11.5 13.25 14 13 12 11 13 13 12 12.25 13 15.5 15 14 13];
. M9 P0 X$ ~  ?, o. S9 ]/ Asubplot(1,2,1);
plot(A1,B1,'.-');
hold on;
0 H& I, Z+ M+ R" C. k5 tp2=polyfit(A1,B1,2)
f2=polyval(p2,A1);
plot(A1,B1,'.',A1,f2,'-');
hold on;
/ E6 Q" f  r  Z' wtitle('白内障病人情况');
" A& E9 e( S! M8 d0 P2 L" V$ ]4 zxlabel('白内障病人门诊时间距7月12日的时间差');ylabel('白内障病人平均
2 t" B0 {3 z& W5 v1 @/ z等待时间');
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];
6 Q' l+ I* F3 Wsubplot(1,2,2);
3 L" G: Y3 S. h+ J7 D" M( {plot(A2,B2,'.-');
1 O' z9 W2 m7 R% _; Zhold on;
p3=polyfit(A2,B2,2)
  X6 U7 w# e& d$ `f3=polyval(p3,A2);
plot(A2,B2,'.',A2,f3,'-');
; F! P, N* F" C2 m7 vhold on;
title('白内障双眼病人情况');
xlabel('白内障双眼病人门诊时间距7月12日的时间差');ylabel('白内障双眼
病人平均等待时间');
( S8 w7 S8 ^( |1 n7 K* Caxis equal;
  \4 s* b! `2 T* q( _. p+ Yfigure(3);
/ }: i0 \  i! Y& l5 f8 k  t$ K# N: K, C2 bA3=[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
% w+ F& x! O; W$ t6 x& W8 g3 E28 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
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 {# _; i3 J7 vsubplot(2,2,1);
plot(A3,B3,'.-');
: X- v9 K- `/ ^& U" j/ thold on;
p4=polyfit(A3,B3,2)
7 P" n. t3 t8 K& [" Gf4=polyval(p4,A3);
4 A6 B. p8 r3 `5 J' m: Aplot(A3,B3,'.',A3,f4,'-');
hold on;
title('视网膜病人情况');
# M2 ~& {/ M( m4 bxlabel('门诊时间距7 月12 日的时间差');ylabel('视网膜病人平均等待时间
+ @. \+ a2 {- b" A* o');
  i* x/ I( U/ P/ ?; u2 b% e- e" yaxis 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];
subplot(2,2,2);
5 m/ j3 |* i' c- H+ X% @( K9 `plot(A4,B4,'.-');
) F4 Q4 T( x$ C: Chold on;
3 U$ x/ r' ^2 k! t# d' ^9 b$ yp5=polyfit(A4,B4,2)
f5=polyval(p5,A4);
1 O5 ?) Z3 S0 Jplot(A4,B4,'.',A4,f5,'-');
hold on;
' c; ?! t2 I  h8 y) Gtitle('青光眼病人情况');
xlabel('门诊时间距7 月12日的时间差');ylabel('视网膜病人平均等待时间
8 o" B/ G2 P  Y5 f6 j7 {2 y');
, v2 Z8 ?1 M% w: }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];
6 v1 q# W( t; G* f/ PB4=[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;
5 a* Y: d- G3 u5 Qp5=polyfit(A4,B4,2)
f5=polyval(p5,A4);
plot(A4,B4,'.',A4,f5,'-');
2 [: {: z% H" `( ?9 Z( _hold on;
, X) y" H$ y+ S8 Ltitle('青光眼病人情况');
xlabel('门诊时间距7 月12 日的时间差');ylabel('青光眼病人平均等待时间
. C) t: C& t. O; D+ ~! C1 Q');
, a, D" q6 o9 ]( Q7 E7 P: y" `8 l- k! I- ^axis equal;
8 c2 @2 M( R9 n6 a/ HA5=[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];
, C* R  M0 x) DB5=[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];
, p: b: u& X2 Q" p" E3 Esubplot(2,2,3);
: c4 g0 G( A( a4 o/ ]+ L  Bplot(A5,B5,'.-');
hold on;
7 D; _8 G) W! P4 }  f1 K( Xp6=polyfit(A5,B5,2)
f6=polyval(p6,A5);
plot(A5,B5,'.',A5,f6,'-');
hold on;
, ~. n: t* J% Y7 e  P/ @/ }title('外伤病人情况');
xlabel('外伤病人门诊时间距离7月12日的时间差');ylabel('外伤病人平均等
  p  C* r7 }& J0 \( s& F# d待时间');
& W* o) i% n& D/ s  h6 u8 laxis equal;
1 V& f8 J, j% U/ t) j# u6 j2 L二、层次分析法的一致性指标判定的matlab 实现:
; c& m% H0 g% y0 @! pfunction H=ff(A)
t=eig(A)
, n: B& n# G! W+ c4 E- M+ Bw=max(t)
+ X8 D4 h8 i- Y& x" Fq=(w-4)/3/.9
& c" J. C" ^. r- t& Z& n  gs=[(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))
  d$ V1 h* b8 Q6 Z! S+ @1 C# |$ E(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))
4 g- k. r( J# H  ^" i(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( K2 V; u/ a) K# T+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)
# G% L4 v7 l2 r: Y1 ^3 p+A(10)+A(11)+A(12)+A(13)+A(14)+A(15)+A(16))]
; O5 B$ v3 B; c: }

linmatsas

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

迷途的羊

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

linmatsas

回复 迷途的羊 的帖子

是自己编的还用别人运行么……不过要是楼主自己编的……小弟就佩服死啦
   

紫辰

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

chenrye121

佩服两个字…………

水木年华zzu

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

wuming1429

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

hhao

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

abc321

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