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

和子

clear
7 |" Z- r# R6 e- p%本程序用于做双高斯拟合，拟合式子为
%yi=r(1).*exp(-((x-r(2))./r(3)).^2) + r(4).*exp(-((x-r(5))./r(6)).^2);
" v: U  t$ o. s9 j" n% I" x) m8 _%采用的方法是高斯-牛顿法
; @2 L8 }4 d9 B" Y9 [. ^%x,y为做双高斯拟合的点，通过下面的式子产生
( |) Q; K! p8 F6 Q- |9 t( l- G) Wx=[0:.2:10];y=exp(-(x-5).^2)+.5*exp(-(x-3).^2)+.1*randn(size(x));
%假定r初始值为1~6
r=1:6;
0 C! {: w3 {- q$ ~r=r';
y_size=size(x);
' E) `' V* t  O& px_size=size(y);
if x_size(1)==1
x=x';
end
if y_size(1)==1
y=y';
end
4 A- R* c1 L5 P% @yi=[];   
/ h1 {, `* ?3 u" X4 hR_square=0;
. ?$ t7 q/ q1 O! l- f: ZB=zeros(length(r),1);  
SSE=10000000;
4 d' I5 d) E( P$ L6 wwhile 1
1 d& h1 T# b# `" Z/ |. Zk=1;
$ Z2 z' f+ y+ K: R%控制下系数增量的步长
5 [$ Z& C3 S( _4 v# R4 O; sfor j=1:7
  T$ y9 Z/ d$ t2 }    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);
5 ?1 }# q0 y; S: X/ \7 v  y    RSSE=SSE;
    yy=y-yi;
    SSE=sum(yy.^2);
8 q% l  I0 i2 I; b0 _    if RSSE>=SSE
        break;
0 k0 m5 @& g) W% g  [- Y: U. G. I, [    else
        k=0.5^j;
        r=r1;
    end
end
SST=sum((y-mean(y)).^2);
) |/ W; d0 G9 J9 h- `2 uR_square=1-SSE/SST;
%R_square为确定系数与拟合优度有关
4 ^. s( Y6 N6 E% _: Cif R_square>0.9
( N5 t5 F& `) e     break;
7 z0 O  G& N( @! L* jend
%下面的算式是对原式做泰勒展开后省略二阶以上导数得到的，具体可参看高斯牛顿法过程
D_a1=exp(-(r(2) - x).^2./r(3).^2);
& u4 @/ `4 k% R* A, }D_b1=-(r(1).*exp(-(r(2) - x).^2./r(3).^2).*(2.*r(2) - 2.*x))./r(3).^2;
6 D+ `* P! E& S2 O5 k) {D_c1=(2.*r(1).*exp(-(r(2) - x).^2./r(3).^2).*(r(2) - x).^2)./r(3).^3;
, q( |, L5 }1 Z, L4 B- q  u) KD_a2=exp(-(r(5) - x).^2./r(6).^2);
1 g( z. R* k: x" ?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;
D=[D_a1 D_b1 D_c1 D_a2 D_b2 D_c2];
B=D\yy;
end


得到的结果不好，运行慢，而且很快出现
. W% p3 H3 `& `. s  `Warning: Rank deficient, rank = 1, tol =  4.079239e-17.
> In shiyan_shuanggaosi at 53
哪位大神有好的思路指点下我

: k8 K3 T4 b6 _' A