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

wuming1429

一、使用曲线拟合模拟病人等待时间的模拟matlab程序实现：
[A]=xlsread('C:\MATLAB7\work\问题3数据5');
* @: d' e+ x5 F9 F7 m, hlen=length(A);
  i+ ?- I. [0 p' I  mi=1:len;
2 n% m) E- G  e1 S! sj=A(:,5);
figure(1);
% g+ _( ~7 v0 iplot(i,j,'.-');
% y  e) Y- r% D; I: f6 {9 khold on;
i=i';
6 b  F6 a* f, @  b4 ?& fp1=polyfit(i,j,2)
0 B  C2 [8 `/ I; o' [7 kf1=polyval(p1,i);
plot(i,j,'.',i,f1,'-');
  ^6 o% V% H  ?. n' `% R0 ehold on;
title('病人平均等待时间与病人到达时间的关系');
xlabel('病人门诊时间距离7 月12 日的时间差');ylabel('病人平均等待时间
5 p2 L' T. ], R  m( x');
7 e0 K* {* B$ C0 Naxis equal;
# G! r4 m: w* a$ l# N$ n) dfigure(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
7 m* h9 q; S8 k& d9 q, A  S/ L9 \; r11.5 13.25 14 13 12 11 13 13 12 12.25 13 15.5 15 14 13];
subplot(1,2,1);
# x9 z9 O! B" B7 w9 Lplot(A1,B1,'.-');
hold on;
p2=polyfit(A1,B1,2)
4 W7 ?4 i$ w, t9 _6 p# ^/ l+ if2=polyval(p2,A1);
( \" p4 h  g+ J! x( N! H5 eplot(A1,B1,'.',A1,f2,'-');
hold on;
title('白内障病人情况');
xlabel('白内障病人门诊时间距7月12日的时间差');ylabel('白内障病人平均
/ c- |9 _9 c* ]* L等待时间');
  I6 T' b, u# s. S9 Saxis equal;
; h4 |1 J* N* @, vA2=[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
% i2 Y% K& J+ L; k8 v: r5 Q& f31 32 33 34 35 36 37 ];
0 C4 W! o3 g' u# fB2=[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
# v" ?: Z$ m  f) Y2 ^% u! @2 ?12 12 13.25 14 14 13 12 12 12.4 13.33 13 13 13 12];
: ?6 H! x: v1 }) S( _subplot(1,2,2);
plot(A2,B2,'.-');
" R& d, B( Y1 U6 hhold on;
p3=polyfit(A2,B2,2)
6 H4 A& V' s/ U& A9 a& wf3=polyval(p3,A2);
plot(A2,B2,'.',A2,f3,'-');
hold on;
4 D3 Z/ I3 e: T& Y  Ftitle('白内障双眼病人情况');
4 p' L  m6 f# Z6 t8 W, J! Exlabel('白内障双眼病人门诊时间距7月12日的时间差');ylabel('白内障双眼
病人平均等待时间');
axis equal;
) Z1 d9 [/ Q0 c4 jfigure(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
/ q7 q, T% x! N28 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];
: H* @- Y$ h4 D1 z" {subplot(2,2,1);
5 Z. r! |7 v4 y+ F8 Pplot(A3,B3,'.-');
hold on;
+ D8 d4 ~4 t% W' Z& Yp4=polyfit(A3,B3,2)
5 Q  Z3 Z: N) h2 d; d4 v3 R* {4 _7 If4=polyval(p4,A3);
plot(A3,B3,'.',A3,f4,'-');
6 D* |7 p) d! shold on;
7 i, O' z; U6 b4 W( _) ]3 y9 Ztitle('视网膜病人情况');
7 u/ ?, N/ |% |" P9 Axlabel('门诊时间距7 月12 日的时间差');ylabel('视网膜病人平均等待时间
( D- o: j0 l2 u4 z" R$ [');
axis equal;
, W1 X/ d1 B( i3 iA4=[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
$ O4 d- }/ R' o# C2 B0 I6 e* l9 P, {( ?37];
- O4 o% M1 I' vB4=[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];
7 ~* {! N: l0 R1 usubplot(2,2,2);
, t; g8 R4 L5 g/ Vplot(A4,B4,'.-');
: z8 W, k; v- T# Uhold on;
$ q: W; g; z8 U/ T/ Ip5=polyfit(A4,B4,2)
f5=polyval(p5,A4);
plot(A4,B4,'.',A4,f5,'-');
- d6 R1 p7 @6 Z* [, fhold on;
" A- C& f- E& R& K8 b, C. utitle('青光眼病人情况');
xlabel('门诊时间距7 月12日的时间差');ylabel('视网膜病人平均等待时间
9 J( f; u9 ^1 n: ]( l');
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];
0 J' M4 |5 }1 ?/ i1 ]) rB4=[12 13 12 12 12 12 13.5 13 13.5 14 13 12 12 11 10.25 12 11 11 15 15
3 r/ x+ r8 J  t# f13 12 12 13 13 12];
/ H4 q- Z/ k' x* D! Esubplot(2,2,2);
plot(A4,B4,'.-');
hold on;
p5=polyfit(A4,B4,2)
f5=polyval(p5,A4);
plot(A4,B4,'.',A4,f5,'-');
( d+ f5 r9 j$ P! \+ Chold on;
title('青光眼病人情况');
: P' x6 c* b/ e3 t+ B* nxlabel('门诊时间距7 月12 日的时间差');ylabel('青光眼病人平均等待时间
1 _8 i3 U0 G/ l$ R');
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];
! G0 `! h/ N1 d# [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];
- w0 _7 f' C+ ]: S: i( `+ m& hsubplot(2,2,3);
5 B  ^7 j4 _  b" ~- e  iplot(A5,B5,'.-');
hold on;
p6=polyfit(A5,B5,2)
3 s9 Y5 M# q+ l: n# xf6=polyval(p6,A5);
plot(A5,B5,'.',A5,f6,'-');
hold on;
) d7 f2 b- U$ U& s; s. Y/ g3 Ytitle('外伤病人情况');
xlabel('外伤病人门诊时间距离7月12日的时间差');ylabel('外伤病人平均等
待时间');
axis equal;
二、层次分析法的一致性指标判定的matlab 实现:
function H=ff(A)
& v5 K3 }: c9 [! J1 Dt=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(
9)+A(10)+A(11)+A(12)+A(13)+A(14)+A(15)+A(16))
1 t5 s0 B4 \' H(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)
9 k- W' U" F7 N- x+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)
, ]6 N( ^+ y$ h8 B: {+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))]
0 M! C) s1 f/ d& n7 X, x

linmatsas

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

迷途的羊

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

linmatsas

回复 迷途的羊 的帖子

是自己编的还用别人运行么……不过要是楼主自己编的……小弟就佩服死啦
0 V4 f7 N9 L+ o0 V8 J2 C- A( z   

紫辰

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

chenrye121

佩服两个字…………

水木年华zzu

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

wuming1429

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

hhao

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

abc321

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