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

wuming1429

一、使用曲线拟合模拟病人等待时间的模拟matlab程序实现：
[A]=xlsread('C:\MATLAB7\work\问题3数据5');
2 h8 L$ b7 \9 \. z! j* Ulen=length(A);
i=1:len;
( ~5 F7 K/ a8 p( e; }j=A(:,5);
! k0 k* y. [- `* yfigure(1);
plot(i,j,'.-');
! N) ^7 W2 A& I$ K1 Ihold on;
" w1 I* M+ D8 k, V; li=i';
p1=polyfit(i,j,2)
$ x3 }4 F5 O! c4 ]6 s+ Kf1=polyval(p1,i);
plot(i,j,'.',i,f1,'-');
hold on;
2 |# b4 R& i: e- F* Wtitle('病人平均等待时间与病人到达时间的关系');
- ^( c, K( P3 m  Y) Rxlabel('病人门诊时间距离7 月12 日的时间差');ylabel('病人平均等待时间
5 \% [0 u6 e: {');
( s* i7 L5 v. A3 h; ?6 n9 a$ S% xaxis equal;
' {! ]& z+ Z% b- G# ~figure(2);
  _7 s% Y) t8 s7 Y, M* g& jA1=[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
- r& Q8 T+ }4 _3 Q11.5 13.25 14 13 12 11 13 13 12 12.25 13 15.5 15 14 13];
subplot(1,2,1);
4 v; ]8 A) i. x  xplot(A1,B1,'.-');
% p, }& }- I1 C& G2 `9 ?8 bhold on;
- _6 y$ Z0 F8 c/ X2 a2 ]p2=polyfit(A1,B1,2)
f2=polyval(p2,A1);
plot(A1,B1,'.',A1,f2,'-');
& Z2 Y6 |1 S  ~( ~hold on;
" h+ n4 u9 [0 S  X. \0 ^title('白内障病人情况');
3 ~& y6 |! v0 Z& S" _xlabel('白内障病人门诊时间距7月12日的时间差');ylabel('白内障病人平均
等待时间');
! \" n5 H6 ?7 U: T0 M- _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];
subplot(1,2,2);
plot(A2,B2,'.-');
4 d( T: o$ Y) e: R5 F/ [( Fhold on;
% C+ c* t) W+ Z+ P: K" l/ r' ~p3=polyfit(A2,B2,2)
* ^5 v; n7 x/ j2 Nf3=polyval(p3,A2);
7 [* s3 L+ B- W& M8 a! P# Lplot(A2,B2,'.',A2,f3,'-');
  f& K8 u% {. f- b0 }hold on;
# u; L( x! k8 rtitle('白内障双眼病人情况');
" u5 z1 Z* [3 M4 ?  h6 ~xlabel('白内障双眼病人门诊时间距7月12日的时间差');ylabel('白内障双眼
病人平均等待时间');
axis equal;
( s" m6 }$ e& I1 @1 J7 n* U! p. @- Kfigure(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
+ v# h$ u! Q6 }1 e28 29 30 31 32 33 34 35 36 37 38 ];
" K2 n; U0 B) W6 {* f, tB3=[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
% I8 ?' X  N: h% B  u/ X11.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,'.-');
hold on;
9 S& w( j6 Q, N# }! e0 p6 Y; ep4=polyfit(A3,B3,2)
  L: }' l+ [) L/ [; W  Pf4=polyval(p4,A3);
7 L' e' a/ D9 ?9 e$ [* tplot(A3,B3,'.',A3,f4,'-');
- `4 X$ q5 F9 E/ |5 D. Vhold on;
1 k! o$ j$ m' |$ I! |7 M* k8 ptitle('视网膜病人情况');
5 ^* {- m7 }' A5 p7 y# Pxlabel('门诊时间距7 月12 日的时间差');ylabel('视网膜病人平均等待时间
  ~6 X: u5 o8 A% Z5 M, m');
. r! w* B$ `6 I2 H9 A4 u3 Saxis equal;
# A  m0 i5 W6 t" a( L. D' BA4=[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
  y1 h3 b% t  I( B37];
8 F/ |, p$ ?% m. H1 N8 K2 R* hB4=[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];
; J7 r; h$ U+ k, C' lsubplot(2,2,2);
* [* F, N% O; @  {8 ?plot(A4,B4,'.-');
hold on;
8 s+ z* Z# I2 m5 K) s1 k/ ~p5=polyfit(A4,B4,2)
, K) i) A) i. l; ?$ Z2 Ef5=polyval(p5,A4);
* p) [7 [* {$ @9 `" P1 B: w& S5 Uplot(A4,B4,'.',A4,f5,'-');
hold on;
title('青光眼病人情况');
9 [) _* T3 `3 P7 r7 P4 d" [* Rxlabel('门诊时间距7 月12日的时间差');ylabel('视网膜病人平均等待时间
+ P6 Y; k: v% T');
axis equal;
5 D# p: D" c+ _+ u7 n$ {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
* f4 v" \7 h$ b/ V  ]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
! f  L& r7 k7 F13 12 12 13 13 12];
subplot(2,2,2);
9 l( x$ N# E& D: D1 q: Lplot(A4,B4,'.-');
hold on;
. @  I& Q( q6 B  bp5=polyfit(A4,B4,2)
+ M% q- p. O! z; @1 t6 xf5=polyval(p5,A4);
" Y( C  @. X* ?- k) Vplot(A4,B4,'.',A4,f5,'-');
hold on;
title('青光眼病人情况');
xlabel('门诊时间距7 月12 日的时间差');ylabel('青光眼病人平均等待时间
');
" r7 e5 w* Q. w3 r: K; G! @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];
1 d6 K; U0 {  V1 q/ [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];
: B% b7 U5 s5 U. `subplot(2,2,3);
  J0 @/ O/ |7 D# m& j" ~plot(A5,B5,'.-');
* l8 B" T- i' N% e6 s# x) a% lhold on;
  N& b4 S0 @. v! a7 J4 r- op6=polyfit(A5,B5,2)
f6=polyval(p6,A5);
plot(A5,B5,'.',A5,f6,'-');
/ x/ x' T: u) y3 `/ N" Bhold on;
title('外伤病人情况');
. Y/ T4 G  F( P# Rxlabel('外伤病人门诊时间距离7月12日的时间差');ylabel('外伤病人平均等
待时间');
$ Q1 x/ I) a2 D5 T4 k, u; Caxis equal;
% w: _% n3 O7 l1 f1 a- N二、层次分析法的一致性指标判定的matlab 实现:
3 x/ L+ g; h+ ?) ]function H=ff(A)
( E; y" r  H( f# ~t=eig(A)
5 }0 W4 F" |5 M% H/ tw=max(t)
! a! ~6 |- Q+ R; G) M, ]: qq=(w-4)/3/.9
# w: ~. S3 |/ zs=[(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))
  L) u& u2 d$ _(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)
+ f2 f8 t' G8 r# ?4 ~' \& i+A(10)+A(11)+A(12)+A(13)+A(14)+A(15)+A(16))
! ]0 x! `* V* w) d0 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)
+A(10)+A(11)+A(12)+A(13)+A(14)+A(15)+A(16))
0 N7 N+ _/ s3 B  r) b, }(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)
! @$ p: T; M; J' `+A(10)+A(11)+A(12)+A(13)+A(14)+A(15)+A(16))]

dongdong111

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

ilava

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

qiuhanty

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

hhao

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

Hyacinth

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

liguoli

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

hab0458

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

星际丑男

强烈支持。楼主万岁

casper

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