| 查看: 523 | 回复: 1 | ||
[求助]
用matlab拟合积分方程- 急用！坐等！非常感谢！
|
方程如下， 自变量是T，因变量是Cp。Tt 和ΔH 是拟合常数。请用 β(T) 对T作图，并根据以上方程拟合得出Tt 和ΔH的值。（Tend是数据的最后一个T值。） 以上是对于曲线只有一个峰值的情况，对于两个峰的情况，方程有一点小改动，如下，  作图要求跟第一种情况一样。因为本人对于matlab知之甚少，哪怕是一点小改动对于我来说都是座高山，所以只能劳烦高手了。附件里有单峰和双峰的数据举例。第一列是T值，第二列是Cp值。 另外以下是我曾经找一个朋友帮忙做的程序，当时对方程的理解错误，把Cp完全积分出来作为Y值拟合的，结果在Tt转折点附近总是积分不好，程序和拟合结果如下，（相应的数据也在附件里面，data.txt)。希望能有参考价值。 clc clear all close all d = load('data.txt', '-ASCII'); XX = d(:,1); YY = d(:,2); %idx = (XX>280); idx = (XX>304); X = XX(idx) Y = YY(idx); R = 8.314; f = @(a,x)a(1)*(1-exp( a(2)/R * (1./x - 1/a(3))) ) ; % hyperbolic-like %f = @(a,x)a(1) * (1-(1-exp( a(2)/R * (1./x - 1/a(3))) ) ./ (1 + exp( a(2)/R * (1./x-1/a(3))) )) ; a0 = [1400 220000 304] [a1,r] = nlinfit(X, Y, f, a0); A = a1(1) B = a1(2) T = a1(3) mse = sum(r.^2)/(length(r)-length(a1)) plot(XX,YY,'k', XX, f(a1,XX), 'r') ylabel('$\Delta \hat{H}(T) \quad / \quad kCal \cdot mol^{-1} \cdot K^{-1}$','Interpreter','LaTex') xlabel('$T \quad / \quad K$','interpreter','latex') legend('data', 'fitted param') axis([XX(1), XX(end), YY(1), YY(end)])  |
回复此楼» 本帖附件资源列表
-
欢迎监督和反馈：小木虫仅提供交流平台，不对该内容负责。
本内容由用户自主发布，如果其内容涉及到知识产权问题，其责任在于用户本人，如对版权有异议，请联系邮箱：xiaomuchong@tal.com - 附件 1 : 双峰.txt
- 附件 2 : 单峰.txt
- 附件 3 : data.txt
2014-08-21 02:26:38, 12.11 K
2014-08-21 03:16:56, 6.1 K
2014-08-21 03:17:38, 9.51 K
» 猜你喜欢
散金币祈福
已经有186人回复
-
求助，如何提取ELK的rt-TDDFT在某一时刻的自旋密度分布和电子密度分布
已经有1人回复
-
物理学I论文润色/翻译怎么收费?
已经有251人回复
-
澳门大学IAPME李海峰教授常年招收博士研究生
已经有0人回复
-
澳门大学IAPME李海峰教授常年招收博士研究生
已经有0人回复
-
澳门大学IAPME李海峰教授常年招收博士研究生
已经有0人回复
-
深圳大学应用技术学院招聘凝聚态物理博士后
已经有112人回复
-
深圳大学2027级光电信息工程/物理学研究生推免
已经有69人回复
-
求助中医院临床中医中标的青年标书学习，谢谢
已经有4人回复
-
浙大物院表面物理课题组长期招收普博及直博学生
已经有257人回复
-
【征稿】Springer Nature热电专题 Discover Materials / Discover Chemistry
已经有39人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 最小二乘法的matlab 实现
已经有4人回复
- 求高版本1stopt跑一个程序
已经有7人回复
- 求数学大神用matlab帮我解决一个问题。不胜感激
已经有6人回复
- 求助用matlab拟合动力学方程
已经有9人回复
- 用origin8.0英文版或7.5汉化版对Langmuir进行数据的拟合，如何操作。本人就一组数据。
已经有3人回复
- mathcad怎样用四阶龙格库塔法求解二阶微分方程？
已经有8人回复
- 回归分析，急用！
已经有9人回复
- 求前辈帮忙用1stOpt你个一下数据，写论文急用，先谢过
已经有14人回复
- 如何用软件绘制米氏方程曲线啊？？急求
已经有3人回复
- MATLAB微分方程参数拟合问题，求大神
已经有7人回复
- 求助 用matlab拟合三元线性回归方程及分析
已经有9人回复
- 请问matlab拟合二元方程的问题
已经有3人回复
- 1stopt或matlab如何求解以下的非线性方程并拟合出相应曲线？
已经有13人回复
- matlab 拟合反应动力学参数结果很差。大家帮忙看一下
已经有14人回复
- matlab微分方程组参数拟合，以周为单位求解，汇总后以年为单位进行数值比较
已经有4人回复
- matlab拟合方程参数时初值的选择
已经有15人回复
- matlab 求指点 动力学方程拟合过程中导数的获取
已经有10人回复
- 【公式拟合问题】如何拟合下面这样的一个公式，多谢大牛。
已经有6人回复
- 急切求助用matlab曲线拟合色散方程
已经有10人回复
- 【求助】使用Matlab拟合反应动力学方程问题
已经有7人回复
2楼2014-08-21 22:29:33
20