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

登录注册

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

24小时热门版块排行榜

北京石油化工学院2026年研究生招生接收调剂公告

返回列表

lf0612

木虫 (小有名气)

应助: 0 (幼儿园)
金币: 1912.8
散金: 100
帖子: 289
在线: 120.1小时
虫号: 964589
注册: 2010-03-08
专业: 计算机软件

[求助] 呼叫版主，在线紧急求助，关于matlab中微分方程组参数拟合得问题！已有2人参与

参考版主“月只蓝”以前帖子中的代码，写了如下关于微分方程组中参数拟合得代码，进行了修改，结果有问题，因为是新手，代码的含义不全懂，加之时间急着要，所以特来求助，往各位高手多多帮助，在此先谢谢了！

function k1k2k3
format long
clear all
clc
tspan = [0 1 2 3 4 5 6 7 8 9 10 ];
x0=[1 1 1 1 ];
k0 = [1 1 1 1 1 1 1 1 1 1 1 ]*1e-4;
lb = [0 0 0 0 0 0 0 0 0 0 0];
ub = [1 1 1 1 1 1 1 1 1 1 1 ]*1e5;

data=[0.012,0.073,0.089,0.092,0.316,0.498,1.810,2.314,0.679,0.063,0.044;
0.005,0.084,0.093,0.270,0.406,0.642,0.814,1.026,0.422,0.053,0.014;
0.037,0.638,0.969,2.057,20.419,39.655,69.107,34.099,5.657,0.697,0.346;
0.008,0.224,0.234,0.347,1.860,2.854,3.786,1.883,0.601,0.157,0.079;
]';
size(data);
yexp = data;

[k,resnorm,residual,exitflag,output,lambda,jacobian] =lsqnonlin(@ObjFunc,k0,lb,ub,[],tspan,x0,yexp);
ci = nlparci(k,residual,jacobian);
fprintf('\n\n使用函数lsqnonlin()估计得到的参数值为:\n')
fprintf('\tk1 = %.16f \n',k(1))
fprintf('\tk2 = %.16f \n',k(2))
fprintf('\tk3 = %.16f \n',k(3))
fprintf('\tk4 = %.16f \n',k(4))
fprintf('\tk5 = %.16f \n',k(5))
fprintf('\tk5 = %.16f \n',k(6))
fprintf('\tk5 = %.16f \n',k(7))
fprintf('\tk5 = %.16f \n',k(8))
fprintf('\tk5 = %.16f \n',k(9))

ts=0:1:max(tspan);

[ts ys]=ode45(@KineticsEqs,ts,x0,[],k);
[ttt XXsim] = ode45(@KineticsEqs,tspan,x0,[],k);

figure(1),
plot(ts,ys(:,1),'b',tspan(2:end),yexp(:,1),'or'),legend('计算值','实验值','Location','best');

figure(2),
plot(ts,ys(:,2),'b',tspan(2:end),yexp(:,2),'or'),legend('计算值','实验值','Location','best');

figure(3),
plot(ts,ys(:,3),'b',tspan(2:end),yexp(:,3),'or'),legend('计算值','实验值','Location','best');

figure(4),
plot(ts,ys(:,4),'b',tspan(2:end),yexp(:,4),'or'),legend('计算值','实验值','Location','best');

function f = ObjFunc(k,tspan,x0,yexp)          % 目标函数
[t Xsim] = ode45(@KineticsEqs,tspan,x0,[],k);
Xsim1=Xsim(:,1);
Xsim2=Xsim(:,2);
Xsim3=Xsim(:,3);
Xsim4=Xsim(:,4);
ysim(:,1) = Xsim1(2:end);
ysim(:,2) = Xsim2(2:end);
ysim(:,3) = Xsim3(2:end);
ysim(:,4) = Xsim4(2:end);
size(ysim(:,1));
size(ysim(:,2));
size(ysim(:,3));
size(ysim(:,4));
size(yexp(:,1));
size(yexp(:,2));
size(yexp(:,3));
size(yexp(:,4));
f = [(ysim(:,1)-yexp(:,1)) (ysim(:,2)-yexp(:,2)) (ysim(:,3)-yexp(:,3))  (ysim(:,4)-yexp(:,4)) ];

function dEdt = KineticsEqs(t,C,k)
% ODE模型方程 %g(1)=k(4);g(2)=k(5);g(3)=k(6);m(1)=k(7);m(2)=k(8);m(3)=k(9);
CA=C(1);CB=C(2);CC=C(3);CD=C(4);
%dCAdt = k(1)*CA*(CB/2-1)-k(2)*CB+k(3)*CC;
%dCBdt=-g(1)*CA+g(2)*CB.*(1-CB/5)+g(3)*CC.*(1-CC/10);
%dCCdt=m(1)*CA.*(CA/2-1)-m(2)*CB-m(3)*CC;
dCAdt= -k(1)*CA-k(2)*CA;
dCBdt=k(1)*CA-k(3)*CB-k(4)*CB;
dCCdt=k(3)*CB-k(5)*CC;
dCDdt=k(2)*CA+k(4)*CB+k(5)*CC;
dEdt=[dCAdt;dCBdt;dCCdt;dCDdt];

在此感谢！

回复此楼

» 猜你喜欢

山西大同大学物理学硕士研究生招收调剂生已经有0人回复
[调剂信息]211智能人工感知方向国家青年特聘专家课题组招收调剂研究生已经有0人回复
物理学I论文润色/翻译怎么收费? 已经有50人回复
［调剂信息］211智能人工感知方向国家青年特聘专家课题组招收调剂研究生已经有0人回复
[调剂信息]211智能人工感知方向国家青年特聘专家课题组招收调剂研究生已经有0人回复
0702一志愿吉大B区求调剂本科期间发表一篇Sci 已经有2人回复
070200求调剂，一志愿某211，288分已经有18人回复
法国博士后职位已经有0人回复
重庆交大26年硕士生招生拟调剂通知已出！欢迎加入机器视觉与3D光学成像课题组。已经有0人回复

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

matlab里拟合后，知道y坐标求x坐标已经有4人回复
Matlab中ode45出现大误差，关于用ODE45解偏微分方程的疑问，有待解决，大家共同讨论！已经有4人回复
]紧急求助复介电常数的模型拟合-Matlab or 1stOpt 已经有9人回复
求助使用matlab对一方程进行拟合，自己拟合总是出错！已经有12人回复
matlab 拟合方法已经有4人回复
matlab中参数改变后怎样求得曲线？？？已经有6人回复
求助，matlab/ 1stOpt 自定义微分方程拟合实验数据来求方程中的参数已经有9人回复
使用lsqnonlin函数优化动力学参数，总是得不到合理的结果已经有6人回复
求助，matlab/ 1stOpt 自定义微分方程拟合实验数据来求方程中的参数已经有4人回复
求助，matlab自定义微分方程拟合实验数据来求方程中的参数已经有4人回复
三维图谱中，用matlab如何实现去除一部分数据并实现曲面拟合已经有10人回复
紧急求助，利用Matlab对实验数据进行拟合求解参数。已经有27人回复
微分方程组参数拟合的问题（多参数动力系统参数估计）已经有11人回复
微分方程组参数拟合的问题已经有11人回复
关于matlab微分、及非线性拟合问题已经有5人回复
MATLAB微分方程参数拟合问题，求大神已经有7人回复
求高版本1stopt，解决一个问题，微分方程组参数拟合已经有3人回复
求助拟合参数的置信区间问题已经有6人回复
关于matlab线性规划背包问题的解法求救呀~~方程知道但不知道如何解答···· 已经有4人回复
matlab微分方程组参数拟合，以周为单位求解，汇总后以年为单位进行数值比较已经有4人回复
matlab-常微分方程参数估计已经有12人回复
lsqnonlin函数拟合微分方程组参数拟合问题已经有10人回复
【求助】微分方程拟合求参数已经有7人回复

1楼 2014-06-24 21:08:21

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

回帖支持 ( 显示支持度最高的前 50 名 )

月只蓝

主管区长 (职业作家)

专家经验: +1059
计算强帖: 8
应助: 1712 (讲师)
贵宾: 8.888
金币: 68128.7
散金: 1938
红花: 443
沙发: 4
帖子: 4373
在线: 3291.4小时
虫号: 1122189
注册: 2010-10-14
专业: 宇宙学
管辖: 计算模拟区

【答案】应助回帖

今晚比较忙，今天才着手做，不知道来得及不。程序已经调通

CODE:

function k1k2k3k4k5

format long

clear all

clc

tspan = [0 1 2 3 4 5 6 7 8 9 10];

x0=[0.012 0.005 0.037 0.008];

%k0 = [-0.5 -0.8 -0.8 -1.5 2];  

k0 = [-0.87 -0.9 48.11 -47.80 -0.142];

lb = -[1 1 1 1 1]*1e9;

ub = [1 1 1 1 1]*1e5;

data=[0.073,0.089,0.092,0.316,0.498,1.810,2.314,0.679,0.063,0.044;

    0.084,0.093,0.270,0.406,0.642,0.814,1.026,0.422,0.053,0.014;

    0.638,0.969,2.057,20.419,39.655,69.107,34.099,5.657,0.697,0.346;

   0.224,0.234,0.347,1.860,2.854,3.786,1.883,0.601,0.157,0.079;

   ]';

size(data);

yexp = data;

[k,resnorm,residual,exitflag,output,lambda,jacobian] =lsqnonlin(@ObjFunc,k0,lb,ub,[],tspan,x0,yexp);       

ci = nlparci(k,residual,jacobian);

fprintf('\n\n使用函数lsqnonlin()估计得到的参数值为:\n')

fprintf('\tk1 = %.16f \n',k(1))

fprintf('\tk2 = %.16f \n',k(2))

fprintf('\tk3 = %.16f \n',k(3))

fprintf('\tk4 = %.16f \n',k(4))

fprintf('\tk5 = %.16f \n',k(5))

ts=0:1:max(tspan);

[ts ys]=ode45(@KineticsEqs,ts,x0,[],k);

[ttt XXsim] = ode45(@KineticsEqs,tspan,x0,[],k); 

figure(1),

plot(ts,ys(:,1),'b',tspan(2:end),yexp(:,1),'or'),legend('计算值','实验值','Location','best');

figure(2),

plot(ts,ys(:,2),'b',tspan(2:end),yexp(:,2),'or'),legend('计算值','实验值','Location','best');

figure(3),

plot(ts,ys(:,3),'b',tspan(2:end),yexp(:,3),'or'),legend('计算值','实验值','Location','best');

figure(4),

plot(ts,ys(:,4),'b',tspan(2:end),yexp(:,4),'or'),legend('计算值','实验值','Location','best');

function f = ObjFunc(k,tspan,x0,yexp)           % 目标函数

[t Xsim] = ode45(@KineticsEqs,tspan,x0,[],k);

Xsim1=Xsim(:,1);

Xsim2=Xsim(:,2);

Xsim3=Xsim(:,3);

Xsim4=Xsim(:,4);

ysim(:,1) = Xsim1(2:end);

ysim(:,2) = Xsim2(2:end);

ysim(:,3) = Xsim3(2:end);

ysim(:,4) = Xsim4(2:end);

size(ysim(:,1));

size(ysim(:,2));

size(ysim(:,3));

size(ysim(:,4));

size(yexp(:,1));

size(yexp(:,2));

size(yexp(:,3));

size(yexp(:,4));

f = [(ysim(:,1)-yexp(:,1)) (ysim(:,2)-yexp(:,2)) (ysim(:,3)-yexp(:,3))  (ysim(:,4)-yexp(:,4)) ];

function dEdt = KineticsEqs(t,C,k)    

% ODE模型方程 %g(1)=k(4);g(2)=k(5);g(3)=k(6);m(1)=k(7);m(2)=k(8);m(3)=k(9);

CA=C(1);CB=C(2);CC=C(3);CD=C(4);

%dCAdt = k(1)*CA*(CB/2-1)-k(2)*CB+k(3)*CC;

%dCBdt=-g(1)*CA+g(2)*CB.*(1-CB/5)+g(3)*CC.*(1-CC/10);

%dCCdt=m(1)*CA.*(CA/2-1)-m(2)*CB-m(3)*CC;

dCAdt= -k(1)*CA-k(2)*CA;

dCBdt=k(1)*CA-k(3)*CB-k(4)*CB;

dCCdt=k(3)*CB-k(5)*CC;

dCDdt=k(2)*CA+k(4)*CB+k(5)*CC;

dEdt=[dCAdt;dCBdt;dCCdt;dCDdt];

参数比较多，你再调调初值k0吧。
多参数的拟合，还是用1stopt软件为好，在此推荐专家dingd。

赞一下(1人)

回复此楼

MATLAB、MS小问题、普通问题请发帖求助！时间精力有限，恕不接受无偿私信求助。

5楼2014-06-25 09:55:08

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

普通回帖

lf0612

木虫 (小有名气)

应助: 0 (幼儿园)
金币: 1912.8
散金: 100
帖子: 289
在线: 120.1小时
虫号: 964589
注册: 2010-03-08
专业: 计算机软件

原始的微分方程组如下：
dS/dt=-k1*S-k2*S
dE/dt=k1*S-k3*E-k4*E
dI/dt=k3*E-k5*I
dR/dt=k2*S+-k4*E+k5*I

赞一下

回复此楼

2楼2014-06-24 21:17:46

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

月只蓝

主管区长 (职业作家)

专家经验: +1059
计算强帖: 8
应助: 1712 (讲师)
贵宾: 8.888
金币: 68128.7
散金: 1938
红花: 443
沙发: 4
帖子: 4373
在线: 3291.4小时
虫号: 1122189
注册: 2010-10-14
专业: 宇宙学
管辖: 计算模拟区

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与，应助指数 +1
lf0612: 金币+10, ★有帮助, 5 2014-07-04 16:29:05

如果你一定要用MATLAB的话，我可以帮你把程序调通。
不过对于常微分方程拟合，建议你找一下dingd专家，他有高版本的1stopt软件，计算高效强大。也许MATLAB编了半天程序，拟合结果没有1stopt花十分钟计算的结果好。

赞一下

回复此楼

MATLAB、MS小问题、普通问题请发帖求助！时间精力有限，恕不接受无偿私信求助。

3楼2014-06-24 21:38:38

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

lf0612

木虫 (小有名气)

应助: 0 (幼儿园)
金币: 1912.8
散金: 100
帖子: 289
在线: 120.1小时
虫号: 964589
注册: 2010-03-08
专业: 计算机软件

引用回帖:

3楼: Originally posted by 月只蓝 at 2014-06-24 21:38:38
如果你一定要用MATLAB的话，我可以帮你把程序调通。
不过对于常微分方程拟合，建议你找一下dingd专家，他有高版本的1stopt软件，计算高效强大。也许MATLAB编了半天程序，拟合结果没有1stopt花十分钟计算的结果好。

谢谢你了，帮忙把我把程序调通吧，我对参数的准确度要求不是那么的高，不要差的太多就好了！

赞一下

回复此楼

4楼2014-06-24 21:47:56

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

dingd

铁杆木虫 (职业作家)

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

【答案】应助回帖

感谢参与，应助指数 +1

公式数据对不对啊？效果很差的。

赞一下

回复此楼

6楼2014-06-25 15:49:35

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

lf0612

木虫 (小有名气)

应助: 0 (幼儿园)
金币: 1912.8
散金: 100
帖子: 289
在线: 120.1小时
虫号: 964589
注册: 2010-03-08
专业: 计算机软件

引用回帖:

6楼: Originally posted by dingd at 2014-06-25 15:49:35
公式数据对不对啊？效果很差的。

谢谢了，数据有点问题，不过已经来不及了，所以最后手工试了几个参数，找了个最好的带进去了！再次感谢！

赞一下

回复此楼

7楼2014-07-04 16:31:45

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

lf0612

木虫 (小有名气)

应助: 0 (幼儿园)
金币: 1912.8
散金: 100
帖子: 289
在线: 120.1小时
虫号: 964589
注册: 2010-03-08
专业: 计算机软件

引用回帖:

5楼: Originally posted by 月只蓝 at 2014-06-25 09:55:08
今晚比较忙，今天才着手做，不知道来得及不。程序已经调通

function k1k2k3k4k5
format long
clear all
clc
tspan = ;
x0=;
%k0 = ;
k0 = ;
lb = -*1e9;
ub = *1e5;

data=';
size(data);
yexp ...

谢谢了，数据有点问题，不过已经来不及了，所以最后手工试了几个参数，找了个最好的带进去了！再次感谢！

赞一下

回复此楼

8楼2014-07-04 16:31:57

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

lf0612

木虫 (小有名气)

应助: 0 (幼儿园)
金币: 1912.8
散金: 100
帖子: 289
在线: 120.1小时
虫号: 964589
注册: 2010-03-08
专业: 计算机软件

上网不方便，所以现在才回复二位，不好意思，谢谢了！

赞一下

回复此楼

9楼2014-07-04 16:32:28

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

小牧的影子

铁虫 (小有名气)

应助: 1 (幼儿园)
金币: 462.4
散金: 25
帖子: 181
在线: 38.1小时
虫号: 2650509
注册: 2013-09-12
性别: GG
专业: 光学

不错不错，学习了

回复此楼

学习，学习是立足之本。

10楼2015-06-04 21:47:45

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

相关版块跳转我要订阅楼主 lf0612 的主题更新

返回列表

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

[考研] 一志愿武汉理工大学080200机械工程308分，求调剂 +4	终不似从前 2026-04-05	4/200	2026-04-06 11:46 by 考研学校招点人
[考研] 319求调剂 +3	handrui 2026-04-05	3/150	2026-04-06 09:33 by jp9609
[考研] 0703求调剂383分 +9	W55j 2026-04-03	9/450	2026-04-06 06:50 by houyaoxu
[考研] 求调剂求调剂 +8	121. 2026-04-02	8/400	2026-04-05 20:15 by lys0704
[考研] 调剂 +3	好好读书。 2026-04-02	3/150	2026-04-05 13:02 by arrow8852
[考研] 271分求调剂学校 +12	zph158488！ 2026-04-02	13/650	2026-04-05 10:13 by lqwchd
[考研] 311分 22408 求调剂 +3	bing_bot 2026-04-03	3/150	2026-04-05 00:43 by chongya
[考研] 288环境专硕,求调材料方向 +13	lllllos 2026-04-04	14/700	2026-04-04 23:34 by lqwchd
[考研] 材料383求调剂 +5	郭阳阳阳成 2026-04-04	5/250	2026-04-04 19:06 by dongzh2009
[考研] 363求调剂 +7	zh096 2026-04-04	7/350	2026-04-04 17:13 by dongzh2009
[考研] 本科211，专业085404，293分请求调剂 +5	莲菜就是藕吧 2026-04-04	5/250	2026-04-04 14:08 by 这是一个无聊的�
[考研] 总分328生物与医药考数学求调剂 +7	aaadim 2026-04-02	9/450	2026-04-03 22:53 by syh9288
[考研] 求调剂 +8	akdhjs 2026-04-03	8/400	2026-04-03 18:17 by 戴维ING
[考研] 372分材料与化工（085600）一志愿湖南大学求调剂 +3	蓝笺片 2026-04-03	4/200	2026-04-03 17:58 by Jimmyandyou
[考研] 一志愿华东理工大学，080500学硕，317分，求调剂 +13	s1145 2026-03-31	15/750	2026-04-03 11:44 by msi123
[考研] 302求调剂 +9	zyx上岸！ 2026-04-02	9/450	2026-04-02 23:07 by 马儿快快地跑
[考研] 求调剂 +7	Aniyaio 2026-04-02	7/350	2026-04-02 16:42 by zzsw+
[考研] 材料化工340求调剂 +5	jhx777 2026-03-30	5/250	2026-04-02 12:45 by smileboy2006
[考研] 085600 一志愿9 总分351 求调剂学校 +7	czhcz 2026-03-31	9/450	2026-04-01 19:24 by 无际的草原
[考研] 考研调剂求助 +7	13287130938 2026-03-31	7/350	2026-03-31 16:39 by 690616278

24小时热门版块排行榜

[求助] 呼叫版主，在线紧急求助，关于matlab中微分方程组参数拟合得问题！ 已有2人参与

» 猜你喜欢

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

【答案】应助回帖

【答案】应助回帖

【答案】应助回帖

[求助] 呼叫版主，在线紧急求助，关于matlab中微分方程组参数拟合得问题！已有2人参与