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

和子

clear
%本程序用于做双高斯拟合，拟合式子为
%yi=r(1).*exp(-((x-r(2))./r(3)).^2) + r(4).*exp(-((x-r(5))./r(6)).^2);
, X* {0 Y( `8 e4 C' a! C% ?%采用的方法是高斯-牛顿法
%x,y为做双高斯拟合的点，通过下面的式子产生
0 F0 V7 m2 A  F" A. Wx=[0:.2:10];y=exp(-(x-5).^2)+.5*exp(-(x-3).^2)+.1*randn(size(x));
%假定r初始值为1~6
) }( n' h* D& T* _r=1:6;
! T: N' A8 b8 G' Or=r';
y_size=size(x);
$ t. p" s; F! ~x_size=size(y);
if x_size(1)==1
x=x';
end
if y_size(1)==1
y=y';
end
8 L$ {6 l; H. a* w' o3 J& z2 |yi=[];   
R_square=0;
B=zeros(length(r),1);  
SSE=10000000;
2 t$ S8 g+ c0 Z- ^while 1
k=1;
%控制下系数增量的步长
0 [! T0 x5 m- N5 C3 d# Sfor 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);
! F, w/ P" H: N: T$ M    RSSE=SSE;
8 \- k# k1 ~; Q# a5 g3 j% o4 a    yy=y-yi;
9 S8 S6 W& X% l( H    SSE=sum(yy.^2);
    if RSSE>=SSE
        break;
    else
0 \( s* n1 P  I* v3 |        k=0.5^j;
# h% p, q3 \! }3 F" k2 S: _        r=r1;
    end
3 X/ P% ~+ h0 g% S% w& N7 Gend
SST=sum((y-mean(y)).^2);
8 ^+ \0 ?( s5 L1 nR_square=1-SSE/SST;
%R_square为确定系数与拟合优度有关
if R_square>0.9
  {6 q; S* n# f1 v     break;
end
%下面的算式是对原式做泰勒展开后省略二阶以上导数得到的，具体可参看高斯牛顿法过程
9 P1 h! g& A: d. gD_a1=exp(-(r(2) - x).^2./r(3).^2);
# Q' F; X9 X3 p, eD_b1=-(r(1).*exp(-(r(2) - x).^2./r(3).^2).*(2.*r(2) - 2.*x))./r(3).^2;
. z1 F* o" d1 b6 bD_c1=(2.*r(1).*exp(-(r(2) - x).^2./r(3).^2).*(r(2) - x).^2)./r(3).^3;
+ c  R% @( a4 _D_a2=exp(-(r(5) - x).^2./r(6).^2);
) h+ {( b$ v% a" E0 G. dD_b2=-(r(4).*exp(-(r(5) - x).^2./r(6).^2).*(2.*r(5) - 2.*x))./r(6).^2;
4 g0 b) _. A# S5 AD_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];
' {/ l8 R2 `% o% z& KB=D\yy;
end
' K- P4 `" b2 ~0 U! S- \6 A
& z2 m9 b2 {# R7 j) X* u# s' {
/ h! @8 b9 e5 H4 w; L得到的结果不好，运行慢，而且很快出现
Warning: Rank deficient, rank = 1, tol =  4.079239e-17.
> In shiyan_shuanggaosi at 53
哪位大神有好的思路指点下我