24小时热门版块排行榜    

查看: 811  |  回复: 3

山鹰HPU

新虫 (初入文坛)

[求助] MATLAB多远函数拟合 已有1人参与

现有以下数据要进行拟合,采用关系式:
y=5.2/(5.2+u)*exp(57.2-11*u)*erf((45*R-u*x)/(2*(2.1*R*x)^0.5))+ 5.2/(5.2+u)*exp(57.2-11*u)*erf((45*R+u*x)/(2*(2.1*R*x)^0.5))+
6.4/(R*k)*exp(111.4-k*x)*erf((45*R+5.2*x)/(4.2*R*x)
其中 u=(27.04+8.4*R*k)^0.5
erf(x)是MATLAB里的内置余误差函数。
[x y]=[1        0
2        2.81
3        4.22
4        2.81
5        3.52
5.67        4.92
6.17        2.81
6.67        4.22
7.17        4.92
7.67        5.62
8.17        4.92
8.67        6.33
9.17        5.62
9.67        5.62
10.17        7.03
10.67        6.33
11.17        7.03
11.67        7.73
12.17        4.92
14.17        9.84
15.17        9.84
16.17        14.05
17.17        15.46
18.17        16.86
19.17        16.16
20.17        17.56
21.17        17.56
22.17        16.16
23.17        16.86
24.17        16.16
25.17        16.86
26.17        15.46
27.17        16.86
28.17        18.27
29.17        18.27
30.17        16.16
31.17        12.65
32.17        14.05
33.17        14.75
34.17        14.75
35.17        14.75
36.17        15.46
37.17        13.35
38.17        14.05
];
拟合之后把代码,参数,图发过来,谢谢!
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐,对您同样有帮助:

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

hyh4402031

金虫 (小有名气)

【答案】应助回帖

感谢参与,应助指数 +1
都没有变量是什么怎么帮你,nlinfit试试,多元变量写成x(1),x(2)……
2楼2014-11-07 19:21:33
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

山鹰HPU

新虫 (初入文坛)

引用回帖:
2楼: Originally posted by hyh4402031 at 2014-11-07 19:21:33
都没有变量是什么怎么帮你,nlinfit试试,多元变量写成x(1),x(2)……

已经给出x和y的一组数据,关系式里只有R和k是未知量,不能拟合出参数R和K吗,本人对MATLAB不是太懂,还请多多帮忙,谢谢!
3楼2014-11-07 21:09:11
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

hyh4402031

金虫 (小有名气)

★ ★ ★
费姚永芬: 金币+3, 谢谢 交流 2014-11-08 12:01:36
引用回帖:
3楼: Originally posted by 山鹰HPU at 2014-11-07 21:09:11
已经给出x和y的一组数据,关系式里只有R和k是未知量,不能拟合出参数R和K吗,本人对MATLAB不是太懂,还请多多帮忙,谢谢!...

可以的,但这不收多元拟合,这叫一元非线性拟合。1 使用工具箱的话,可以使用cftool,里面可以自定义函数。2 使用命令栏的话,也可以使用fittypep定义函数,再使用fit拟合。
比nlinfit的好处是,可以限定参数的范围
4楼2014-11-08 10:06:45
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 山鹰HPU 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[硕博家园] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +3 f1dmPx051StA 2026-08-03 7/350 2026-08-04 09:31 by v1TJC4qJauAB
[博后之家] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +3 jRl3mE6ddZGq 2026-08-03 5/250 2026-08-04 09:27 by v1TJC4qJauAB
[教师之家] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +4 jRl3mE6ddZGq 2026-08-03 10/500 2026-08-04 09:26 by v1TJC4qJauAB
[博后之家] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +4 cu9Nq1xK233Z 2026-08-03 9/450 2026-08-04 09:17 by v1TJC4qJauAB
[教师之家] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +5 cu9Nq1xK233Z 2026-08-03 10/500 2026-08-04 09:15 by v1TJC4qJauAB
[基金申请] 系统今天又提示维护了,估计离放榜不远了 +12 winnerche 2026-07-29 16/800 2026-08-04 08:18 by Equinoxhua
[有机交流] 一个有机合成实验室都需要哪些设备? 50+3 kf2781974 2026-07-31 10/500 2026-08-04 08:05 by 88817753
[找工作] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +3 cu9Nq1xK233Z 2026-08-03 9/450 2026-08-04 07:37 by k7OM8YghWbkC
[硕博家园] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +4 cu9Nq1xK233Z 2026-08-03 13/650 2026-08-04 07:30 by k7OM8YghWbkC
[教师之家] 基础研究怎么拉横向,学校到款任务越来越多,难以完成 拉横向,都有哪些途径啊 +9 锦衣卫寒战 2026-07-28 13/650 2026-08-04 07:21 by Ermito
[硕博家园] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +3 jRl3mE6ddZGq 2026-08-03 6/300 2026-08-04 04:47 by k7OM8YghWbkC
[公派出国] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +3 cu9Nq1xK233Z 2026-08-03 7/350 2026-08-04 04:46 by k7OM8YghWbkC
[考研] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +3 cu9Nq1xK233Z 2026-08-03 8/400 2026-08-04 03:21 by k7OM8YghWbkC
[基金申请] 什么时候能放榜呀? +3 Jacob678 2026-08-03 3/150 2026-08-03 16:14 by gltch
[基金申请] 面上再次挂了,太难了,躺也躺不了,倦也卷不过,小学校之殇! +19 低垂的野花 2026-07-31 25/1250 2026-08-03 09:25 by gy116024
[基金申请] 2026年国自然面上资助率 +16 布布和一二 2026-07-30 22/1100 2026-08-03 09:20 by gy116024
[基金申请] 微信指数没变化,科研之友没阅读 +15 wangze12014 2026-07-28 19/950 2026-08-02 20:05 by 蔡棒棒菂
[考博] 2027年申博 50+3 射雕英雄胜 2026-07-30 3/150 2026-08-02 09:36 by lfy8008
[高分子] HXDI做水性聚氨酯乳液,是不是特别容易出渣 15+3 yuyusuv 2026-07-29 3/150 2026-07-31 09:21 by huizingga
[基金申请] 你们的时间戳变了吗 +3 archvillain 2026-07-30 4/200 2026-07-30 18:53 by levinzhwen
信息提示
请填处理意见