 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 查看: 2175 | 回復: 15 | |||
[求助]
最優(yōu)化方法求教�。�！ 已有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 因為方程之間耦聯(lián)，不能得到解析解。用bfgs法優(yōu)化求解參數時，導函數不知道怎么求解。本人小白，沒有最優(yōu)化基礎，希望大神能指導一下，怎么求dc/k，或者有沒有比較好的不用導函數的優(yōu)化方法推薦？ 發(fā)自小木蟲Android客戶端 |
回復此樓» 收錄本帖的淘帖專輯推薦
 優(yōu)化 |
» 猜你喜歡
2026年第五屆電氣、電子與信息工程國際會議（ISEEIE 2026）
已經有0人回復
-
哈爾濱理工大學物理系招收物理學考研調劑
已經有5人回復
-
物理學I論文潤色/翻譯怎么收費?
已經有79人回復
-
0702一志愿吉大B區(qū)求調劑
已經有5人回復
-
求調劑
已經有0人回復
-
0702一志愿吉大B區(qū)求調劑有論文
已經有0人回復
-
請問還有沒有用Latex寫文章的小伙伴們？
已經有0人回復
-
光學工程學碩調劑信息
已經有30人回復
-
歡迎加入課題組
已經有0人回復
-
散金幣，求好運，祝面上順利！
已經有62人回復
-
華東師范大學芯片設計徐瓏真實水平咋樣
已經有1人回復
» 本主題相關價值貼推薦，對您同樣有幫助:
- Materials Studio 7.0 計算PGA(一種高分子)密度，最終目的期望得到好的PGA分子結構
已經有2人回復
- 多年的醫(yī)學課題設計得出一些自己的切身體會
已經有1人回復
- 易錯PCR效果
已經有7人回復
- 催化版上2008-04-25前 沸石分子篩，ZSM-系、SAPO-系 共466個話題集錦
已經有110人回復
huab1984666
新蟲 (著名寫手)
- 應助: 9 (幼兒園)
- 金幣: 241.5
- 散金: 2094
- 紅花: 12
- 帖子: 1320
- 在線: 210.3小時
- 蟲號: 516828
- 注冊: 2008-03-03
- 性別: GG
- 專業(yè): 膠體與界面化學
【答案】應助回帖
★ ★ ★ ★ ★
感謝參與，應助指數 +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
- 專業(yè): 膠體與界面化學
3樓2017-06-27 21:02:38
4樓2017-06-27 21:57:42
huab1984666
新蟲 (著名寫手)
- 應助: 9 (幼兒園)
- 金幣: 241.5
- 散金: 2094
- 紅花: 12
- 帖子: 1320
- 在線: 210.3小時
- 蟲號: 516828
- 注冊: 2008-03-03
- 性別: GG
- 專業(yè): 膠體與界面化學
5樓2017-06-27 22:34:50
huab1984666
新蟲 (著名寫手)
- 應助: 9 (幼兒園)
- 金幣: 241.5
- 散金: 2094
- 紅花: 12
- 帖子: 1320
- 在線: 210.3小時
- 蟲號: 516828
- 注冊: 2008-03-03
- 性別: GG
- 專業(yè): 膠體與界面化學
6樓2017-06-27 22:36:05
7樓2017-06-28 09:26:25
huab1984666
新蟲 (著名寫手)
- 應助: 9 (幼兒園)
- 金幣: 241.5
- 散金: 2094
- 紅花: 12
- 帖子: 1320
- 在線: 210.3小時
- 蟲號: 516828
- 注冊: 2008-03-03
- 性別: GG
- 專業(yè): 膠體與界面化學
8樓2017-06-28 10:42:47
9樓2017-06-28 10:55:20
huab1984666
新蟲 (著名寫手)
- 應助: 9 (幼兒園)
- 金幣: 241.5
- 散金: 2094
- 紅花: 12
- 帖子: 1320
- 在線: 210.3小時
- 蟲號: 516828
- 注冊: 2008-03-03
- 性別: GG
- 專業(yè): 膠體與界面化學
10樓2017-06-28 10:59:30
相關版塊跳轉
我要訂閱樓主 vvaa 的主題更新
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 0817 化學工程 299分求調劑 有科研經歷 有二區(qū)文章 +14 | rare12345 2026-03-18 | 14/700 |
|
|---|---|---|---|---|
|
[考研] 能源材料化學課題組招收碩士研究生8-10名 +4 | 脫穎而出 2026-03-16 | 12/600 |
|
|
|
[考研]
 271材料工程求調劑
+6
|
.6lL 2026-03-18 | 6/300 |
|
|
|
[考研] 0703化學調劑 +5 | pupcoco 2026-03-17 | 8/400 |
|
|
|
[考研] 332求調劑 +3 | ydfyh 2026-03-17 | 3/150 |
|
|
|
[考研] 0817調劑 +3 | 沒有答案_ 2026-03-14 | 3/150 |
|
|
|
[考研] 311求調劑 +4 | 冬十三 2026-03-18 | 4/200 |
|
|
|
[考研] 【同濟軟件】軟件（085405）考研求調劑 +3 | 2026eternal 2026-03-18 | 3/150 |
|
|
|
[考研] 311求調劑 +11 | 冬十三 2026-03-15 | 12/600 |
|
|
|
[考研] 312求調劑 +8 | 陌宸希 2026-03-16 | 9/450 |
|
|
|
[考研] 0703化學求調劑 總分331 +3 | ZY-05 2026-03-13 | 3/150 |
|
|
|
[考博] 26博士申請 +3 | 1042136743 2026-03-17 | 3/150 |
|
|
|
[考研] 326求調劑 +5 | 上岸的小葡 2026-03-15 | 6/300 |
|
|
|
[考研] 材料專碩326求調劑 +6 | 墨煜姒莘 2026-03-15 | 7/350 |
|
|
|
[考博] 26申博 +4 | 八6八68 2026-03-16 | 4/200 |
|
|
|
[考研] 一志愿南京大學，080500材料科學與工程，調劑 +4 | Jy? 2026-03-16 | 4/200 |
|
|
|
[考研] 283求調劑 +3 | 聽風就是雨； 2026-03-16 | 3/150 |
|
|
|
[考研] 0856求調劑 +3 | 劉夢微 2026-03-15 | 3/150 |
|
|
|
[考研] 本科南京大學一志愿川大藥學327 +3 | 麥田耕者 2026-03-14 | 3/150 |
|
|
|
[考研] 26調劑/材料科學與工程/總分295/求收留 +9 | 2026調劑俠 2026-03-12 | 9/450 |
|