matlab中关于nlinfit函数初始值的选定

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

送您一个全额奖学金名额~ !

立即领取

感谢您参与论坛问题回答

经管之家送您两个论坛币！

+2 论坛币

本人是matlab的菜鸟啊，问大神们一个问题啊，做非线性拟合时选取的初值不同，系数拟合的结果往往也不同，到底什么样的结果才是最好的呢？本人想在真实值附近找一个区间，通过一个较小的步长来穷举，最后选残差最小的那组系数。。。无奈程序老是报错，求大神们看看啊！！！

主程序：
A=[33.35                      1854 34.78                      1857 36.22                      1948 37.66                      201739.09                      209240.53                      2197 41.76                      2608 42.81                      2937 43.9                      311944.94                      3654];process=A(:,1);  %自变量price=A(:,2);  %因变量aa=0.0001:0.0001:0.5;  %穷举第一个系数bb=0.001:0.001:0.5;  %穷举第二个系数cc=1701:1900;  %穷举第三个系数rr=[];  %设定一个三维空矩阵放残差low=[1e-008,1e-008,1500];  %系数beta的下线up=[2,2,2500];  %系数beta的上限for i=1:5000;    for m=1:500;        for k=1:200;beta0=[aa(1,i),bb(1,m),cc(1,k)]';  %初值[beta,r]=nlinfit(process,price,@myfun,beta0);  %调用myfun函数if beta<low | beta>up;  %系数返回值不符合要求的话，残差就定义为0    flag=0;else    flag=r'*r;  %计算残差end;rr(i,m,k)=flag;  %得出三维残差阵        end;    end;end;rmin=nanmin(nanmin(nanmin(rr)))   %找出最小的残差[ii,mm,kk]=find(rr==rmin)  %找出最小残差的位置
函数：function yy=myfun(beta,x)a=beta(1);b=beta(2);c=beta(3);yy=a*exp(b*x)+c;
此程序一直报错：??? Error using ==> nlinfit>checkFunVals at 356MODELFUN has returned Inf or NaN values.
Error in ==> nlinfit>LMfit at 303    if funValCheck && ~isfinite(sse),    checkFunVals(r); end
Error in ==> nlinfit at 166    [beta,J,lsiter,cause] = LMfit(X,y,    model,beta,options,verbose,maxiter);
Error in ==> fangjia at 24[beta,r]=nlinfit(process,price,@myfun,beta0);
我看不懂啊，怎样才能在函数不收敛的时候跳过去接着计算下一个啊，大神们快帮我看看！！！

扫码加我 拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

分享0 收藏0 回帖

关键词：nlinfit MATLAB matla atlab Nlin matlab nlinfit 初值设置 穷举

第一次发帖，编辑的时候明明是分行的，怎么发出来就全堆一起了？

再发一遍：
A=[33.35    1854
34.78      1857
36.22      1948
37.66      2017
39.09       2092
40.53       2197
41.76      2608
42.81      2937
43.9       3119
44.94        3654
];
process=A(:,1);%自变量
price=A(:,2); %因变量
aa=0.0001:0.0001:0.5;  %穷举第一个系数
bb=0.001:0.001:0.5; %穷举第二个系数
cc=1701:1900; %穷举第三个系数
rr=[]; %设定一个三维空矩阵放残差
low=[1e-008,1e-008,1500];  %系数beta的下线
up=[2,2,2500]; %系数beta的上限
for i=1:5000;
    for m=1:500;
        for k=1:200;
beta0=[aa(1,i),bb(1,m),cc(1,k)]';  %初值
[beta,r]=nlinfit(process,price,@myfun,beta0); %调用myfun函数
if beta<low | beta>up; %系数返回值不符合要求的话，残差就定义为0
    flag=0;
else
    flag=r'*r;  %计算残差的平方和
end;
rr(i,m,k)=flag;%得出三维残差阵
        end;
    end;
end;
rmin=nanmin(nanmin(nanmin(rr))) %找出最小的残差
[ii,mm,kk]=find(rr==rmin) %找出最小残差的位置

函数：

function yy=myfun(beta,x)
a=beta(1);
b=beta(2);
c=beta(3);
yy=a*exp(b*x)+c;

报错：
??? Error using ==> nlinfit>checkFunVals at 356
MODELFUN has returned Inf or NaN values.

Error in ==> nlinfit>LMfit at 303
    if funValCheck && ~isfinite(sse),
    checkFunVals(r); end

Error in ==> nlinfit at 166
    [beta,J,lsiter,cause] = LMfit(X,y,
    model,beta,options,verbose,maxiter);

Error in ==> fangjia at 24
[beta,r]=nlinfit(process,price,@myfun,beta0);

我看不懂啊，怎样才能在函数不收敛的时候跳过去接着计算下一个啊，大神们快帮我看看！！！

楼主  最后解决了吗？能否把代码贴出来呀