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

wuming1429

一、使用曲线拟合模拟病人等待时间的模拟matlab程序实现：
" m7 ]# s7 K; v4 U[A]=xlsread('C:\MATLAB7\work\问题3数据5');
len=length(A);
i=1:len;
j=A(:,5);
figure(1);
plot(i,j,'.-');
0 Y" p3 L$ c! l7 |7 M  vhold on;
i=i';
! O- _: U4 S) M% _p1=polyfit(i,j,2)
f1=polyval(p1,i);
8 u7 O$ a( ?+ M$ u9 I% a" B8 D$ {plot(i,j,'.',i,f1,'-');
hold on;
title('病人平均等待时间与病人到达时间的关系');
8 T* h$ [% Z% {/ P+ e7 P' Fxlabel('病人门诊时间距离7 月12 日的时间差');ylabel('病人平均等待时间
');
axis equal;
7 Y. L$ t! _* Y7 @% @figure(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
$ L" ]# m, J9 i+ n9 g4 u! Z+ z. L32 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);
; L7 z" y0 @3 D: v  H8 Q& Rplot(A1,B1,'.-');
hold on;
p2=polyfit(A1,B1,2)
; ?8 _2 V" @8 }* ~1 _3 v- wf2=polyval(p2,A1);
0 P; ^# M& }/ O( R4 L; O1 Wplot(A1,B1,'.',A1,f2,'-');
! K) P- |" A; u  u, l$ Zhold on;
title('白内障病人情况');
8 O  Y, v. U' S2 r* ^! @# z% Q# zxlabel('白内障病人门诊时间距7月12日的时间差');ylabel('白内障病人平均
- V" }& Q2 T5 V7 e& u/ Q等待时间');
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
% `& }2 k; U1 f1 \. t2 U2 ~31 32 33 34 35 36 37 ];
. Y- {! K) c* D' c$ j8 HB2=[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];
8 k- Y9 ], q* O( A+ h7 Psubplot(1,2,2);
plot(A2,B2,'.-');
* o) Q# Y. k8 B1 N+ khold on;
p3=polyfit(A2,B2,2)
f3=polyval(p3,A2);
plot(A2,B2,'.',A2,f3,'-');
hold on;
title('白内障双眼病人情况');
xlabel('白内障双眼病人门诊时间距7月12日的时间差');ylabel('白内障双眼
病人平均等待时间');
axis equal;
& W. B: b. u2 c9 Q# R# J5 Ifigure(3);
" w, Q' ^# V+ @. B6 x  g+ n" {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
28 29 30 31 32 33 34 35 36 37 38 ];
# q; S4 z5 J) _2 H# XB3=[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];
subplot(2,2,1);
! Y) F7 ?8 M) T. j7 v5 b# dplot(A3,B3,'.-');
hold on;
) ]1 w, g. V" O: K2 j0 Rp4=polyfit(A3,B3,2)
% \6 C: ]- E- B- B' A$ f, hf4=polyval(p4,A3);
plot(A3,B3,'.',A3,f4,'-');
( N9 m8 N8 [* q" ~! {hold on;
* V& v% G0 I' {6 V- Vtitle('视网膜病人情况');
- q: k) i3 Q  R* \xlabel('门诊时间距7 月12 日的时间差');ylabel('视网膜病人平均等待时间
');
5 G" B, v8 [; r4 ?% faxis equal;
' `  w0 }5 {. b+ L! x5 c0 f( L9 rA4=[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
, v" \5 ^- q: t# l: \) @/ r37];
B4=[12 13 12 12 12 12 13.5 13 13.5 14 13 12 12 11 10.25 12 11 11 15 15
" O$ v" q5 S) n3 y$ a; s5 [4 T# X13 12 12 13 13 12];
subplot(2,2,2);
1 ~7 m1 }0 {% c+ \8 i+ m8 M1 {plot(A4,B4,'.-');
hold on;
p5=polyfit(A4,B4,2)
f5=polyval(p5,A4);
4 i+ n( X6 F% @5 Qplot(A4,B4,'.',A4,f5,'-');
( Q7 u4 }- p+ U. L( U  V4 Hhold on;
title('青光眼病人情况');
' ~5 T' z3 D  ~xlabel('门诊时间距7 月12日的时间差');ylabel('视网膜病人平均等待时间
');
axis equal;
2 j4 b5 n2 G+ MA4=[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];
3 A1 z1 [0 G9 |/ A2 x4 ^B4=[12 13 12 12 12 12 13.5 13 13.5 14 13 12 12 11 10.25 12 11 11 15 15
$ F, P( B8 K7 x3 X# }; Z4 H0 ?13 12 12 13 13 12];
' {1 @8 X2 b# l( D; C3 s1 v$ p9 Msubplot(2,2,2);
  D  b' u. t: ~% f5 [7 Xplot(A4,B4,'.-');
hold on;
p5=polyfit(A4,B4,2)
! W# [* B, F2 qf5=polyval(p5,A4);
plot(A4,B4,'.',A4,f5,'-');
hold on;
9 V% d# G0 \5 F& t& btitle('青光眼病人情况');
9 f3 v1 ^/ i, W9 e7 @7 Ixlabel('门诊时间距7 月12 日的时间差');ylabel('青光眼病人平均等待时间
- }5 k( Y+ q8 R0 {');
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];
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];
/ z" f  I; v' T  ]subplot(2,2,3);
6 g1 T4 i, W9 U2 Kplot(A5,B5,'.-');
. p+ F3 V7 }8 r6 z2 W! Xhold on;
p6=polyfit(A5,B5,2)
: l5 x  {+ K( n  ]& i& Uf6=polyval(p6,A5);
plot(A5,B5,'.',A5,f6,'-');
. [/ ]4 l& p% S% `( I( Uhold on;
3 B: @' z4 ]& r( C2 m1 V0 Rtitle('外伤病人情况');
xlabel('外伤病人门诊时间距离7月12日的时间差');ylabel('外伤病人平均等
待时间');
axis equal;
二、层次分析法的一致性指标判定的matlab 实现:
; K# k* t$ ~2 J9 g; ~+ {1 w' f) D: w+ _function H=ff(A)
( x: _5 U2 G9 Vt=eig(A)
" C* F/ E$ d# f0 F6 `  h1 pw=max(t)
5 d. K$ e: K# O8 c& `1 vq=(w-4)/3/.9
+ e, j+ g! o8 l3 E) @) cs=[(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(
: \8 D1 F% _+ L3 z8 |' T% ^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))
1 }' T+ z+ z% ^5 @# g7 u1 a(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)
1 w' R9 S* H! n+A(10)+A(11)+A(12)+A(13)+A(14)+A(15)+A(16))
: c% G( V: ~& z3 K  b0 h' u! c3 l" W(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))]
6 ]8 \% `  g5 X0 y* [" w+ ~& g- v

linmatsas

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

迷途的羊

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

linmatsas

回复 迷途的羊 的帖子
" A& v% U, e. d
# q0 q+ L& e7 U# E& E1 X$ f' w是自己编的还用别人运行么……不过要是楼主自己编的……小弟就佩服死啦
! j8 q3 i" ]9 }+ K" K   

紫辰

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

chenrye121

佩服两个字…………

水木年华zzu

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

wuming1429

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

hhao

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

abc321

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