24小时热门版块排行榜

返回列表

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

heshuangrong

新虫 (小有名气)

应助: 2 (幼儿园)
金币: 784.2
散金: 25
帖子: 96
在线: 93.5小时
虫号: 1498593
注册: 2011-11-18
专业: 化学反应工程

[求助] MATLAB拟合动力学问题已有1人参与

程序是书上的，店Debug--Run没法运行，求教各位大侠

function KineticDataFit
clear all
clc

% 动力学数据
t = [0 20  40  60  120  180  300];
CA = [10  8  6  5  3  2  1];

% 用最小二乘样条拟合法计算微分dCA/dt--使用不经过实验点的B样条插值函数
knots = 3;
K = 3;    % 三次B样条
sp = spap2(knots,K,t,CA)
sp = spap2(newknt(sp),K,t,CA);
pp = fnder(sp)       % 计算B样条函数的导函数
dCAdt = fnval(pp,t)    % 计算t处的导函数值

% 绘制图形
ti = linspace(t(1),t(end),200);
CAi = fnval(sp,ti)
plot(t,CA,'ro',ti,CAi,'b-')
xlabel('t')
ylabel('C_A')
figure
fnplt(pp)
% dCAdti = fnval(pp,ti)
% plot(ti,dCAdti,'-')
xlabel('t')
ylabel('dC/dt')

% 线性拟合
rA = dCAdt;
y = log(-rA);
x = log(CA);
p = polyfit(x,y,1);
k = exp(p(2))
n = p(1)

回复此楼

» 猜你喜欢

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

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

向大神求教一个用matlab拟合动力学参数的问题，求指导已经有5人回复
给我大神求教一个用matlab拟合动力学参数的问题，求指导已经有0人回复
Matlab拟合动力学参数问题，求指导！已经有9人回复
matlab 拟合多参数动力学方程求程序已经有4人回复
MATLAB拟合动力学问题已经有11人回复
matlab拟合动力学参数，非常谢谢~ 已经有0人回复
求助matlab拟合动力学K值已经有4人回复
matlab拟合吸附动力学已经有0人回复
使用MATLAB拟合动力学参数问题已经有18人回复
拜求！matlab拟合动力学参数已经有14人回复
求助用matlab拟合动力学方程已经有9人回复
求助用matlab拟合动力学方程已经有13人回复
matlab 拟合反应动力学参数已经有8人回复
matlab拟合酶促反应机制已经有7人回复
求大家帮帮忙看看matlab拟合动力学方程的问题，已经走投无路了，求大家帮帮忙啊！已经有19人回复
matlab 拟合反应动力学参数结果很差。大家帮忙看一下已经有14人回复
用matlab 拟合动力学级数已经有2人回复
matlab拟合拟均相动力学参数已经有6人回复
matlab拟合动力学参数已经有4人回复
【求助】使用Matlab拟合反应动力学方程问题已经有7人回复

1楼 2015-09-01 14:23:46

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

heshuangrong

新虫 (小有名气)

应助: 2 (幼儿园)
金币: 784.2
散金: 25
帖子: 96
在线: 93.5小时
虫号: 1498593
注册: 2011-11-18
专业: 化学反应工程

引用回帖:

4楼: Originally posted by lvgao116 at 2015-09-01 17:01:02
spap2是一个matlab自带函数，你运行后matlab给你的提示是没有定义该函数或变量，说明你的matlab版本中没有这个函数。ps：我的matlab版本是R2008a，7.6.0....

我的是6.5,去下个高点的版本看看

发自小木虫Android客户端