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

和子

clear
3 r7 s6 x5 d2 ?$ f2 e9 W) b1 s' P%本程序用于做双高斯拟合，拟合式子为
: W( t/ e8 S5 N% t) S%yi=r(1).*exp(-((x-r(2))./r(3)).^2) + r(4).*exp(-((x-r(5))./r(6)).^2);
%采用的方法是高斯-牛顿法
" x5 d; b( L. p$ e%x,y为做双高斯拟合的点，通过下面的式子产生
x=[0:.2:10];y=exp(-(x-5).^2)+.5*exp(-(x-3).^2)+.1*randn(size(x));
; |: c2 M  G6 B/ \9 e* ^4 W%假定r初始值为1~6
  R3 X4 Q. s' Mr=1:6;
( M: w' c- I4 i0 n5 ]8 j6 br=r';
y_size=size(x);
: Y/ ^5 O5 p2 b/ `/ q# xx_size=size(y);
if x_size(1)==1
x=x';
1 X. p% T6 ?' a. k0 j4 a6 aend
if y_size(1)==1
8 ]$ E# k% m9 X# V) I) t+ My=y';
end
yi=[];   
2 |2 [2 J# {; u6 B: a0 |& YR_square=0;
B=zeros(length(r),1);  
SSE=10000000;
. J; b! Y( B" q8 Ewhile 1
k=1;
%控制下系数增量的步长
! _' e1 C) ]) A2 ^for 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);
    RSSE=SSE;
    yy=y-yi;
    SSE=sum(yy.^2);
/ R; o' T3 N# e2 o) c; n    if RSSE>=SSE
        break;
    else
        k=0.5^j;
        r=r1;
% \! m- H; j* \6 h* w  X    end
end
SST=sum((y-mean(y)).^2);
) q3 j5 H. E- f9 TR_square=1-SSE/SST;
%R_square为确定系数与拟合优度有关
4 t7 B7 A) `. q( ?4 a: |  u3 nif R_square>0.9
     break;
end
5 n5 a8 W7 @& @7 _' Y# Z' q, X%下面的算式是对原式做泰勒展开后省略二阶以上导数得到的，具体可参看高斯牛顿法过程
D_a1=exp(-(r(2) - x).^2./r(3).^2);
- i! P# ^0 d9 l" i5 x  b  Q* Y- yD_b1=-(r(1).*exp(-(r(2) - x).^2./r(3).^2).*(2.*r(2) - 2.*x))./r(3).^2;
5 x; t5 [1 Q3 XD_c1=(2.*r(1).*exp(-(r(2) - x).^2./r(3).^2).*(r(2) - x).^2)./r(3).^3;
) h" y3 u0 d3 |! sD_a2=exp(-(r(5) - x).^2./r(6).^2);
D_b2=-(r(4).*exp(-(r(5) - x).^2./r(6).^2).*(2.*r(5) - 2.*x))./r(6).^2;
D_c2=(2.*r(4).*exp(-(r(5) - x).^2./r(6).^2).*(r(5) - x).^2)./r(6).^3;
) @( A& l6 a+ ~# Q% A. X5 YD=[D_a1 D_b1 D_c1 D_a2 D_b2 D_c2];
. Q" P- }3 |# H* KB=D\yy;
4 Y/ [4 h$ W$ _end
9 z! ^, {& `7 `
9 `& l4 y9 V. `' e* U4 i
6 w4 Z+ j9 x2 S1 b得到的结果不好，运行慢，而且很快出现
/ W0 b1 J5 M6 I1 F8 |% WWarning: Rank deficient, rank = 1, tol =  4.079239e-17.
# p( L4 W) U5 o, k& \5 B9 b1 o> In shiyan_shuanggaosi at 53
哪位大神有好的思路指点下我
4 s5 K" [5 b8 k5 N
# T% S$ l0 R9 K+ n+ s