24小时热门版块排行榜    

查看: 2543  |  回复: 18
【悬赏金币】回答本帖问题,作者小小莓邮将赠送您 5 个金币
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

小小莓邮

铁虫 (小有名气)

[求助] matlab拟合参数显示初始点是一个局部最小值已有5人参与

小弟最近刚学习用matlab拟合参数,但编写完运行显示Initial point is a local minimum.得到的k值为k0,求给位大神给指点一下怎么解决这种问题!急求,新人啥也不懂啊!先在此谢过啦!
function zixie
clear all;
clc;
data=[0 8.32 0;
    1 7.37 0.211;
    2 6.74 0.738;
    3 5.92 1.106;
    4 5.68 1.296;
    5 5.01 1.583;
    6 6.09 1.364;
    8 5.35 1.84;
    10 4.99 2];
x0=[8.32 0];
k0=[10 10 10 10 10 ];
lb=[0 0 0 0 0];
ub=[+inf +inf +inf +inf +inf];
yexp=data(:,2:3);
[k,resnorm,residual,exitflag,output,lambda,jacobin]=lsqnonlin(@objFunc1,k0,lb,ub,[],x0,yexp)
fprintf('\tk1=%.4\n',k(1)),
fprintf('\tk2=%.4\n',k(2)),
fprintf('\tk3=%.4\n',k(3)),
fprintf('\tk4=%.4\n',k(4)),
fprintf('\tk5=%.4\n',k(5))
function f=objFunc1(k,x0,yexp)
tspan=[0 1 2 3 4 5 6 8 10];
[t x]=ode45(@funceqs,tspan,x0,[],k);
y(:,1)=x(:,1);
y(:,2)=x(:,2);
f1=y(:,1)-yexp(:,1);
f2=y(:,2)-yexp(:,2);
f=[f1;f2];
function dxdt=funceqs(t,x,k)
dx1dt=-k(1)*k(2)*x(1)*sqrt(k(3)*x(2))/(k(1)+k(2)*x(1))-k(4)*k(5)*x(2);
dx2dt=k(1)*k(2)*x(1)*sqrt(k(3)*x(2))/(k(1)+k(2)*x(1))+k(5)*x(2);
dxdt=[dx1dt;dx2dt];

Initial point is a local minimum.


Optimization completed because the size of the gradient at the initial point
is less than the default value of the function tolerance.


<stopping criteria details>
k =

    10   10   10    10    10
Optimization completed: The final point is the initial point.
The first-order optimality measure, 0.000000e+00, is less than
options.TolFun = 1.000000e-06.

Optimization Metric                                   Options
relative first-order optimality =   0.00e+00         TolFun =   1e-06 (default)

>>
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐,对您同样有帮助:

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

月只蓝

主管区长 (职业作家)

【答案】应助回帖


感谢参与,应助指数 +1
fegg7502: 金币+1, 3ks 2014-02-27 08:43:13
k0=[10 10 10 10 10 ] 换成其他数值试试。
MATLAB、MS小问题、普通问题请发帖求助!时间精力有限,恕不接受无偿私信求助。
2楼2014-02-26 22:02:02
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

月只蓝

主管区长 (职业作家)

【答案】应助回帖

引用回帖:
7楼: Originally posted by 小小莓邮 at 2014-02-27 09:49:39
也就是说这造成这种现象是初值的原因是吧!有没有什么方法可以先初步估计一下初值,然后再以这个值为初值进一步估计参数?...

先调用遗传算法预估一个结果,以此结果作为初值,再调用lsqnonlin函数计算。
MATLAB、MS小问题、普通问题请发帖求助!时间精力有限,恕不接受无偿私信求助。
11楼2014-02-27 10:10:24
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 小小莓邮 的主题更新
不应助 确定回帖应助 (注意:应助才可能被奖励,但不允许灌水,必须填写15个字符以上)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[论文投稿] 关于论文第一通讯作者的问题? +12 winsaint 2024-06-26 12/600 2024-06-29 14:15 by 杨18354098226
[基金申请] 感觉地区基金太难了,申请了多年不中 +24 123为wqw 2024-06-27 31/1550 2024-06-29 14:14 by Siuhong Yiu
[考博] 30岁 想读非全博,求建议 +7 八千里路云雨泥 2024-06-27 19/950 2024-06-29 11:24 by 八千里路云雨泥
[基金申请] 我劝基金委要善良+一个大胆的设想 +13 我是大王 2024-06-28 14/700 2024-06-29 09:32 by Higraduate
[基金申请] 2024上会名单已出 +5 dream200606 2024-06-28 7/350 2024-06-29 00:11 by 忆念7
[找工作] offer选择 +13 羡渔樵 2024-06-26 13/650 2024-06-28 23:59 by 高敖曹
[药学] 大佬们来讨论下,为什么分子里加入了有毒性的官能团分子反而没毒了 +4 HF111001 2024-06-24 5/250 2024-06-28 22:45 by 渣男2019
[教师之家] 高校辞职,要求赔偿,这到底合不合理 +25 传动_海神 2024-06-23 37/1850 2024-06-28 12:39 by caiqingchong
[基金申请] 评审专家会不会很在意申请人的单位啊 +11 lancet0903 2024-06-24 14/700 2024-06-28 09:45 by xingmoon0869
[硕博家园] 联培博士文章第一单位署名问题交流~ +12 橙成成c 2024-06-23 25/1250 2024-06-28 07:14 by 逆行的路人
[论文投稿] 为何全世界都用JCR分区,咱们要弄个中科院分区?意义何在? +6 意得辑_editage 2024-06-25 7/350 2024-06-27 20:59 by keyaner23
[基金申请] 国基在研影响申请结果吗 +14 WOWO159357 2024-06-26 23/1150 2024-06-27 20:27 by wspglt
[基金申请] 基金评审时,评审专家还回去看申请人代表作之外的文章吗? +12 孤单12站 2024-06-26 18/900 2024-06-27 18:30 by 孤单12站
[论文投稿] Nano energy 投稿一个月with editor 到今天状态没有变过 5+3 甜甜甜。 2024-06-26 3/150 2024-06-27 18:00 by nono2009
[有机交流] 大佬们,打的核磁氢谱与chemdraw预测的有偏差 +7 来了个 2024-06-26 8/400 2024-06-27 14:38 by hptianyan
[基金申请] 博后面上今天下午会公布吗?有无消息? +26 hajkdfdf 2024-06-24 45/2250 2024-06-27 09:19 by kyukitu
[硕博家园] 【45岁以上博士】柔性人才引进项目,有补贴 +9 Dreamsummit 2024-06-24 18/900 2024-06-26 12:28 by hujm159
[教师之家] 神奇的中医 +8 水冰月月野兔 2024-06-24 10/500 2024-06-26 12:10 by lyfbangong
[有机交流] 对苯乙烯磺酰氯的合成机理 10+3 该死的科研 2024-06-24 5/250 2024-06-25 17:30 by 王学士
[公派出国] 博士csc联培会看重第一学历学校层次吗 +4 也就这样 2024-06-23 4/200 2024-06-24 08:18 by 晓目崇
信息提示
请填处理意见