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

wuming1429

一、使用曲线拟合模拟病人等待时间的模拟matlab程序实现：
8 G0 A/ h% @8 I[A]=xlsread('C:\MATLAB7\work\问题3数据5');
len=length(A);
i=1:len;
j=A(:,5);
figure(1);
+ A3 M( {. z0 Q8 s) Z' Q5 splot(i,j,'.-');
hold on;
i=i';
p1=polyfit(i,j,2)
/ F* L6 s- ?- t/ m/ A  @6 M3 W% {. mf1=polyval(p1,i);
3 o3 k8 z  b/ a; Tplot(i,j,'.',i,f1,'-');
hold on;
title('病人平均等待时间与病人到达时间的关系');
' Q: L% O0 X# yxlabel('病人门诊时间距离7 月12 日的时间差');ylabel('病人平均等待时间
');
axis equal;
figure(2);
) R' _+ c. T2 d' ]' A. `8 O8 g' pA1=[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
) ^1 p# M! k/ \$ {: V2 L3 {7 x6 w32 34 35 36 37 38 39 40 41 42];
3 K: z& J7 J6 i2 z' \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
4 Q, o6 w- i/ W11.5 13.25 14 13 12 11 13 13 12 12.25 13 15.5 15 14 13];
. g0 f: K+ o7 |subplot(1,2,1);
1 a% M1 L6 z: A/ x4 K% lplot(A1,B1,'.-');
hold on;
p2=polyfit(A1,B1,2)
6 i6 J3 V3 S* s+ lf2=polyval(p2,A1);
plot(A1,B1,'.',A1,f2,'-');
  r# R( D+ K4 n2 fhold on;
title('白内障病人情况');
xlabel('白内障病人门诊时间距7月12日的时间差');ylabel('白内障病人平均
: b' L8 a! m/ ?7 ~2 X' ^等待时间');
0 b' j0 K8 |  N7 ~2 Baxis equal;
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
2 j: f7 x, w. t) Z) o31 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 s7 J  \. p& |6 s12 12 13.25 14 14 13 12 12 12.4 13.33 13 13 13 12];
subplot(1,2,2);
plot(A2,B2,'.-');
+ A5 g% v' z, [! s- z) B% _hold on;
( N9 `% S" s8 pp3=polyfit(A2,B2,2)
8 B  @! r3 n0 l9 O: Uf3=polyval(p3,A2);
plot(A2,B2,'.',A2,f3,'-');
& \& J7 G" @$ R5 ?. _4 mhold on;
, R- w1 V, u( b; F' A$ ntitle('白内障双眼病人情况');
xlabel('白内障双眼病人门诊时间距7月12日的时间差');ylabel('白内障双眼
病人平均等待时间');
: [9 I! p$ d0 N1 [/ j1 j1 waxis 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
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
  V8 y! H4 Z5 h# Z11.67 12 12.25 12 11 11.75 13 14 14 13 12 12 12.33 13 13 13 12.5 12 12];
$ M9 ~4 p9 r" @/ {subplot(2,2,1);
plot(A3,B3,'.-');
hold on;
p4=polyfit(A3,B3,2)
( B2 L5 Q- W6 F0 [' ef4=polyval(p4,A3);
plot(A3,B3,'.',A3,f4,'-');
hold on;
title('视网膜病人情况');
/ x$ d1 \; g& I" P& lxlabel('门诊时间距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
8 X7 J! a1 w+ |9 c, V37];
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);
" R2 n$ p! i  pplot(A4,B4,'.-');
hold on;
$ }. J5 [4 @0 Y1 ^% `p5=polyfit(A4,B4,2)
# [% @" r8 H3 j; d9 g' _! x' Kf5=polyval(p5,A4);
plot(A4,B4,'.',A4,f5,'-');
3 K! W* q/ l) O0 U4 Phold on;
title('青光眼病人情况');
' ]+ P9 ^7 \# axlabel('门诊时间距7 月12日的时间差');ylabel('视网膜病人平均等待时间
3 X3 W% ^& }' d( W4 F');
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
37];
/ @9 y: t. N8 R3 C) S# BB4=[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);
8 t6 v. n1 F' t/ p3 L' [3 F: jplot(A4,B4,'.-');
hold on;
; g# }1 R5 w5 B+ g  R7 L* d9 J- pp5=polyfit(A4,B4,2)
f5=polyval(p5,A4);
plot(A4,B4,'.',A4,f5,'-');
; V4 W- G9 ^6 r' m( phold on;
8 o3 o- X+ _/ E6 P- Q* z$ A: }; ltitle('青光眼病人情况');
xlabel('门诊时间距7 月12 日的时间差');ylabel('青光眼病人平均等待时间
/ J' x( j: X3 J7 U. {');
7 Z  \9 a# v/ }axis equal;
9 m% G/ B1 n- K7 H1 |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];
1 \2 r+ G/ m' p' ZB5=[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];
# R* [3 o: Z- Y" Tsubplot(2,2,3);
plot(A5,B5,'.-');
hold on;
p6=polyfit(A5,B5,2)
' \* R' h8 j  of6=polyval(p6,A5);
* e* ?0 o+ L+ C& t$ |plot(A5,B5,'.',A5,f6,'-');
6 y8 h8 l) z! v6 q7 rhold on;
3 B) O7 q" X* Jtitle('外伤病人情况');
4 e+ {  _! A; q; l; vxlabel('外伤病人门诊时间距离7月12日的时间差');ylabel('外伤病人平均等
' [& R* z2 ~; v0 F1 c2 p8 i; }* W待时间');
axis equal;
二、层次分析法的一致性指标判定的matlab 实现:
function H=ff(A)
9 J- y# ?5 l4 @! U& Q4 I% Ft=eig(A)
  h4 k" K+ J% X6 Mw=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(
0 Z4 s* v: i% v# e9)+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)
+A(10)+A(11)+A(12)+A(13)+A(14)+A(15)+A(16))
( j& l0 s3 n4 u- w(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))
(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))]

linmatsas

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

迷途的羊

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

linmatsas

回复 迷途的羊 的帖子

是自己编的还用别人运行么……不过要是楼主自己编的……小弟就佩服死啦
6 ^9 A8 y" ]- X7 K. K- Q   

紫辰

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

chenrye121

佩服两个字…………

水木年华zzu

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

wuming1429

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

hhao

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

abc321

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