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

和子

clear
%本程序用于做双高斯拟合，拟合式子为
2 Z7 q% _6 w6 O  R%yi=r(1).*exp(-((x-r(2))./r(3)).^2) + r(4).*exp(-((x-r(5))./r(6)).^2);
%采用的方法是高斯-牛顿法
%x,y为做双高斯拟合的点，通过下面的式子产生
! j/ [) r  D) @; Jx=[0:.2:10];y=exp(-(x-5).^2)+.5*exp(-(x-3).^2)+.1*randn(size(x));
%假定r初始值为1~6
% r7 g2 d; g! K: [8 c# q2 x# ^r=1:6;
r=r';
' {6 q( A, Y# R# e$ D3 e, [5 Gy_size=size(x);
8 s& T7 F6 K6 \" N& fx_size=size(y);
if x_size(1)==1
x=x';
end
" t% }2 X: s% L3 d! G+ N- Xif y_size(1)==1
y=y';
end
yi=[];   
# r8 Y9 r$ _4 l& M( m% a1 z9 k% f" GR_square=0;
B=zeros(length(r),1);  
SSE=10000000;
( j* \) R: s% _2 T9 D& i8 vwhile 1
1 G# r! a7 d6 Z" n' y, f3 Qk=1;
' J8 K# d" i/ j! a0 ~! ^! n% n%控制下系数增量的步长
for j=1:7
    r1=r;
1 o- L# i( y; W! n2 f5 G    r=r+k.*B;
' _, n- x- c# G# t* b3 `    yi=r(1).*exp(-((x-r(2))./r(3)).^2) + r(4).*exp(-((x-r(5))./r(6)).^2);
4 g" k- k5 w  K( J    RSSE=SSE;
    yy=y-yi;
/ ^+ ~2 [: R" Q2 Y. E/ M. o( d4 [    SSE=sum(yy.^2);
    if RSSE>=SSE
        break;
    else
        k=0.5^j;
        r=r1;
1 ~: M( d8 _; m! X6 o3 G+ n1 k    end
) O  H$ \& V7 ~) h  v" Tend
- Z% |5 J4 Q  i$ k6 `7 ?3 C2 c9 NSST=sum((y-mean(y)).^2);
7 F# ]3 u; Q6 b2 Z7 R- Y& \R_square=1-SSE/SST;
%R_square为确定系数与拟合优度有关
if R_square>0.9
     break;
end
%下面的算式是对原式做泰勒展开后省略二阶以上导数得到的，具体可参看高斯牛顿法过程
D_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;
, ?2 O1 }! o% z( }D_a2=exp(-(r(5) - x).^2./r(6).^2);
! C' y" h% ^  GD_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# P, x" C" A+ k( p) YD=[D_a1 D_b1 D_c1 D_a2 D_b2 D_c2];
B=D\yy;
6 S2 _* J& t" X7 p: z, Y$ L) Oend
( }, g# A# |' U0 A; x

得到的结果不好，运行慢，而且很快出现
( J5 l( y' y, o: ^$ DWarning: Rank deficient, rank = 1, tol =  4.079239e-17.
> In shiyan_shuanggaosi at 53
哪位大神有好的思路指点下我
# H; |0 f) t, [/ S( ^