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

和子

clear
( j/ \- o# x2 r, P6 D$ X3 H7 H%本程序用于做双高斯拟合，拟合式子为
1 p2 A  |1 B+ n  }6 s%yi=r(1).*exp(-((x-r(2))./r(3)).^2) + r(4).*exp(-((x-r(5))./r(6)).^2);
7 U( e- Z+ l3 ^, [1 K5 X; x8 U! @2 Z%采用的方法是高斯-牛顿法
%x,y为做双高斯拟合的点，通过下面的式子产生
2 s; Z0 m8 ^6 ?# r% Rx=[0:.2:10];y=exp(-(x-5).^2)+.5*exp(-(x-3).^2)+.1*randn(size(x));
. Y2 U( v+ Q; N%假定r初始值为1~6
r=1:6;
r=r';
y_size=size(x);
x_size=size(y);
if x_size(1)==1
x=x';
end
if y_size(1)==1
- I3 J# l  Z7 t. p) Cy=y';
( v  w: w6 y. Bend
' W! X" x9 P, R. e1 myi=[];   
* ?. c: }7 u1 {8 C: w- PR_square=0;
B=zeros(length(r),1);  
1 n- F$ D, w5 q* f. OSSE=10000000;
while 1
k=1;
%控制下系数增量的步长
for j=1:7
* c0 m2 u: U! S9 `. ?  C/ ]    r1=r;
5 d! o# N5 j2 `8 s0 K: s: H, ^    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);
    if RSSE>=SSE
0 C2 q0 G4 W. x9 D$ i        break;
2 n' K: `( ?# s$ a5 m$ n    else
        k=0.5^j;
) Y; L& G, k" z        r=r1;
5 d6 t& b! Y  ]% p    end
8 K2 ]; x) m" w+ F' g" @! gend
. z# o/ H1 y0 o8 k1 K7 E' i: BSST=sum((y-mean(y)).^2);
R_square=1-SSE/SST;
%R_square为确定系数与拟合优度有关
if R_square>0.9
     break;
end
" R" {2 v0 @( h( F: I2 I%下面的算式是对原式做泰勒展开后省略二阶以上导数得到的，具体可参看高斯牛顿法过程
& z# r. K; T$ J4 w- e/ C& yD_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;
# p/ J% n( x3 q. |1 ^D_a2=exp(-(r(5) - x).^2./r(6).^2);
# z( }! Z3 F2 H/ eD_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;
D=[D_a1 D_b1 D_c1 D_a2 D_b2 D_c2];
3 U$ }1 h; i/ p( vB=D\yy;
end


; |6 h; M. w# a. ]1 S/ m% }. v$ k2 V得到的结果不好，运行慢，而且很快出现
+ Z; o/ J7 u+ H, M2 E) cWarning: Rank deficient, rank = 1, tol =  4.079239e-17.
1 E6 J* I3 \6 P( a3 @> In shiyan_shuanggaosi at 53
哪位大神有好的思路指点下我
2 j! W, V" R$ E0 x- k