24小时热门版块排行榜    

查看: 1469  |  回复: 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。同时画出原图和拟合图形。

Matlab自定义拟合函数编程问题
1.gif
回复此楼
我就是我
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

dingd

铁杆木虫 (职业作家)

【答案】应助回帖

不是两组x, y数值吗?对应不同的P, Y, M, t, S, ρ, R,T,u值

如果只有这一组,1stOpt试试:
CODE:
Constant P=101325, Y1=0.0112, M=0.032, t1=4/60, S=9.08, p1=2857.14, R=8.314,T=973,u=815.29;
Function y= exp(k*P*Y1*M*t1/S)/( exp(k*P*Y1*M*t1/S)+exp(k*p1*R*T*x/u)-1);
Data;
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];

均方差(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
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 10 个回答

dingd

铁杆木虫 (职业作家)

【答案】应助回帖

感谢参与,应助指数 +1
两组值在哪儿?公式最好先给出文本格式
2楼2015-11-25 21:00:47
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

liuliping3

银虫 (小有名气)

引用回帖:
2楼: Originally posted by dingd at 2015-11-25 21:00:47
两组值在哪儿?公式最好先给出文本格式

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
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

liuliping3

银虫 (小有名气)

引用回帖:
2楼: Originally posted by dingd at 2015-11-25 21:00:47
两组值在哪儿?公式最好先给出文本格式

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
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 不应该看fileCode +6 且听虎啸 2026-08-12 8/400 2026-08-13 12:07 by 小木虫6752397
[基金申请] 重要来源:本周末出结果 +9 瞬息宇宙 2026-08-12 9/450 2026-08-13 10:12 by Vivilian
[基金申请] 我的国基提前知道中了,可是同事的操作让我实在接受不了,怎么会有这样的人 +9 家与远方 2026-08-10 14/700 2026-08-13 09:41 by 臭臭不臭01
[基金申请] 结合人工智能,周易传统文化,filecode打分制来了,3分以上希望很大。 +3 Tide man 2026-08-12 4/200 2026-08-13 08:35 by ZJTJZ
[基金申请] 2019年青年基金涵评意见,大家看看几个A,几个B? +11 Tide man 2026-08-11 11/550 2026-08-13 07:35 by 撸猫猫
[基金申请] 好奇怪的filecode +5 布布和一二 2026-08-08 6/300 2026-08-12 16:11 by 云上清扬
[基金申请] filecode +13 documentary 2026-08-10 15/750 2026-08-12 15:56 by 云上清扬
[论文投稿] 职称评审,求友友推荐见刊最快的期刊 +4 工厂打螺丝 2026-08-08 4/200 2026-08-12 09:18 by hansi2025
[基金申请] 有时候,自然基金真的不能太认真 (我的申报经验) +3 majunge000 2026-08-11 4/200 2026-08-11 20:13 by lch2012
[基金申请] 帮忙看看fileCode +7 wwncly 2026-08-10 13/650 2026-08-11 19:36 by 冰心玉壶晴
[基金申请] 确定了,国自然21号放榜 +6 布布和一二 2026-08-10 7/350 2026-08-10 19:15 by 2000zf36392
[基金申请] 静等基金结果 +5 gjjjzhong 2026-08-10 16/800 2026-08-10 17:19 by Tide man
[基金申请] 国自然结果 +4 Vierhys 2026-08-10 8/400 2026-08-10 15:06 by Vierhys
[基金申请] 2026国自然放榜时间 +9 布布和一二 2026-08-08 9/450 2026-08-10 11:22 by xxxx2020
[基金申请] filecode与中标关系的预测 +5 布布和一二 2026-08-07 5/250 2026-08-09 16:15 by 袁向阳007
[基金申请] 关于filecode,很负责任的告诉大家 +6 爱看书的可乐 2026-08-08 7/350 2026-08-08 22:13 by a_niu
[基金申请] 国基金的申报应该改成非等额制,评价高的钱多评价低的钱少,但是增加资助率 +7 a089 2026-08-07 7/350 2026-08-08 18:05 by gltch
[基金申请] 关于filecode +4 布布和一二 2026-08-07 7/350 2026-08-07 22:55 by zhanghaozhu
[基金申请] 固定端突然变了,今天 +6 archvillain 2026-08-06 10/500 2026-08-07 16:03 by 医学老男孩
[基金申请] filecode变化情况 +6 布布和一二 2026-08-07 22/1100 2026-08-07 14:45 by 且听虎啸
信息提示
请填处理意见