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

wuming1429

一、使用曲线拟合模拟病人等待时间的模拟matlab程序实现：
" A7 P8 y) q8 t9 F' I[A]=xlsread('C:\MATLAB7\work\问题3数据5');
- G! W6 o* p( n* n8 A6 rlen=length(A);
i=1:len;
; ~% q( v$ ]. V2 {j=A(:,5);
figure(1);
plot(i,j,'.-');
hold on;
' e+ C7 `& ^9 ii=i';
p1=polyfit(i,j,2)
! ^8 G- c  q+ e  [* Gf1=polyval(p1,i);
5 m- M$ O' O0 I, [plot(i,j,'.',i,f1,'-');
hold on;
title('病人平均等待时间与病人到达时间的关系');
xlabel('病人门诊时间距离7 月12 日的时间差');ylabel('病人平均等待时间
: O5 B* s; D* R( j  x1 o');
axis equal;
3 y1 F% S& c( O' Ffigure(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
% o! O* U  `: B11.5 13.25 14 13 12 11 13 13 12 12.25 13 15.5 15 14 13];
& T* R' W4 e4 v9 |3 D$ k9 zsubplot(1,2,1);
3 W* R+ x9 o/ }" J: oplot(A1,B1,'.-');
6 \' e% c7 a3 b0 f4 uhold on;
1 j) N* c9 }: t8 e! t9 I: Dp2=polyfit(A1,B1,2)
3 P5 s/ Y: n3 T& p5 If2=polyval(p2,A1);
! @+ O  D3 W% Eplot(A1,B1,'.',A1,f2,'-');
; U/ n3 m; L$ \1 ~. c7 ?hold on;
title('白内障病人情况');
xlabel('白内障病人门诊时间距7月12日的时间差');ylabel('白内障病人平均
等待时间');
& Q5 j; I5 u, H# {  Vaxis equal;
! k: T( f+ K6 X- n( IA2=[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
' ?' w1 x& C. ?) R6 Z1 x31 32 33 34 35 36 37 ];
2 K9 M; R. M+ @) h  _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
12 12 13.25 14 14 13 12 12 12.4 13.33 13 13 13 12];
subplot(1,2,2);
5 M' e( x2 \8 s# c; ?  X/ A# Yplot(A2,B2,'.-');
9 D* p$ A' z8 A7 }8 E: Hhold on;
; a1 G9 T' Q; K+ O! Z+ Xp3=polyfit(A2,B2,2)
( ^6 F7 ~6 A7 Mf3=polyval(p3,A2);
% h- R) H# ~! Q* Z# ~1 ?plot(A2,B2,'.',A2,f3,'-');
hold on;
title('白内障双眼病人情况');
xlabel('白内障双眼病人门诊时间距7月12日的时间差');ylabel('白内障双眼
病人平均等待时间');
* b% ?( u2 ?0 j( \4 N: ^9 h) T6 m7 `axis 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
' V: s  T5 M* f: f+ g" L28 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
8 m. ~3 I2 z4 m2 O1 F* x4 k11.67 12 12.25 12 11 11.75 13 14 14 13 12 12 12.33 13 13 13 12.5 12 12];
/ `$ i  ?, S; T4 h2 X) o) b# ksubplot(2,2,1);
plot(A3,B3,'.-');
1 E! t" S+ c! h! e( B% uhold on;
2 G+ j$ o% o4 ]$ F. Tp4=polyfit(A3,B3,2)
" P1 Q1 x5 g1 tf4=polyval(p4,A3);
! w& E* Y( F8 O  B' iplot(A3,B3,'.',A3,f4,'-');
- f5 Y  E' S% y8 m2 k! u2 [hold on;
title('视网膜病人情况');
xlabel('门诊时间距7 月12 日的时间差');ylabel('视网膜病人平均等待时间
: h  G: Y( M9 {$ F');
* v7 T' [0 E' X2 Oaxis 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];
, }9 h1 ^$ V7 s4 AB4=[12 13 12 12 12 12 13.5 13 13.5 14 13 12 12 11 10.25 12 11 11 15 15
4 ^/ h: T8 C1 m% D13 12 12 13 13 12];
subplot(2,2,2);
plot(A4,B4,'.-');
hold on;
p5=polyfit(A4,B4,2)
; w' b* U/ S! i/ @0 jf5=polyval(p5,A4);
plot(A4,B4,'.',A4,f5,'-');
: g8 Q" n, [% B! Q: A& }0 |hold on;
0 y: Q* ^/ C0 w+ o4 p( X( Htitle('青光眼病人情况');
2 `( Q) N' R+ hxlabel('门诊时间距7 月12日的时间差');ylabel('视网膜病人平均等待时间
1 K  J. ~  D5 K. j7 {');
1 L4 [2 k# v& ^. oaxis 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];
; w; O% S, K4 ]+ H3 |( `5 RB4=[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];
0 [2 E0 q9 i0 I" j0 a* msubplot(2,2,2);
5 D" g5 x4 Y2 [- z% \% Xplot(A4,B4,'.-');
% {# _; P9 ^/ [9 x6 jhold on;
& F; I6 F4 l; ip5=polyfit(A4,B4,2)
f5=polyval(p5,A4);
plot(A4,B4,'.',A4,f5,'-');
hold on;
title('青光眼病人情况');
- P3 _  w( o6 d" n; {* A/ Q) V% Hxlabel('门诊时间距7 月12 日的时间差');ylabel('青光眼病人平均等待时间
');
3 I9 q; [+ Y7 }; V7 |axis equal;
' q) N# K+ a7 D+ ~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];
# W2 g7 [) @, m+ t( \6 Jsubplot(2,2,3);
plot(A5,B5,'.-');
0 \1 f! g8 s- X7 Y$ j3 c3 ~$ P0 Y- Khold on;
# H" O. |3 Y3 k6 _p6=polyfit(A5,B5,2)
f6=polyval(p6,A5);
3 ?2 m& ^( r# A3 d* \( C  d' x! }plot(A5,B5,'.',A5,f6,'-');
hold on;
# I* q% a# F7 T  U4 t' @title('外伤病人情况');
xlabel('外伤病人门诊时间距离7月12日的时间差');ylabel('外伤病人平均等
$ Y, A8 i! q7 N0 W" J4 j待时间');
axis equal;
! n* N# i5 r7 s+ k) g# Z" G# W二、层次分析法的一致性指标判定的matlab 实现:
1 L6 M3 u7 |. Y6 v" x' xfunction H=ff(A)
* i$ _9 S4 z" I& O, ]t=eig(A)
w=max(t)
q=(w-4)/3/.9
# z2 y: f) Q2 n" O" j" Q5 r$ G- X& {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(
2 H2 Y+ s0 y& ?$ z% L  S$ m7 Z3 ^) A$ z9)+A(10)+A(11)+A(12)+A(13)+A(14)+A(15)+A(16))
+ W- `! |7 O& q" A+ c5 U; \6 U  g(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)
+ n$ |8 Y1 t% l% @' C9 x+A(10)+A(11)+A(12)+A(13)+A(14)+A(15)+A(16))
; r/ Q4 b  G6 Y0 Y' e2 ^(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))]
! L2 F' {* r9 ]. L2 K4 @. t

linmatsas

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

迷途的羊

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

linmatsas

回复 迷途的羊 的帖子
2 J9 J! R1 t3 D
. U! n; U; r3 A0 T5 P& U是自己编的还用别人运行么……不过要是楼主自己编的……小弟就佩服死啦
   

紫辰

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

chenrye121

佩服两个字…………

水木年华zzu

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

wuming1429

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

hhao

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

abc321

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