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

wuming1429

一、使用曲线拟合模拟病人等待时间的模拟matlab程序实现：
- ]* `! L* u, f* O7 l1 L[A]=xlsread('C:\MATLAB7\work\问题3数据5');
len=length(A);
, s9 n" L, z) }' n1 }$ ?i=1:len;
  \: S- F, a) B6 ]7 ~j=A(:,5);
figure(1);
( V- q, P8 J- |: p/ [' n! N- Nplot(i,j,'.-');
& I; V7 z/ I0 V2 D" Bhold on;
7 v. U( ?8 n: Hi=i';
8 E$ U# ]1 d$ |3 cp1=polyfit(i,j,2)
. u; y2 C  r  J* }4 S* \f1=polyval(p1,i);
plot(i,j,'.',i,f1,'-');
9 N" C: A+ |0 bhold on;
6 ?# o8 R# S' D7 X+ G3 V. W2 Ctitle('病人平均等待时间与病人到达时间的关系');
4 |6 B5 G9 N% @4 ]$ ^" qxlabel('病人门诊时间距离7 月12 日的时间差');ylabel('病人平均等待时间
');
; i* h' S/ x# d4 daxis equal;
figure(2);
9 j4 X; y3 a' [1 P* o& M1 {. ]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];
7 n0 e/ L! F' [4 \1 x  r# S$ KB1=[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];
0 D% [+ Y& j6 C4 m1 f. C' ?1 `subplot(1,2,1);
5 y% G, l, A4 B, f4 w% fplot(A1,B1,'.-');
hold on;
p2=polyfit(A1,B1,2)
- ~1 U7 C. s! `- j+ B  S# If2=polyval(p2,A1);
  Q/ v' p" Y) {* uplot(A1,B1,'.',A1,f2,'-');
hold on;
) w% e4 J. h3 ?, R; Ltitle('白内障病人情况');
xlabel('白内障病人门诊时间距7月12日的时间差');ylabel('白内障病人平均
等待时间');
7 P2 C8 s- v# |  @# Aaxis equal;
) g: B0 K7 l$ y! Q* QA2=[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
# L" d* ?/ v. n) F31 32 33 34 35 36 37 ];
B2=[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
, h. y  K5 q& s( o12 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;
. g- b$ b  R' }! y/ e; q; Yp3=polyfit(A2,B2,2)
f3=polyval(p3,A2);
plot(A2,B2,'.',A2,f3,'-');
hold on;
title('白内障双眼病人情况');
/ ~' X% Z+ U& x0 n0 lxlabel('白内障双眼病人门诊时间距7月12日的时间差');ylabel('白内障双眼
病人平均等待时间');
axis equal;
% E) X, S3 L; t4 T: e; z8 N/ rfigure(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
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
11.67 12 12.25 12 11 11.75 13 14 14 13 12 12 12.33 13 13 13 12.5 12 12];
/ S$ t" z# w: ^+ O  E% C9 osubplot(2,2,1);
plot(A3,B3,'.-');
/ N$ C; a; B% A% k1 m* j# f8 T& F+ @) ^hold on;
p4=polyfit(A3,B3,2)
! g8 `! o/ B: ^8 r$ O+ q/ hf4=polyval(p4,A3);
- o& Q- U/ B4 H( @& q, l, Cplot(A3,B3,'.',A3,f4,'-');
0 v; i) J; H# e( O, v8 W( shold on;
" e$ B4 I6 G$ h' Z- otitle('视网膜病人情况');
6 N' @  M( {; Axlabel('门诊时间距7 月12 日的时间差');ylabel('视网膜病人平均等待时间
5 |6 d) @. o& x% K8 ?');
" y6 {0 t) R- q0 ?) aaxis equal;
* I* e- R0 j  e) t; N% r) P' g( L5 JA4=[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];
& g7 ?/ B/ J" d, WB4=[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);
plot(A4,B4,'.-');
hold on;
9 k$ l2 `$ K. I$ a$ w: Gp5=polyfit(A4,B4,2)
7 ^% i6 ^* V; U4 j. n8 ?( Kf5=polyval(p5,A4);
& |! |5 G+ |2 u# wplot(A4,B4,'.',A4,f5,'-');
hold on;
title('青光眼病人情况');
/ _3 m" ?$ z0 q6 L& m3 t0 zxlabel('门诊时间距7 月12日的时间差');ylabel('视网膜病人平均等待时间
');
0 p$ V: v7 a, k$ M# D+ Kaxis 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
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
  ~; j& `) B7 u/ [13 12 12 13 13 12];
subplot(2,2,2);
  b5 J6 S. N( F& ^7 |) eplot(A4,B4,'.-');
hold on;
- `5 y( e3 F) l6 Z% z4 Yp5=polyfit(A4,B4,2)
8 w& y' N, I* y6 c1 Z& if5=polyval(p5,A4);
plot(A4,B4,'.',A4,f5,'-');
hold on;
* K- u% A* J4 V  U0 b9 [title('青光眼病人情况');
1 g; M" c2 H& R1 B" k" rxlabel('门诊时间距7 月12 日的时间差');ylabel('青光眼病人平均等待时间
');
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];
subplot(2,2,3);
plot(A5,B5,'.-');
8 {/ l4 p- ]2 c, M2 _7 Ghold on;
p6=polyfit(A5,B5,2)
8 q# ^! R1 |; nf6=polyval(p6,A5);
  I, L/ ]+ s( W! n4 K: T5 q" B/ Pplot(A5,B5,'.',A5,f6,'-');
( w" z+ {2 h% E; ?. F0 _2 hhold on;
title('外伤病人情况');
3 d, Z% e# q3 O( \8 N4 ?xlabel('外伤病人门诊时间距离7月12日的时间差');ylabel('外伤病人平均等
' h* ~4 A% C# W) J' @待时间');
axis equal;
二、层次分析法的一致性指标判定的matlab 实现:
3 e8 F7 `7 v0 q  yfunction H=ff(A)
2 U9 W6 i* a0 M# P3 R: Q/ nt=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))
; s# i+ Z" H, R  F" T(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)
+A(10)+A(11)+A(12)+A(13)+A(14)+A(15)+A(16))
  l0 b2 W0 |: O  y( g(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)
0 y# t& }& k; n: ^$ d' A, }+A(10)+A(11)+A(12)+A(13)+A(14)+A(15)+A(16))]

linmatsas

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

迷途的羊

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

linmatsas

回复 迷途的羊 的帖子

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

紫辰

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

chenrye121

佩服两个字…………

水木年华zzu

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

wuming1429

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

hhao

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

abc321

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