| 查看: 2363 | 回复: 6 | |||
[求助]
matlab拟合拟均相动力学参数
|
|
模型如附件所示，要求里面的k+,k-,m,n,p，我自己编了个程序，但是一直有问题 M文件 function f=fortyfour_D(c,t,rD) t=[5,10,15,20,30,45,60,80,100,120,150,180,210,240,300,360]; CD=[749.16 ,994.49 ,1382.00,1589.18,1909.81,2222.86,2525.68,2721.90,2896.67,2951.42,2960.67,2999.60,3017.97,3048.91,3072.27,3063.85]; knots = 3;K = 3; %三次B样条 sp=spap2(knots,K,t,CD); sp=spap2(newknt(sp),K,t,CD); pp=fnder(sp); % 计算B样条函数的导函数 dCDdt=fnval(pp,t); % 计算t处的导函数值 Wcat=85.6312 %催化剂浓度 rD=dCDdt./Wcat; f=c(1)*CA.^c(2)*CB.^c(3)-c(4)*CD.^c(5)))-rD; 命令如下 c0=[1,1,1,1,1]; for i=1:100 c=lsqnonlin('fortyfour_D ',c0); c0=c; end c % 绘制图形 ti = linspace(t(1),t(end),200); %y=linspace（a，b，n）生成一个行向量，该向量将a与b之间平分为n个点，包含端点a和b。 CDi = fnval(sp,ti); %计算样条函数区间内任意一点的值 plot(t,CD,'ro',ti,CDi,'b-')，xlabel('t'),ylabel('C_D') legend('实验数据','拟合曲线') c=[1,1,1,1,1];for i=1:100;c=lsqnonlin('fortyfour_D ',c);c;end CD=[749.16 ,994.49 ,1382.00,1589.18,1909.81,2222.86,2525.68,2721.90,2896.67,2951.42,2960.67,2999.60,3017.97,3048.91,3072.27,3063.85]; CA = [2732.43 2468.85 2223.97 2023.64 1696.41 1343.42 1054.34 902.20 756.91 717.86 703.53 631.00 652.41 649.53 624.92 744.95]; CB= [3127.04 2722.48 2414.51 2267.55 1927.70 1500.17 1285.07 1111.91 924.27 886.32 830.86 821.82 801.97 780.39 794.41 783.73]; rD=[0.5631,0.5282,0.4934,0.4585,0.3888,0.2842,0.1796,0.0482,0.0404,0.0326,0.0209,0.0092,0.0046,0.0040,0.0028,0.0015]; >> c0=[1,1,1,1,1]; for i=1:1000 c=lsqcurvefit('ffD ',c0,CD,rD); c0=c; end c 请大家帮忙看看我程序哪里有问题[ Last edited by lxyy on 2011-12-30 at 16:05 ] |
回复此楼» 本帖附件资源列表
-
欢迎监督和反馈：小木虫仅提供交流平台，不对该内容负责。
本内容由用户自主发布，如果其内容涉及到知识产权问题，其责任在于用户本人，如对版权有异议，请联系邮箱：xiaomuchong@tal.com - 附件 1 : clip_image002.gif
2011-12-30 16:05:36, 503 bytes
» 猜你喜欢
【新加坡】纳米电子器件项目组有“联合培养博士生”名额
已经有0人回复
-
中德博士后2026项目5月30日截止- 国内博士后入站、德国两年全职
已经有0人回复
-
物理学I论文润色/翻译怎么收费?
已经有197人回复
-
面向全球申请人的中德博士后2026项目5月30日截止- 国内博后入站德国两年全职
已经有0人回复
-
【新加坡】纳米电子器件项目组有“联合培养博士生”名额
已经有2人回复
-
重庆交大材料学院26年硕士招生复试实施细则已出。欢迎加入光子学微结构与器件课题组。
已经有0人回复
-
北工大宋晓艳教授团队招收2026年9月入学博士【学博和专博名额均有】！
已经有1人回复
-
天津城建大学理学院化学、环境工程、材料物理等方向接收10人
已经有0人回复
-
天津城建大学理学院化学、环境工程、材料物理等方向接收调剂，10人
已经有0人回复
-
计算与实验科学工程国际会议(8月7-11日,香港) 二维材料分会S5与声学换能器分会S8
已经有0人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- matlab拟合动力学参数
已经有4人回复
- 非线性拟合估计方程参数（matlab，spss，origin都可以）
已经有4人回复
- matlab非线性参数拟合问题
已经有7人回复
- 请教matlab拟合
已经有9人回复
- 求助：化学反应动力学数据拟合得到Arrenius方程参数（A，belta，Ea）
已经有4人回复
- 动力学数据拟合
已经有5人回复
- 动力学方程拟合 参数求解
已经有4人回复
- 【求助】用matlab最优化方法进行参数拟合
已经有17人回复
- 求助matlab拟合参数的问题
已经有2人回复
- 【求助】使用MATLAB拟合出按照公式的系数
已经有19人回复
- 【求助】使用Matlab拟合反应动力学方程问题
已经有7人回复
dingd
铁杆木虫 (职业作家)
- 计算强帖: 4
- 应助: 1641 (讲师)
- 金币: 15037.3
- 散金: 101
- 红花: 234
- 帖子: 3410
- 在线: 1223.7小时
- 虫号: 291104
- 注册: 2006-10-28
★
dbb627(金币+1): 欢迎交流 2011-12-31 08:41:48
dbb627(金币+1): 欢迎交流 2011-12-31 08:41:48
|
不知理解的对否，用1stOpt试下： 均方差(RMSE): 0.0187338627500751 残差平方和(SSE): 0.00561532181661841 相关系数(R): 0.996169068786095 相关系数之平方(R^2): 0.992352813606156 决定系数(DC): 0.99235277918993 参数 最佳估算 ---------- ------------- k1 14.3465921423275 m 0.0595653946350665 n -0.452413851540345 k2 1.79201130280021E-8 p 2.22793339858939  |
2楼2011-12-30 22:02:01
3楼2011-12-30 22:39:49
dingd
铁杆木虫 (职业作家)
- 计算强帖: 4
- 应助: 1641 (讲师)
- 金币: 15037.3
- 散金: 101
- 红花: 234
- 帖子: 3410
- 在线: 1223.7小时
- 虫号: 291104
- 注册: 2006-10-28
4楼2011-12-31 08:18:50
|
M文件 function f=q(c,t,CD,rD) t=[0 0.08333 0.166667 0.25 0.5 0.75 1 1.5 2 2.666667 3.333333 4 5]; CD=[0 0.53259 1.154323 1.621033 2.644088 3.163488 3.461446 3.710468 3.784974 3.746928 3.761759 3.713404 3.643825]; rD=[0.0710291014584411 0.0700187830917774 0.0655437729479848 0.0576582625377526 0.0329826401194250 0.0166083513736405 0.00912074378901541 0.00278932296049053 0.000202497370868135 0 0 0 0]； f= c(1)*((4.6026-CD).^c(2)*(5.064335-CD).^c(3)-c(4)*CD.^c(5))-rD; 命令如下 c0=[0.1,0.1,0.1,0.1,0.1]; for i=1:50 c=lsqnonlin('q ',c0); c0=c; end c 我把程序简化了下，还是一直运行不了 |
5楼2011-12-31 11:56:24
bao2007pin
铜虫 (正式写手)
- 应助: 2 (幼儿园)
- 金币: 2047.6
- 散金: 207
- 帖子: 307
- 在线: 170.5小时
- 虫号: 1119550
- 注册: 2010-10-11
- 性别: GG
- 专业: 有机合成
6楼2011-12-31 12:43:36
7楼2011-12-31 14:00:47
