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

和子

clear
%本程序用于做双高斯拟合，拟合式子为
* S4 r8 ~0 t( k%yi=r(1).*exp(-((x-r(2))./r(3)).^2) + r(4).*exp(-((x-r(5))./r(6)).^2);
%采用的方法是高斯-牛顿法
%x,y为做双高斯拟合的点，通过下面的式子产生
x=[0:.2:10];y=exp(-(x-5).^2)+.5*exp(-(x-3).^2)+.1*randn(size(x));
%假定r初始值为1~6
r=1:6;
' r* D7 i" ^- U, W5 `1 y/ h; Mr=r';
y_size=size(x);
" I" ~( [# u+ j4 jx_size=size(y);
if x_size(1)==1
x=x';
end
4 X5 H; X1 z2 \6 k7 Q; V: |if y_size(1)==1
y=y';
end
6 W4 `/ c; E0 H# cyi=[];   
R_square=0;
! ~4 R: B8 p, YB=zeros(length(r),1);  
SSE=10000000;
while 1
/ C* B: ^, ^1 V6 a; E/ U' yk=1;
%控制下系数增量的步长
for j=1:7
    r1=r;
% p! b3 i3 g; x) B' ~# y* n    r=r+k.*B;
    yi=r(1).*exp(-((x-r(2))./r(3)).^2) + r(4).*exp(-((x-r(5))./r(6)).^2);
, H4 Z+ B/ x) J) L    RSSE=SSE;
# b1 Z  m( ?# y* f6 R" J$ M    yy=y-yi;
    SSE=sum(yy.^2);
0 {3 i' S" C8 ]8 t) U    if RSSE>=SSE
        break;
7 v4 V" L7 h' o0 H0 z* w    else
5 n" N% F+ |% l: |- O        k=0.5^j;
        r=r1;
    end
) }6 b8 f# b+ ]4 Fend
SST=sum((y-mean(y)).^2);
R_square=1-SSE/SST;
%R_square为确定系数与拟合优度有关
1 }% B7 O' l8 K% Kif R_square>0.9
$ g# K- {, g& ~8 ~     break;
. T0 j; k3 O  U; H1 wend
5 g/ L" l( c/ w$ g, k%下面的算式是对原式做泰勒展开后省略二阶以上导数得到的，具体可参看高斯牛顿法过程
  `1 z" I& }+ y% u; ?5 g  eD_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;
2 c; u7 _- J  r! G+ R1 X( ^D_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);
D_b2=-(r(4).*exp(-(r(5) - x).^2./r(6).^2).*(2.*r(5) - 2.*x))./r(6).^2;
( k, i, \# L, `8 [D_c2=(2.*r(4).*exp(-(r(5) - x).^2./r(6).^2).*(r(5) - x).^2)./r(6).^3;
8 }. T8 g! F  @1 v/ dD=[D_a1 D_b1 D_c1 D_a2 D_b2 D_c2];
B=D\yy;
end

/ Q  g& o- k& K6 t4 J
- b) V' Q' p4 e5 q1 b得到的结果不好，运行慢，而且很快出现
Warning: Rank deficient, rank = 1, tol =  4.079239e-17.
: E( ]7 C8 G# G7 t8 K> In shiyan_shuanggaosi at 53
哪位大神有好的思路指点下我