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

登录注册

应《网络安全法》要求，自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];

在此感谢！

回复此楼

» 猜你喜欢

欢迎调剂到赣南师范大学智能制造与未来能源学院已经有27人回复
重庆交大26年硕士生招生拟调剂通知已出！欢迎加入光子学微结构与器件课题组。已经有0人回复
物理学I论文润色/翻译怎么收费? 已经有189人回复
UJN物理学专业调剂已经有12人回复
0702一志愿吉大B区求调剂本科期间发表一篇Sci 已经有0人回复
北京纳米能源与系统研究所王中林院士/曹南颖研究员课题组招收2026硕士研究生1名已经有0人回复
山西大同大学物理学硕士研究生招收调剂生已经有0人回复
[调剂信息]211智能人工感知方向国家青年特聘专家课题组招收调剂研究生已经有0人回复
［调剂信息］211智能人工感知方向国家青年特聘专家课题组招收调剂研究生已经有0人回复
[调剂信息]211智能人工感知方向国家青年特聘专家课题组招收调剂研究生已经有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的回帖

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的回帖

查看全部 13 个回答

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的回帖

月只蓝

主管区长 (职业作家)

专家经验: +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的回帖

查看全部 13 个回答

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

[考研] 一志愿a区211，085601-307分求调剂 +11	党嘉豪 2026-03-31	24/1200	2026-04-02 21:54 by yuq
[考研] 一志愿华中农业071010，总分320求调剂 +6	困困困困坤坤 2026-04-02	6/300	2026-04-02 21:28 by dongzh2009
[考研] 土木304求调剂 +4	兔突突突， 2026-04-02	5/250	2026-04-02 21:16 by 兔突突突，
[考研] 环境科学与工程334分求调剂 +7	王一一依依 2026-03-30	9/450	2026-04-02 21:15 by 1104338198
[考研] 一志愿哈尔滨工业大学材料与化工方向336分 +14	辰沐5211314 2026-03-26	14/700	2026-04-02 21:12 by 1104338198
[考研] 一志愿北京科技大学材料学硕328分求调剂 +6	1段时间 2026-03-31	7/350	2026-04-02 13:57 by 3041
[考研] 298求B区调剂 +4	zzz，，r 2026-04-02	5/250	2026-04-02 12:17 by 土木硕士招生
[考研] 329求调剂 +8	miaodesi 2026-04-02	10/500	2026-04-02 10:24 by 不吃魚的貓
[考研] 食品学硕362求调剂 +3	xuanxianxian 2026-04-01	3/150	2026-04-01 21:05 by 啊李999
[考研] 【求调剂】085601材料工程专硕 \| 总分272 \| +10	脚滑的守法公民 2026-03-27	10/500	2026-04-01 17:23 by pies112
[考研] 285求调剂 +5	FZAC123 2026-03-30	5/250	2026-04-01 15:50 by 韩雨涵
[考研] 材料与化工调剂一志愿大连海事085600，349 +9	吃的不少 2026-03-30	9/450	2026-04-01 11:24 by wangjy2002
[考研] 一志愿西交大080500材料学硕349 +6	jqx1258 2026-03-31	7/350	2026-03-31 21:08 by yuq
[考研] 复试调剂 +7	双马尾痞老板2 2026-03-31	7/350	2026-03-31 19:49 by Dyhoer
[考研] 求调剂 +8	11ggg 2026-03-30	8/400	2026-03-31 13:56 by nanaliuyun
[考研] 调剂求院校招收 +7	鹤鲸鸽 2026-03-28	7/350	2026-03-31 11:21 by oooqiao
[考研] 279求调剂 +12	j的立方 2026-03-29	12/600	2026-03-30 20:30 by dick_runner
[考研] 一志愿双一流机械285分求调剂 +4	幸运的三木 2026-03-29	5/250	2026-03-29 14:49 by Miko19
[考研] 材料与化工（0856）304求B区调剂 +8	邱gl 2026-03-27	8/400	2026-03-28 12:42 by 唐沐儿
[考研] 081200-314 +3	LILIQQ 2026-03-27	4/200	2026-03-28 09:41 by 保护地球你我做�

24小时热门版块排行榜

lf0612

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

» 猜你喜欢

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

lf0612

lf0612

月只蓝

【答案】应助回帖

月只蓝

【答案】应助回帖

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