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

wuming1429

一、使用曲线拟合模拟病人等待时间的模拟matlab程序实现：
  t; c- T7 M! ~1 W8 o9 S[A]=xlsread('C:\MATLAB7\work\问题3数据5');
len=length(A);
$ H% q  j( r2 N2 E, \+ k5 oi=1:len;
9 Z6 z1 S4 u. v% v  zj=A(:,5);
figure(1);
plot(i,j,'.-');
% y0 m1 l2 ]8 q% Hhold on;
$ D& X6 v& h$ D: j) Z; L* Ci=i';
8 D. r$ t1 _, U% lp1=polyfit(i,j,2)
- p2 {1 y- x2 a% X) c) w4 B, W8 hf1=polyval(p1,i);
. P( K; j' q5 \; gplot(i,j,'.',i,f1,'-');
hold on;
+ \9 y5 `1 s! j4 ?; |9 @9 N) ~title('病人平均等待时间与病人到达时间的关系');
1 Q( K7 D- Z/ K3 N* |$ O: I0 |5 ?6 Wxlabel('病人门诊时间距离7 月12 日的时间差');ylabel('病人平均等待时间
0 |* i6 {9 p: a' h" k');
axis equal;
+ h( L. e/ A& w2 efigure(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
, z" m- u# W$ ^' f7 J% i! j32 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);
* F7 B/ `. x) t7 u# R1 yplot(A1,B1,'.-');
hold on;
7 b" i6 |4 @  k  kp2=polyfit(A1,B1,2)
1 J8 Z" t, [2 vf2=polyval(p2,A1);
plot(A1,B1,'.',A1,f2,'-');
( y! Z, _8 J3 i: A' Fhold on;
; \6 o+ w& {& x- o7 n* X2 Ftitle('白内障病人情况');
1 D# K7 w/ o, V3 W8 n6 i9 qxlabel('白内障病人门诊时间距7月12日的时间差');ylabel('白内障病人平均
: s5 O* ^% t# R8 g, O等待时间');
axis equal;
3 E# T/ J, d1 b$ v+ J1 B  CA2=[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 ];
& N, v( S9 Z3 \/ N5 vB2=[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];
, _# S5 F' Z( x$ X& {subplot(1,2,2);
plot(A2,B2,'.-');
' |2 ~. U+ Z* {/ [! g! ?hold on;
  P1 I5 Z- q# |7 t! Up3=polyfit(A2,B2,2)
f3=polyval(p3,A2);
plot(A2,B2,'.',A2,f3,'-');
hold on;
2 x+ R1 R4 y3 {8 C( rtitle('白内障双眼病人情况');
xlabel('白内障双眼病人门诊时间距7月12日的时间差');ylabel('白内障双眼
病人平均等待时间');
+ p# N8 R& m  B( ]8 C) F9 Y3 b5 Iaxis 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
11.67 12 12.25 12 11 11.75 13 14 14 13 12 12 12.33 13 13 13 12.5 12 12];
1 u% E6 y- I  ?/ N! }subplot(2,2,1);
plot(A3,B3,'.-');
hold on;
p4=polyfit(A3,B3,2)
f4=polyval(p4,A3);
plot(A3,B3,'.',A3,f4,'-');
: s5 O& [0 j' {! hhold on;
% y4 }6 r9 @. Y) P0 K- }' ctitle('视网膜病人情况');
+ U- H7 R% ?' Z2 e9 q; ~6 Y/ txlabel('门诊时间距7 月12 日的时间差');ylabel('视网膜病人平均等待时间
');
" b0 A* V3 B, l4 Waxis equal;
) A2 \' o; P  K. v1 e: @6 fA4=[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];
# p7 C& N8 B* Xsubplot(2,2,2);
3 t( I" @6 Z/ m; x# gplot(A4,B4,'.-');
hold on;
! U7 f+ z" ]; N; I: `p5=polyfit(A4,B4,2)
f5=polyval(p5,A4);
plot(A4,B4,'.',A4,f5,'-');
- m2 Z3 r9 S& Rhold on;
7 Z. ^% k- v2 s- Ztitle('青光眼病人情况');
3 @3 K$ E' @! {2 ?% G7 L6 K5 cxlabel('门诊时间距7 月12日的时间差');ylabel('视网膜病人平均等待时间
');
axis equal;
: ~; d) E* j' ?4 X0 p" yA4=[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);
, O" _6 g$ Q* r7 _! b( p% Dplot(A4,B4,'.-');
7 |2 I# T7 L3 R. O2 M- phold on;
p5=polyfit(A4,B4,2)
f5=polyval(p5,A4);
plot(A4,B4,'.',A4,f5,'-');
hold on;
4 `& r4 `( u. b' o( [( n. m1 @) ~; |title('青光眼病人情况');
. A7 w, C8 k; }+ ^1 }2 x9 vxlabel('门诊时间距7 月12 日的时间差');ylabel('青光眼病人平均等待时间
9 i1 T4 B! ~4 ]' G4 }1 T8 d');
axis equal;
7 g/ F$ l$ C+ `, cA5=[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];
8 r$ @0 z! H2 x8 J. x2 JB5=[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];
( L- i# A% m( E- ?; m/ Y; w' I) Msubplot(2,2,3);
1 d5 K: Z$ B( q. c4 q7 wplot(A5,B5,'.-');
2 t+ o1 R9 s- ?& }hold on;
/ h5 Q1 d  Z) k4 M7 b2 qp6=polyfit(A5,B5,2)
/ A" {6 ^) E: }0 T0 y8 m# if6=polyval(p6,A5);
$ _( o$ N: l; m0 Yplot(A5,B5,'.',A5,f6,'-');
! I5 g/ @& t* u7 D. E1 E' m! jhold on;
title('外伤病人情况');
xlabel('外伤病人门诊时间距离7月12日的时间差');ylabel('外伤病人平均等
待时间');
$ O/ y; H' Y: A! }1 V; Q% Xaxis equal;
4 L: n% ?5 k' u. ]+ f二、层次分析法的一致性指标判定的matlab 实现:
  F+ Z; J* o9 @4 R  x1 qfunction H=ff(A)
t=eig(A)
, g& c' I; t. Y+ ~, ^, Y( V: @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(
9)+A(10)+A(11)+A(12)+A(13)+A(14)+A(15)+A(16))
4 g0 m1 U. ?$ t  z( l. k(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))
& U4 H  v* M' c% C! X" Q2 w( n(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)
! f9 \2 {) X. P+A(10)+A(11)+A(12)+A(13)+A(14)+A(15)+A(16))
7 F* P/ M* j/ |9 _: X(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)
+ U0 A) G. f3 [7 l" c) O0 h+A(10)+A(11)+A(12)+A(13)+A(14)+A(15)+A(16))]

dongdong111

其实代码质量不高。。。。明显可视化设置没到位。然后后面的AHP写的也很乱还有bug

ilava

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

qiuhanty

声明一下：本人看贴和回贴的规则，好贴必看，精华贴必回。

hhao

我来了~~~~~~~~~ 闪人~~~~~~~~~~~~~~~~  

Hyacinth

我来了~~~~~~~~~ 闪人~~~~~~~~~~~~~~~~  

liguoli

楼主的帖子实在是写得太好了。可是我立刻想到，这么好的帖子，倘若别人看不到，那么不是浪费楼主的心血吗？经过痛苦的思想斗争，我终于下定决心，牺牲小我，奉献大我。我要拿出这帖子奉献给世人赏阅，我要把这个帖子一直往上顶，往上顶！顶到所有人都看到为止！  

hab0458

楼主那种裂纸欲出的大手笔，竟使我忍不住一次次的翻开楼主的帖子……   

星际丑男

强烈支持。楼主万岁

casper

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