自己编写了一个关于灰度的程序，但是有点问题。。。

无聊2012

代码如下：

1. clear
   6 _. u/ Z9 L& @8 l
2. syms a b;+ o7 G9 R( y7 @, H5 N2 n0 @
   
3. c=[a b]';
   ) M1 X/ e6 |& {. r* j: E, J
4. A=[89677,99215,109655,120333,135823,159878,182321,209407,246619,300670 ];0 O2 U! T\" F& C1 X, c
   
5. B=cumsum(A);1 R: e) T5 |' h: N; A
   
6. n=length(A);9 y7 ^; `: `8 Y% x
   
7. for i=1:(n-1)
   & U$ H1 x. P. C
8.       C(i)=(B(i)+B(i+1))/2;
   & m- K9 i4 G& x8 V' S. o
9. end$ k2 S% b1 I% k0 b: s+ r
   
10. %计算待定参数的值7 _  r4 k  d2 _: T) u
    
11. D=A;D(1)=[];
    4 b) r) p5 F$ ^! D3 k; L* w& `
12. D=D';! Z! ^- r3 I. ^9 t+ p  g1 M
    
13. E=[-C;ones(1,n-1)];
    / w1 `3 p* L  ]3 v( V
14. c=inv(E*E')*E*D;' v7 _/ u9 l' n$ D4 s
    
15. c=c';
    . D# q4 o& J* k5 y
16. a=c(1);b=c(2);
    2 t. ]8 ^2 I5 \% V) L4 {
17. %预测后续的数据6 d\" Z- `9 }/ R0 K* C, r) F
    
18. F=[];F(1)=A(1);/ U2 I0 Q: e, s; H2 ^! E; T
    
19. for i=2;(n+10)\" D3 h6 N\" B5 y4 a1 s4 ]! q4 P
    
20.    F(i)=(A(1)-b/a)/exp(a*(i-1))+b/a;
    ) b# v3 S8 @0 x, R' x' l7 s
21. end
    - H0 k3 d- Z; s# l: O
22. G=[];G(1)=A(1);3 K, I5 W8 H5 I\" g: y
    
23. for i=2:(n+10)
    % c2 v+ L3 @  [( t
24.     G(i)=F(i)-F(i-1);    %得到预测出来的数据
    : j: N' B1 h: R  }. f& Y
25. end
    0 l1 O7 z9 X2 x( a; `  x
26. t1=1999:2008;
    - s+ b$ h5 g. U7 u7 b, M
27. t2=1999:2018;
    - }7 K3 f+ C; @8 V8 e* V
28. G;
    . {. r' \/ o) E5 x% n- s& t
29. plot(t1,A,'o',t2,G);
    % I3 |) d7 e  a3 U
30. ' q; i1 G9 Y& n# v' P
    
31. 运行结果是这样的：curvefun% k& W8 f1 _+ s# A% B
    
32. 
      J( ?8 V  N. {. b
33. ans =; {7 Q& P  J8 B7 T% w
    
34. - }4 Z# ^) D( `( h0 h) N
    
35.     20% a\" u! |6 `* h
    
36. $ m\" T7 a: b! Z1 V/ t6 q
    
37. Attempted to access F(3); index out of bounds because numel(F)=2.. I9 m$ w% e7 W7 L) a( m, z, |  G
    
38. 
    $ h/ E* V1 s& v2 R7 ]
39. Error in curvefun (line 24)
      p8 w1 f( m5 c9 U
40.     G(i)=F(i)-F(i-1);    %得到预测出来的数据

复制代码