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

wuming1429

一、使用曲线拟合模拟病人等待时间的模拟matlab程序实现：
[A]=xlsread('C:\MATLAB7\work\问题3数据5');
len=length(A);
i=1:len;
# M8 j) M9 g+ sj=A(:,5);
figure(1);
plot(i,j,'.-');
1 h: e+ U& x* Z3 b! l) S4 ghold on;
i=i';
+ z* V# I: h9 ]! ~$ T+ J6 D7 Kp1=polyfit(i,j,2)
+ D/ T4 h* @  p2 xf1=polyval(p1,i);
5 V2 L7 w' ]2 L4 {+ t; p& Vplot(i,j,'.',i,f1,'-');
6 c3 z' J3 {: F- ]+ h9 S3 `hold on;
title('病人平均等待时间与病人到达时间的关系');
xlabel('病人门诊时间距离7 月12 日的时间差');ylabel('病人平均等待时间
');
axis 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];
# T) |2 e) H% n! VB1=[12 12 12 12 12 12.5 13 13 13 13 13.4 14 14 13 12.5 11 11 12 12 12 11
, q. m) |  T% H: C) z" M. G11.5 13.25 14 13 12 11 13 13 12 12.25 13 15.5 15 14 13];
* s9 O" e/ ]7 d0 @subplot(1,2,1);
plot(A1,B1,'.-');
hold on;
1 F# ?+ ~' J& Q8 o. C% s- V* pp2=polyfit(A1,B1,2)
f2=polyval(p2,A1);
plot(A1,B1,'.',A1,f2,'-');
hold on;
% @1 W2 C' t. {title('白内障病人情况');
' C: i. N2 b. @# Oxlabel('白内障病人门诊时间距7月12日的时间差');ylabel('白内障病人平均
等待时间');
0 c; g/ i! k9 e# saxis equal;
, b2 }- K* H+ x: W& A4 A+ B0 _+ R! BA2=[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
, a; j7 r4 }  A3 U. L- q- |) Q31 32 33 34 35 36 37 ];
# t6 P; ^2 |' |1 w+ V) 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
12 12 13.25 14 14 13 12 12 12.4 13.33 13 13 13 12];
  K  c8 y. a# {  I8 z, asubplot(1,2,2);
$ n+ q+ D2 L1 L: F0 M# s: @plot(A2,B2,'.-');
! a2 |% F0 t! U4 W$ K, G6 Vhold on;
p3=polyfit(A2,B2,2)
' _# {2 Z, z; g2 Q+ Jf3=polyval(p3,A2);
. F& {. W# C) r6 [0 s1 E2 Tplot(A2,B2,'.',A2,f3,'-');
hold on;
3 Y9 t, B" S; j3 T) K1 \; p3 l$ Ptitle('白内障双眼病人情况');
- t' p$ ^% t3 O' E( c: m0 `- mxlabel('白内障双眼病人门诊时间距7月12日的时间差');ylabel('白内障双眼
病人平均等待时间');
axis equal;
# a% P7 V( H/ I# E5 l. j2 `& Nfigure(3);
' f9 j2 q" P! Z& q% l# YA3=[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 ];
+ n/ r- Z' f' iB3=[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
/ F; x( P) ~) K& m9 q11.67 12 12.25 12 11 11.75 13 14 14 13 12 12 12.33 13 13 13 12.5 12 12];
* R3 o2 ?' E1 F6 j4 tsubplot(2,2,1);
: ?0 V2 B; q7 M9 ?1 `4 Eplot(A3,B3,'.-');
* J' k" M4 c, F+ yhold on;
p4=polyfit(A3,B3,2)
+ B4 K* e8 t$ U9 ?f4=polyval(p4,A3);
. p/ E" y9 @4 z% s7 F3 fplot(A3,B3,'.',A3,f4,'-');
4 j4 b4 L8 h7 X+ M, B7 Khold on;
title('视网膜病人情况');
xlabel('门诊时间距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
1 Y7 g! s1 a( |3 ?9 ?% V! 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
' }# e, m( a  |13 12 12 13 13 12];
subplot(2,2,2);
plot(A4,B4,'.-');
hold on;
p5=polyfit(A4,B4,2)
& |7 R  y( _' U- ~: H8 K" Xf5=polyval(p5,A4);
, o2 R% h; H# y. }7 Wplot(A4,B4,'.',A4,f5,'-');
hold on;
3 I! W( F- }' U7 W  ftitle('青光眼病人情况');
xlabel('门诊时间距7 月12日的时间差');ylabel('视网膜病人平均等待时间
5 v! V3 R9 b0 r');
3 O  H& W2 n& R$ Y# Qaxis equal;
* E" j8 ?( _. m& H: Y0 HA4=[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];
( U9 K2 L, l( c8 i) ^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);
5 U) o/ D7 P; U! u2 ?plot(A4,B4,'.-');
hold on;
p5=polyfit(A4,B4,2)
f5=polyval(p5,A4);
plot(A4,B4,'.',A4,f5,'-');
3 t- t0 P$ q" q  t4 L& Rhold on;
. _6 O2 ^7 _4 Y9 N6 P# _8 Ytitle('青光眼病人情况');
xlabel('门诊时间距7 月12 日的时间差');ylabel('青光眼病人平均等待时间
');
axis equal;
6 j" v$ I9 ?' [6 V' e) bA5=[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,'.-');
hold on;
p6=polyfit(A5,B5,2)
f6=polyval(p6,A5);
4 X% {7 v& R# b* w( A/ Splot(A5,B5,'.',A5,f6,'-');
hold on;
title('外伤病人情况');
xlabel('外伤病人门诊时间距离7月12日的时间差');ylabel('外伤病人平均等
待时间');
axis equal;
二、层次分析法的一致性指标判定的matlab 实现:
  t  Y/ J8 M: \function H=ff(A)
$ I1 o% z& Y6 a% u( s$ G: B2 pt=eig(A)
w=max(t)
2 @  o: ~0 O- r4 j* Xq=(w-4)/3/.9
; ?* e5 ~1 b, r" ks=[(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(
! l" R5 G' p' J9)+A(10)+A(11)+A(12)+A(13)+A(14)+A(15)+A(16))
' M" h1 \, E. k; E3 C* k& U/ g4 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))
' V) F1 T, F& h(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)
) A9 B) J0 P, j7 r: r; r+A(10)+A(11)+A(12)+A(13)+A(14)+A(15)+A(16))
) o' q3 L  N* L(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)
4 A  I& k0 O+ K# J: S& i2 M4 L+A(10)+A(11)+A(12)+A(13)+A(14)+A(15)+A(16))]

linmatsas

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

迷途的羊

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

linmatsas

回复 迷途的羊 的帖子
( a  A( h& N, I+ z; j8 C
+ y$ p. R1 B6 i/ [是自己编的还用别人运行么……不过要是楼主自己编的……小弟就佩服死啦
/ T+ i2 h  u* {; Q   

紫辰

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

chenrye121

佩服两个字…………

水木年华zzu

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

wuming1429

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

hhao

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

abc321

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