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

登录注册

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

24小时热门版块排行榜

返回列表

2011207156

金虫 (小有名气)

应助: 0 (幼儿园)
金币: 1335.7
红花: 2
帖子: 194
在线: 153.8小时
虫号: 1829583
注册: 2012-05-22
专业: 化学反应工程

[求助] matlab求模型参数有错误，求指导已有2人参与

以下的程序是我照葫芦画瓢改的，不知道哪里错了。。。。。。本人刚开始接触matlab，求大神指导啊！

function KineticsEst1_Int % 用积分法分析得到kh，b，和kd
clear all; clc
global nm
t=[31 32 33 34 35 36 37 38 39 40 41 42 43 44 45];
nm=[20.41220 22.73032 24.55624 25.88679 26.69019 27.16334 27.45866 27.64919 27.57615 27.70635 27.77303 27.86195 27,91911 27.95721 27.98897]';

% 非线性拟合
beta0=[0.0001 2 0.0001]; tspan=[31 32 33 34 35 36 37 38 39 40 41 42 43 44 45]; n0=20.4122;
[beta,resnorm,resid,exitflag,output,lambda,jacobian] = ...
lsqnonlin(@OptObjFunc,beta0,[],[],[],tspan,n0,nm)
ci=nlparci(beta,resid,jacobian)

% 拟合效果图（实验与拟合的比较）
[t4plot n4plot]=ode45(@KineticsEqs,[tspan(1) tspan(end)],n0,[],beta);
plot(t,nm,'bo',t4plot,n4plot,'k-')
legend('Exp','Model'), xlabel('时间t，s'), ylabel('累积量n,mol')

%残差关于拟合值的残差图
[t nc]=ode45(@KineticsEqs,tspan,nO,[],beta);
figure, plot(nc,resid,'*'), xlabel('累积值（mol)'), ylabel('残差R(mol)'), refline(0,0)

%参数辨识结果
fprintf('\n\nEstimated Parameters:\n'), fprintf('\tkh=%.4f±%.4f\n',beta(1),ci(1,2)-beta(1))
fprintf('\tb = %.2f ± %.2f\n',beta(2),ci(2,2)-beta(2))
fprintf('\tkd = %.2f ± %.2f\n',beta(3),ci(3,2)-beta(3))
% ------------------------------------------------------------------
function f = OptObjFunc(beta,tspan,n0,nm)
[t nc] = ode45(@KineticsEqs,tspan,n0,[],beta);
f = nc - nm;
% ------------------------------------------------------------------
function dndt = KineticsEqs(t,n,beta)
dndt = 0.001414/(1/(4.5727*10^4*beta(1)*(1-0.02119*n)^beta(2))+1/beta(3)/45727); % kh= beta(1), b= beta(2),kd=beta(3)

报错：
??? Error using ==> minus
Matrix dimensions must agree.

Error in ==> KineticsEst1_Int4>OptObjFunc at 29
f = nc - nm;

Error in ==> lsqnonlin at 203
initVals.F = feval(funfcn{3},xCurrent,varargin{:});

Error in ==> KineticsEst1_Int4 at 9
[beta,resnorm,resid,exitflag,output,lambda,jacobian] = ...

Caused by:
Failure in initial user-supplied objective function
evaluation. LSQNONLIN cannot continue.

回复此楼

» 猜你喜欢

比利时鲁文大学固体物理系招收CSC博士已经有0人回复
🌟 比利时新鲁汶大学（UCLouvain）诚邀CSC博士加入Pascal Gehring教授团队已经有0人回复
物理学I论文润色/翻译怎么收费? 已经有179人回复
第一性原理计算方向2026级博士申请 PRB*1，四级484 已经有1人回复
求助VISSIM破解版软件已经有0人回复
求2026年在台湾举行的物理和材料领域国际学术会议信息已经有0人回复
求国际会议网站已经有1人回复
求取一些关于纳米材料和纳米技术相关的英文PPT。已经有0人回复
【复旦大学】二维材料方向招收2026年博士研究生1名已经有0人回复
北京纳米能源与系统研究所王中林院士/曹南颖研究员课题组2026级硕/博/博后招生已经有10人回复

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

空间状态模型在MATLAB里的表示？已经有10人回复
matlab动力学模型编程求助已经有6人回复
]紧急求助复介电常数的模型拟合-Matlab or 1stOpt 已经有9人回复
关于求极值的问题已经有4人回复
matlab非线性拟合模型参数，大神你在哪！已经有73人回复
毕不了业了，发帖求助Origin8.0多元非线性拟合的问题已经有9人回复
求解这类方程用什么软件编程较好···· 已经有8人回复
求一个MATLAB的汽车模型已经有4人回复
仿真算例中的被控对象是系统的数学模型吗已经有11人回复
重金请教-如何用MAtlab对力学模型进行参数识别已经有12人回复
MATLAB关于lsqcurvefit程序几个错误，求大神指导一下已经有6人回复
matlab求解非线性方程组，错误提示怎么解决已经有5人回复
matlab拟合模型参数已经有5人回复
二元幂函数参数求解已经有10人回复
重金求助Matlab图像加噪声问题已经有8人回复
Matlab中RBF建模，仿真是遇到的问题已经有6人回复
【求助】求MATLAB DLCA模型程序已经有11人回复

1楼2014-06-30 11:03:08

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

spiderone

木虫 (小有名气)

应助: 2 (幼儿园)
金币: 3472.1
散金: 12
红花: 1
帖子: 157
在线: 108.4小时
虫号: 2520847
注册: 2013-06-25
性别: GG
专业: 化学反应工程

【答案】应助回帖

感谢参与，应助指数 +1

矩阵维数不对，nc和nm的个数应该不对，可以仔细看看。

赞一下(1人)

回复此楼

2楼2014-06-30 16:47:09

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

月只蓝

主管区长 (职业作家)

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

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与，应助指数 +1
2011207156: 金币+20, ★★★★★最佳答案, 非常谢谢哦~ 2014-07-01 16:01:45

CODE:

function KineticsEst1_Int_modified_by_Yuezhilan

clear all;clc

format long

tspan=[31 32 33 34 35 36 37 38 39 40 41 42 43 44 45]-31; 

yexp=[22.73032 24.55624 25.88679 26.69019 27.16334 27.45866 27.64919 27.57615 27.70635 27.77303 27.86195 27.91911 27.95721 27.98897]';

beta0=[100 10 0.1 ];

y0=20.4122;

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

ub=[1e9 1e9 1e2];   

yy=[y0 yexp'];

k0=beta0;

[k,resnorm,residual,exitflag,output,lambda,jacobian] = ...

    lsqnonlin(@ObjFunc,k0,lb,ub,[],tspan,y0,yexp);

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

fprintf('\t待拟合参数 k1 = %.6f\n',k(1))

fprintf('\t待拟合参数 k2 = %.6f\n',k(2))

fprintf('\t待拟合参数 k3 = %.6f\n',k(3))

fprintf('  残差平方和= %.6f\n\n',resnorm)

ts=0:1:max(tspan);

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

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

y=XXsim(2:end);

xexp=yexp;

R2=1-sum((xexp-y).^2)./sum((xexp-mean(y)).^2);

fprintf('\n\t决定系数R-Square = %.6f',R2);

figure(1)

plot(ts,ys,'b',tspan,yy,'or'),legend('计算值','实验值','Location','best');

yr=y-yexp;

figure(2)

plot(tspan(2:end),yr,'r*',[-1 15],[0 0]),axis([-1 15 -0.5 0.5]);

figure(3)

plot(yexp,y,'ro',[21 29],[21 29],'b-');

%---------------------------------------------------------

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

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

ysim = Xsim(2:end);

size(ysim);

size(yexp);

f=ysim-yexp;

 %----------------------------------------------------------

function dydt = KineticsEqs(t,y,k)

beta(1)=k(1);

beta(2)=k(2);

beta(3)=k(3);

dydt = 0.001414/(1/(4.5727*10^4*beta(1)*(1-0.02119*y)^beta(2))+1/beta(3)/45727);

计算结果显示：
使用函数lsqnonlin()估计得到的参数值为:
待拟合参数 k1 = 73143442.357330
待拟合参数 k2 = 27.497659
待拟合参数 k3 = 0.035271
残差平方和= 0.241151

决定系数R-Square = 0.992227

附图1.png

赞一下(3人)

回复此楼

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

3楼2014-06-30 20:05:14

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

2011207156

金虫 (小有名气)

应助: 0 (幼儿园)
金币: 1335.7
红花: 2
帖子: 194
在线: 153.8小时
虫号: 1829583
注册: 2012-05-22
专业: 化学反应工程

引用回帖:

3楼: Originally posted by 月只蓝 at 2014-06-30 20:05:14
function KineticsEst1_Int_modified_by_Yuezhilan
clear all;clc
format long
tspan=-31;
yexp=';

beta0=;
y0=20.4122;
lb=*1e6;
ub=;
yy=;
k0=beta0;
= ...
lsqnonlin(@ObjFunc,k0,lb,u ...

真心感谢大神~

回复此楼

4楼2014-07-01 16:02:35

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