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

和子

clear
  U' c2 p* t* O3 H+ O5 `/ M+ |%本程序用于做双高斯拟合，拟合式子为
8 H! x  t, f/ b/ G%yi=r(1).*exp(-((x-r(2))./r(3)).^2) + r(4).*exp(-((x-r(5))./r(6)).^2);
%采用的方法是高斯-牛顿法
* Q* m3 D5 ~0 M4 N: }) a- E# l%x,y为做双高斯拟合的点，通过下面的式子产生
x=[0:.2:10];y=exp(-(x-5).^2)+.5*exp(-(x-3).^2)+.1*randn(size(x));
( ?; n" u( }6 J%假定r初始值为1~6
r=1:6;
! ~# `; H% V, O7 F" P) Sr=r';
6 a2 K9 W! k$ S! \y_size=size(x);
x_size=size(y);
/ ?( p( e; `) f. T7 Jif x_size(1)==1
  I  _/ H2 w: I/ Rx=x';
1 U' z# S( o2 B4 c+ }, vend
, _& X7 ?% W$ }. J  Jif y_size(1)==1
2 B) _! J2 s. f( ]y=y';
end
8 K9 d/ g) u1 t- O8 u4 Myi=[];   
R_square=0;
2 P$ E& q) b6 qB=zeros(length(r),1);  
SSE=10000000;
while 1
k=1;
! A. `# }3 g; j) t%控制下系数增量的步长
for j=1:7
    r1=r;
) M' B0 B) u* @( d    r=r+k.*B;
    yi=r(1).*exp(-((x-r(2))./r(3)).^2) + r(4).*exp(-((x-r(5))./r(6)).^2);
7 h0 H& G# Z% w1 T* n    RSSE=SSE;
    yy=y-yi;
& b# E. g# E' [  T    SSE=sum(yy.^2);
$ ^8 u' W* ^5 i+ u( w9 R/ U    if RSSE>=SSE
        break;
    else
4 w% r* X% K( g' l6 P5 B        k=0.5^j;
! s, d* u" J2 q, m1 l* C- a6 V7 p        r=r1;
    end
end
1 n& T/ {# b; v9 x$ RSST=sum((y-mean(y)).^2);
: H8 i! U- V' x* S1 ?R_square=1-SSE/SST;
' i3 Y+ [" M  R%R_square为确定系数与拟合优度有关
if R_square>0.9
     break;
end
%下面的算式是对原式做泰勒展开后省略二阶以上导数得到的，具体可参看高斯牛顿法过程
D_a1=exp(-(r(2) - x).^2./r(3).^2);
  ]: x6 T6 Y( w3 R$ C4 n0 f' g- GD_b1=-(r(1).*exp(-(r(2) - x).^2./r(3).^2).*(2.*r(2) - 2.*x))./r(3).^2;
; U. ~* {, i1 O4 J4 PD_c1=(2.*r(1).*exp(-(r(2) - x).^2./r(3).^2).*(r(2) - x).^2)./r(3).^3;
D_a2=exp(-(r(5) - x).^2./r(6).^2);
8 g1 ^' v) r9 B* I* \0 }D_b2=-(r(4).*exp(-(r(5) - x).^2./r(6).^2).*(2.*r(5) - 2.*x))./r(6).^2;
: L6 B/ q. |5 \: y! bD_c2=(2.*r(4).*exp(-(r(5) - x).^2./r(6).^2).*(r(5) - x).^2)./r(6).^3;
D=[D_a1 D_b1 D_c1 D_a2 D_b2 D_c2];
) A9 d+ O4 b$ V7 d& w; s3 V6 Q: s3 rB=D\yy;
# |7 p+ g3 R6 w7 ^, `end
8 y' X; _5 Q5 A
" c/ t' U- @& F: {% Q" K5 P1 L
得到的结果不好，运行慢，而且很快出现
Warning: Rank deficient, rank = 1, tol =  4.079239e-17.
# R- t+ k3 u, g6 F7 f4 J. _# @> In shiyan_shuanggaosi at 53
4 |  S" I. p. H# C- v+ Y哪位大神有好的思路指点下我