 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 5 | 1/1 | 返回列表 |
| 查看: 2190 | 回復(fù): 15 | |||
| 當(dāng)前只顯示滿足指定條件的回帖，點擊這里查看本話題的所有回帖 | |||
vvaa新蟲 (初入文壇)
|
[求助]
最優(yōu)化方法求教�。�！ 已有1人參與
|
||
|
我的課題要做一個動力學(xué)模型，目前打算用集總的方法，建立一個微分方程組，比如 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)化求解參數(shù)時，導(dǎo)函數(shù)不知道怎么求解。本人小白，沒有最優(yōu)化基礎(chǔ)，希望大神能指導(dǎo)一下，怎么求dc/k，或者有沒有比較好的不用導(dǎo)函數(shù)的優(yōu)化方法推薦？ 發(fā)自小木蟲Android客戶端 |
回復(fù)此樓» 收錄本帖的淘帖專輯推薦
 優(yōu)化 |
» 猜你喜歡
中德博士后2026項目5月30日截止- 國內(nèi)博士后入站、德國兩年全職
已經(jīng)有0人回復(fù)
-
面向全球申請人的中德博士后2026項目5月30日截止- 國內(nèi)博后入站德國兩年全職
已經(jīng)有0人回復(fù)
-
物理學(xué)I論文潤色/翻譯怎么收費?
已經(jīng)有120人回復(fù)
-
【新加坡】納米電子器件項目組有“聯(lián)合培養(yǎng)博士生”名額
已經(jīng)有2人回復(fù)
-
重慶交大材料學(xué)院26年碩士招生復(fù)試實施細則已出。歡迎加入光子學(xué)微結(jié)構(gòu)與器件課題組。
已經(jīng)有0人回復(fù)
-
北工大宋曉艷教授團隊招收2026年9月入學(xué)博士【學(xué)博和專博名額均有】！
已經(jīng)有1人回復(fù)
-
天津城建大學(xué)理學(xué)院化學(xué)、環(huán)境工程、材料物理等方向接收10人
已經(jīng)有0人回復(fù)
-
天津城建大學(xué)理學(xué)院化學(xué)、環(huán)境工程、材料物理等方向接收調(diào)劑，10人
已經(jīng)有0人回復(fù)
-
計算與實驗科學(xué)工程國際會議(8月7-11日,香港) 二維材料分會S5與聲學(xué)換能器分會S8
已經(jīng)有0人回復(fù)
-
0702一志愿吉大B區(qū)求調(diào)劑有論文
已經(jīng)有2人回復(fù)
» 本主題相關(guān)價值貼推薦，對您同樣有幫助:
- Materials Studio 7.0 計算PGA(一種高分子)密度，最終目的期望得到好的PGA分子結(jié)構(gòu)
已經(jīng)有2人回復(fù)
- 多年的醫(yī)學(xué)課題設(shè)計得出一些自己的切身體會
已經(jīng)有1人回復(fù)
- 易錯PCR效果
已經(jīng)有7人回復(fù)
- 催化版上2008-04-25前 沸石分子篩，ZSM-系、SAPO-系 共466個話題集錦
已經(jīng)有110人回復(fù)
vvaa
新蟲 (初入文壇)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 293.6
- 帖子: 30
- 在線: 12.6小時
- 蟲號: 3912957
- 注冊: 2015-06-07
- 專業(yè): 化學(xué)反應(yīng)工程
9樓2017-06-28 10:55:20
huab1984666
新蟲 (著名寫手)
- 應(yīng)助: 9 (幼兒園)
- 金幣: 241.5
- 散金: 2094
- 紅花: 12
- 帖子: 1320
- 在線: 210.3小時
- 蟲號: 516828
- 注冊: 2008-03-03
- 性別: GG
- 專業(yè): 膠體與界面化學(xué)
【答案】應(yīng)助回帖
★ ★ ★ ★ ★
感謝參與，應(yīng)助指數(shù) +1
vvaa: 金幣+5, ★★★很有幫助 2017-06-28 09:27:15
感謝參與，應(yīng)助指數(shù) +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
新蟲 (著名寫手)
- 應(yīng)助: 9 (幼兒園)
- 金幣: 241.5
- 散金: 2094
- 紅花: 12
- 帖子: 1320
- 在線: 210.3小時
- 蟲號: 516828
- 注冊: 2008-03-03
- 性別: GG
- 專業(yè): 膠體與界面化學(xué)
3樓2017-06-27 21:02:38
vvaa
新蟲 (初入文壇)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 293.6
- 帖子: 30
- 在線: 12.6小時
- 蟲號: 3912957
- 注冊: 2015-06-07
- 專業(yè): 化學(xué)反應(yīng)工程
4樓2017-06-27 21:57:42
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 一志愿北京化工大學(xué)材料與化工（085600）296求調(diào)劑 +22 | 稻妻小編 2026-03-26 | 22/1100 |
|
|---|---|---|---|---|
|
[考研] 材料科學(xué)與工程 317求調(diào)劑 +5 | JKSOIID 2026-03-26 | 5/250 |
|
|
|
[考研] 一志愿華東師范大學(xué)有機化學(xué)專業(yè)，初試351分，復(fù)試被刷求調(diào)劑! +6 | 真名有冰 2026-03-29 | 7/350 |
|
|
|
[考研] 化學(xué)工程085602 305分求調(diào)劑 +25 | RichLi_ 2026-03-25 | 25/1250 |
|
|
|
[考研] 一志愿南開大學(xué)0710生物學(xué)359求調(diào)劑 +5 | 兔兔兔111223314 2026-03-29 | 5/250 |
|
|
|
[考研] 290求調(diào)劑 +3 | dfffsar 2026-03-29 | 3/150 |
|
|
|
[考研] 343求調(diào)劑 +6 | 愛羈絆 2026-03-29 | 6/300 |
|
|
|
[考研] 生物學(xué)學(xué)碩，一志愿湖南大學(xué)，初試成績338 +6 | YYYYYNNNNN 2026-03-26 | 7/350 |
|
|
|
[考研] 083000學(xué)碩274求調(diào)劑 +8 | Li李魚 2026-03-26 | 8/400 |
|
|
|
[考研] 312，生物學(xué)求調(diào)劑 +3 | 小譯同學(xué)abc 2026-03-28 | 3/150 |
|
|
|
[考研] 266分，求材料冶金能源化工等調(diào)劑 +7 | 哇呼哼呼哼 2026-03-27 | 9/450 |
|
|
|
[考研] 藥學(xué)105500求調(diào)劑 +3 | Ssun。。 2026-03-28 | 3/150 |
|
|
|
[考研] 081200-314 +3 | LILIQQ 2026-03-27 | 4/200 |
|
|
|
[考研] 070300化學(xué)求調(diào)劑 +4 | 起個名咋這么難 2026-03-27 | 4/200 |
|
|
|
[考研] 336材料求調(diào)劑 +7 | 陳瀅瑩 2026-03-26 | 9/450 |
|
|
|
[考研] 中國科學(xué)院深圳先進技術(shù)研究院-光纖傳感課題組招生-中國科學(xué)院大學(xué)、深圳理工大學(xué)聯(lián)培 +5 | YangTyu1 2026-03-26 | 5/250 |
|
|
|
[考研] 機械學(xué)碩310分，數(shù)一英一，一志愿211本科雙非找調(diào)劑信息 +3 | @357 2026-03-25 | 3/150 |
|
|
|
[考研] 080500求調(diào)劑 +3 | zzzzfan 2026-03-24 | 3/150 |
|
|
|
[考研] 277分求調(diào)劑，跨調(diào)材料 +3 | 考研調(diào)劑lxh 2026-03-24 | 3/150 |
|
|
|
[考研]
 085404電子信息284分求調(diào)劑
+4
|
13659058978 2026-03-24 | 4/200 |
|