24小时热门版块排行榜    

查看: 2542  |  回复: 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的回帖

pdl9527

专家顾问 (小有名气)

【答案】应助回帖


感谢参与,应助指数 +1
fegg7502: 金币+1, 鼓励交流 2014-02-27 08:43:41
先用原方程线性化求得K作为初值,然后再用最小二乘拟合~~
5楼2014-02-27 00:00:32
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 小小莓邮 的主题更新
不应助 确定回帖应助 (注意:应助才可能被奖励,但不允许灌水,必须填写15个字符以上)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 青拔答辩通知出来了 +5 5011fudan 2024-06-28 5/250 2024-06-29 15:59 by Agela1129
[论文投稿] 关于论文第一通讯作者的问题? +12 winsaint 2024-06-26 12/600 2024-06-29 14:15 by 杨18354098226
[论文投稿] 求大佬们推荐一下与生物材料、骨修复主题相关见刊快的期刊 +3 莫咯老头 2024-06-24 4/200 2024-06-29 14:11 by 杨18354098226
[基金申请] 科研之友最近访客那出现来自北京的未知用户,什么情况? +11 lihangshi 2024-06-28 15/750 2024-06-29 13:40 by THK1912
[基金申请] 能看出是否上会了吗 +11 articlefan 2024-06-23 17/850 2024-06-29 13:28 by articlefan
[访问学者] 蹲! 希望今天能有大红奖状 +15 326lhpqk 2024-06-28 19/950 2024-06-29 11:47 by wenxixiaoer
[访问学者] 蹲CSC结果,会延迟到七月初吗? +9 326lhpqk 2024-06-27 11/550 2024-06-29 11:36 by 加菲的树袋熊
[基金申请] 博后BC就没人问么?存在感好低啊! +5 qade544 2024-06-28 5/250 2024-06-29 00:20 by 666中
[基金申请] 出结果了 +25 wuyue526 2024-06-28 31/1550 2024-06-29 00:16 by wjl_135
[药学] 大佬们来讨论下,为什么分子里加入了有毒性的官能团分子反而没毒了 +4 HF111001 2024-06-24 5/250 2024-06-28 22:45 by 渣男2019
[基金申请] 恳请大家帮我看看如何提高 +10 drwuwenhan 2024-06-28 13/650 2024-06-28 22:14 by 好人与坏人
[考博] 25申请博士 +3 lxcqwer 2024-06-27 3/150 2024-06-28 12:28 by highxixi
[基金申请] 博后特别资助与面上时间 +23 alex5387 2024-06-27 31/1550 2024-06-28 09:56 by 心上之田
[论文投稿] 为何全世界都用JCR分区,咱们要弄个中科院分区?意义何在? +6 意得辑_editage 2024-06-25 7/350 2024-06-27 20:59 by keyaner23
[有机交流] 大佬们,打的核磁氢谱与chemdraw预测的有偏差 +7 来了个 2024-06-26 8/400 2024-06-27 14:38 by hptianyan
[有机交流] 有机物的酸酐如何准确表征 10+3 方酱 2024-06-24 7/350 2024-06-26 17:53 by 宁静远行
[基金申请] 基金申请书名称有变化 +5 xuel2011 2024-06-25 7/350 2024-06-26 00:07 by 老虎当猫养
[考博] 申博好难 +6 自强不息a?a 2024-06-24 6/300 2024-06-25 23:02 by 考研吒儿
[硕博家园] 数据不好 +5 Hetai 2024-06-23 7/350 2024-06-25 12:37 by 1591099
[金属] 寻找钛合金热压缩代做 +4 liuyang358 2024-06-23 4/200 2024-06-24 13:50 by 搬砖狗不放弃
信息提示
请填处理意见