| 查看: 1340 | 回复: 4 | ||
l1003785517新虫 (小有名气)
|
[求助]
Matlab 用多变量线性回归方法估计动力学参数 已有1人参与
|
|
各位虫友好！ 我在用多变量线性回归方法估计动力学参数的时候写了如下代码（根据书上例题进行修改的），但是程序运行不了，说是目标函数返回了未定义参数。 有哪位帮忙看一下是哪里写错了吗？我自己一直没检查出来哪里错了，谢谢大家！ function KineticsEst % 动力学参数估计 clear all clc C_DC=[6.31317 6.29412 6.27507 6.26237 6.24713 6.23316 6.221095 6.210935]; C_TT=[0.029302 0.04459 0.059878 0.069433 0.080899 0.091091 0.099372 0.107016]; C_NP=[0.007644 0.011466 0.015288 0.018473 0.022295 0.026117 0.029939 0.032487]; r1=[0.0574 0.0407 0.0241 0.0278 0.0315 0.0268 0.022 0.0172]; % 用多变量线性回归方法估计动力学参数 R = C_DC/r1; y = R; X = [ones(size(y)) C_DC C_TT C_NP]; b=X\y; % 或b = X\y k1= 1/b(1); k2 = b(2)*k1; k3 = b(3)*k1; k4 = b(4)*k1; % 用lsqnonlin()--求解非线性最小二乘法(非线性数据拟合)问题 beta0 = [k1 k2 k3 k4]; lb = [0 0 0 0]; ub = [+inf +inf +inf +inf]; [beta,resnorm,resid,exitflag,output,lambda,jacobian] = ... lsqnonlin(@ObjFun,beta0,lb,ub,[], C_DC,C_TT,C_NP,r1); ci = nlparci(beta,resid,jacobian); % 模型适定性判别 Ne = length(r1); Np = length(beta); [rho2, F] = rho2_F(k1, r1, resnorm, Ne, Np); % 残差关于拟合值的残差图 rc = RateEqs(beta,C_DC,C_TT,C_NP); plot(rc,resid,'*') xlabel('反应速率拟合值, mol/(min g催化剂)') ylabel('残差R, mol/(min g催化剂)') refline(0,0) % 参数辨识结果 fprintf('\n\nEstimated Parameters:\n') fprintf('\tk1 = %.4f ± %.4f\n',beta(1),ci(1,2)-beta(1)) fprintf('\tk3 = %.4f ± %.4f\n',beta(2),ci(2,2)-beta(2)) fprintf('\tk4 = %.4f ± %.4f\n',beta(3),ci(3,2)-beta(3)) fprintf('\tk5 = %.4f ± %.4f\n',beta(4),ci(4,2)-beta(4)) fprintf(' 决定性指标ρ^2: %.3f\n',rho2) fprintf(' F比: %.3f\n\n',F) % ------------------------------------------------------------------ function f = ObjFun(beta,C_DC,C_TT,C_NP,r1) rc =RateEqs(beta,C_DC,C_TT,C_NP); f = r1 - rc; % ------------------------------------------------------------------ function rc =RateEqs(beta,C_DC,C_TT,C_NP) % Rate equation rc = beta(1)*C_DC./(1+beta(2)*C_DC+beta(3)*C_TT+beta(4)*C_NP); |
回复此楼» 猜你喜欢
瑞典林雪平大学博士后招聘｜PEC Water Splitting 方向
已经有130人回复
-
PRB投稿，Acknowledgment sent to author状态十天了
已经有3人回复
-
物理学I论文润色/翻译怎么收费?
已经有78人回复
-
PVA溶解
已经有2人回复
-
基金申请
已经有44人回复
-
CSC与新西兰维多利亚大学PhD奖学金项目
已经有0人回复
-
新西兰Robinson研究所 招聘CSC公派访问人员
已经有0人回复
-
帮我的英语口语老师找学生
已经有0人回复
-
什么时候开奖？
已经有7人回复
-
散金币祈福
已经有94人回复
hzlhm
至尊木虫 (著名写手)
- 应助: 387 (硕士)
- 金币: 17854.5
- 红花: 53
- 帖子: 2880
- 在线: 606.9小时
- 虫号: 1517335
- 注册: 2011-11-30
- 性别: GG
- 专业: 常微分方程与动力系统
2楼2018-07-11 19:52:37
3楼2018-07-13 06:46:24
l1003785517
新虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 11.4
- 散金: 30
- 帖子: 59
- 在线: 22.1小时
- 虫号: 4326948
- 注册: 2016-01-01
- 专业: 能源化工
4楼2018-07-15 14:41:43
l1003785517
新虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 11.4
- 散金: 30
- 帖子: 59
- 在线: 22.1小时
- 虫号: 4326948
- 注册: 2016-01-01
- 专业: 能源化工
5楼2018-07-15 14:42:05
40