 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 2 | 1/1 | 返回列表 |
| 查看: 835 | 回復(fù): 1 | |||
| 【懸賞金幣】回答本帖問題，作者szy242424將贈(zèng)送您 15 個(gè)金幣 | |||
szy242424新蟲 (初入文壇)
|
[求助]
1stop能模擬優(yōu)化12個(gè)參數(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的計(jì)算值','C1的實(shí)驗(yàn)值','C2的計(jì)算值','C2的實(shí)驗(yàn)值','C3的計(jì)算值','C3的實(shí)驗(yàn)值','C4的計(jì)算值','C4的實(shí)驗(yàn)值') 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);%第一個(gè)微分方程，等號(hào)前面是C(1)對(duì)t的微分 dC2dt =k(5)*dC1dt+k(6)*C(1);%第二個(gè)微分方程，等號(hào)前面是C(2)對(duì)t的微分 dC3dt =k(7)*dC1dt+k(8)*C(1);%第三個(gè)微分方程，等號(hào)前面是C(3)對(duì)t的微分 dC4dt =-(1/k(9)*dC1dt+k(10)*C(1)+1/k(11)*dC2dt+1/k(12)*dC3dt);%第四個(gè)微分方程，等號(hào)前面是C(4)對(duì)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的取值點(diǎn)，x0為微分方程組的初始值，k為微分方程的系數(shù)，返回t為微分方程組解的取值點(diǎn)，Xsim為微分方程組的解 Xsim1=Xsim(:,1);%提取Xsim的第一列 Xsim2=Xsim(:,2);%提取Xsim的第二列 Xsim3=Xsim(:,3);%提取Xsim的第三列 Xsim4=Xsim(:,4);%提取Xsim的第四列 ysim(:,1) = Xsim1(2:end);%微分方程的第一個(gè)變量，從第二個(gè)解到最后一個(gè)解賦值給ysim的第一列 ysim(:,2) = Xsim2(2:end);%微分方程的第二個(gè)變量，從第二個(gè)解到最后一個(gè)解賦值給ysim的第二列 ysim(:,3) = Xsim3(2:end);%微分方程的第三個(gè)變量，從第二個(gè)解到最后一個(gè)解賦值給ysim的第三列 ysim(:,4) = Xsim4(2:end);%微分方程的第四個(gè)變量，從第二個(gè)解到最后一個(gè)解賦值給ysim的第四列 f = [(ysim(:,1)-yexp(:,1)) (ysim(:,2)-yexp(:,2)) (ysim(:,3)-yexp(:,3)) (ysim(:,4)-yexp(:,4))];%形成目標(biāo)優(yōu)化函數(shù) |
回復(fù)此樓» 猜你喜歡
幫我的英語口語老師找學(xué)生
已經(jīng)有4人回復(fù)
-
計(jì)算三階力常數(shù)用thirdorder擴(kuò)胞時(shí)參數(shù)如何選擇？臨近原子數(shù)怎么取值
已經(jīng)有7人回復(fù)
-
物理學(xué)I論文潤(rùn)色/翻譯怎么收費(fèi)?
已經(jīng)有214人回復(fù)
-
2026年循環(huán)經(jīng)濟(jì)功能材料國(guó)際會(huì)議（ICFMCE 2026）
已經(jīng)有0人回復(fù)
-
2026年第五屆電氣、電子與信息工程國(guó)際會(huì)議（ISEEIE 2026）
已經(jīng)有0人回復(fù)
-
哈爾濱理工大學(xué)物理系招收物理學(xué)考研調(diào)劑
已經(jīng)有1人回復(fù)
-
0702一志愿吉大B區(qū)求調(diào)劑
已經(jīng)有2人回復(fù)
-
求調(diào)劑
已經(jīng)有0人回復(fù)
-
0702一志愿吉大B區(qū)求調(diào)劑有論文
已經(jīng)有0人回復(fù)
獨(dú)孤神宇
版主 (知名作家)
- 應(yīng)助: 490 (碩士)
- 貴賓: 0.008
- 金幣: 31016.3
- 散金: 802
- 紅花: 122
- 沙發(fā): 1
- 帖子: 5600
- 在線: 856.5小時(shí)
- 蟲號(hào): 3522474
- 注冊(cè): 2014-11-06
- 性別: GG
- 專業(yè): 機(jī)械動(dòng)力學(xué)
- 管轄: 計(jì)算模擬
2樓2020-12-08 19:45:28
相關(guān)版塊跳轉(zhuǎn)
我要訂閱樓主 szy242424 的主題更新
| 2 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 308求調(diào)劑 +3 | 是Lupa啊 2026-03-11 | 3/150 |
|
|---|---|---|---|---|
|
[考研] 接受26屆調(diào)劑生 +22 | 豬豬豬毛 2026-03-06 | 22/1100 |
|
|
|
[考研] 298求調(diào)劑 +3 | Vv呀！ 2026-03-10 | 3/150 |
|
|
|
[考研] 306求調(diào)劑 +3 | 唐薏薏 2026-03-09 | 3/150 |
|
|
|
[考研] 2026考研求調(diào)劑-材料類-本科211一志愿985-初試301分 +10 | 蟲友233 2026-03-07 | 10/500 |
|
|
|
[考研] 0817學(xué)碩華東區(qū)求調(diào)劑 +3 | 30660438 2026-03-08 | 3/150 |
|
|
|
[考研] 求調(diào)劑！ +3 | 朔朔話 2026-03-09 | 3/150 |
|
|
|
[考研] 327分求調(diào)劑086 +4 | 西紅柿?小帥 2026-03-09 | 7/350 |
|
|
|
[考博] 求材料讀博院校 +9 | yanglei131 2026-03-08 | 9/450 |
|
|
|
[考研] 085600材料與化工 326 求調(diào)劑 +4 | 熱愛生活ing 2026-03-09 | 4/200 |
|
|
|
[考研] 0832食品科學(xué)與工程293調(diào)劑 +3 | 東東不東 2026-03-07 | 3/150 |
|
|
|
[考研] 考研調(diào)劑，一志愿山東大學(xué)材料與化工，328分，政治51 +5 | 關(guān)你西紅柿929 2026-03-08 | 6/300 |
|
|
|
[考研] 0817找調(diào)劑 +6 | kk扛 2026-03-08 | 6/300 |
|
|
|
[教師之家] 交大前校長(zhǎng)王樹國(guó)：現(xiàn)在最先進(jìn)的科技并不在大學(xué)實(shí)驗(yàn)室，而是在企業(yè)研究院 +4 | zju2000 2026-03-08 | 6/300 |
|
|
|
[考研] 081700學(xué)碩一志愿北京化工大學(xué)數(shù)二英一過六級(jí)有競(jìng)賽求調(diào)劑 +5 | galaxary 2026-03-07 | 7/350 |
|
|
|
[考研]
 085600材料與化工 292分 求調(diào)劑
+6
|
程晴之 2026-03-06 | 6/300 |
|
|
|
[考研] 哈爾濱理工大學(xué)2026年研究生調(diào)劑，材料科學(xué)與化學(xué)工程學(xué)院研究生調(diào)劑 +3 | xinliu866 2026-03-06 | 3/150 |
|
|
|
[考研] 材料與化工354調(diào)劑 +4 | Lucy-xiao 2026-03-06 | 7/350 |
|
|
|
[考研] 求調(diào)劑，學(xué)校研究所都可以，材料與化工267分 +6 | wmx1 2026-03-05 | 6/300 |
|
|
|
[考研] 304求調(diào)劑 +3 | 曼殊2266 2026-03-04 | 3/150 |
|