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

wuming1429

一、使用曲线拟合模拟病人等待时间的模拟matlab程序实现：
[A]=xlsread('C:\MATLAB7\work\问题3数据5');
! @/ {( h! y1 P7 H, u/ @len=length(A);
& B% Y3 ~7 k; ?i=1:len;
j=A(:,5);
6 }( H- ~8 i- A* `figure(1);
; A: s' l& Z) ?plot(i,j,'.-');
hold on;
$ _1 H2 y9 P- |4 j* L6 Vi=i';
p1=polyfit(i,j,2)
8 V" x& i5 n7 d% X( Z2 d. ^3 w  Ff1=polyval(p1,i);
7 {* d$ {6 C' U% R* w/ Bplot(i,j,'.',i,f1,'-');
hold on;
title('病人平均等待时间与病人到达时间的关系');
6 B. L) e6 I9 b" S+ e0 z* c* A/ I9 _xlabel('病人门诊时间距离7 月12 日的时间差');ylabel('病人平均等待时间
');
axis equal;
+ r+ N3 e$ `. Y0 w8 I, Kfigure(2);
2 a1 c' ~6 Q! a+ u% EA1=[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
/ {! V0 P& S$ m32 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);
% R7 |5 F/ |. U' L3 v4 ?plot(A1,B1,'.-');
& {  y: M6 m2 q. ~hold on;
9 J, U' [+ {/ zp2=polyfit(A1,B1,2)
9 ?9 v. R/ v% p, o2 I- mf2=polyval(p2,A1);
1 b! @4 Y- I5 C! p0 T* w. qplot(A1,B1,'.',A1,f2,'-');
. s& [$ ]! W1 f# V5 m9 h( n: t1 H7 Ohold on;
, B6 E1 w, w# L! p, t# wtitle('白内障病人情况');
$ _- D2 Q5 [5 \8 Yxlabel('白内障病人门诊时间距7月12日的时间差');ylabel('白内障病人平均
  m# f% R3 v& E4 h* }4 Q9 [等待时间');
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
/ X; R" g& A4 z; E! H6 }5 j31 32 33 34 35 36 37 ];
% t. \- c% a' k$ X: |- O* G" p. [$ TB2=[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];
- {7 C4 Z' X6 Nsubplot(1,2,2);
; [) P! ]) G' T+ A5 Y: i" Pplot(A2,B2,'.-');
4 P# A1 v5 c9 A; [7 Ghold on;
/ M, l# p( Y+ |# T3 [p3=polyfit(A2,B2,2)
f3=polyval(p3,A2);
" Q0 [) q/ @9 l& O! W' ?9 x" t4 Gplot(A2,B2,'.',A2,f3,'-');
hold on;
: d; \* u8 `7 |$ Btitle('白内障双眼病人情况');
& q6 x+ x5 j7 y+ d# bxlabel('白内障双眼病人门诊时间距7月12日的时间差');ylabel('白内障双眼
1 u$ r7 L9 c8 \" H. t, p7 \6 s; x/ l病人平均等待时间');
axis equal;
+ J3 i4 q+ k5 ^$ J7 w) tfigure(3);
' F$ E% ?% e8 F5 a3 g/ n- fA3=[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 ];
8 \0 y! n% o  ?+ j* z8 T' YB3=[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
' q) J: v; u; B8 b* ~; B11.67 12 12.25 12 11 11.75 13 14 14 13 12 12 12.33 13 13 13 12.5 12 12];
3 j- |' m1 g2 t. H7 z# S) V  F- vsubplot(2,2,1);
plot(A3,B3,'.-');
hold on;
: ^6 `/ s# V: H( ?p4=polyfit(A3,B3,2)
f4=polyval(p4,A3);
plot(A3,B3,'.',A3,f4,'-');
hold on;
; `, n5 x# d) {9 \title('视网膜病人情况');
xlabel('门诊时间距7 月12 日的时间差');ylabel('视网膜病人平均等待时间
');
! b7 D! a9 R% Y- `( V( t" xaxis equal;
6 H$ ?' L3 ?0 g0 A0 H# x0 qA4=[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];
5 A# e9 _$ |; y6 V1 E5 S8 w+ j0 OB4=[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];
! \! T& |3 w7 ysubplot(2,2,2);
3 U7 Z# I" N  X- U7 x% gplot(A4,B4,'.-');
0 X  f# Q0 H6 G4 j1 u) O6 W0 K6 U; ahold on;
7 `  C( U$ a" e+ [; Z, ap5=polyfit(A4,B4,2)
7 O/ q2 q! W* {8 \& P# O- t, Uf5=polyval(p5,A4);
plot(A4,B4,'.',A4,f5,'-');
( c) z+ B! ?5 y' E5 m3 w  \0 Qhold on;
title('青光眼病人情况');
xlabel('门诊时间距7 月12日的时间差');ylabel('视网膜病人平均等待时间
');
6 O& W  ?8 Z  Faxis 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
4 ]4 H. [; o# b0 D2 D% ~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
6 ^6 S5 |" |1 a  O' N; F4 N13 12 12 13 13 12];
subplot(2,2,2);
0 y+ _) \! K  Z  |3 K6 D) @plot(A4,B4,'.-');
hold on;
p5=polyfit(A4,B4,2)
5 k, n. o( z: |" f  E; Vf5=polyval(p5,A4);
plot(A4,B4,'.',A4,f5,'-');
hold on;
: h& b6 g9 D6 B% ]title('青光眼病人情况');
xlabel('门诊时间距7 月12 日的时间差');ylabel('青光眼病人平均等待时间
9 r" [4 J$ D+ G9 Y* r1 w8 `; e$ r');
8 B9 W# F5 U" g) baxis 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];
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];
" p+ W5 [6 J' j9 w1 h' H3 Vsubplot(2,2,3);
plot(A5,B5,'.-');
hold on;
# q& [/ b% n+ t: j8 `5 ?p6=polyfit(A5,B5,2)
f6=polyval(p6,A5);
plot(A5,B5,'.',A5,f6,'-');
6 c! M! P/ l$ z. rhold on;
# A, Y/ y) E5 N% T# x6 l; }title('外伤病人情况');
xlabel('外伤病人门诊时间距离7月12日的时间差');ylabel('外伤病人平均等
3 T; n3 m1 v* `' ]9 ^- V( ]' C4 Y1 f待时间');
axis equal;
2 b5 P' m2 @) Q4 e% ^8 g二、层次分析法的一致性指标判定的matlab 实现:
' r, ^2 [' a( ^( Q& R1 z: hfunction H=ff(A)
t=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(
2 ~/ A$ n, h* `$ Z$ K9 c: ~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)
+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)
3 [. h1 w- Y' i$ |/ L) Z7 P& K- `! ]+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))]

linmatsas

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

迷途的羊

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

linmatsas

回复 迷途的羊 的帖子
) n' d, T* ^% }1 {$ x5 D. \
是自己编的还用别人运行么……不过要是楼主自己编的……小弟就佩服死啦
; W& P* {" d- ?- m( n# ]+ t& c   

紫辰

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

chenrye121

佩服两个字…………

水木年华zzu

[A]=xlsread('C:\MATLAB7\work\问题3数据5');
0 B1 l8 h: [6 a/ t" N. [: x' B在不同的地方运行上面一句要改成本地数据保存的路径的

wuming1429

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

hhao

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

abc321

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