| 查看: 1079 | 回复: 6 | ||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | ||
heshuangrong新虫 (小有名气)
|
[求助]
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) |
回复此楼» 猜你喜欢
重庆交大26年硕士生招生拟调剂通知已出！欢迎加入机器视觉与3D光学成像课题组。
已经有0人回复
-
**
已经有1人回复
-
物理学I论文润色/翻译怎么收费?
已经有126人回复
-
0702一志愿吉大B区求调剂 本科期间发表一篇Sci
已经有3人回复
-
基底STO,薄膜SRO，XRD里面的振荡，是laue震荡还是kiessig振荡? 怎么判断？
已经有2人回复
-
B区学生调剂-兰州交通大学材料科学与工程学院
已经有8人回复
-
山西大同大学物理学专业还有调剂名额，欢迎调剂！
已经有12人回复
-
桂林理工大学物理学专业招收调剂
已经有11人回复
-
VASP 的一组 GPU / CPU 基准测试记录
已经有0人回复
-
美国顶级物理期刊《应用物理快报》（APL）的编辑欺骗和歧视作者及AIP的官僚主义傲慢
已经有1人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 向大神求教一个用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人回复
heshuangrong
新虫 (小有名气)
- 应助: 2 (幼儿园)
- 金币: 784.2
- 散金: 25
- 帖子: 96
- 在线: 93.5小时
- 虫号: 1498593
- 注册: 2011-11-18
- 专业: 化学反应工程
送红花一朵 
3楼2015-09-01 15:49:15
lvgao116
铁杆木虫 (著名写手)
- 应助: 4 (幼儿园)
- 金币: 20624.4
- 红花: 13
- 沙发: 1
- 帖子: 1691
- 在线: 153小时
- 虫号: 907551
- 注册: 2009-11-19
- 性别: GG
- 专业: 应用地球物理学
【答案】应助回帖
感谢参与，应助指数 +1
| 不要输入带有function语句的那一行。或者用%注释。我已经试过没问题。 |
» 本帖已获得的红花（最新10朵）
2楼2015-09-01 15:21:47
lvgao116
铁杆木虫 (著名写手)
- 应助: 4 (幼儿园)
- 金币: 20624.4
- 红花: 13
- 沙发: 1
- 帖子: 1691
- 在线: 153小时
- 虫号: 907551
- 注册: 2009-11-19
- 性别: GG
- 专业: 应用地球物理学
4楼2015-09-01 17:01:02
lvgao116
铁杆木虫 (著名写手)
- 应助: 4 (幼儿园)
- 金币: 20624.4
- 红花: 13
- 沙发: 1
- 帖子: 1691
- 在线: 153小时
- 虫号: 907551
- 注册: 2009-11-19
- 性别: GG
- 专业: 应用地球物理学
【答案】应助回帖
|
sp = form: 'B-' knots: [0 0 0 50 120 300 300 300] coefs: [10.0579 6.9123 3.5802 1.3462 0.9975] number: 5 order: 3 dim: 1 pp = form: 'B-' knots: [0 0 78.8178 182.1889 300 300] coefs: [-0.1187 -0.0312 -0.0104 -0.0061] number: 4 order: 2 dim: 1 dCAdt = -0.1187 -0.0965 -0.0743 -0.0521 -0.0229 -0.0109 -0.0061 CAi = Columns 1 through 11 10.0307 9.8531 9.6780 9.5054 9.3353 9.1678 9.0027 8.8402 8.6803 8.5228 8.3679 Columns 12 through 22 8.2155 8.0656 7.9182 7.7734 7.6310 7.4913 7.3540 7.2192 7.0870 6.9573 6.8301 Columns 23 through 33 6.7054 6.5833 6.4637 6.3466 6.2320 6.1200 6.0105 5.9035 5.7990 5.6970 5.5976 Columns 34 through 44 5.5007 5.4063 5.3144 5.2251 5.1382 5.0539 4.9722 4.8929 4.8162 4.7420 4.6703 Columns 45 through 55 4.6011 4.5345 4.4703 4.4088 4.3497 4.2931 4.2391 4.1876 4.1386 4.0916 4.0452 Columns 56 through 66 3.9992 3.9536 3.9085 3.8639 3.8197 3.7760 3.7328 3.6900 3.6476 3.6057 3.5643 Columns 67 through 77 3.5233 3.4828 3.4427 3.4031 3.3640 3.3253 3.2870 3.2492 3.2119 3.1750 3.1386 Columns 78 through 88 3.1027 3.0672 3.0321 2.9975 2.9634 2.9297 2.8965 2.8638 2.8315 2.7996 2.7682 Columns 89 through 99 2.7373 2.7068 2.6768 2.6472 2.6181 2.5895 2.5613 2.5335 2.5062 2.4794 2.4530 Columns 100 through 110 2.4271 2.4017 2.3767 2.3521 2.3280 2.3044 2.2812 2.2585 2.2363 2.2144 2.1931 Columns 111 through 121 2.1722 2.1518 2.1318 2.1123 2.0932 2.0746 2.0564 2.0387 2.0215 2.0047 1.9884 Columns 122 through 132 1.9725 1.9568 1.9413 1.9258 1.9104 1.8950 1.8798 1.8646 1.8496 1.8346 1.8197 Columns 133 through 143 1.8048 1.7901 1.7754 1.7608 1.7463 1.7319 1.7176 1.7033 1.6892 1.6751 1.6611 Columns 144 through 154 1.6472 1.6333 1.6196 1.6059 1.5923 1.5788 1.5654 1.5521 1.5388 1.5256 1.5126 Columns 155 through 165 1.4996 1.4866 1.4738 1.4610 1.4484 1.4358 1.4233 1.4108 1.3985 1.3862 1.3740 Columns 166 through 176 1.3619 1.3499 1.3380 1.3262 1.3144 1.3027 1.2911 1.2796 1.2682 1.2568 1.2456 Columns 177 through 187 1.2344 1.2233 1.2122 1.2013 1.1905 1.1797 1.1690 1.1584 1.1479 1.1374 1.1271 Columns 188 through 198 1.1168 1.1066 1.0965 1.0865 1.0765 1.0667 1.0569 1.0472 1.0376 1.0281 1.0186 Columns 199 through 200 1.0093 1.0000 k = 0.0053 n = 1.3854 |
5楼2015-09-01 17:01:44
20