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

wuming1429

一、使用曲线拟合模拟病人等待时间的模拟matlab程序实现：
[A]=xlsread('C:\MATLAB7\work\问题3数据5');
len=length(A);
i=1:len;
/ ]7 o: i* [# z7 l( _j=A(:,5);
- s6 F# t- u/ y! C. \" ?/ Kfigure(1);
plot(i,j,'.-');
hold on;
i=i';
p1=polyfit(i,j,2)
f1=polyval(p1,i);
" ]$ p  H6 J7 W7 _8 q- s7 u& wplot(i,j,'.',i,f1,'-');
/ ?0 p9 m* @* z  Qhold on;
title('病人平均等待时间与病人到达时间的关系');
0 f. {4 k+ v. d# Q& y( k, M3 R# i; Uxlabel('病人门诊时间距离7 月12 日的时间差');ylabel('病人平均等待时间
: i; m& V$ X( F/ D');
6 g% c9 ]) N/ U0 m+ b! Xaxis equal;
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
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
11.5 13.25 14 13 12 11 13 13 12 12.25 13 15.5 15 14 13];
subplot(1,2,1);
plot(A1,B1,'.-');
3 K/ o3 f8 L1 z) x3 Z! vhold on;
p2=polyfit(A1,B1,2)
f2=polyval(p2,A1);
- L. o4 M- g1 K. \' {plot(A1,B1,'.',A1,f2,'-');
: m3 E7 d" j+ e, Q+ g- a* Uhold on;
7 u. J2 ]& X6 ?" r  _$ L. ?title('白内障病人情况');
9 O# T5 f2 w; _7 Yxlabel('白内障病人门诊时间距7月12日的时间差');ylabel('白内障病人平均
9 j- T% y6 u/ K2 V9 y等待时间');
" b  s$ N. p- |- ^1 Z1 Qaxis equal;
8 E8 D( y- X  [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 ];
. |% P  D' s) d" ~, |" qB2=[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
+ {! b; a& g1 N2 n12 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;
p3=polyfit(A2,B2,2)
f3=polyval(p3,A2);
9 y" w/ q/ ?! nplot(A2,B2,'.',A2,f3,'-');
hold on;
' q5 A1 ~3 g0 R/ Vtitle('白内障双眼病人情况');
xlabel('白内障双眼病人门诊时间距7月12日的时间差');ylabel('白内障双眼
8 h: z0 Q4 u# }# v. ~; A) i- f病人平均等待时间');
axis equal;
2 o/ x# Y" n, p0 ?; X% ^, ifigure(3);
0 P, z9 `7 u! l3 v2 t5 I; i2 l( [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 ];
  ^" l, q# k! ]* P% r2 rB3=[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];
9 n6 P: T/ q' W: B  \subplot(2,2,1);
5 V$ U( e" l- O, [( Splot(A3,B3,'.-');
! _# a2 F" H) t0 O( f7 ahold on;
p4=polyfit(A3,B3,2)
- d( o4 f% B. j5 Vf4=polyval(p4,A3);
plot(A3,B3,'.',A3,f4,'-');
! Y- u/ p: y: f5 n' o1 }! Xhold on;
8 L+ a( {  Y) jtitle('视网膜病人情况');
xlabel('门诊时间距7 月12 日的时间差');ylabel('视网膜病人平均等待时间
7 _( t5 G+ `0 s: x- {; e7 W8 S');
axis equal;
, u: r. T. Q; U- c' ~+ ^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
* p, V( C  k) `, z' g37];
B4=[12 13 12 12 12 12 13.5 13 13.5 14 13 12 12 11 10.25 12 11 11 15 15
: y* F' X; U6 K4 X& ?13 12 12 13 13 12];
subplot(2,2,2);
plot(A4,B4,'.-');
hold on;
, e6 u; U$ M" O* Tp5=polyfit(A4,B4,2)
f5=polyval(p5,A4);
' S3 o; u9 b/ A( K7 N% d* }' b1 xplot(A4,B4,'.',A4,f5,'-');
4 D4 |/ q; \% H: n) Lhold on;
title('青光眼病人情况');
, v: x0 l0 P# G# w6 Bxlabel('门诊时间距7 月12日的时间差');ylabel('视网膜病人平均等待时间
');
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
+ }" u2 [2 `, E: F2 h& o- r% q37];
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];
7 H* {5 {, E: c2 d& B, T/ `subplot(2,2,2);
plot(A4,B4,'.-');
hold on;
p5=polyfit(A4,B4,2)
; X( }) V9 Q% U' ]$ vf5=polyval(p5,A4);
plot(A4,B4,'.',A4,f5,'-');
5 E3 c1 R+ Q9 {3 B& J7 w/ Fhold on;
( ^" G- I) n3 r- Z$ y8 W6 i2 K; `title('青光眼病人情况');
9 i5 B' C% r# m0 \3 L# Yxlabel('门诊时间距7 月12 日的时间差');ylabel('青光眼病人平均等待时间
$ t4 x+ f1 I/ x$ U; t: _" \4 v  w');
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];
: b8 R. r6 ]2 A' m' KB5=[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];
2 i% L( \) }: K( zsubplot(2,2,3);
6 o# H/ L' Y/ `7 K' g8 Vplot(A5,B5,'.-');
  L4 `: n. W2 Thold on;
% j: T$ d0 K' e0 D! e8 U' dp6=polyfit(A5,B5,2)
! V2 W& I5 Z; P" ff6=polyval(p6,A5);
plot(A5,B5,'.',A5,f6,'-');
hold on;
title('外伤病人情况');
3 J1 B1 N& S: U0 `! Q! |xlabel('外伤病人门诊时间距离7月12日的时间差');ylabel('外伤病人平均等
' J- j4 h4 g1 D" I9 w0 C8 z待时间');
axis equal;
二、层次分析法的一致性指标判定的matlab 实现:
function H=ff(A)
t=eig(A)
) o) s& R! }$ ?w=max(t)
q=(w-4)/3/.9
  N/ u$ ~* f1 J" j, xs=[(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(
  g8 J% O2 b" Y6 G$ B( s2 z; ]8 v9)+A(10)+A(11)+A(12)+A(13)+A(14)+A(15)+A(16))
$ W2 E6 s3 G  I" x3 d7 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)
+A(10)+A(11)+A(12)+A(13)+A(14)+A(15)+A(16))
/ o% `4 c7 ~3 }- k(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))
! M$ ^8 S$ c, e(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))]
- ^' f; ?3 j8 ^' B1 e" m# k

linmatsas

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

迷途的羊

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

linmatsas

回复 迷途的羊 的帖子

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

紫辰

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

chenrye121

佩服两个字…………

水木年华zzu

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

wuming1429

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

hhao

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

abc321

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