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

wuming1429

一、使用曲线拟合模拟病人等待时间的模拟matlab程序实现：
1 ?# ^7 M" I. i7 [- l[A]=xlsread('C:\MATLAB7\work\问题3数据5');
" A# N- z: ?% x7 llen=length(A);
i=1:len;
0 B6 U0 {. z' s! p) ^j=A(:,5);
3 J0 g- A- U, Dfigure(1);
plot(i,j,'.-');
" d' W* f9 T7 \/ `hold on;
) ]4 l% P* {' v7 G6 Ei=i';
3 A- `8 g' I3 ]9 V  M5 e* w) {p1=polyfit(i,j,2)
f1=polyval(p1,i);
) \3 e" O( x/ }  c; e2 Gplot(i,j,'.',i,f1,'-');
hold on;
) [" l+ Q9 B# @$ Utitle('病人平均等待时间与病人到达时间的关系');
xlabel('病人门诊时间距离7 月12 日的时间差');ylabel('病人平均等待时间
');
  u8 B" |# k& U' Faxis equal;
; I# ~0 y% `3 L: b" [6 ?/ L" hfigure(2);
! s  H& U' @; ?7 }2 |5 [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
9 l) a1 k, ^* ^& d# \+ n32 34 35 36 37 38 39 40 41 42];
% M" o- _! |5 r1 i1 D. g; eB1=[12 12 12 12 12 12.5 13 13 13 13 13.4 14 14 13 12.5 11 11 12 12 12 11
/ X6 C# N* {' X4 }11.5 13.25 14 13 12 11 13 13 12 12.25 13 15.5 15 14 13];
subplot(1,2,1);
  l; p1 R! Y5 F3 q( y& l& Eplot(A1,B1,'.-');
hold on;
p2=polyfit(A1,B1,2)
' ^3 c1 \. u* j4 n7 df2=polyval(p2,A1);
: P6 F( h' S* v. \( e) u- ^* aplot(A1,B1,'.',A1,f2,'-');
2 w$ a6 Q5 m; ^4 u* b0 @( mhold on;
title('白内障病人情况');
& K+ {3 C: ?1 R1 e' N0 h$ dxlabel('白内障病人门诊时间距7月12日的时间差');ylabel('白内障病人平均
等待时间');
axis equal;
7 z' P& x  o  X5 u7 j: @1 |. ]7 eA2=[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
7 ?( \) l* Q5 a31 32 33 34 35 36 37 ];
9 K" L, ?* S, w1 r' E/ j  F6 u7 mB2=[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
4 `. G3 @$ u3 D: r12 12 13.25 14 14 13 12 12 12.4 13.33 13 13 13 12];
subplot(1,2,2);
plot(A2,B2,'.-');
hold on;
& v( Z& P& n* Z# d& o2 _* up3=polyfit(A2,B2,2)
f3=polyval(p3,A2);
plot(A2,B2,'.',A2,f3,'-');
hold on;
1 l0 t* O/ _. c4 T$ B" i  [+ ]title('白内障双眼病人情况');
6 S) O, W0 u! y  G/ _xlabel('白内障双眼病人门诊时间距7月12日的时间差');ylabel('白内障双眼
病人平均等待时间');
axis equal;
+ O- m, g4 O; Y: M8 n2 h% w+ N6 |- Lfigure(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
& m& T# T  [" p, \1 ^( {28 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
9 Q7 a+ Q' q9 X! w: V11.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);
$ G3 Q+ B- ]6 a% v( N# |6 Hplot(A3,B3,'.-');
hold on;
p4=polyfit(A3,B3,2)
f4=polyval(p4,A3);
plot(A3,B3,'.',A3,f4,'-');
hold on;
title('视网膜病人情况');
% O5 m1 a/ y; r/ b" [xlabel('门诊时间距7 月12 日的时间差');ylabel('视网膜病人平均等待时间
');
axis equal;
" S" \1 l- M. C: p, 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
0 C9 J% g; T# k4 a3 w4 D7 v37];
B4=[12 13 12 12 12 12 13.5 13 13.5 14 13 12 12 11 10.25 12 11 11 15 15
9 F6 K' G+ X3 l4 f13 12 12 13 13 12];
subplot(2,2,2);
plot(A4,B4,'.-');
- ?& S7 Y$ U3 D1 [  lhold on;
p5=polyfit(A4,B4,2)
f5=polyval(p5,A4);
" _7 ]' U& t. b* Q2 y$ ~( lplot(A4,B4,'.',A4,f5,'-');
1 C2 ?- D  @. `$ ^+ I& @hold on;
( p/ ~7 L3 S+ ]title('青光眼病人情况');
xlabel('门诊时间距7 月12日的时间差');ylabel('视网膜病人平均等待时间
& {$ h/ l6 u2 M( ?/ I; R5 G! j');
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
4 B) i  T4 K/ X5 r) a3 O2 X37];
B4=[12 13 12 12 12 12 13.5 13 13.5 14 13 12 12 11 10.25 12 11 11 15 15
, r* m& y8 c+ y  A9 O13 12 12 13 13 12];
subplot(2,2,2);
1 p% Q( D! M7 E6 n; Q' L2 {% cplot(A4,B4,'.-');
- d4 `* n8 j2 o! w; b4 nhold on;
p5=polyfit(A4,B4,2)
f5=polyval(p5,A4);
plot(A4,B4,'.',A4,f5,'-');
hold on;
" Y4 m5 s% e/ A5 I# L+ @1 ]title('青光眼病人情况');
* P- S* |/ ]5 Y0 c9 q5 cxlabel('门诊时间距7 月12 日的时间差');ylabel('青光眼病人平均等待时间
');
axis equal;
. V/ ]" Z7 z% B" V) r; FA5=[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];
6 C8 a/ c/ p' k2 f4 Y+ `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];
subplot(2,2,3);
plot(A5,B5,'.-');
' `; Q# x; ?$ K5 F( ]" Nhold on;
p6=polyfit(A5,B5,2)
* n0 l' @. X7 ]& Z: v5 g; nf6=polyval(p6,A5);
; h! ~# y; i6 J- mplot(A5,B5,'.',A5,f6,'-');
hold on;
4 N* ~1 P1 M( C6 \% j: Otitle('外伤病人情况');
xlabel('外伤病人门诊时间距离7月12日的时间差');ylabel('外伤病人平均等
6 z3 r. B* r5 ?4 _1 j0 G- Y& W待时间');
) ?% H+ P( r; w; s: U4 Yaxis equal;
二、层次分析法的一致性指标判定的matlab 实现:
2 D- Q! ~% R$ Q* afunction H=ff(A)
; S4 `' o( J; P/ Kt=eig(A)
& K) @; h/ h6 l% zw=max(t)
! A  p; T4 [( J7 @6 b! Vq=(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))
+ d$ Z( ~1 k+ C3 |. C(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)
& h2 E/ M/ p% p, j2 ?+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))
6 e; \3 e  t, g% z(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))]
4 `2 `/ e3 K$ P% C  N

linmatsas

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

迷途的羊

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

linmatsas

回复 迷途的羊 的帖子

5 k' N! ?" a5 G4 G6 R是自己编的还用别人运行么……不过要是楼主自己编的……小弟就佩服死啦
7 ]) K) F8 P1 {+ O6 C: ^  \   

紫辰

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

chenrye121

佩服两个字…………

水木年华zzu

[A]=xlsread('C:\MATLAB7\work\问题3数据5');
4 }% k) x, M* c在不同的地方运行上面一句要改成本地数据保存的路径的

wuming1429

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

hhao

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

abc321

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