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

wuming1429

一、使用曲线拟合模拟病人等待时间的模拟matlab程序实现：
[A]=xlsread('C:\MATLAB7\work\问题3数据5');
len=length(A);
i=1:len;
j=A(:,5);
. l; P) \" B5 ^5 y, R' Rfigure(1);
plot(i,j,'.-');
0 l# y: H* X# ^4 v0 E+ u; a, L( @1 Whold on;
; L  S. a9 _9 B! R) ?( xi=i';
8 K0 O( S/ }* ]p1=polyfit(i,j,2)
f1=polyval(p1,i);
plot(i,j,'.',i,f1,'-');
hold on;
% Y! @4 m1 N* C; Q/ _7 Ltitle('病人平均等待时间与病人到达时间的关系');
0 ^6 J' H+ g, ?2 Z) m& |# s' q4 Hxlabel('病人门诊时间距离7 月12 日的时间差');ylabel('病人平均等待时间
');
axis equal;
! ?% G, S$ \+ O( |1 dfigure(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];
' i4 g3 v. K+ d- ^/ M8 X7 xB1=[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];
; y% o' Q# i  s1 G: ysubplot(1,2,1);
plot(A1,B1,'.-');
hold on;
p2=polyfit(A1,B1,2)
f2=polyval(p2,A1);
plot(A1,B1,'.',A1,f2,'-');
hold on;
2 K9 ~5 L& Y2 i+ m" Ititle('白内障病人情况');
. \, M1 G1 ?8 k% J3 r5 T& Dxlabel('白内障病人门诊时间距7月12日的时间差');ylabel('白内障病人平均
- U+ X% s! R: _; s! f9 j4 ?等待时间');
, ~. X5 o/ @" g0 M" i7 Q; \+ h, }; aaxis 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
31 32 33 34 35 36 37 ];
, N& Z- b5 o/ G4 O$ V/ PB2=[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
* X& H! y: ~5 X4 M1 z* N! F& V* d12 12 13.25 14 14 13 12 12 12.4 13.33 13 13 13 12];
subplot(1,2,2);
" t! F+ ~7 L# C$ l& M" Wplot(A2,B2,'.-');
+ L2 I& q+ G' U3 j1 ^7 V$ f6 }4 `hold on;
& }% D. T3 [) z9 ]. Cp3=polyfit(A2,B2,2)
f3=polyval(p3,A2);
, k4 u& R) N" }. N  D- U( Jplot(A2,B2,'.',A2,f3,'-');
hold on;
title('白内障双眼病人情况');
  C7 J( e8 q3 g- J/ I6 v0 y4 e- L7 G/ _xlabel('白内障双眼病人门诊时间距7月12日的时间差');ylabel('白内障双眼
病人平均等待时间');
axis equal;
! y6 W" T: k( `% n, Pfigure(3);
* ]' y. Y* H& ^  MA3=[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
$ a) K2 L8 f! [9 w28 29 30 31 32 33 34 35 36 37 38 ];
7 W6 o0 E- w& w6 F- fB3=[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
0 b- `. b* V# g7 Z11.67 12 12.25 12 11 11.75 13 14 14 13 12 12 12.33 13 13 13 12.5 12 12];
subplot(2,2,1);
plot(A3,B3,'.-');
hold on;
p4=polyfit(A3,B3,2)
f4=polyval(p4,A3);
9 {% U' [3 x( X+ Eplot(A3,B3,'.',A3,f4,'-');
hold on;
; K0 f, n9 l6 t$ O% G- C, P4 b- J+ `title('视网膜病人情况');
xlabel('门诊时间距7 月12 日的时间差');ylabel('视网膜病人平均等待时间
');
axis equal;
. b/ p% _- n- H8 w1 _' t2 eA4=[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
! B) p  P, ?% |2 b1 G& {6 F37];
B4=[12 13 12 12 12 12 13.5 13 13.5 14 13 12 12 11 10.25 12 11 11 15 15
' {1 @/ J1 T0 ?! R6 F13 12 12 13 13 12];
% i: R$ `) _2 R* N! q2 csubplot(2,2,2);
/ H2 S2 V$ ^1 u/ }' L, F) a6 Splot(A4,B4,'.-');
hold on;
$ \, S: o4 A' M' E- yp5=polyfit(A4,B4,2)
5 t7 K* R9 [! }5 Xf5=polyval(p5,A4);
! ~: {/ V) C- U& [/ Q# C2 D! jplot(A4,B4,'.',A4,f5,'-');
hold on;
title('青光眼病人情况');
xlabel('门诊时间距7 月12日的时间差');ylabel('视网膜病人平均等待时间
- m6 N# b+ z( T% C; G! G');
% O; H' }7 P# H9 k! l+ Paxis 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
. v" e7 {* [) |, P5 _37];
) m) s* G9 u" Z9 ^B4=[12 13 12 12 12 12 13.5 13 13.5 14 13 12 12 11 10.25 12 11 11 15 15
8 y5 e% `/ }, _13 12 12 13 13 12];
subplot(2,2,2);
" ^2 O/ e$ j3 {/ S0 @8 }5 Splot(A4,B4,'.-');
hold on;
0 @; h$ \( P5 ^/ X: Q" Tp5=polyfit(A4,B4,2)
f5=polyval(p5,A4);
; Z8 ]5 v7 d5 }: R) S. }2 W7 Lplot(A4,B4,'.',A4,f5,'-');
4 {0 `  c* `/ hhold on;
5 L3 k1 t. q, e7 {. Stitle('青光眼病人情况');
: t  h5 {" l0 R2 Y0 }6 V" [xlabel('门诊时间距7 月12 日的时间差');ylabel('青光眼病人平均等待时间
: f' b* C+ q, F. l');
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];
* g0 U0 Q2 O  N- [0 Rsubplot(2,2,3);
plot(A5,B5,'.-');
hold on;
( }: K, z$ u& D$ q, v/ F7 J- j) wp6=polyfit(A5,B5,2)
  G5 g! X0 k! }0 E8 ff6=polyval(p6,A5);
6 |. \# Z3 _' s; }plot(A5,B5,'.',A5,f6,'-');
) v* x: b' u. P& Y0 Chold on;
: s1 @4 u/ _( E8 f" ntitle('外伤病人情况');
% B% e" q% }* Y3 X0 f6 Zxlabel('外伤病人门诊时间距离7月12日的时间差');ylabel('外伤病人平均等
, R9 q* Q; u" y( j5 t) _待时间');
% z' ]" w  L) V' _  Waxis equal;
# N2 @5 n7 f/ X8 P0 c- h二、层次分析法的一致性指标判定的matlab 实现:
function H=ff(A)
t=eig(A)
# v/ p/ r- S7 ^, C9 T. `* k( Vw=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))
- o& C6 x. E2 ?' E(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))
' t- `) s* R/ R7 c(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)
2 u# z9 i' q' W+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))]
9 A( S5 u" {1 J+ x( o7 j  k% N

linmatsas

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

迷途的羊

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

linmatsas

回复 迷途的羊 的帖子

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

紫辰

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

chenrye121

佩服两个字…………

水木年华zzu

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

wuming1429

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

hhao

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

abc321

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