| 查看: 1268 | 回复: 4 | ||||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | ||||
[求助]
ln(1+r*(1-x)/x)=rh(1/T-1/484.898) 知道一组T和对应的x值,需要拟合得到r和h的值 已有1人参与
|
||||
|
原公式为ln(1+r*(1-x)/x)=rh(1/T-1/484.898)，也即x=r/(exp(r*h*(1/T-1/484.898))+r-1)，一组对应的T和x值分别为 T X 293.15 0.000126 298.15 0.00016 303.15 0.000205 308.15 0.000253 313.15 0.000315 318.15 0.000385 323.15 0.000466 这个应该可以用Matlab程序编辑计算，输入T和相应的x值，计算得到r和h的值， 这组数据 拟合得到的r值应该很小，在0.008附近，h值很大，在350000附近。 拟合得到的r和h值，代入T得到的相应x计算值与x实际值偏差小于2%。 急求，跪求 希望大神指点。。。 |
回复此楼» 猜你喜欢
» 本主题相关价值贴推荐，对您同样有帮助:
求助，拟合曲线上，如固定X值，求Y值，Y值的误差怎么得到？
已经有3人回复
- 一组数据想用y=x来拟合，求取r2，该用什么软件？怎么做？
已经有6人回复
- 关于origin数据拟合方程后，求一个y值对应的x值
已经有7人回复
- 如何评价一组数值模拟结果与y=x的拟合程度
已经有24人回复
- 求助大侠帮忙，有关matlab模拟数据
已经有29人回复
3楼2015-05-29 21:43:51
feixiaolin
荣誉版主 (文坛精英)
-
专家经验: +518 - 应助: 942 (博后)
- 贵宾: 1.275
- 金币: 3880
- 散金: 58785
- 红花: 532
- 沙发: 11
- 帖子: 24215
- 在线: 2601.8小时
- 虫号: 2139575
- 注册: 2012-11-21
- 专业: 光学信息获取与处理
- 管辖: 数学
★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
zyq19911118(pippi6代发): 金币+30 2015-05-29 21:51:36
zyq19911118(pippi6代发): 金币+30 2015-05-29 21:51:36
|
参数 最佳估算 ---------- ------------- r 0.0337841231633771 h 122839.17638377 ====== 结果输出 ===== No 实测值x 计算值x 1 0.000126 0.0001256 2 0.00016 0.0001594 3 0.000205 0.0002008 4 0.000253 0.0002511 5 0.000315 0.0003119 6 0.000385 0.0003849 7 0.000466 0.0004722 -------------------- 你那个区间限定似乎不妥，误差超出。 |
2楼2015-05-29 17:09:51
wurongjun
专家顾问 (职业作家)
-
专家经验: +831 - 数学EPI: 9
- 应助: 791 (博后)
- 贵宾: 0.308
- 金币: 24609
- 散金: 310
- 红花: 75
- 帖子: 3004
- 在线: 881.4小时
- 虫号: 1368482
- 注册: 2011-08-14
- 性别: GG
- 专业: 计算数学与科学工程计算
- 管辖: 数学
【答案】应助回帖
★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与，应助指数 +1
zyq19911118(feixiaolin代发): 金币+40 2015-05-30 11:05:17
感谢参与，应助指数 +1
zyq19911118(feixiaolin代发): 金币+40 2015-05-30 11:05:17
|
程序: function nihrh Data=[ 293.15 0.000126 298.15 0.00016 303.15 0.000205 308.15 0.000253 313.15 0.000315 318.15 0.000385 323.15 0.000466]; a0=[0.008,350000]%初值 b=nlinfit(Data(:,1),Data(:,2),@frh,a0); X=frh(b,Data(:,1)); disp(['最优值r=',num2str(b(1),6),',最优值h='num2str(b(2),6)]) disp(['拟合值:']) disp([num2str(X,6)]) function x=frh(a,T) x=a(1)./(exp(a(1)*a(2)*(1./T-1/484.898))+a(1)-1); 运行: >> nihrh a0 = 1.0e+005 * 0.0000 3.5000 最优值r=0.00840754,最优值h=350000 拟合值: 0.000161801 0.000192137 0.000227013 0.000266957 0.000312548 0.000364426 0.000423294 |
4楼2015-05-30 10:43:45
feixiaolin
荣誉版主 (文坛精英)
-
专家经验: +518 - 应助: 942 (博后)
- 贵宾: 1.275
- 金币: 3880
- 散金: 58785
- 红花: 532
- 沙发: 11
- 帖子: 24215
- 在线: 2601.8小时
- 虫号: 2139575
- 注册: 2012-11-21
- 专业: 光学信息获取与处理
- 管辖: 数学
5楼2015-05-30 11:03:11
