24小时热门版块排行榜    

查看: 744  |  回复: 2
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

瓦兹河畔之夜

木虫 (正式写手)

[求助] 由于目标函数值本来就很小造成的无法最优化问题 已有2人参与

各位前辈好,最进求一个最优化(拟合)问题时遇到的问题。
目标函数是理想值和实际数学模型之差,但由于本身实际值就不超过1,实际值也是1以下的(因为是拟合累计分布函数),所以造成了目标函数值本身很小,不能很好地拟合数学模型。
不管给出什么初值,如果用最小二乘的方法,lsqnonlin函数就会在初值点找到最小值,原因是满足了函数误差的要求。如果用遗传算法,因为目标函数值本身很小的问题,也无法得到理想的结果。而且收敛原因也是满足了函数误差的要求。
请问这类问题怎么解决,如何得到真正的最小值?
谢谢,祝春节愉快!

贴上function
function f=modifiedweibull(x)
alpha=x(1);
beta=x(2);
lambda=x(3);
p=[0.91 1.05 1.26 1.67 1.96 2.23 2.16 2.75]*0.01;
n=[2 2 2 2 3 3 3 3];
t=[800 1000 1200 1500 1750 2000 2150 2300];

omega=n.*t./sum(n.*t);
omega=sqrt(omega);
f=sqrt(2)*omega.*(1-exp(-alpha*t.^beta.*exp(lambda.*t))-p);
f=sum(f);%这个是为了遗传算法的,如果是最小二乘就没有这句

优化参数的下界是[0 1 0]
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

baobiao007

木虫 (职业作家)

中国特色

【答案】应助回帖

对计算精度要求高的情况下,二进制编码的遗传算法满足不了要求,可以考虑实数编码的遗传算法
我同意叔本华的观点,人们投身艺术和科学领域的强烈愿望之一就是逃离痛苦、残酷和枯燥无味的现实生活,逃离自己飘忽不定的七情六欲的桎梏。--爱因斯坦
3楼2014-03-05 09:38:58
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 3 个回答

onesupeng

金虫 (职业作家)

【答案】应助回帖

感谢参与,应助指数 +1
可不可以乘以一个因子,相当于尺度化已下目标函数,等结果出来,你再换回去
长期招收博士生,参见http://fsl-unsw.com
2楼2014-02-03 12:25:13
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 帮忙看看fileCode +4 wwncly 2026-08-10 8/400 2026-08-10 22:07 by demonstreets
[基金申请] 小木虫上这么多卖论文的,真有人买论文么?感觉没必要啊 +3 Tide man 2026-08-10 3/150 2026-08-10 21:48 by GZUZM
[基金申请] filecode +8 documentary 2026-08-10 8/400 2026-08-10 19:20 by loufangrui
[基金申请] 确定了,国自然21号放榜 +6 布布和一二 2026-08-10 7/350 2026-08-10 19:15 by 2000zf36392
[基金申请] 综述论文作为代表作会不会影响评审专家的印象分? +10 yufeiwaner 2026-08-09 11/550 2026-08-10 18:47 by yufeiwaner
[基金申请] 应该是93bebmhtak前后十一个字符比较关键 +14 Lanmanbaby 2026-08-09 22/1100 2026-08-10 14:57 by wwncly
[基金申请] 8月时间戳变的,举个手。玩一下,释放压力 +11 archvillain 2026-08-04 13/650 2026-08-10 14:41 by seekfind
[基金申请] 据悉今年马上要出结果了 +7 瞬息宇宙 2026-08-10 8/400 2026-08-10 12:42 by Vivilian
[基金申请] 2026国自然放榜时间 +9 布布和一二 2026-08-08 9/450 2026-08-10 11:22 by xxxx2020
[基金申请] 这样的filecode谁见过 +11 布布和一二 2026-08-08 22/1100 2026-08-10 11:10 by wmfsnow
[基金申请] 面上项目filecode邪修 +5 西山十月 2026-08-09 7/350 2026-08-10 07:32 by 仁砚薪传
[基金申请] 关于豆爷回答的JTJC与%2F数量 +5 yang182083 2026-08-06 7/350 2026-08-08 18:29 by zhanghaozhu
[基金申请] 国基金的申报应该改成非等额制,评价高的钱多评价低的钱少,但是增加资助率 +7 a089 2026-08-07 7/350 2026-08-08 18:05 by gltch
[基金申请] 基金中了 +14 laoda193707 2026-08-06 14/700 2026-08-08 00:23 by 实验小白ha
[基金申请] 化学口download_prp&fileCode的固定段好像这几天一直没变,有变的大神么? +3 Tide man 2026-08-07 4/200 2026-08-07 22:39 by Tide man
[基金申请] 大家散了吧,后缀研究没有意义,别浪费时间了,过好目前的每一天,不要焦虑 +5 Tide man 2026-08-06 7/350 2026-08-07 13:11 by 医学老男孩
[基金申请] filecode +14 等待解的谜 2026-08-06 19/950 2026-08-07 12:20 by wlwhappy
[基金申请] filecode +8 布布和一二 2026-08-06 11/550 2026-08-06 20:41 by tangpu318
[基金申请] 求各位大神看下 100+6 hpkpkpkp 2026-08-05 33/1650 2026-08-06 14:49 by zhiyanjiang
[基金申请] 纯娱乐,不喜欢勿喷 +7 Tide man 2026-08-04 10/500 2026-08-04 15:10 by loufangrui
信息提示
请填处理意见