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

wuming1429

一、使用曲线拟合模拟病人等待时间的模拟matlab程序实现：
( R$ f% w2 K2 L[A]=xlsread('C:\MATLAB7\work\问题3数据5');
7 i+ {! x. Y2 o: T/ y* r* G8 ]( x% Wlen=length(A);
0 o9 Q7 B/ j, Q* S7 bi=1:len;
j=A(:,5);
+ K: D8 K) j9 g9 Q& \& l4 {figure(1);
plot(i,j,'.-');
: a/ k0 d7 c' s7 z/ ~' J# m0 Phold on;
' a1 t+ l* L" I# j0 ei=i';
# N  y, v/ J% ~; N8 V, |p1=polyfit(i,j,2)
f1=polyval(p1,i);
# y9 W0 W, G( H) ~0 [& Dplot(i,j,'.',i,f1,'-');
hold on;
# s3 I8 A8 ?4 utitle('病人平均等待时间与病人到达时间的关系');
( [- h- c. R! a) Y! Z9 e- mxlabel('病人门诊时间距离7 月12 日的时间差');ylabel('病人平均等待时间
');
axis equal;
- U4 S3 C! k' jfigure(2);
+ p- p7 E, `5 i  p3 c2 ZA1=[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
, S3 G, u; b2 D$ ?11.5 13.25 14 13 12 11 13 13 12 12.25 13 15.5 15 14 13];
subplot(1,2,1);
( P0 o/ t; a! qplot(A1,B1,'.-');
1 T4 S, O' m' |! [  H, vhold on;
0 O/ W: R: N0 L5 s, O- A$ Qp2=polyfit(A1,B1,2)
f2=polyval(p2,A1);
8 I1 \$ _" S& G5 oplot(A1,B1,'.',A1,f2,'-');
hold on;
title('白内障病人情况');
xlabel('白内障病人门诊时间距7月12日的时间差');ylabel('白内障病人平均
1 [% [$ Q, W. d- m1 {9 C等待时间');
' A+ d8 E9 L2 s1 [1 q  R# daxis equal;
1 z4 k4 t+ y4 Y$ Q' ~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
0 R1 n  g, m8 a* [( z' q/ O4 l31 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
4 N8 F, O: @9 {$ b* g' @! V; i12 12 13.25 14 14 13 12 12 12.4 13.33 13 13 13 12];
' W$ X6 y1 D' o3 w" E0 nsubplot(1,2,2);
6 e: Z: E5 V5 y0 T3 \plot(A2,B2,'.-');
hold on;
2 D# ~7 r" Q$ \p3=polyfit(A2,B2,2)
' p" R* }, {) N4 }# I, Cf3=polyval(p3,A2);
1 M1 y1 Z6 l  {plot(A2,B2,'.',A2,f3,'-');
, W& H; _7 P# Vhold on;
title('白内障双眼病人情况');
xlabel('白内障双眼病人门诊时间距7月12日的时间差');ylabel('白内障双眼
; w, W, M$ ]( X. p病人平均等待时间');
4 y3 [$ b: P) V( q) c: Baxis equal;
figure(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
' n8 R. r1 Q9 V$ ?/ j: o28 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];
% M5 d. x: s( K) F: lsubplot(2,2,1);
plot(A3,B3,'.-');
hold on;
7 n8 W: h7 }  h- k! j/ H# R& V+ H; \p4=polyfit(A3,B3,2)
f4=polyval(p4,A3);
plot(A3,B3,'.',A3,f4,'-');
+ Q( O3 Q! n. h6 W/ U, `hold on;
title('视网膜病人情况');
xlabel('门诊时间距7 月12 日的时间差');ylabel('视网膜病人平均等待时间
9 M% E2 [. l  z' a" Q% ^');
axis equal;
  I; K$ U$ p3 G3 W: T. [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
/ _; Y( }' ~& D) ]  b+ _* M37];
' ?* ~4 I) \( o7 w7 u+ hB4=[12 13 12 12 12 12 13.5 13 13.5 14 13 12 12 11 10.25 12 11 11 15 15
6 B+ t) r( u, |% H/ ]) [# S- M13 12 12 13 13 12];
subplot(2,2,2);
plot(A4,B4,'.-');
hold on;
" {5 \2 U) m" `, }# o7 l7 q# Y5 Yp5=polyfit(A4,B4,2)
f5=polyval(p5,A4);
plot(A4,B4,'.',A4,f5,'-');
hold on;
title('青光眼病人情况');
) X: R( m5 I  n/ Hxlabel('门诊时间距7 月12日的时间差');ylabel('视网膜病人平均等待时间
7 C5 A: U7 }' q* ?; q');
) P8 {) M6 c" W# t! baxis equal;
% k- T- t3 J' ]8 R3 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];
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];
+ T7 J( }! F( j  c1 wsubplot(2,2,2);
/ d" H9 R% r, r9 |- E! z8 i! cplot(A4,B4,'.-');
hold on;
/ q7 m8 B# t4 s9 Q' a7 ?+ P$ ?p5=polyfit(A4,B4,2)
f5=polyval(p5,A4);
5 f  m" V9 L  W' G& l) B* Qplot(A4,B4,'.',A4,f5,'-');
, f& \; S+ K1 d3 b9 n4 S  Ehold on;
title('青光眼病人情况');
xlabel('门诊时间距7 月12 日的时间差');ylabel('青光眼病人平均等待时间
');
3 i% U4 f+ z/ m, W$ P; ?. baxis 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];
& q8 {. q) g" h7 t2 CB5=[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];
+ W* ?0 L& D7 z0 o# d$ n/ xsubplot(2,2,3);
( c. v- g, f4 ^9 Jplot(A5,B5,'.-');
0 P# ]2 }6 ]7 f8 r/ Q1 Mhold on;
p6=polyfit(A5,B5,2)
3 g2 @( c  b, `( ^) Y7 @  D$ X8 hf6=polyval(p6,A5);
plot(A5,B5,'.',A5,f6,'-');
% o6 L3 f( t3 @6 X& W) Thold on;
title('外伤病人情况');
. p9 L: P1 |/ C1 Z5 R; ]xlabel('外伤病人门诊时间距离7月12日的时间差');ylabel('外伤病人平均等
6 W: P/ u, |! w" R1 i5 Z% y) x& X, t6 z待时间');
axis equal;
. U# M0 D+ H+ o  Q7 {二、层次分析法的一致性指标判定的matlab 实现:
  _; o" ^- C, T4 f2 W9 `: f7 ^; Kfunction H=ff(A)
t=eig(A)
w=max(t)
q=(w-4)/3/.9
# T. M# D% g2 f" Z" Bs=[(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))
(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)
! [' l0 E2 u8 c" L. f4 ~2 [% |+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)
/ z5 L3 [% N$ O& F, S" \# o: d4 t8 h+A(10)+A(11)+A(12)+A(13)+A(14)+A(15)+A(16))
4 v1 n4 C# F) ~8 z% D  I(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)
- Z# Q* j& R; T. x! \7 V+A(10)+A(11)+A(12)+A(13)+A(14)+A(15)+A(16))]
+ a0 L/ D3 d! P' d

linmatsas

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

迷途的羊

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

linmatsas

回复 迷途的羊 的帖子
* N% Q) [( j- X
是自己编的还用别人运行么……不过要是楼主自己编的……小弟就佩服死啦
   

紫辰

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

chenrye121

佩服两个字…………

水木年华zzu

[A]=xlsread('C:\MATLAB7\work\问题3数据5');
. U( v- X6 K( O6 g+ m+ t6 p3 K在不同的地方运行上面一句要改成本地数据保存的路径的

wuming1429

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

hhao

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

abc321

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