| 查看: 2371 | 回复: 4 | |||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | |||
tdwcn木虫 (正式写手)
金虫
|
[交流]
【求助】matlab函数拟合，拟合结果并保存，方便调用 已有1人参与
|
||
|
现有一组数据，如下面所示，想编写一个子函数：函数参数为x,y，函数的返回值为x,y通过最小二乘法拟合（或者其他的拟合方法）得到的拟合方程。 并把这个方程保存下来，方便下次可以直接代入x的不同值（如5.3,6.8之类的数据），得到对应的y值。 请各位高手帮忙。 x=[0 5 10 15 20 25 30 35 40 45 50 55 60 65 70 75 80 85 90 95 100]; y=[0 8.8 18 27.6 35.6 42.2 48.1 54.6 60.1 63.4 68.6 73.1 76.9 80.2 84 88.4 91.7 93.2 95.6 98.1 100]; 函数模型为：f(a,x)=a(1)*x.^a(2).*exp(a(3)*x); |
回复此楼» 猜你喜欢
售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急
已经有6人回复
-
售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急
已经有8人回复
-
售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急
已经有8人回复
-
售SCI一区T0P文章，我:8O.55.1.O.54,科目全,可伽急
已经有10人回复
-
2026国自然函评费到账
已经有11人回复
-
售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急
已经有3人回复
-
售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急
已经有3人回复
-
93BebMhtakh前后11位开头都是大写
已经有8人回复
-
售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急
已经有5人回复
-
范进中举一文的中心思想
已经有4人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- matlab拟合方程参数时初值的选择
已经有15人回复
- matlab拟合求值
已经有17人回复
- 在matlab中拟合曲面的表达式
已经有29人回复
- matlab拟合拟均相动力学参数
已经有6人回复
- 需要用matlab拟合数据的外行人急需帮助
已经有10人回复
- matlab的简单函数拟合
已经有6人回复
- 求助matlab隐函数参数拟合
已经有16人回复
- 【求助】使用MATLAB拟合出按照公式的系数
已经有19人回复
- 【求助】急求关于matlab做非线性拟合的函数
已经有8人回复
- 如何用MATLAB拟合指数函数
已经有3人回复
- 【其他】自定义的二元函数拟合
已经有6人回复
tdwcn
木虫 (正式写手)
金虫
- 应助: 0 (幼儿园)
- 金币: 4351.2
- 散金: 2
- 红花: 2
- 帖子: 479
- 在线: 130.6小时
- 虫号: 669188
- 注册: 2008-12-06
- 性别: GG
- 专业: 光学
5楼2010-10-21 10:35:33
tdwcn
木虫 (正式写手)
金虫
- 应助: 0 (幼儿园)
- 金币: 4351.2
- 散金: 2
- 红花: 2
- 帖子: 479
- 在线: 130.6小时
- 虫号: 669188
- 注册: 2008-12-06
- 性别: GG
- 专业: 光学
2楼2010-08-15 14:22:11
★
余泽成(金币+1):谢谢参与应助！ 2010-08-17 19:59:34
tdwcn(金币+1): 2010-10-21 10:34:49
余泽成(金币+1):谢谢参与应助！ 2010-08-17 19:59:34
tdwcn(金币+1): 2010-10-21 10:34:49
|
本帖内容被屏蔽 |
3楼2010-08-17 15:08:44
urion
铁杆木虫 (职业作家)
- 应助: 8 (幼儿园)
- 贵宾: 0.361
- 金币: 12933.2
- 散金: 55
- 红花: 2
- 帖子: 3619
- 在线: 199.6小时
- 虫号: 437304
- 注册: 2007-08-26
- 专业: 外国考古
tdwcn(金币+1): 2010-10-21 10:34:52
|
参考一些代码: function nonlinefit clc;clear; t=[0 4 8 40]; y=[20.09 64.52 85.83 126.75]; beta=nlinfit(t,y,@myfunc,[1 1 1 1]) a=beta(1) k1=beta(2) k2=beta(3) m=beta(4) tt=0:1:40 yy=a+k1*exp(m*tt)+k2*exp(-m*tt) plot(t,y,'o',tt,yy) function y1=myfunc(const,t) a=const(1); k1=const(2); k2=const(3); m=const(4); y1=a+k1*exp(m*t)+k2*exp(-m*t); |
4楼2010-10-19 20:40:05
