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

wuming1429

一、使用曲线拟合模拟病人等待时间的模拟matlab程序实现：
[A]=xlsread('C:\MATLAB7\work\问题3数据5');
: {8 s$ F# `' Q- l! Ilen=length(A);
i=1:len;
: [, J1 o* D5 C9 _- T! f/ C. w# Cj=A(:,5);
6 a6 ~" ^. e7 n! ]; E  F' ^1 x$ afigure(1);
plot(i,j,'.-');
hold on;
i=i';
7 v; u! w+ G/ o' F+ op1=polyfit(i,j,2)
f1=polyval(p1,i);
0 C+ h( R1 n) {- l; o+ tplot(i,j,'.',i,f1,'-');
hold on;
" E* N+ r8 A5 ?8 V$ c* C0 Gtitle('病人平均等待时间与病人到达时间的关系');
5 f2 c" ]: C+ W* w8 K0 E( Hxlabel('病人门诊时间距离7 月12 日的时间差');ylabel('病人平均等待时间
  E" C% y$ W, `- ]');
axis equal;
6 Y5 Z3 J, y( @figure(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];
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);
  q9 \) [4 @3 v* I9 Jplot(A1,B1,'.-');
% N1 ?8 v. e0 Q# H: \# S$ Phold on;
( d+ a4 f: O5 @; V1 cp2=polyfit(A1,B1,2)
f2=polyval(p2,A1);
9 p1 R" e& m+ I1 F8 f( Wplot(A1,B1,'.',A1,f2,'-');
' @. [+ q* `; {3 T/ K: zhold on;
2 C% g; P& F+ [% [# E" ]; S8 otitle('白内障病人情况');
xlabel('白内障病人门诊时间距7月12日的时间差');ylabel('白内障病人平均
, J: M2 }% N$ a$ N* R  ]* i等待时间');
7 |! W1 ~& s* P8 waxis 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 ];
7 H% r" m! ]8 C( Y$ X9 HB2=[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
/ K; q; A  T0 A( ~& b* f, Y6 R. h12 12 13.25 14 14 13 12 12 12.4 13.33 13 13 13 12];
9 y7 j5 f, g$ ~2 o/ Jsubplot(1,2,2);
plot(A2,B2,'.-');
/ F5 k7 N6 U' H  x" V, M% zhold on;
' y2 P1 l' h8 M2 q3 hp3=polyfit(A2,B2,2)
f3=polyval(p3,A2);
3 X# s/ c  ]6 p; Fplot(A2,B2,'.',A2,f3,'-');
8 I. P- _% ~' k+ yhold on;
8 i  C( E& {- Rtitle('白内障双眼病人情况');
( c* g" I" T5 J# ]xlabel('白内障双眼病人门诊时间距7月12日的时间差');ylabel('白内障双眼
病人平均等待时间');
1 J' P! H5 c2 `+ Taxis equal;
; \: |/ g5 R+ }figure(3);
7 s  w' d' G7 U3 e% h) [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 ];
: M  o  n& f, D; g) W" m3 W  ?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 D6 f7 H. v! J& G% V: Bsubplot(2,2,1);
! z7 y4 M) y$ S6 [plot(A3,B3,'.-');
* L! N) n  f4 W# z6 Xhold on;
p4=polyfit(A3,B3,2)
f4=polyval(p4,A3);
1 ~( ]5 F& g! ?/ Zplot(A3,B3,'.',A3,f4,'-');
! D' h# d3 }2 o. M) Dhold on;
title('视网膜病人情况');
xlabel('门诊时间距7 月12 日的时间差');ylabel('视网膜病人平均等待时间
  b9 `- V2 G! @8 V');
& g; H7 m4 _: _axis equal;
( h% G) x. M9 w) r; b' 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
8 S0 q$ @: a' U$ W2 K5 f37];
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];
- h- `! j$ d$ W" [8 ?1 m% gsubplot(2,2,2);
plot(A4,B4,'.-');
1 A  I4 @3 y; I3 k3 Rhold on;
p5=polyfit(A4,B4,2)
9 |6 V& q2 l8 |$ F% b# c5 x8 s2 [f5=polyval(p5,A4);
! u8 x' }+ w- i- l# r5 Gplot(A4,B4,'.',A4,f5,'-');
hold on;
title('青光眼病人情况');
xlabel('门诊时间距7 月12日的时间差');ylabel('视网膜病人平均等待时间
');
7 |2 S4 _8 M$ h; f6 ?( |& {. Daxis 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
- P0 T. `. R: F" S37];
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];
- c8 |! Y) U/ Wsubplot(2,2,2);
plot(A4,B4,'.-');
hold on;
3 E( t1 t  s, X. ~5 sp5=polyfit(A4,B4,2)
6 L4 e4 @- ^4 C7 d1 ~f5=polyval(p5,A4);
3 |( }5 {% W. [/ R8 pplot(A4,B4,'.',A4,f5,'-');
3 ^8 Y( o( G8 t% K: fhold on;
title('青光眼病人情况');
xlabel('门诊时间距7 月12 日的时间差');ylabel('青光眼病人平均等待时间
');
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];
. i2 v% |( s% J0 d' bB5=[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];
1 ?+ B% R( D( m+ e2 N) f% msubplot(2,2,3);
! C, ?2 j" A' K8 ?( l6 t) u' I1 |plot(A5,B5,'.-');
/ x& a5 b- j& N7 p: f2 ghold on;
p6=polyfit(A5,B5,2)
1 }* ^" x6 `- Z9 F  cf6=polyval(p6,A5);
plot(A5,B5,'.',A5,f6,'-');
  M2 v0 G1 q* V; Fhold on;
title('外伤病人情况');
xlabel('外伤病人门诊时间距离7月12日的时间差');ylabel('外伤病人平均等
6 K$ Z% j3 a2 c1 B& _* a. l待时间');
axis equal;
5 y! V" b5 ~8 ^二、层次分析法的一致性指标判定的matlab 实现:
" e, o, U/ Y+ R4 W$ q$ tfunction H=ff(A)
t=eig(A)
w=max(t)
9 ?0 [5 l" u2 dq=(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))
(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)
( }0 W( c. B" ^% e+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))]
) A; K; C6 z4 h- f' j+ o$ ?6 H

dongdong111

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

ilava

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

qiuhanty

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

hhao

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

Hyacinth

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

liguoli

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

hab0458

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

星际丑男

强烈支持。楼主万岁

casper

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