| 查看: 2460 | 回复: 15 | |||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | |||
[求助]
最优化方法求教！！！ 已有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客户端 |
回复此楼» 收录本帖的淘帖专辑推荐
 优化 |
» 猜你喜欢
PRB投稿，Acknowledgment sent to author状态十天了
已经有3人回复
-
PVA溶解
已经有2人回复
-
物理学I论文润色/翻译怎么收费?
已经有70人回复
-
基金申请
已经有44人回复
-
CSC与新西兰维多利亚大学PhD奖学金项目
已经有0人回复
-
新西兰Robinson研究所 招聘CSC公派访问人员
已经有0人回复
-
帮我的英语口语老师找学生
已经有0人回复
-
什么时候开奖？
已经有13人回复
-
散金币祈福
已经有94人回复
-
青基已中
已经有1人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- Materials Studio 7.0 计算PGA(一种高分子)密度，最终目的期望得到好的PGA分子结构
已经有2人回复
- 多年的医学课题设计得出一些自己的切身体会
已经有1人回复
- 易错PCR效果
已经有7人回复
- 催化版上2008-04-25前 沸石分子筛，ZSM-系、SAPO-系 共466个话题集锦
已经有110人回复
7楼2017-06-28 09:26:25
huab1984666
新虫 (著名写手)
- 应助: 9 (幼儿园)
- 金币: 241.5
- 散金: 2094
- 红花: 12
- 帖子: 1320
- 在线: 210.3小时
- 虫号: 516828
- 注册: 2008-03-03
- 性别: GG
- 专业: 胶体与界面化学
【答案】应助回帖
★ ★ ★ ★ ★
感谢参与，应助指数 +1
vvaa: 金币+5, ★★★很有帮助 2017-06-28 09:27:15
感谢参与，应助指数 +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-'); |
2楼2017-06-27 20:53:54
huab1984666
新虫 (著名写手)
- 应助: 9 (幼儿园)
- 金币: 241.5
- 散金: 2094
- 红花: 12
- 帖子: 1320
- 在线: 210.3小时
- 虫号: 516828
- 注册: 2008-03-03
- 性别: GG
- 专业: 胶体与界面化学
3楼2017-06-27 21:02:38
4楼2017-06-27 21:57:42
