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

无聊2012

代码如下：

1. clear & c# P$ l  E  [- M9 ^+ Z
   
2. syms a b;
   % t* f' m: ?/ Q2 l. [! Y( D
3. c=[a b]';
   2 A3 ~: v0 @8 y6 m1 n+ v
4. A=[89677,99215,109655,120333,135823,159878,182321,209407,246619,300670 ];
   / a8 E2 [% F7 r# n$ ~
5. B=cumsum(A);! Y0 A8 p3 w7 [4 f; O, ]
   
6. n=length(A);
   ; T3 H+ d; c/ X2 x8 S0 N' w
7. for i=1:(n-1)
   2 ?0 y5 K5 [$ g3 \/ ^1 X8 }
8.       C(i)=(B(i)+B(i+1))/2;
   ( J, q( Q7 D3 N2 t3 Z, l* Q% e  z
9. end& g$ r; p  K# a4 m
   
10. %计算待定参数的值  G* Z& X& V+ G/ A
    
11. D=A;D(1)=[];  j# V8 n6 ]: N. l5 y6 S
    
12. D=D';
    8 B  _- H& L. h1 y
13. E=[-C;ones(1,n-1)];. }: z) ^: m- ^, O3 c
    
14. c=inv(E*E')*E*D;
    9 N! G+ t# @4 ?6 }2 }. \7 R
15. c=c';
    7 W$ ]/ u\" l( W0 I8 J
16. a=c(1);b=c(2);
    % y1 I' ^/ g4 R9 b2 K
17. %预测后续的数据
    ( D' [6 [: v& P0 Y\" J4 _
18. F=[];F(1)=A(1);* Y\" S8 u0 u5 }( M% u
    
19. for i=2;(n+10)
    7 o6 c' ?\" m' T: e( s
20.    F(i)=(A(1)-b/a)/exp(a*(i-1))+b/a;
    & e, w; G8 F) H2 x5 @
21. end; V2 _$ A) l( ?% g; E
    
22. G=[];G(1)=A(1);
    1 W$ r! o# g+ A6 E* h
23. for i=2:(n+10)
    - I, C1 ^\" U\" v# O0 V- a9 s9 m& i0 ~
24.     G(i)=F(i)-F(i-1);    %得到预测出来的数据
    3 N# a\" k1 x! e; w/ F/ e+ L; P
25. end
    ( T- v: P' U1 s6 G0 _. p
26. t1=1999:2008;
    0 }* `5 r9 e/ f* Z/ g/ F; B
27. t2=1999:2018;
    4 n) @1 `1 ?8 u
28. G;( X: h0 h& n: R
    
29. plot(t1,A,'o',t2,G);% k: ]/ ]! h4 X& v2 p
    
30. 2 D; C\" I. O; j! z9 `2 H
    
31. 运行结果是这样的：curvefun1 }# n8 k- ]- k; H/ G
    
32. 
    * M% f8 \4 H\" m! T* a# }. z$ V
33. ans =+ J. f2 N: G: ^9 i
    
34. 
    5 n, |$ `3 f( s
35.     20$ j4 P1 y1 ]8 a( u+ h3 X
    
36. 
    $ R+ s. l\" w\" c; C* S; a9 p& W3 r
37. Attempted to access F(3); index out of bounds because numel(F)=2.& R/ t4 f# v- K  {3 _2 m6 U
    
38. ; J/ h0 N6 T: Z2 x\" E2 u' _
    
39. Error in curvefun (line 24)
    ; B5 I5 C0 y* R
40.     G(i)=F(i)-F(i-1);    %得到预测出来的数据

复制代码