24小时热门版块排行榜    

查看: 983  |  回复: 4

cg2016

新虫 (初入文坛)

[求助] 使用matlab进行参数拟合可以在搜索的时候设置待求参数的大小顺序吗? 已有2人参与

在使用matlab做参数拟合时,使用了遗传算法(ga)和levenberg-marquardt算法,在ga算法里面可以设置待求参数的大小顺序吗?比如我要求的是5个k,那么我想要让k1>k2>k3>k4>k5,这该怎么实现?
回复此楼

» 猜你喜欢

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

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

吟风舞墨

银虫 (小有名气)

最好是贴出流程图
贴出部分程序代码
2楼2016-12-08 18:20:00
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

cg2016

新虫 (初入文坛)

引用回帖:
2楼: Originally posted by 吟风舞墨 at 2016-12-08 18:20:00
最好是贴出流程图
贴出部分程序代码

lb = zeros(1,11);   % 参数下限
ub = ones(1,11)*1e+8;                % 参数上限
%遗传算法全局优化
[k, fval, exitflag, output, population, scores] = ga(@ObjFunc4ga,nvars,[],[],[],[],lb,ub,[],Options);
  fprintf('\n使用函数ga()估计得到的参数值为:\n')
     for i=1:11
         fprintf('\t k %d = %.6f\n',i,k(i));
     end
  fprintf('  The sum of the squares is: %.8f\n\n',fval)
%其中有些k的大小顺序是一定的k(2)<=k(4)<=k(6)<=k(8)
if k(2)<=k(4)
    return %这里该怎么返回重新进行全局优化
else
    continue;
end
%上述结果达到要求后进行下一步局部优化
  [k_1,resnorm,residual,exitflag,output,lambda,jacobian] =...
                 lsqnonlin(@ObjFunc4LNL,k,lb,ub,Options,ExpData);
            
%其余的就是微分方程组
3楼2016-12-08 19:19:35
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

月只蓝

主管区长 (职业作家)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与,应助指数 +1
cg2016: 金币+10, ★★★★★最佳答案 2016-12-09 09:39:53
在构造目标函数的时候,使用if语句,诸如
if k1>k2>k3>k4>k5 (这里的写法只是示例,实际代码需要用到“与”运算语句)
obj=...
esle
obj=+inf;
这种类似于罚函数法,因为+inf不可能是最优值,所以搜索的结果将满足k1>k2>k3>k4>k5的约束。
MATLAB、MS小问题、普通问题请发帖求助!时间精力有限,恕不接受无偿私信求助。
4楼2016-12-08 20:51:26
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

dingd

铁杆木虫 (职业作家)

【答案】应助回帖

感谢参与,应助指数 +1
加约束,变为有约束的拟和来求解即可。

[ 发自手机版 http://muchong.com/3g ]
5楼2016-12-08 23:13:46
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 cg2016 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 国自然面上复盘~欢迎讨论 (金币+5) +11 晴天加油 2026-08-26 12/600 2026-08-28 01:22 by li09465290@
[基金申请] 面上合作单位盖章 +5 ssyjh 2026-08-27 5/250 2026-08-27 20:50 by gdfollow
[教师之家] 导师吐槽:我怎么摊上了这么个极品研究生! +7 苏东坡二世 2026-08-23 7/350 2026-08-27 18:14 by 瞬息宇宙
[基金申请] 为什么 国际(地区)合作与交流项目 没有放榜? 10+3 majunge000 2026-08-26 11/550 2026-08-27 08:42 by 北京莱茵编辑
[基金申请] 我不理解! +15 Edward_pc 2026-08-26 23/1150 2026-08-26 20:34 by zzuzxg
[基金申请] 2026年的国家社科基金项目通讯评审的新规则与新动向、新挑战 +7 process2012 2026-08-23 10/500 2026-08-26 19:23 by hmhminy
[基金申请] 范进中举一文的中心思想 +9 炎黄贵胄 2026-08-22 10/500 2026-08-26 15:40 by semaglutide
[基金申请] 2026年8月25日国自然放榜前突然收到列入评审专家邮件,有关系吗? +25 木水思豆 2026-08-25 28/1400 2026-08-26 14:53 by draco1987
[基金申请] 能否退出参与的面上项目解除限项 +23 koalala 2026-08-24 26/1300 2026-08-26 14:29 by 宝贝虫子
[基金申请] 为什么国自然不能直接公布 +4 bjdxyxy 2026-08-26 4/200 2026-08-26 13:12 by qingmu1201
[基金申请] 国际合作可查了,中了面上 (EPI+1)(金币+50) +18 Ldrop2023 2026-08-26 18/900 2026-08-26 11:15 by cmrandy
[基金申请] 系统进不去 +4 yanglien 2026-08-26 5/250 2026-08-26 11:10 by wenfengw83
[基金申请] 国合现在查不到了吗? +10 chengyan1220 2026-08-24 20/1000 2026-08-26 08:57 by peasantsprig
[基金申请] 牛来!米来!面来! +8 beefly 2026-08-26 8/400 2026-08-26 08:37 by xuzhipiao
[基金申请] 在坚冰还盖着北海的时候,我看到了怒放的梅花。 (金币+10) +6 ziyangfang 2026-08-25 9/450 2026-08-25 20:26 by huagongfeihu
[基金申请] 明天应该可查了!? +6 chengyan1220 2026-08-23 6/300 2026-08-25 19:45 by zfd97
[基金申请] 某些机构,以效率低为荣,以效率低作为存在感 +9 yuleib84 2026-08-25 10/500 2026-08-25 17:14 by alexon
[基金申请] 没有任何消息-是不是就凉了 +9 图啦图啦 2026-08-24 10/500 2026-08-25 11:59 by 南海小哥
[基金申请] 人气不行了 +11 fansofjerry 2026-08-21 11/550 2026-08-25 11:04 by 孤独的英雄6
[教师之家] 跳槽后在研项目怎么办? +5 简单化xn 2026-08-22 10/500 2026-08-23 12:38 by 简单化xn
信息提示
请填处理意见