| 查看: 1783 | 回复: 5 | |||
rhinoxixi金虫 (正式写手)
|
[交流]
【求助】递推公式的参数拟合已有2人参与
|
|
自己建立的一个模型，是一个递推公式： Yn=a+b*(Xn)^c+d*Yn-1 要拟合参数，除了自己编程之外，有没有什么软件可以做到的？ matlab行吗？how？ |
回复此楼» 猜你喜欢
第一性原理计算方向2026级博士申请 PRB*1，四级484
已经有1人回复
-
求助VISSIM破解版软件
已经有0人回复
-
物理学I论文润色/翻译怎么收费?
已经有298人回复
-
求2026年在台湾举行的物理和材料领域国际学术会议信息
已经有0人回复
-
求国际会议网站
已经有1人回复
-
求取一些关于纳米材料和纳米技术相关的英文PPT。
已经有0人回复
-
【复旦大学】二维材料方向招收2026年博士研究生1名
已经有0人回复
-
北京纳米能源与系统研究所 王中林院士/曹南颖研究员课题组2026级硕/博/博后招生
已经有10人回复
-
荷兰Utrecht University超快太赫兹光谱王海教授课题招收2026 CSC博士生
已经有17人回复
-
反铁磁体中的磁性切换：两种不同的机制已成功可视化
已经有0人回复
» 本主题相关价值贴推荐，对您同样有帮助:
根据数据，如何得到拟合公式
已经有7人回复- 为什么同样的数据，origin 和excel做成的拟合公式不同？
已经有13人回复
- 【致H兄】数列题的万能多项式拟合求法
已经有28人回复
- 拟合模型的好坏
已经有4人回复
- 【公式拟合问题】如何拟合下面这样的一个公式，多谢大牛。
已经有6人回复
- 请问各位origin8版里面的非线性拟合中，选择的LangmuirEXT1公式是什么意思啊
已经有14人回复
- 实验数据的数学模型构建、曲线拟合方面问题求助
已经有6人回复
- Origin参数拟合，自定义公式编译通不过
已经有1人回复
- origin 8.0自定义公式拟合对公式长度有限制吗？
已经有7人回复
- 【求助】使用MATLAB拟合出按照公式的系数
已经有19人回复
- 【求助】Arcgis地统计学半变异函数模型拟合
已经有7人回复
- 如何对阿伦尼乌斯公式进行曲线拟合从而得到活化能和指前因子
已经有14人回复
- 【交流】有没有药动学的高手：不同给药方式拟合出来的房室模型不一样，大家遇到过没？
已经有7人回复
- origin中非线性拟合中logistic模型的疑问？
已经有6人回复
2楼2009-04-11 16:41:17
rhinoxixi
金虫 (正式写手)
- 应助: 0 (幼儿园)
- 金币: 2499.1
- 散金: 15
- 帖子: 434
- 在线: 177.6小时
- 虫号: 476811
- 注册: 2007-12-13
- 专业: 污染物行为过程及其环境效
谢谢热心的fspdlh！
|
具体说明如下 ： n Xn Yn 1 0.82900 0.12500 2 0.81600 0.14600 3 0.64300 0.24300 4 0.57000 0.24700 5 0.47200 0.26800 6 0.36500 0.28400 7 0.24500 0.30700 8 0.11800 0.27100 9 0.37200 0.30000 10 0.55800 0.32400 Y0=0 Y1=a+b*X1^c+d*Y0 Y2=a+b*X2^c+d*Y1 ....... 我要拟合出a、b、c、d 其实我知道编程用最小二乘法可以求出来，但是我不是很懂编程，所以问有没有软件算。如果有人fortran很厉害，也可以帮我编个程序来拟合参数，谢谢啦！ |
3楼2009-04-11 22:31:11
★ ★ ★ ★ ★
rhinoxixi(金币+5,VIP+0):你实在是太太太好人了，像天使一样～谢谢你！ 4-12 00:46
rhinoxixi(金币+5,VIP+0):你实在是太太太好人了，像天使一样～谢谢你！ 4-12 00:46
|
用MATLAB做的，beta为结果，不知是否符合要求，可能换别的初值会更好，楼主可以多试试。 function [beta r]=datafit(beta) close all hold on load data.txt X=data(:,2); Y=data(:,3); beta=nlinfit(X,Y,'fun',beta); Y1=fun(beta,X); r=corrcoef(Y,Y1); beta=beta'; r=r(1,2); function Y=fun(beta,X) a=beta(1); b=beta(2); c=beta(3); d=beta(4); n=length(X); Y(1)=a+b*X(1)^c; for i=2:n Y(i)=a+b*X(i)^c+d*Y(i-1); end Y=Y'; >> [beta r]=datafit([0 1 0 1]) Warning: Iteration limit exceeded. Returning results from final iteration. > In nlinfit at 220 In datafit at 7 beta = -4.4118 4.5129 0.0035 0.7410 r = 0.9689 >> [ Last edited by fspdlh on 2009-4-11 at 23:29 ] |
4楼2009-04-11 23:26:06
dingd
铁杆木虫 (职业作家)
- 计算强帖: 4
- 应助: 1641 (讲师)
- 金币: 15037.3
- 散金: 101
- 红花: 234
- 帖子: 3410
- 在线: 1223.5小时
- 虫号: 291104
- 注册: 2006-10-28
★ ★ ★
小木虫(金币+0.5):给个红包，谢谢回帖交流
dongdong3881(金币+2): 多谢指导 2011-02-21 15:17:50
小木虫(金币+0.5):给个红包，谢谢回帖交流
dongdong3881(金币+2): 多谢指导 2011-02-21 15:17:50
|
用1stOpt处理，结果更好更稳定，也不需要猜初值： 均方差(RMSE): 0.0154322799500792 残差平方和(SSE): 0.00238155264457617 相关系数(R): 0.971104117667667 相关系数之平方(R^2): 0.943043207351097 参数 最佳估算 -------------------- ------------- a 0.100663466550055 b -2.12112269614125E-24 c -23.8930576866545 d 0.69639916154168 |
5楼2011-02-21 09:57:14
zt_chem
木虫 (正式写手)
- 应助: 12 (小学生)
- 金币: 2147.8
- 红花: 6
- 帖子: 872
- 在线: 981.6小时
- 虫号: 286920
- 注册: 2006-10-15
- 性别: GG
- 专业: 催化化学
6楼2012-01-06 18:42:10