24小时热门版块排行榜

返回列表

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

l1003785517

新虫 (小有名气)

应助: 0 (幼儿园)
金币: 11.4
散金: 30
帖子: 59
在线: 22.1小时
虫号: 4326948
注册: 2016-01-01
专业: 能源化工

[求助] 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);

回复此楼

» 猜你喜欢

【复旦大学】二维材料方向招收2026年博士研究生1名已经有0人回复
北京纳米能源与系统研究所王中林院士/曹南颖研究员课题组2026级硕/博/博后招生已经有10人回复
物理学I论文润色/翻译怎么收费? 已经有283人回复
荷兰Utrecht University超快太赫兹光谱王海教授课题招收2026 CSC博士生已经有23人回复
反铁磁体中的磁性切换：两种不同的机制已成功可视化已经有0人回复
求标准粉末衍射卡号 ICDD 01-076-1802 已经有0人回复
新西兰Robinson研究所招收全奖PhD 已经有0人回复
石墨烯转移--二氧化硅衬底石墨烯已经有0人回复

1楼 2018-07-09 17:05:38

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

Mxx99

新虫 (初入文坛)

应助: 0 (幼儿园)
金币: 23
帖子: 4
在线: 1.1小时
虫号: 9108945
注册: 2018-06-18
性别: MM
专业: 化学反应工程

引用回帖:

1楼: Originally posted by l1003785517 at 2018-07-09 17:05:38
各位虫友好！
我在用多变量线性回归方法估计动力学参数的时候写了如下代码（根据书上例题进行修改的），但是程序运行不了，说是目标函数返回了未定义参数。
有哪位帮忙看一下是哪里写错了吗？我自己一直没 ...

你在哪本书里看的例题，可以发一下，我也想学习这个

发自小木虫IOS客户端