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

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

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

返回列表

當(dāng)前只顯示滿足指定條件的回帖，點(diǎn)擊這里查看本話題的所有回帖

For_study

金蟲(chóng) (小有名氣)

木蟲(chóng)

應(yīng)助: 0 (幼兒園)
金幣: 1089.8
散金: 506
帖子: 174
在線: 108.8小時(shí)
蟲(chóng)號(hào): 2085809
注冊(cè): 2012-10-25
性別: MM
專(zhuān)業(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)

% 動(dòng)力學(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)

% 動(dòng)力學(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

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

回復(fù)此樓

» 猜你喜歡

哈爾濱理工大學(xué)物理系招收物理學(xué)考研調(diào)劑已經(jīng)有7人回復(fù)
0702一志愿吉大B區(qū)求調(diào)劑已經(jīng)有5人回復(fù)
物理學(xué)I論文潤(rùn)色/翻譯怎么收費(fèi)? 已經(jīng)有93人回復(fù)
求調(diào)劑已經(jīng)有0人回復(fù)
0702一志愿吉大B區(qū)求調(diào)劑有論文已經(jīng)有0人回復(fù)
請(qǐng)問(wèn)還有沒(méi)有用Latex寫(xiě)文章的小伙伴們？已經(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à)值貼推薦，對(duì)您同樣有幫助:

遺傳算法程序已經(jīng)有0人回復(fù)
福州大學(xué)2012年獲國(guó)家自然科學(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年度國(guó)家杰出青年科學(xué)基金建議資助項(xiàng)目申請(qǐ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年度國(guó)家杰出青年科學(xué)基金建議資助項(xiàng)目申請(qǐng)人名單的通告已經(jīng)有7人回復(fù)
2011年度國(guó)家杰出青年科學(xué)基金名單已公布，所以8月19日公布面上青年地區(qū)基金已經(jīng)有50人回復(fù)
關(guān)于公布2011年度國(guó)家杰出青年科學(xué)基金建議資助項(xiàng)目申請(qǐng)人名單的通告已經(jīng)有21人回復(fù)
關(guān)于公布2011年度國(guó)家杰出青年科學(xué)基金建議資助項(xiàng)目申請(qǐ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的回帖

dingd

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

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

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

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

換1stOpt試試，微分方程擬合簡(jiǎn)單好用的多。

贊一下

回復(fù)此樓

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

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

查看全部 9 個(gè)回答

For_study

金蟲(chóng) (小有名氣)

木蟲(chóng)

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

沒(méi)人嗎。。。。。。。。。。？

回復(fù)此樓

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

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

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

For_study

金蟲(chóng) (小有名氣)

木蟲(chóng)

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

自己再頂一下。。。。。

回復(fù)此樓

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

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

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

For_study

金蟲(chóng) (小有名氣)

木蟲(chóng)

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

引用回帖:

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

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

贊一下

回復(fù)此樓

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

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

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

查看全部 9 個(gè)回答

普通表情龍兔虎貓

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

[考研] 本人考085602 化學(xué)工程專(zhuān)碩 +21	不知道叫什么！ 2026-03-15	23/1150	2026-03-22 00:12 by BruceLiu320
[考研] 一志愿北京化工大學(xué)070300 學(xué)碩336求調(diào)劑 +3	vv迷 2026-03-21	4/200	2026-03-21 23:05 by f19980501
[考研] 326求調(diào)劑 +5	諾貝爾化學(xué)獎(jiǎng)覬?/a> 2026-03-15	8/400	2026-03-21 19:33 by ColorlessPI
[基金申請(qǐng)] 學(xué)校已經(jīng)提交到NSFC，還能修改嗎？ 40+4	babangida 2026-03-19	9/450	2026-03-21 16:12 by babangida
[考研] 材料與化工（0856）304求 B區(qū) 調(diào)劑 +3	邱gl 2026-03-21	3/150	2026-03-21 13:47 by lature00
[考研] 279分求調(diào)劑一志愿211 +14	chaojifeixia 2026-03-19	15/750	2026-03-21 13:24 by zhukairuo
[考研] 一志愿中國(guó)石油大學(xué)（華東）本科齊魯工業(yè)大學(xué) +3	石能偉 2026-03-17	3/150	2026-03-21 02:22 by JourneyLucky
[考研] 085700資源與環(huán)境308求調(diào)劑 +12	墨墨漠 2026-03-18	13/650	2026-03-21 01:42 by JourneyLucky
[考研] 271材料工程求調(diào)劑 +8	.6lL 2026-03-18	8/400	2026-03-21 00:58 by JourneyLucky
[考研] 323求調(diào)劑 +3	洼小桶 2026-03-18	3/150	2026-03-20 22:54 by JourneyLucky
[考研] 329求調(diào)劑 +9	想上學(xué)吖吖 2026-03-19	9/450	2026-03-20 22:01 by luoyongfeng
[考研] A區(qū)線材料學(xué)調(diào)劑 +5	周周無(wú)極 2026-03-20	5/250	2026-03-20 21:33 by laoshidan
[考研] 一志愿西南交通專(zhuān)碩材料355 本科雙非求調(diào)劑 +5	西南交通專(zhuān)材355 2026-03-19	5/250	2026-03-20 21:10 by JourneyLucky
[考研] 085410人工智能專(zhuān)碩317求調(diào)劑（0854都可以） +4	xbxudjdn 2026-03-18	4/200	2026-03-20 09:07 by 不168
[考研] 求調(diào)劑 +3	暗涌afhb 2026-03-16	3/150	2026-03-20 00:28 by 河南大學(xué)校友
[考研] 293求調(diào)劑 +11	zjl的號(hào) 2026-03-16	16/800	2026-03-18 08:10 by zhukairuo
[考研] 326求調(diào)劑 +5	上岸的小葡 2026-03-15	6/300	2026-03-17 17:26 by ruiyingmiao
[考研] 材料工程專(zhuān)碩274一志愿211求調(diào)劑 +6	薛云鵬 2026-03-15	6/300	2026-03-17 11:05 by 學(xué)員h26Tkc
[考研] 275求調(diào)劑 +4	太陽(yáng)花天天開(kāi)心 2026-03-16	4/200	2026-03-17 10:53 by 功夫瘋狂
[考研] 070303 總分349求調(diào)劑 +3	LJY9966 2026-03-15	5/250	2026-03-16 14:24 by xwxstudy