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

wuming1429

一、使用曲线拟合模拟病人等待时间的模拟matlab程序实现：
8 W* m: p6 ~3 U7 x6 g3 w" u[A]=xlsread('C:\MATLAB7\work\问题3数据5');
# p0 O# {+ P# A" w1 g8 i4 slen=length(A);
i=1:len;
j=A(:,5);
figure(1);
plot(i,j,'.-');
hold on;
i=i';
6 g: \2 ^4 R% y* e  c" i: Q9 `/ }p1=polyfit(i,j,2)
3 U. i) o6 e+ y( A' j; `( f& Qf1=polyval(p1,i);
& S4 E  P. j4 i4 j5 Qplot(i,j,'.',i,f1,'-');
" `5 \+ s$ w8 b. Z- Whold on;
* g& @- E8 y$ ktitle('病人平均等待时间与病人到达时间的关系');
xlabel('病人门诊时间距离7 月12 日的时间差');ylabel('病人平均等待时间
');
axis equal;
( ^% a9 z0 ]' q8 a$ |- R- v+ j" r8 Kfigure(2);
7 L4 i$ y- e4 L; v) p9 KA1=[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];
3 w5 P, ]5 g9 ^  s2 E, W6 wB1=[12 12 12 12 12 12.5 13 13 13 13 13.4 14 14 13 12.5 11 11 12 12 12 11
" Y+ h3 i6 Z3 }! J4 X5 p11.5 13.25 14 13 12 11 13 13 12 12.25 13 15.5 15 14 13];
subplot(1,2,1);
plot(A1,B1,'.-');
hold on;
$ B1 F+ |# v7 _; T' ^p2=polyfit(A1,B1,2)
f2=polyval(p2,A1);
plot(A1,B1,'.',A1,f2,'-');
hold on;
1 ^7 r6 _" {1 @title('白内障病人情况');
xlabel('白内障病人门诊时间距7月12日的时间差');ylabel('白内障病人平均
等待时间');
/ z! ~' u" A- C$ m; Daxis equal;
) w4 ?' E) \4 i2 H. Q( L/ F3 gA2=[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
# v: h. F; ?/ Y4 ]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
) x* j# H3 F; g" [12 12 13.25 14 14 13 12 12 12.4 13.33 13 13 13 12];
8 W  L. D' h- H* x; r( C# I+ }subplot(1,2,2);
7 \$ h$ d" t! z' s1 R9 l- M% h2 uplot(A2,B2,'.-');
' I& N" m$ H- V  A2 |! Jhold on;
0 h+ B5 ?) g2 x& G- Z6 y& dp3=polyfit(A2,B2,2)
f3=polyval(p3,A2);
plot(A2,B2,'.',A2,f3,'-');
: i# R6 x* z$ j* K% i: Uhold on;
0 F& p- ?- k; l( A. v2 B% ~title('白内障双眼病人情况');
xlabel('白内障双眼病人门诊时间距7月12日的时间差');ylabel('白内障双眼
病人平均等待时间');
( y/ t. K' w. D8 l& A( yaxis equal;
figure(3);
) V% O; c6 a+ `3 s( tA3=[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 ];
& v3 _% K8 c2 a. x6 C+ EB3=[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
% ^6 \! J8 z2 Q: A11.67 12 12.25 12 11 11.75 13 14 14 13 12 12 12.33 13 13 13 12.5 12 12];
" ~3 U9 L/ A( B$ [( Asubplot(2,2,1);
& a* E* ~& u/ n& Eplot(A3,B3,'.-');
, y0 s# }0 u4 K  F/ E' v, ~7 Bhold on;
p4=polyfit(A3,B3,2)
9 X4 {& p; c( _9 @& D! ~f4=polyval(p4,A3);
% Z+ M/ u8 v  O- M' ~- b# i0 Q, Nplot(A3,B3,'.',A3,f4,'-');
+ J' N8 o9 m! P3 j& V# ?9 qhold on;
' a0 k( p; G/ Jtitle('视网膜病人情况');
xlabel('门诊时间距7 月12 日的时间差');ylabel('视网膜病人平均等待时间
4 g1 r9 l/ p# t' V6 d');
axis equal;
3 }% `8 ?9 x8 z0 v% ?% a. z1 A4 ^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
% S: z, C' z/ }- _2 {13 12 12 13 13 12];
4 P* ~" R9 X# S) S2 V$ C3 q2 ^subplot(2,2,2);
plot(A4,B4,'.-');
0 \/ K" y/ x8 {6 Ohold on;
7 p9 t9 M4 F' k( Sp5=polyfit(A4,B4,2)
9 }( A# u9 u8 n0 p) Rf5=polyval(p5,A4);
plot(A4,B4,'.',A4,f5,'-');
" L' a# |1 r8 f+ I$ O( z% nhold on;
title('青光眼病人情况');
  R! @( x+ |; [1 Fxlabel('门诊时间距7 月12日的时间差');ylabel('视网膜病人平均等待时间
');
axis equal;
  Q. K+ [8 l3 a7 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
37];
0 }) D/ @. P8 M, J9 b' [* x  O9 [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];
8 g5 B7 Q* y* Qsubplot(2,2,2);
plot(A4,B4,'.-');
hold on;
p5=polyfit(A4,B4,2)
. y$ I' L% `: u0 i' df5=polyval(p5,A4);
plot(A4,B4,'.',A4,f5,'-');
hold 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];
+ e3 Q- ?8 `4 b+ 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];
subplot(2,2,3);
plot(A5,B5,'.-');
hold on;
p6=polyfit(A5,B5,2)
9 \4 P7 U4 m& ]5 L1 |" K- M4 cf6=polyval(p6,A5);
9 `9 ?4 e' e5 d$ ]6 g% K( Nplot(A5,B5,'.',A5,f6,'-');
hold on;
title('外伤病人情况');
xlabel('外伤病人门诊时间距离7月12日的时间差');ylabel('外伤病人平均等
待时间');
axis equal;
二、层次分析法的一致性指标判定的matlab 实现:
function H=ff(A)
9 ]4 R- U5 ?  p8 ^  N" Qt=eig(A)
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(
+ W4 O! O+ Z  {# R/ b9)+A(10)+A(11)+A(12)+A(13)+A(14)+A(15)+A(16))
! J0 @. x8 H- m1 i7 A* ~(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)
! {/ X6 e9 K: i2 b- A/ m+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)
! W% H' u7 p& r8 ?. W+A(10)+A(11)+A(12)+A(13)+A(14)+A(15)+A(16))]

dongdong111

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

ilava

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

qiuhanty

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

hhao

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

Hyacinth

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

liguoli

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

hab0458

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

星际丑男

强烈支持。楼主万岁

casper

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