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

wuming1429

一、使用曲线拟合模拟病人等待时间的模拟matlab程序实现：
1 ?; [" _9 l+ q) ~+ @[A]=xlsread('C:\MATLAB7\work\问题3数据5');
/ a2 H, F4 w! ]! clen=length(A);
: o4 f& s6 q% li=1:len;
j=A(:,5);
' l8 H# q9 N2 y$ e1 D3 N% I* Pfigure(1);
plot(i,j,'.-');
hold on;
5 ^% o6 ^1 B  Ti=i';
: X9 q( U4 A' u" {1 i2 y; A% \p1=polyfit(i,j,2)
f1=polyval(p1,i);
' h* S3 [* l" m* d( Z5 Kplot(i,j,'.',i,f1,'-');
hold on;
title('病人平均等待时间与病人到达时间的关系');
+ N$ t  N- [2 V7 l! Bxlabel('病人门诊时间距离7 月12 日的时间差');ylabel('病人平均等待时间
');
7 x! n: A7 ]( n& n$ Saxis equal;
figure(2);
4 T( G$ U% Q6 v# B' QA1=[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
- {" w$ ~2 }1 @3 p32 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);
* j! B' I0 d& i6 `- kplot(A1,B1,'.-');
hold on;
0 U7 D9 |7 Y/ ?5 p, J- A3 Up2=polyfit(A1,B1,2)
; B( W3 W* ]/ w" i. `f2=polyval(p2,A1);
9 r7 d: E0 R4 K- Q$ }4 P1 W6 splot(A1,B1,'.',A1,f2,'-');
hold on;
title('白内障病人情况');
xlabel('白内障病人门诊时间距7月12日的时间差');ylabel('白内障病人平均
8 p7 m- n+ E" s, p等待时间');
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 ];
+ W4 T/ L, T/ v( L1 e$ B0 I+ EB2=[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];
9 R$ ~2 H0 e, R. A; T, gsubplot(1,2,2);
plot(A2,B2,'.-');
hold on;
! i4 l. Y% Z4 q" m% `: J0 n9 _p3=polyfit(A2,B2,2)
5 |8 l  \. ^7 x6 h2 z$ c/ @f3=polyval(p3,A2);
% t7 _0 l% {/ e8 y8 \, }plot(A2,B2,'.',A2,f3,'-');
0 ?+ Y) b: t. u; x% d8 Xhold on;
6 N3 ^# _1 K, d; Ctitle('白内障双眼病人情况');
xlabel('白内障双眼病人门诊时间距7月12日的时间差');ylabel('白内障双眼
+ z" z) c7 a% W4 Z病人平均等待时间');
+ B( h, ]2 A6 x! `axis equal;
( J, V# [3 d* A. q3 v% Ofigure(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
9 d2 o) `7 \2 |28 29 30 31 32 33 34 35 36 37 38 ];
# b$ n- y1 j: e# J4 k8 q6 vB3=[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
+ X& Y5 q0 R2 x2 h2 [11.67 12 12.25 12 11 11.75 13 14 14 13 12 12 12.33 13 13 13 12.5 12 12];
" F" f; z! G8 }1 e' Fsubplot(2,2,1);
plot(A3,B3,'.-');
* t* C" p- l2 w- shold on;
p4=polyfit(A3,B3,2)
f4=polyval(p4,A3);
plot(A3,B3,'.',A3,f4,'-');
. \# z9 n- s" ~" a# ehold on;
title('视网膜病人情况');
xlabel('门诊时间距7 月12 日的时间差');ylabel('视网膜病人平均等待时间
  g* i: B" A& e$ p  ?5 k* c');
& d0 ?9 f3 c7 l% b: v) naxis equal;
& S0 _" T8 |0 `8 s: X5 F& T# aA4=[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
13 12 12 13 13 12];
subplot(2,2,2);
" u& n% @: E: w0 r; dplot(A4,B4,'.-');
hold on;
p5=polyfit(A4,B4,2)
# F- C2 N+ O: R$ S# u+ q) v% \" B* Xf5=polyval(p5,A4);
plot(A4,B4,'.',A4,f5,'-');
, P! c; u' r7 b( @hold on;
6 @/ U/ Y& ]; @. f. ytitle('青光眼病人情况');
' L/ O6 }& r) T* }" Gxlabel('门诊时间距7 月12日的时间差');ylabel('视网膜病人平均等待时间
');
- C4 G8 I* k4 q. Q, Uaxis 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
: l! u3 q  X5 r' R6 B4 m37];
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];
6 O7 r- M3 N1 usubplot(2,2,2);
  w7 m: }' }6 `+ N5 ?3 zplot(A4,B4,'.-');
% p$ o! \2 Y8 ^# Y2 f* _hold on;
p5=polyfit(A4,B4,2)
f5=polyval(p5,A4);
% p6 X% q! `* j! Bplot(A4,B4,'.',A4,f5,'-');
( b  ]  K3 K8 W% \: e1 fhold on;
& j0 x+ d6 E% ^& ~! O5 I! Mtitle('青光眼病人情况');
# r2 c. y% C  S! U: _% [1 Vxlabel('门诊时间距7 月12 日的时间差');ylabel('青光眼病人平均等待时间
');
axis equal;
  u  M/ Y' B; M( z% ]# IA5=[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];
subplot(2,2,3);
' k" c7 g$ z8 @" hplot(A5,B5,'.-');
hold on;
$ q, b8 l+ z5 U7 D: cp6=polyfit(A5,B5,2)
f6=polyval(p6,A5);
% j( ~1 y# k- Tplot(A5,B5,'.',A5,f6,'-');
hold on;
title('外伤病人情况');
xlabel('外伤病人门诊时间距离7月12日的时间差');ylabel('外伤病人平均等
待时间');
axis equal;
7 ^$ P, u' s' P二、层次分析法的一致性指标判定的matlab 实现:
, f( C0 t: a/ l( r# w( wfunction H=ff(A)
' y( ?" ~  l# _/ w! X# lt=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))
& F- E) M0 l( B$ ]/ L' M(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)
; d1 g# N) J  z9 ~$ d+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))
( p6 A8 }+ b7 x7 ]8 F8 E5 m(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)
, r; M- Q8 w$ W7 a7 W" }+A(10)+A(11)+A(12)+A(13)+A(14)+A(15)+A(16))]
$ d# Y. I4 f' Y0 N

linmatsas

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

迷途的羊

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

linmatsas

回复 迷途的羊 的帖子

是自己编的还用别人运行么……不过要是楼主自己编的……小弟就佩服死啦
   

紫辰

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

chenrye121

佩服两个字…………

水木年华zzu

[A]=xlsread('C:\MATLAB7\work\问题3数据5');
6 b/ E. K. C: y' ?+ ?2 \; \, y  |0 ?( G在不同的地方运行上面一句要改成本地数据保存的路径的

wuming1429

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

hhao

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

abc321

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