 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 2 | 1/1 | 返回列表 |
| 查看: 833 | 回復(fù): 1 | ||
| 【懸賞金幣】回答本帖問題，作者szy242424將贈送您 15 個金幣 | ||
szy242424新蟲 (初入文壇)
|
[求助]
1stop能模擬優(yōu)化12個參數(shù)嗎，下面是我的matlab的程序，怎么用1stop優(yōu)化 已有1人參與
|
|
|
function kk format long e k0=[230,0.3,50,0.3,1.2,0.001,0.1,0.001,0.1,0.001,0.1,0.1]; lb=[230,0.3,50,0.3,1.2,0.001,0.1,0.001,0.1,0.001,0.1,0.1]; ub =[300,1,100,3,10.4, 1.2,2.4,1.2,0.8, 0.1,4.5,4.5]; data=... [0 34.2775 0 0 0.562205354 6 33.9805 0 0 0.618817941 12 31.941 0 0 0.797936454 18 28.739 0.9935 0 1.554008384 24 26.3835 1.3815 0 2.141789344 30 23.8955 1.5745 1.4795 2.543955264 36 21.334 1.908 1.7625 3.017273616 42 19.128 2.9885 2.038 3.295696176 60 11.2065 4.693 3.262 3.274041088 66 8.9005 5.1295 3.581 3.4070652 72 6.4395 5.6395 4.041 3.753546608 84 1.496 6.6735 4.71 3.595773824]; x0=data(1,2:end); tspan=data(:,1); yexp = [data(2:end,2) data(2:end,3) data(2:end,4) data(2:end,5)]; [k,resnorm,residual,nlparci,jacobian] =lsqnonlin(@ObjFunc,k0,lb,ub,[],tspan,x0,yexp); ts=[0,84]; [ts ,ys] = ode23s(@KineticsEqs,ts,x0,[],k); yy = [data(:,2) data(:,3) data(:,4) data(:,5)]; plot(ts,ys(:,1),'b',tspan,yy(:,1),'bo'); hold on plot(ts,ys(:,2),'r',tspan,yy(:,2),'r*'); plot(ts,ys(:,3),'k',tspan,yy(:,3),'k+'); plot(ts,ys(:,4),'g',tspan,yy(:,4),'g+'); legend('C1的計算值','C1的實驗值','C2的計算值','C2的實驗值','C3的計算值','C3的實驗值','C4的計算值','C4的實驗值') end function dCdt = KineticsEqs(t ,C,k) % ODE模型方程，C為濃度，k為微分方程中的系數(shù)，t為導(dǎo)數(shù) dC1dt =0.19*C(4)/(C(4)+0.78+C(4)^2/300.8)*(1-C(2)/k(1))^k(2)*(1-C(3)/k(3))^k(4)*C(1);%第一個微分方程，等號前面是C(1)對t的微分 dC2dt =k(5)*dC1dt+k(6)*C(1);%第二個微分方程，等號前面是C(2)對t的微分 dC3dt =k(7)*dC1dt+k(8)*C(1);%第三個微分方程，等號前面是C(3)對t的微分 dC4dt =-(1/k(9)*dC1dt+k(10)*C(1)+1/k(11)*dC2dt+1/k(12)*dC3dt);%第四個微分方程，等號前面是C(4)對t的微分 dCdt = [dC1dt; dC2dt;dC3dt;dC4dt];%微分方程組 function f = ObjFunc(k0,tspan,x0,yexp) % 目標(biāo)函數(shù) [t, Xsim] = ode23s(@KineticsEqs,tspan,x0,[],k0);% 求解常微分方程，其中tspan為t的取值點，x0為微分方程組的初始值，k為微分方程的系數(shù)，返回t為微分方程組解的取值點，Xsim為微分方程組的解 Xsim1=Xsim(:,1);%提取Xsim的第一列 Xsim2=Xsim(:,2);%提取Xsim的第二列 Xsim3=Xsim(:,3);%提取Xsim的第三列 Xsim4=Xsim(:,4);%提取Xsim的第四列 ysim(:,1) = Xsim1(2:end);%微分方程的第一個變量，從第二個解到最后一個解賦值給ysim的第一列 ysim(:,2) = Xsim2(2:end);%微分方程的第二個變量，從第二個解到最后一個解賦值給ysim的第二列 ysim(:,3) = Xsim3(2:end);%微分方程的第三個變量，從第二個解到最后一個解賦值給ysim的第三列 ysim(:,4) = Xsim4(2:end);%微分方程的第四個變量，從第二個解到最后一個解賦值給ysim的第四列 f = [(ysim(:,1)-yexp(:,1)) (ysim(:,2)-yexp(:,2)) (ysim(:,3)-yexp(:,3)) (ysim(:,4)-yexp(:,4))];%形成目標(biāo)優(yōu)化函數(shù) |
回復(fù)此樓» 猜你喜歡
德國亥姆霍茲Hereon中心招收兩位醫(yī)用鎂合金腐蝕與LPSO相變方向2026公派博士生
已經(jīng)有0人回復(fù)
-
推薦一款可以AI輔助寫作的Latex編輯器SmartLatexEditor，超級好用，推薦試試
已經(jīng)有11人回復(fù)
-
物理學(xué)I論文潤色/翻譯怎么收費?
已經(jīng)有296人回復(fù)
-
2026-CJ開始申報了
已經(jīng)有1人回復(fù)
-
西安電子科學(xué)大學(xué)杭州研究院劉麗香教授招收智能多模態(tài)傳感器和微型儲能器件方向博士
已經(jīng)有13人回復(fù)
-
重慶交通大學(xué)光子學(xué)微結(jié)構(gòu)與器件課題組招收2026年碩士研究生信息
已經(jīng)有1人回復(fù)
-
一志愿鄭大材料學(xué)碩298分，求調(diào)劑
已經(jīng)有6人回復(fù)
-
尋合作：應(yīng)力腐蝕多尺度模擬
已經(jīng)有2人回復(fù)
-
考研交流
已經(jīng)有0人回復(fù)
獨孤神宇
版主 (知名作家)
- 應(yīng)助: 490 (碩士)
- 貴賓: 0.008
- 金幣: 31016.3
- 散金: 802
- 紅花: 122
- 沙發(fā): 1
- 帖子: 5600
- 在線: 856.5小時
- 蟲號: 3522474
- 注冊: 2014-11-06
- 性別: GG
- 專業(yè): 機械動力學(xué)
- 管轄: 計算模擬
2樓2020-12-08 19:45:28
相關(guān)版塊跳轉(zhuǎn)
我要訂閱樓主 szy242424 的主題更新
| 2 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 085600 材料與化工 298 +11 | 小西笑嘻嘻 2026-03-03 | 11/550 |
|
|---|---|---|---|---|
|
[考研] 085701環(huán)境工程 求調(diào)劑 +5 | xiiiia 2026-03-04 | 5/250 |
|
|
|
[考研] 325求調(diào)劑 +3 | 學(xué)家科 2026-03-04 | 3/150 |
|
|
|
[考研] 復(fù)試調(diào)劑 +4 | 杰尼龜aaa 2026-03-03 | 4/200 |
|
|
|
[考研] 中國林科院林化所（南京）2026年招收化學(xué)/材料/環(huán)境工程等背景碩士研究生3名 +3 | realstar2006 2026-02-27 | 3/150 |
|
|
|
[考研] 266材料化工求調(diào)劑 +3 | 哇塞王帥 2026-03-03 | 3/150 |
|
|
|
[考研] 江蘇省農(nóng)科院招調(diào)劑1名 +5 | Qwertyuop 2026-03-01 | 5/250 |
|
|
|
[基金申請] 沒有青基直接申請面上，感覺自己瘋了 +5 | kevin63t 2026-03-02 | 6/300 |
|
|
|
[考研] 化工專碩348，一志愿985求調(diào)劑 +8 | 弗格個 2026-02-28 | 11/550 |
|
|
|
[考研] 一志愿中科大能動297求調(diào)劑，本科川大 +3 | 邵11 2026-03-03 | 3/150 |
|
|
|
[考研] 267求調(diào)劑 +6 | 釣魚佬as 2026-03-02 | 6/300 |
|
|
|
[考研] 調(diào)劑 +5 | 13853210211 2026-03-02 | 7/350 |
|
|
|
[考研] 0856調(diào)劑 +10 | 劉夢微 2026-02-28 | 10/500 |
|
|
|
[考研] 材料284求調(diào)劑，一志愿鄭州大學(xué)英一數(shù)二專碩 +15 | 想上岸的土撥鼠 2026-02-28 | 15/750 |
|
|
|
[考研] 一志愿東北大學(xué)材料專碩328，求調(diào)劑 +3 | shs1083 2026-03-02 | 3/150 |
|
|
|
[考研] 295求調(diào)劑 +8 | 19171856320 2026-02-28 | 8/400 |
|
|
|
[基金申請] 此成果不能導(dǎo)入原因：元數(shù)據(jù)必填信息不完整，可 進行補充。 +4 | Kittylucky 2026-03-02 | 5/250 |
|
|
|
[考研] 材料類求調(diào)劑 +11 | wana_kiko 2026-02-28 | 14/700 |
|
|
|
[考研] 311求調(diào)劑 +9 | 南迦720 2026-02-28 | 10/500 |
|
|
|
[論文投稿]
求助coordination chemistry reviews 的寫作模板
 10+3
|
ljplijiapeng 2026-02-27 | 4/200 |
|