求熟悉高斯牛顿法的大神帮忙看看我的程序

和子

clear
! ]1 ^: p' T2 N%本程序用于做双高斯拟合，拟合式子为
/ s' w! }% c+ r. p# J+ a, w& I8 _8 |7 E3 c+ {%yi=r(1).*exp(-((x-r(2))./r(3)).^2) + r(4).*exp(-((x-r(5))./r(6)).^2);
%采用的方法是高斯-牛顿法
3 _$ B7 K; H6 r  p0 ^4 B4 E%x,y为做双高斯拟合的点，通过下面的式子产生
. S( r4 q* {$ ?: E) fx=[0:.2:10];y=exp(-(x-5).^2)+.5*exp(-(x-3).^2)+.1*randn(size(x));
! S* W" ^3 `7 |5 R1 c- G%假定r初始值为1~6
r=1:6;
+ i- P4 ~& f& g, M' a: er=r';
y_size=size(x);
x_size=size(y);
if x_size(1)==1
x=x';
end
if y_size(1)==1
: n  `/ u) l1 Y' A: B0 by=y';
* F! O0 Z) J: o# j- eend
yi=[];   
/ c5 a+ Q; Q) ]6 |) d& [& YR_square=0;
B=zeros(length(r),1);  
SSE=10000000;
6 i2 Z& ]5 l7 J. T* Owhile 1
1 x. E" h2 V; ?3 I3 x( i& ck=1;
* D; t, B9 L3 @$ e, Q! ~: b%控制下系数增量的步长
# j5 Z7 C/ B" x" g- Ufor j=1:7
    r1=r;
    r=r+k.*B;
    yi=r(1).*exp(-((x-r(2))./r(3)).^2) + r(4).*exp(-((x-r(5))./r(6)).^2);
& o$ _( o2 H  ]' W8 p: J% f& w    RSSE=SSE;
    yy=y-yi;
    SSE=sum(yy.^2);
  P! V/ e" X* r) ~6 I! V: i" q    if RSSE>=SSE
        break;
* @: c7 v9 ?* @! q    else
        k=0.5^j;
$ a0 ~  y& W0 W; a        r=r1;
    end
end
9 a3 ^6 e2 j0 ]# F5 t1 g! l6 pSST=sum((y-mean(y)).^2);
R_square=1-SSE/SST;
%R_square为确定系数与拟合优度有关
if R_square>0.9
     break;
  `7 X4 Z2 T9 O, B% D7 vend
%下面的算式是对原式做泰勒展开后省略二阶以上导数得到的，具体可参看高斯牛顿法过程
D_a1=exp(-(r(2) - x).^2./r(3).^2);
D_b1=-(r(1).*exp(-(r(2) - x).^2./r(3).^2).*(2.*r(2) - 2.*x))./r(3).^2;
D_c1=(2.*r(1).*exp(-(r(2) - x).^2./r(3).^2).*(r(2) - x).^2)./r(3).^3;
7 n$ i! T& B/ D3 M1 x% i2 J5 sD_a2=exp(-(r(5) - x).^2./r(6).^2);
! V8 \) s+ P5 `4 ?D_b2=-(r(4).*exp(-(r(5) - x).^2./r(6).^2).*(2.*r(5) - 2.*x))./r(6).^2;
$ T5 _; A7 p* ~0 R4 ]7 \2 j9 Z# j6 e* L, D# `D_c2=(2.*r(4).*exp(-(r(5) - x).^2./r(6).^2).*(r(5) - x).^2)./r(6).^3;
: _( ^. I4 e; Y" V$ y4 gD=[D_a1 D_b1 D_c1 D_a2 D_b2 D_c2];
B=D\yy;
end

2 \+ e! o1 Z. S; p  }: _3 }
9 ]* B  Z( h( ]; ~- C" s+ |得到的结果不好，运行慢，而且很快出现
Warning: Rank deficient, rank = 1, tol =  4.079239e-17.
) H* I+ y( u; M9 f> In shiyan_shuanggaosi at 53
哪位大神有好的思路指点下我