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

和子

clear
%本程序用于做双高斯拟合，拟合式子为
%yi=r(1).*exp(-((x-r(2))./r(3)).^2) + r(4).*exp(-((x-r(5))./r(6)).^2);
%采用的方法是高斯-牛顿法
%x,y为做双高斯拟合的点，通过下面的式子产生
, w, K  Z3 X  c: Px=[0:.2:10];y=exp(-(x-5).^2)+.5*exp(-(x-3).^2)+.1*randn(size(x));
%假定r初始值为1~6
r=1:6;
r=r';
2 R$ J5 {8 ]: t# l0 by_size=size(x);
8 V' J$ T  L6 ~1 \x_size=size(y);
if x_size(1)==1
0 o, D5 `; L+ X- S: rx=x';
end
$ G0 _& Q8 E) p' f# R$ z6 nif y_size(1)==1
y=y';
( w# [# |( h3 K9 d" \end
yi=[];   
( A( [- t! B. H- x9 B& H$ tR_square=0;
B=zeros(length(r),1);  
SSE=10000000;
" L- P$ e6 _6 @$ |. twhile 1
) O8 g5 l0 M6 k/ C; T# ^k=1;
* ?% r. O9 K% k8 q% k4 n3 |%控制下系数增量的步长
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);
$ t3 K. [1 k! s: I5 B; }. ~* X    if RSSE>=SSE
# I' ~% L% [  k$ D$ r, {! t# G        break;
4 d, K! ^9 L$ @  G    else
        k=0.5^j;
+ I6 I2 P: S  d1 ?, Z2 f        r=r1;
! r$ @# C2 h7 @    end
end
8 ~" J( o% N( R$ p: B4 w) ASST=sum((y-mean(y)).^2);
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;
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;
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];
2 \- I; W8 B. O8 z/ x8 i/ bB=D\yy;
1 f$ k+ r. X) Z+ v" _! `$ |" i8 yend


$ f- m5 ~0 m; a. T# V$ f  A0 ?得到的结果不好，运行慢，而且很快出现
# E0 N: t/ r! x; y7 g5 JWarning: Rank deficient, rank = 1, tol =  4.079239e-17.
> In shiyan_shuanggaosi at 53
哪位大神有好的思路指点下我
/ j( H( u; R. Z- f/ b4 _4 K
2 K: _* B  h; M8 \; x7 R9 Q