应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

Znn3bq.jpeg
小木虫论坛-学术科研互动平台 » 计算模拟区 » 计算模拟 » 心得、技巧、经验分享 » 最优化方法求教！！！
北京石油化工学院2026年研究生招生接收调剂公告
161/2返回列表12下一页
查看: 2196  |  回复: 15
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

vvaa

新虫 (初入文坛)

[求助] 最优化方法求教！！！ 已有1人参与

我的课题要做一个动力学模型，目前打算用集总的方法，建立一个微分方程组，比如
dc1/dt=-k1c1+k2c2
dc2/dt=k1c1-(k2+k3)c2+k4c3
dc3/dt=k3c2-(k4+k5)c3+k6c4
dc4/dt=k5c3-(k6+k7)c4
dc5/dt=k7c4
因为方程之间耦联，不能得到解析解。用bfgs法优化求解参数时，导函数不知道怎么求解。本人小白，没有最优化基础，希望大神能指导一下，怎么求dc/k，或者有没有比较好的不用导函数的优化方法推荐？

发自小木虫Android客户端
回复此楼

» 收录本帖的淘帖专辑推荐

优化

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
1楼 2017-06-27 18:15:43
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
回帖支持 ( 显示支持度最高的前 50 名 )

huab1984666

新虫 (著名写手)

【答案】应助回帖

★ ★ ★ ★ ★
感谢参与，应助指数 +1
vvaa: 金币+5, ★★★很有帮助 2017-06-28 09:27:15
function dy=myfunx(t,y,Parameters)

%%
k1=Parameters(1);
k2=Parameters(2);
k3=Parameters(3);
k4=Parameters(4);
k5=Parameters(5);
k6=Parameters(6);
k7=Parameters(7);
%%
dy=zeros(5,1);
dy(1)=-k1*y(1)+k2*y(2);
dy(2)=k1*y(1)-(k2+k3)*y(2)+k4*y(2);
dy(3)=k3*y(2)-(k4+k5)*y(3)+k6*y(4);
dy(4)=k5*y(5)-(k6+k7)*y(4);
dy(5)=k7*y(4);
%%%%


function f=Kinetic(beta)
global c;
global t;
Parameters=beta;
t0=t(1);
tn=t(end);
%%
options=odeset('RelTol',1e-4;'AbsTol',[1e-4 1e-4 1e-4 1e-4 1e-5]);
[tx,cx]=ode15s(@myfunx,[t0 tn],[xx 0 0 0 0],options,Parameters);%%xx是t=0的时候，c的值。
c_x=interp1(tx,cx,t,'spline');
c1=c_x(:,1);c2=c_x(:,2);c3=c_x(:,3);c4=c_x(:,4);c5=c_x(:,5);
ce=c1+c2+c3+c4+c5;
N=length(t);
for i=1:1:N
    rx(i)=(ce(i)-c(i)).^2;
end
f=sum(rx);


global c; c=[- - - - -];%实验值
global t; t=[- - - - -];实验值
lb=[0 0 0 0 0 0 0];
ub=[2 2 2 2 2 2 2];
options =psoptimset('Display','Iter','MaxIter',2000,'TolFun',1e-8,'CompleteSearch','on');
beta0=[0.7 0.01 0.01 0.0176 0.02 0.03 0.05];%% k1-7初始值
[beta,fval,exitflag,output] = patternsearch(@Kinetic,beta0,[],[],[],[],lb,ub,options );
%%
Parameters=beta;
t0=t(1);
tn=t(end);
%%
Options=odeset('RelTol',1e-4,'AbsTol',[1e-4 1e-4 1e-4 1e-4 1e-5]);
[tx,cx]=ode15s(@myfunx,[t0 tn],[1 0 0 0 0],Options,Parameters);
c_x=interp1(tx,cx,t,'spline');
c1=c_x(:,1);c2=c_x(:,2);c3=c_x(:,3);c4=c_x(:,4);c5=c_x(:,5);
ce=c1+c2+c3+c4+c5;

plot(t,c,'o');hold on;
plot(t,ce,'r-');
赞一下(1人) 回复此楼
春风又绿江南岸，明月何时照我还。
2楼2017-06-27 20:53:54
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通回帖

huab1984666

新虫 (著名写手)
引用回帖:
2楼: Originally posted by huab1984666 at 2017-06-27 20:53:54
function dy=myfunx(t,y,Parameters)

%%
k1=Parameters(1);
k2=Parameters(2);
k3=Parameters(3);
k4=Parameters(4);
k5=Parameters(5);
k6=Parameters(6);
k7=Parameters(7);
%%
dy=zeros(5,1);
dy(1 ...

global c; c=[- - - - -];%实验值
global t; t=[- - - - -];实验值

c,t都是列向量
赞一下 回复此楼
春风又绿江南岸，明月何时照我还。
3楼2017-06-27 21:02:38
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

vvaa

新虫 (初入文坛)
引用回帖:
3楼: Originally posted by huab1984666 at 2017-06-27 21:02:38
global c; c=;%实验值
global t; t=;实验值

c,t都是列向量...

你好，这个是matlab程序吗？是什么优化方法？没看懂

发自小木虫Android客户端
赞一下 回复此楼
4楼2017-06-27 21:57:42
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

huab1984666

新虫 (著名写手)
直接搜索，patternsearch

发自小木虫Android客户端
赞一下 回复此楼
春风又绿江南岸，明月何时照我还。
5楼2017-06-27 22:34:50
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

huab1984666

新虫 (著名写手)
引用回帖:
4楼: Originally posted by vvaa at 2017-06-27 21:57:42
你好，这个是matlab程序吗？是什么优化方法？没看懂
...

就是你测的浓度c,t是对应的。用你的模型来拟合，得到k1-7未知参数，就是最优化

发自小木虫Android客户端
赞一下 回复此楼
春风又绿江南岸，明月何时照我还。
6楼2017-06-27 22:36:05
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

vvaa

新虫 (初入文坛)
引用回帖:
6楼: Originally posted by huab1984666 at 2017-06-27 22:36:05
就是你测的浓度c,t是对应的。用你的模型来拟合，得到k1-7未知参数，就是最优化
...

嗯嗯，谢谢，我学习学习模式搜索法。还有一个问题，我每个实验条件下要做几个时间点才能求出来7个k值呢？

发自小木虫Android客户端
赞一下 回复此楼
7楼2017-06-28 09:26:25
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

huab1984666

新虫 (著名写手)
引用回帖:
7楼: Originally posted by vvaa at 2017-06-28 09:26:25
嗯嗯，谢谢，我学习学习模式搜索法。还有一个问题，我每个实验条件下要做几个时间点才能求出来7个k值呢？
...

这个看你自己啊。我觉得试验点越多越好
赞一下 回复此楼
春风又绿江南岸，明月何时照我还。
8楼2017-06-28 10:42:47
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

vvaa

新虫 (初入文坛)
引用回帖:
8楼: Originally posted by huab1984666 at 2017-06-28 10:42:47
这个看你自己啊。我觉得试验点越多越好...

实验点太多的话，实验量太大。根据这个微分方程组，你觉得最少要做几个点？我数学不太好，搞不太明白

发自小木虫Android客户端
赞一下 回复此楼
9楼2017-06-28 10:55:20
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

huab1984666

新虫 (著名写手)
引用回帖:
9楼: Originally posted by vvaa at 2017-06-28 10:55:20
实验点太多的话，实验量太大。根据这个微分方程组，你觉得最少要做几个点？我数学不太好，搞不太明白
...

你反应时间多长啊？我觉得30-40个很好啊，作图也好看
赞一下 回复此楼
春风又绿江南岸，明月何时照我还。
10楼2017-06-28 10:59:30
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 vvaa 的主题更新
161/2返回列表12下一页
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 求材料调剂，一志愿郑州大学289分 +16 硕星赴 2026-04-03 16/800 2026-04-07 20:23 by 云民大李老师
[考研] 一志愿南科大生物学297分，求调剂推荐 +8 Y-yyusx 2026-04-06 9/450 2026-04-07 19:38 by biomichael
[考研] 22408 调剂材料 +7 我叫ez 2026-04-06 8/400 2026-04-07 17:12 by 蓝云思雨
[考研] 生物与医药273求调剂 +14 荔题南墙 2026-04-05 15/750 2026-04-07 16:29 by 可口可乐不加冰�
[考研] 化学调剂 +18 艾志恒 2026-04-03 19/950 2026-04-07 16:00 by 起飞的比熊1
[考研] 333求调剂 +6 合乘杨习夕 2026-04-06 6/300 2026-04-07 09:44 by 猪会飞
[考研] 一志愿北交大材料工程总分358求调剂 +10 cs0106 2026-04-05 12/600 2026-04-06 19:41 by 无际的草原
[考研] 318求调剂 +12 ykyhsa 2026-04-05 14/700 2026-04-06 17:46 by fuyu_
[考研] 296求调剂 +3 汪！？！ 2026-04-05 4/200 2026-04-05 20:13 by 啵啵啵0119
[考研] 工科08专硕机械275求调剂 +3 AaAa7420 2026-04-02 3/150 2026-04-05 13:26 by jp9609
[考研] 323分（计算机视觉和大模型项目）能直接上手 +3 chaoxiicy 2026-04-01 3/150 2026-04-05 00:50 by chongya
[考研] 11408,335分，本科211，求调剂，可转专业 +5 鳄梨大鳄鱼 2026-04-03 5/250 2026-04-04 22:49 by chongya
[考研] 341求调剂 +3 洛多罗 2026-04-02 4/200 2026-04-04 21:36 by 智能智慧
[考研] 283求调剂 +4 mcbbc 2026-04-03 5/250 2026-04-04 20:51 by imissbao
[考研] 总分328生物与医药考数学求调剂 +7 aaadim 2026-04-02 9/450 2026-04-03 22:53 by syh9288
[考研] 工科341分调剂 +3 洛多罗 2026-04-03 3/150 2026-04-03 14:20 by 1753564080
[考研] 260求调剂 +3 朱芷琳 2026-04-02 3/150 2026-04-03 08:44 by yulian1987
[考研] 一志愿北京科技大学085601材料工程英一数二初试总分335求调剂 +8 双马尾痞老板2 2026-04-02 9/450 2026-04-02 14:45 by 5896
[考研] 能源动力 调剂 +3 不破不立0 2026-04-02 3/150 2026-04-02 12:46 by ffffjjjj
[考研] 一志愿厦门大学化学工程（专硕）-数二英二406分-求调剂 +5 厦大化工 2026-04-01 5/250 2026-04-02 10:03 by jp9609
信息提示
关闭
请填处理意见
关闭