版块导航: 正在加载中...

登录注册

应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜

>论坛更新日志 (6726)
>考研 (1905)
>导师招生 (1553)
>文献求助 (512)
>考博 (312)
>虫友互识 (300)
>硕博家园 (164)
>博后之家 (126)
>基金申请 (112)
>招聘信息布告栏 (108)
>论文投稿 (99)
>休闲灌水 (96)
>第一性原理 (79)
>教师之家 (68)
>公派出国 (50)
>外文书籍求助 (32)

返回列表

当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

For_study

金虫 (小有名气)

木虫

应助: 0 (幼儿园)
金币: 1089.8
散金: 506
帖子: 174
在线: 108.8小时
虫号: 2085809
注册: 2012-10-25
性别: MM
专业: 环境化工

[求助] 遗传算法优化微分方程中的参数已有1人参与

%拟合微分方程参数
%直接拟合反应的活化能和指前因子
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);

%输出参数
fprintf('\n\n遗传算法的估计数值:\n');
disp(h0);

%构造适应度函数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)

% 动力学微分方程
% 共分为8集总，原料（饱和分SS为1，芳香分SA为2，胶质沥青质SR为3），柴油（D）为4,汽油（G）为5，液化气（LPG）为6，干气（Gas）为7，焦炭（C）为8，
% 由于遗传算法只能返回向量，所以反应常数不能是矩阵，只能是向量，第一个数表示反应物，第二个数表示生成物
%反应常数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常数，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重芳烃失活系数，N碱氮吸附失活系数，deact催化剂结焦失活系数，dens密度，Swh真实重时空速
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)

% 动力学微分方程
% 共分为8集总，原料（饱和分SS为1，芳香分SA为2，胶质沥青质SR为3），柴油（D）为4,汽油（G）为5，液化气（LPG）为6，干气（Gas）为7，焦炭（C）为8，
% 由于遗传算法只能返回向量，所以反应常数不能是矩阵，只能是向量，第一个数表示反应物，第二个数表示生成物
%反应常数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常数，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重芳烃失活系数，N碱氮吸附失活系数，deact催化剂结焦失活系数，dens密度，Swh真实重时空速
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

问题叙述，遗传算法未设定lb和ub时，程序可以正常运行，获得结果，当设定lb和ub后，程序依然可以运行（不会报错），但是只能运算一次优化后的结果（设置断点时发现）。遗传算法好像没有继续优化下去，希望有大神可以帮忙解答一下问题。。。。谢谢。。。。

回复此楼

» 猜你喜欢

光学工程学硕调剂信息已经有32人回复
欢迎加入课题组已经有0人回复
物理学I论文润色/翻译怎么收费? 已经有212人回复
散金币，求好运，祝面上顺利！已经有79人回复
华东师范大学芯片设计徐珑真实水平咋样已经有1人回复
调剂已经有0人回复
【新加坡】纳米电子器件项目组有“联合培养博士生”名额已经有0人回复
中德博士后2026项目5月30日截止- 国内博士后入站、德国两年全职已经有0人回复
面向全球申请人的中德博士后2026项目5月30日截止- 国内博后入站德国两年全职已经有0人回复
【新加坡】纳米电子器件项目组有“联合培养博士生”名额已经有2人回复
重庆交大材料学院26年硕士招生复试实施细则已出。欢迎加入光子学微结构与器件课题组。已经有0人回复

» 本主题相关价值贴推荐，对您同样有帮助:

遗传算法程序已经有0人回复
福州大学2012年获国家自然科学基金资助项目清单（截止8月21日）已经有5人回复
2012年深圳大学资助项目清单已经有7人回复
matlab在科技上的应用已经有17人回复
数学建模资料课件已经有164人回复
2011杰青名单已经有11人回复
2011年度国家杰出青年科学基金建议资助项目申请人名单已经有0人回复
2011年度杰青名单公布已经有7人回复
2011杰青名单出炉已经有26人回复
2011年杰青名单及立项已经有1人回复
关于公布2011年度国家杰出青年科学基金建议资助项目申请人名单的通告已经有7人回复
2011年度国家杰出青年科学基金名单已公布，所以8月19日公布面上青年地区基金已经有50人回复
关于公布2011年度国家杰出青年科学基金建议资助项目申请人名单的通告已经有21人回复
关于公布2011年度国家杰出青年科学基金建议资助项目申请人已经有65人回复
【资源】《化学化工中的数学方法及MATLAB实现》.PDF 已经有79人回复

第一颗纽扣扣错了。。。。

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

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

For_study

金虫 (小有名气)

木虫

应助: 0 (幼儿园)
金币: 1089.8
散金: 506
帖子: 174
在线: 108.8小时
虫号: 2085809
注册: 2012-10-25
性别: MM
专业: 环境化工

引用回帖:

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

笑脸是因为笑脸的代码和我打出来的一样，所以在网页上显示为笑脸，其实笑脸的地方是：）。。。。。。。

赞一下

回复此楼

第一颗纽扣扣错了。。。。

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

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

查看全部 9 个回答

For_study

金虫 (小有名气)

木虫

应助: 0 (幼儿园)
金币: 1089.8
散金: 506
帖子: 174
在线: 108.8小时
虫号: 2085809
注册: 2012-10-25
性别: MM
专业: 环境化工

没人吗。。。。。。。。。。？

回复此楼

第一颗纽扣扣错了。。。。

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

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

For_study

金虫 (小有名气)

木虫

应助: 0 (幼儿园)
金币: 1089.8
散金: 506
帖子: 174
在线: 108.8小时
虫号: 2085809
注册: 2012-10-25
性别: MM
专业: 环境化工

自己再顶一下。。。。。

回复此楼

第一颗纽扣扣错了。。。。

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

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

dingd

铁杆木虫 (职业作家)

计算强帖: 4
应助: 1641 (讲师)
金币: 15037.3
散金: 101
红花: 234
帖子: 3410
在线: 1223.7小时
虫号: 291104
注册: 2006-10-28

【答案】应助回帖

感谢参与，应助指数 +1

换1stOpt试试，微分方程拟合简单好用的多。

赞一下

回复此楼

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

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

查看全部 9 个回答

最具人气热帖推荐 [查看全部]		作者	回/看	最后发表

[考研] 085601求调剂总分293英一数二 +3	钢铁大炮 2026-03-24	3/150	2026-03-24 22:03 by bingxueer79
[考研] 085600材料与化工调剂 +8	A-哆啦Z梦 2026-03-23	13/650	2026-03-24 21:05 by greychen00
[考研] 0854 考研调剂招生了！AI 方向 +5	pk3725069 2026-03-19	17/850	2026-03-24 17:30 by zhouxuan..
[考研] 求调剂一志愿武汉理工大学材料工程（085601） +5	WW.' 2026-03-23	7/350	2026-03-24 14:50 by sprinining
[考研] 321求调剂 +4	Ymlll 2026-03-24	4/200	2026-03-24 14:44 by sprinining
[考研] 293求调剂 +6	加一一九 2026-03-24	6/300	2026-03-24 14:29 by JourneyLucky
[考研] 一志愿北京化工大学材料与化工 264分各科过A区国家线 +3	哈哈157349 2026-03-21	3/150	2026-03-24 14:11 by zhyzzh
[基金申请] 请教下大家 2026年国家基金申请是双盲审吗？ +3	lishucheng1 2026-03-22	5/250	2026-03-24 08:22 by gltch
[考研] 276求调剂。有半年电池和半年高分子实习经历 +9	材料学257求调剂 2026-03-23	10/500	2026-03-24 07:36 by wangy0907
[考研] 求调剂 +7	十三加油 2026-03-21	7/350	2026-03-23 23:48 by 热情沙漠
[考研] 一志愿南京理工大学085701资源与环境302分求调剂 +5	葵梓卫队 2026-03-18	7/350	2026-03-23 16:26 by lingjue
[考研] 350求调剂 +6	weudhdk 2026-03-19	6/300	2026-03-23 15:47 by tangyuan0840221
[考研] 315分，诚求调剂，材料与化工085600 +3	13756423260 2026-03-22	3/150	2026-03-22 20:11 by edmund7
[考研] 一志愿西北大学，070300化学学硕，总分287，双非一本，求调剂。 +3	晨昏线与星海 2026-03-20	3/150	2026-03-22 16:00 by ColorlessPI
[考研] 299求调剂 +5	shxchem 2026-03-20	7/350	2026-03-21 17:09 by ColorlessPI
[考研] 材料与化工（0856）304求 B区调剂 +3	邱gl 2026-03-21	3/150	2026-03-21 13:47 by lature00
[考研] 材料与化工 322求调剂 +4	然11 2026-03-19	4/200	2026-03-20 22:12 by luoyongfeng
[考研] 一志愿西南交通专硕材料355 本科双非求调剂 +5	西南交通专材355 2026-03-19	5/250	2026-03-20 21:10 by JourneyLucky
[考研] 一志愿吉林大学材料学硕321求调剂 +11	Ymlll 2026-03-18	15/750	2026-03-20 19:40 by 丁丁*
[考研] 261求B区调剂，科研经历丰富 +3	牛奶很忙 2026-03-20	4/200	2026-03-20 19:34 by JourneyLucky

24小时热门版块排行榜

For_study

[求助] 遗传算法优化微分方程中的参数 已有1人参与

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:

For_study

For_study

For_study

dingd

【答案】应助回帖

[求助] 遗传算法优化微分方程中的参数已有1人参与