| 查看: 1441 | 回复: 9 | ||
liuliping3银虫 (小有名气)
|
[求助]
Matlab自定义拟合函数编程问题 已有2人参与
|
|
matlab里面做自定义拟合，我不会编码，由于时间紧迫，想让大家帮帮忙，不甚感激 x=[0.5 1 1.5 2 2.5 3 3.5 4 4.5 5 5.5 6 6.5 7 8.7 10.4 12.1 13.8 15.5 17.2 18.9 20.6 22.3 24]; y=[0.96564 0.96454 0.95738 0.9543 0.82423 0.6315 0.51949 0.36035 0.2554 0.20826 0.14967 0.06784 0.05 0.01123 0.00749 0.000220264 0 0 0 0 0 0 0 0]; 用Matlab拟合方程的参数，方程如下： 其中：k为待求的参数， P=101325, Y=0.0112, M=0.032, t=4/60, S=9.08, ρ=2857.14， R=8.314，T=973，u=815.29 是对应于两组x, y的已知值。不同P, Y, M, t, S, ρ， R，T，u下的x, y值不同，x, y值都为已知。现在希望利用这两组数据求参数k。同时画出原图和拟合图形。  1.gif |
回复此楼» 猜你喜欢
关于代码变化问题，想知道的进来
已经有14人回复
-
这样的filecode谁见过
已经有5人回复
-
国基金的申报应该改成非等额制，评价高的钱多评价低的钱少，但是增加资助率
已经有4人回复
-
售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急
已经有4人回复
-
售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急
已经有6人回复
-
fileCode有新解读？
已经有16人回复
-
售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急
已经有5人回复
-
售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急
已经有4人回复
-
售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急
已经有4人回复
-
售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急
已经有3人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- Origin自定义模拟与matlab中cftool工具箱拟合结果比较
已经有4人回复
- 通过Matlab软件自行计算了一个自由曲面拟合成非球面多项式或者zernike多项式的问题
已经有3人回复
- 关于使用origin数据拟合的问题请教。
已经有6人回复
- 【求助】请各位高手告诉我：这种想法能不能用matlab实现？200金币求助
已经有16人回复
dingd
铁杆木虫 (职业作家)
- 应助: 1641 (讲师)
- 金币: 15037.3
- 散金: 101
- 红花: 234
- 帖子: 3410
- 在线: 1223.7小时
- 虫号: 291104
- 注册: 2006-10-28
2楼2015-11-25 21:00:47
liuliping3
银虫 (小有名气)
- 应助: 1 (幼儿园)
- 金币: 60.8
- 散金: 50
- 帖子: 54
- 在线: 36.6小时
- 虫号: 2125053
- 注册: 2012-11-13
- 性别: MM
- 专业: 分离过程
|
x=[0.5 1 1.5 2 2.5 3 3.5 4 4.5 5 5.5 6 6.5 7 8.7 10.4 12.1 13.8 15.5 17.2 18.9 20.6 22.3 24]; y=[0.96564 0.96454 0.95738 0.9543 0.82423 0.6315 0.51949 0.36035 0.2554 0.20826 0.14967 0.06784 0.05 0.01123 0.00749 0.000220264 0 0 0 0 0 0 0 0]; y=exp(k*P*Y*M*t/S)/(exp(k*P*Y*M*t/S)+exp(k*ρ*R*T*x/u)-1); P=101325, Y=0.0112, M=0.032, t=4/60, S=9.08, ρ=2857.14， R=8.314，T=973，u=815.29 |
3楼2015-11-25 21:56:49
liuliping3
银虫 (小有名气)
- 应助: 1 (幼儿园)
- 金币: 60.8
- 散金: 50
- 帖子: 54
- 在线: 36.6小时
- 虫号: 2125053
- 注册: 2012-11-13
- 性别: MM
- 专业: 分离过程
|
x=[0.5 1 1.5 2 2.5 3 3.5 4 4.5 5 5.5 6 6.5 7 8.7 10.4 12.1 13.8 15.5 17.2 18.9 20.6 22.3 24]; y=[0.96564 0.96454 0.95738 0.9543 0.82423 0.6315 0.51949 0.36035 0.2554 0.20826 0.14967 0.06784 0.05 0.01123 0.00749 0.000220264 0 0 0 0 0 0 0 0]; y= exp(k*P*Y*M*t/S)/( exp(k*P*Y*M*t/S)+exp(k*ρ*R*T*x/u)-1); P=101325, Y=0.0112, M=0.032, t=4/60, S=9.08, ρ=2857.14， R=8.314，T=973，u=815.29 |
4楼2015-11-25 21:59:44
dingd
铁杆木虫 (职业作家)
- 应助: 1641 (讲师)
- 金币: 15037.3
- 散金: 101
- 红花: 234
- 帖子: 3410
- 在线: 1223.7小时
- 虫号: 291104
- 注册: 2006-10-28
【答案】应助回帖
|
不是两组x, y数值吗？对应不同的P, Y, M, t, S, ρ， R，T，u值 如果只有这一组，1stOpt试试： 均方差(RMSE):0.13879469945569 残差平方和(SSE):0.462335246327889 相关系数(R): 0.958172545145573 相关系数之平方(R^2): 0.918094626270745 确定系数(DC): 0.860178141633635 卡方系数(Chi-Square): 1.00096945974885 参数 最佳估算 ---------- ------------- k 8.31057566523805E-6 |
5楼2015-11-25 22:46:38
liuliping3
银虫 (小有名气)
- 应助: 1 (幼儿园)
- 金币: 60.8
- 散金: 50
- 帖子: 54
- 在线: 36.6小时
- 虫号: 2125053
- 注册: 2012-11-13
- 性别: MM
- 专业: 分离过程
6楼2015-11-26 09:09:52
liuliping3
银虫 (小有名气)
- 应助: 1 (幼儿园)
- 金币: 60.8
- 散金: 50
- 帖子: 54
- 在线: 36.6小时
- 虫号: 2125053
- 注册: 2012-11-13
- 性别: MM
- 专业: 分离过程
7楼2015-11-26 09:23:17
dingd
铁杆木虫 (职业作家)
- 应助: 1641 (讲师)
- 金币: 15037.3
- 散金: 101
- 红花: 234
- 帖子: 3410
- 在线: 1223.7小时
- 虫号: 291104
- 注册: 2006-10-28
8楼2015-11-26 11:12:15
9楼2015-11-26 14:28:57
liuliping3
银虫 (小有名气)
- 应助: 1 (幼儿园)
- 金币: 60.8
- 散金: 50
- 帖子: 54
- 在线: 36.6小时
- 虫号: 2125053
- 注册: 2012-11-13
- 性别: MM
- 专业: 分离过程
10楼2015-11-26 15:01:50
20