版塊導(dǎo)航: 正在加載中...

登錄注冊

應(yīng)《網(wǎng)絡(luò)安全法》要求，自2017年10月1日起，未進(jìn)行實(shí)名認(rèn)證將不得使用互聯(lián)網(wǎng)跟帖服務(wù)。為保障您的帳號能夠正常使用，請盡快對帳號進(jìn)行手機(jī)號驗(yàn)證，感謝您的理解與支持！

24小時(shí)熱門版塊排行榜

返回列表

For_study

金蟲 (小有名氣)

木蟲

應(yīng)助: 0 (幼兒園)
金幣: 1089.8
散金: 506
帖子: 174
在線: 108.8小時(shí)
蟲號: 2085809
注冊: 2012-10-25
性別: MM
專業(yè): 環(huán)境化工

[求助] 遺傳算法優(yōu)化微分方程中的參數(shù) 已有1人參與

%擬合微分方程參數(shù)
%直接擬合反應(yīng)的活化能和指前因子
function mynewtry4
clear;
clc;

%遺傳算法
lb1=5*ones(1,22); lb2=-100000*ones(1,22); lb=[lb1,lb2]; ub1=100*ones(1,22); ub2=-1000*ones(1,22); ub=[ub1,ub2];
options = gaoptimset('Generations',1000,'StallGenLimit',50,...
'StallTimeLimit',50,'TolFun',1e-12,'TolCon',1e-12,'MutationFcn',@mutationadaptfeasible);
[h0,fval,exitflag,reason,output,final_pop]=ga(@my_funtest1,44,options);

%輸出參數(shù)
fprintf('\n\n遺傳算法的估計(jì)數(shù)值:\n');
disp(h0);

%構(gòu)造適應(yīng)度函數(shù)my_funtest
function yfit=my_funtest1(h)
tspan1=0:0.05:0.5;
y01=[0.561,0.298,0.141,0 0 0 0 0];
[t,ycal1]=ode45(@myfun1,tspan1,y01,[],h);

y02=ycal1(11,

;
tspan2=0.5:0.05:1;
[t,ycal2]=ode45(@myfun2,tspan2,y02,[],h);
ycalculate=ycal2(11,

yreal2=[0.0587 0 0 0.1988 0.3889 0.2309 0.0397 0.083];
for i=4:8
ff(i)=(ycalculate(i)-yreal2(i))^2;
end
yfit=sum(ff)
end

function dy = myfun1(t,y,h)

% 動力學(xué)微分方程
% 共分為8集總，原料（飽和分SS為1，芳香分SA為2，膠質(zhì)瀝青質(zhì)SR為3），柴油（D）為4,汽油（G）為5，液化氣（LPG）為6，干氣（Gas）為7，焦炭（C）為8，
% 由于遺傳算法只能返回向量，所以反應(yīng)常數(shù)不能是矩陣，只能是向量，第一個(gè)數(shù)表示反應(yīng)物，第二個(gè)數(shù)表示生成物
%反應(yīng)常數(shù)k（1，4）=k（1），k（1，5）=k（2），k（1，6）=k（3），k（1，7）=k（4），k（1，8）=k（5），k（2，4）=k（6），
%  k（2，5）=k（7），k（2，6）=k（8），k（2，7）=k（9），k（2，8）=k（10），k（3，4）=k（11），k（3，5）=k（12），
%  k（3，6）=k（13），k（3，7）=k（14），k（3，8）=k（15），k（4，5）=k（16），k（4，6）=k（17），k（4，7）=k（18），
%  k（4，8）=k（19），k（5，6）=k（20），k（5，7）=k（21），k（5，8）=k（22）
%
%
A=1;N=1;
deact=1;
dens=13.357;  %單位是kg/m^3
Swh=90; %單位s
Tem=788; R=8.314;

k=zeros(1,22);
k(1)=h(1)*exp(-h(23)/(R*Tem));  % R常數(shù)，8.314，Tem溫度，h表示活化能23-44，J/mol 指前因子1-22 kg/(m^3.s)
k(2)=h(2)*exp(-h(24)/(R*Tem));
k(3)=h(3)*exp(-h(25)/(R*Tem));
k(4)=h(4)*exp(-h(26)/(R*Tem));
k(5)=h(5)*exp(-h(27)/(R*Tem));
k(6)=h(6)*exp(-h(28)/(R*Tem));
k(7)=h(7)*exp(-h(29)/(R*Tem));
k(8)=h(8)*exp(-h(30)/(R*Tem));
k(9)=h(9)*exp(-h(31)/(R*Tem));
k(10)=h(10)*exp(-h(32)/(R*Tem));
k(11)=h(11)*exp(-h(33)/(R*Tem));
k(12)=h(12)*exp(-h(34)/(R*Tem));
k(13)=h(13)*exp(-h(35)/(R*Tem));

k(14)=h(14)*exp(-h(36)/(R*Tem));
k(15)=h(15)*exp(-h(37)/(R*Tem));
k(16)=h(16)*exp(-h(38)/(R*Tem));
k(17)=h(17)*exp(-h(39)/(R*Tem));
k(18)=h(18)*exp(-h(40)/(R*Tem));
k(19)=h(19)*exp(-h(41)/(R*Tem));

k(20)=h(20)*exp(-h(42)/(R*Tem));
k(21)=h(21)*exp(-h(43)/(R*Tem));
k(22)=h(22)*exp(-h(44)/(R*Tem));

dy(1)=-(k(1)+k(2)+k(3)+k(4)+k(5))*y(1)*A*N*deact*dens/Swh;  % A重芳烴失活系數(shù)，N堿氮吸附失活系數(shù)，deact催化劑結(jié)焦失活系數(shù)，dens密度，Swh真實(shí)重時(shí)空速
dy(2)=-(k(6)+k(7)+k(8)+k(9)+k(10))*y(2)*A*N*deact*dens/Swh;
dy(3)=-(k(11)+k(12)+k(13)+k(14)+k(15))*y(3)*A*N*deact*dens/Swh;
dy(4)=(k(1)*y(1)+k(6)*y(2)+k(11)*y(3)-(k(16)+k(17)+k(18)+k(19))*y(4))*A*N*deact*dens/Swh;
dy(5)=(k(2)*y(1)+k(7)*y(2)+k(12)*y(3)+k(16)*y(4)-(k(20)+k(21)+k(22))*y(5))*A*N*deact*dens/Swh;
dy(6)=(k(3)*y(1)+k(8)*y(2)+k(13)*y(3)+k(17)*y(4)+k(20)*y(5))*A*N*deact*dens/Swh;
dy(7)=(k(4)*y(1)+k(9)*y(2)+k(14)*y(3)+k(18)*y(4)+k(21)*y(5))*A*N*deact*dens/Swh;
dy(8)=(k(5)*y(1)+k(10)*y(2)+k(15)*y(3)+k(19)*y(4)+k(22)*y(5))*A*N*deact*dens/Swh;
dy=dy';

end

function dy = myfun2(t,y,h)

% 動力學(xué)微分方程
% 共分為8集總，原料（飽和分SS為1，芳香分SA為2，膠質(zhì)瀝青質(zhì)SR為3），柴油（D）為4,汽油（G）為5，液化氣（LPG）為6，干氣（Gas）為7，焦炭（C）為8，
% 由于遺傳算法只能返回向量，所以反應(yīng)常數(shù)不能是矩陣，只能是向量，第一個(gè)數(shù)表示反應(yīng)物，第二個(gè)數(shù)表示生成物
%反應(yīng)常數(shù)k（1，4）=k（1），k（1，5）=k（2），k（1，6）=k（3），k（1，7）=k（4），k（1，8）=k（5），k（2，4）=k（6），
%  k（2，5）=k（7），k（2，6）=k（8），k（2，7）=k（9），k（2，8）=k（10），k（3，4）=k（11），k（3，5）=k（12），
%  k（3，6）=k（13），k（3，7）=k（14），k（3，8）=k（15），k（4，5）=k（16），k（4，6）=k（17），k（4，7）=k（18），
%  k（4，8）=k（19），k（5，6）=k（20），k（5，7）=k（21），k（5，8）=k（22）
%
%
A=1;N=1;
deact=1;
dens=13.357;Swh=90;
Tem=761; R=8.314;

k=zeros(1,22);
k(1)=h(1)*exp(-h(23)/(R*Tem));  % R常數(shù)，8.314，Tem溫度，h表示活化能1-22，之前因子23-44
k(2)=h(2)*exp(-h(24)/(R*Tem));
k(3)=h(3)*exp(-h(25)/(R*Tem));
k(4)=h(4)*exp(-h(26)/(R*Tem));
k(5)=h(5)*exp(-h(27)/(R*Tem));
k(6)=h(6)*exp(-h(28)/(R*Tem));
k(7)=h(7)*exp(-h(29)/(R*Tem));
k(8)=h(8)*exp(-h(30)/(R*Tem));
k(9)=h(9)*exp(-h(31)/(R*Tem));
k(10)=h(10)*exp(-h(32)/(R*Tem));
k(11)=h(11)*exp(-h(33)/(R*Tem));
k(12)=h(12)*exp(-h(34)/(R*Tem));
k(13)=h(13)*exp(-h(35)/(R*Tem));

k(14)=h(14)*exp(-h(36)/(R*Tem));
k(15)=h(15)*exp(-h(37)/(R*Tem));
k(16)=h(16)*exp(-h(38)/(R*Tem));
k(17)=h(17)*exp(-h(39)/(R*Tem));
k(18)=h(18)*exp(-h(40)/(R*Tem));
k(19)=h(19)*exp(-h(41)/(R*Tem));

k(20)=h(20)*exp(-h(42)/(R*Tem));
k(21)=h(21)*exp(-h(43)/(R*Tem));
k(22)=h(22)*exp(-h(44)/(R*Tem));

dy(1)=-(k(1)+k(2)+k(3)+k(4)+k(5))*y(1)*A*N*deact*dens/Swh;  % A重芳烴失活系數(shù)，N堿氮吸附失活系數(shù)，deact催化劑結(jié)焦失活系數(shù)，dens密度，Swh真實(shí)重時(shí)空速
dy(2)=-(k(6)+k(7)+k(8)+k(9)+k(10))*y(2)*A*N*deact*dens/Swh;
dy(3)=-(k(11)+k(12)+k(13)+k(14)+k(15))*y(3)*A*N*deact*dens/Swh;
dy(4)=(k(1)*y(1)+k(6)*y(2)+k(11)*y(3)-(k(16)+k(17)+k(18)+k(19))*y(4))*A*N*deact*dens/Swh;
dy(5)=(k(2)*y(1)+k(7)*y(2)+k(12)*y(3)+k(16)*y(4)-(k(20)+k(21)+k(22))*y(5))*A*N*deact*dens/Swh;
dy(6)=(k(3)*y(1)+k(8)*y(2)+k(13)*y(3)+k(17)*y(4)+k(20)*y(5))*A*N*deact*dens/Swh;
dy(7)=(k(4)*y(1)+k(9)*y(2)+k(14)*y(3)+k(18)*y(4)+k(21)*y(5))*A*N*deact*dens/Swh;
dy(8)=(k(5)*y(1)+k(10)*y(2)+k(15)*y(3)+k(19)*y(4)+k(22)*y(5))*A*N*deact*dens/Swh;
dy=dy';

end
end

問題敘述，遺傳算法未設(shè)定lb和ub時(shí)，程序可以正常運(yùn)行，獲得結(jié)果，當(dāng)設(shè)定lb和ub后，程序依然可以運(yùn)行（不會報(bào)錯(cuò)），但是只能運(yùn)算一次優(yōu)化后的結(jié)果（設(shè)置斷點(diǎn)時(shí)發(fā)現(xiàn)）。遺傳算法好像沒有繼續(xù)優(yōu)化下去，希望有大神可以幫忙解答一下問題。。。。謝謝。。。。

回復(fù)此樓

» 猜你喜歡

哈爾濱理工大學(xué)物理系招收物理學(xué)考研調(diào)劑已經(jīng)有7人回復(fù)
0702一志愿吉大B區(qū)求調(diào)劑已經(jīng)有5人回復(fù)
物理學(xué)I論文潤色/翻譯怎么收費(fèi)? 已經(jīng)有88人回復(fù)
求調(diào)劑已經(jīng)有0人回復(fù)
0702一志愿吉大B區(qū)求調(diào)劑有論文已經(jīng)有0人回復(fù)
請問還有沒有用Latex寫文章的小伙伴們？已經(jīng)有0人回復(fù)
光學(xué)工程學(xué)碩調(diào)劑信息已經(jīng)有31人回復(fù)
歡迎加入課題組已經(jīng)有0人回復(fù)
散金幣，求好運(yùn)，祝面上順利！已經(jīng)有68人回復(fù)
華東師范大學(xué)芯片設(shè)計(jì)徐瓏真實(shí)水平咋樣已經(jīng)有1人回復(fù)
調(diào)劑已經(jīng)有0人回復(fù)

» 本主題相關(guān)價(jià)值貼推薦，對您同樣有幫助:

遺傳算法程序已經(jīng)有0人回復(fù)
福州大學(xué)2012年獲國家自然科學(xué)基金資助項(xiàng)目清單（截止8月21日）已經(jīng)有5人回復(fù)
2012年深圳大學(xué)資助項(xiàng)目清單已經(jīng)有7人回復(fù)
matlab在科技上的應(yīng)用已經(jīng)有17人回復(fù)
數(shù)學(xué)建模資料課件已經(jīng)有164人回復(fù)
2011杰青名單已經(jīng)有11人回復(fù)
2011年度國家杰出青年科學(xué)基金建議資助項(xiàng)目申請人名單已經(jīng)有0人回復(fù)
2011年度杰青名單公布已經(jīng)有7人回復(fù)
2011杰青名單出爐已經(jīng)有26人回復(fù)
2011年杰青名單及立項(xiàng) 已經(jīng)有1人回復(fù)
關(guān)于公布2011年度國家杰出青年科學(xué)基金建議資助項(xiàng)目申請人名單的通告已經(jīng)有7人回復(fù)
2011年度國家杰出青年科學(xué)基金名單已公布，所以8月19日公布面上青年地區(qū)基金已經(jīng)有50人回復(fù)
關(guān)于公布2011年度國家杰出青年科學(xué)基金建議資助項(xiàng)目申請人名單的通告已經(jīng)有21人回復(fù)
關(guān)于公布2011年度國家杰出青年科學(xué)基金建議資助項(xiàng)目申請人已經(jīng)有65人回復(fù)
【資源】《化學(xué)化工中的數(shù)學(xué)方法及MATLAB實(shí)現(xiàn)》.PDF 已經(jīng)有79人回復(fù)

第一顆紐扣扣錯(cuò)了。。。。

1樓 2015-12-08 20:52:11

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

For_study

金蟲 (小有名氣)

木蟲

應(yīng)助: 0 (幼兒園)
金幣: 1089.8
散金: 506
帖子: 174
在線: 108.8小時(shí)
蟲號: 2085809
注冊: 2012-10-25
性別: MM
專業(yè): 環(huán)境化工

沒人嗎。。。。。。。。。。？

回復(fù)此樓

第一顆紐扣扣錯(cuò)了。。。。

2樓2015-12-09 14:14:41

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

For_study

金蟲 (小有名氣)

木蟲

應(yīng)助: 0 (幼兒園)
金幣: 1089.8
散金: 506
帖子: 174
在線: 108.8小時(shí)
蟲號: 2085809
注冊: 2012-10-25
性別: MM
專業(yè): 環(huán)境化工

自己再頂一下。。。。。

回復(fù)此樓

第一顆紐扣扣錯(cuò)了。。。。

3樓2015-12-09 14:15:00

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

dingd

鐵桿木蟲 (職業(yè)作家)

計(jì)算強(qiáng)帖: 4
應(yīng)助: 1641 (講師)
金幣: 15037.3
散金: 101
紅花: 234
帖子: 3410
在線: 1223.7小時(shí)
蟲號: 291104
注冊: 2006-10-28

【答案】應(yīng)助回帖

感謝參與，應(yīng)助指數(shù) +1

換1stOpt試試，微分方程擬合簡單好用的多。

贊一下

回復(fù)此樓

4樓2015-12-09 16:06:39

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

yufeishu

銀蟲 (小有名氣)

應(yīng)助: 3 (幼兒園)
金幣: 467.5
帖子: 149
在線: 39.2小時(shí)
蟲號: 3222635
注冊: 2014-05-21
性別: GG
專業(yè): 熱物性與熱物理測試技術(shù)

然而并看不懂

回復(fù)此樓

5樓2015-12-10 10:00:17

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

ybkooo

至尊木蟲 (著名寫手)

應(yīng)助: 32 (小學(xué)生)
金幣: 11360
散金: 178
紅花: 15
帖子: 2310
在線: 543.5小時(shí)
蟲號: 2710519
注冊: 2013-10-09
性別: GG
專業(yè): 控制理論與方法

好多笑臉.

回復(fù)此樓

abcd

6樓2015-12-10 21:36:38

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

For_study

金蟲 (小有名氣)

木蟲

應(yīng)助: 0 (幼兒園)
金幣: 1089.8
散金: 506
帖子: 174
在線: 108.8小時(shí)
蟲號: 2085809
注冊: 2012-10-25
性別: MM
專業(yè): 環(huán)境化工

引用回帖:

6樓: Originally posted by ybkooo at 2015-12-10 21:36:38
好多笑臉.

笑臉是因?yàn)樾δ樀拇a和我打出來的一樣，所以在網(wǎng)頁上顯示為笑臉，其實(shí)笑臉的地方是：）。。。。。。。

贊一下

回復(fù)此樓

第一顆紐扣扣錯(cuò)了。。。。

7樓2015-12-11 15:04:45

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

For_study

金蟲 (小有名氣)

木蟲

應(yīng)助: 0 (幼兒園)
金幣: 1089.8
散金: 506
帖子: 174
在線: 108.8小時(shí)
蟲號: 2085809
注冊: 2012-10-25
性別: MM
專業(yè): 環(huán)境化工

引用回帖:

4樓: Originally posted by dingd at 2015-12-09 16:06:39
換1stOpt試試，微分方程擬合簡單好用的多。

參數(shù)太多，估計(jì)只有付費(fèi)版的才能解決問題。。。。

贊一下

回復(fù)此樓

第一顆紐扣扣錯(cuò)了。。。。

8樓2015-12-11 15:05:32

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

For_study

金蟲 (小有名氣)

木蟲

應(yīng)助: 0 (幼兒園)
金幣: 1089.8
散金: 506
帖子: 174
在線: 108.8小時(shí)
蟲號: 2085809
注冊: 2012-10-25
性別: MM
專業(yè): 環(huán)境化工

自己再給自己加把油。。。。

贊一下

回復(fù)此樓

第一顆紐扣扣錯(cuò)了。。。。

9樓2015-12-11 15:06:01

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

相關(guān)版塊跳轉(zhuǎn) 我要訂閱樓主 For_study 的主題更新

返回列表

最具人氣熱帖推薦 [查看全部]		作者	回/看	最后發(fā)表

[考研] 初試 317 +7	半拉月丙 2026-03-20	7/350	2026-03-21 22:26 by peike
[考研] 考研化學(xué)學(xué)碩調(diào)劑，一志愿985 +5	張vvvv 2026-03-15	7/350	2026-03-21 19:23 by ColorlessPI
[考研] 求調(diào)劑 +5	十三加油 2026-03-21	5/250	2026-03-21 18:48 by 學(xué)員8dgXkO
[考研] 0703化學(xué)調(diào)劑 +4	妮妮ninicgb 2026-03-21	4/200	2026-03-21 18:39 by 學(xué)員8dgXkO
[考研] 材料學(xué)碩333求調(diào)劑 +3	北道巷 2026-03-18	3/150	2026-03-21 18:17 by 學(xué)員8dgXkO
[考研] 311求調(diào)劑 +3	勇敢的小吳 2026-03-20	3/150	2026-03-21 17:40 by ColorlessPI
[考研] 265求調(diào)劑 +9	梁梁校校 2026-03-17	9/450	2026-03-21 02:17 by JourneyLucky
[考研] 一志愿南京理工大學(xué)085701資源與環(huán)境302分求調(diào)劑 +4	葵梓衛(wèi)隊(duì) 2026-03-18	6/300	2026-03-20 23:02 by JourneyLucky
[考研] 一志愿西安交通大學(xué) 學(xué)碩 354求調(diào)劑211或者雙一流 +3	我想要讀研究生 2026-03-20	3/150	2026-03-20 20:13 by JourneyLucky
[考研] 353求調(diào)劑 +3	拉鉤不許變 2026-03-20	3/150	2026-03-20 19:56 by JourneyLucky
[考研] 求調(diào)劑 +3	@taotao 2026-03-20	3/150	2026-03-20 19:35 by JourneyLucky
[考研] 招收調(diào)劑碩士 +4	lidianxing 2026-03-19	12/600	2026-03-20 12:25 by lidianxing
[考研] 085601材料工程專碩求調(diào)劑 +10	慕寒mio 2026-03-16	10/500	2026-03-19 15:26 by 丁丁*
[考研] 286求調(diào)劑 +6	lemonzzn 2026-03-16	10/500	2026-03-19 14:31 by lemonzzn
[考研] 一志愿985，本科211，0817化學(xué)工程與技術(shù)319求調(diào)劑 +10	Liwangman 2026-03-15	10/500	2026-03-19 10:25 by 無際的草原
[考研] 344求調(diào)劑 +6	knight344 2026-03-16	7/350	2026-03-18 20:13 by walc
[考研] 085601專碩，總分342求調(diào)劑，地區(qū)不限 +5	share_joy 2026-03-16	5/250	2026-03-18 14:48 by haxia
[考研] 283求調(diào)劑 +3	聽風(fēng)就是雨； 2026-03-16	3/150	2026-03-17 07:41 by 熱情沙漠
[考研] 321求調(diào)劑 +5	大米飯！ 2026-03-15	5/250	2026-03-16 16:33 by houyaoxu
[考研] 085601材料工程315分求調(diào)劑 +3	yang_0104 2026-03-15	3/150	2026-03-15 10:58 by peike

24小時(shí)熱門版塊排行榜

[求助] 遺傳算法優(yōu)化微分方程中的參數(shù) 已有1人參與

» 猜你喜歡

» 本主題相關(guān)價(jià)值貼推薦，對您同樣有幫助:

【答案】應(yīng)助回帖

» 本主題相關(guān)價(jià)值貼推薦，對您同樣有幫助: