版块导航: 正在加载中...

登录注册

应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜

返回列表

ponder_zsm

新虫 (初入文坛)

应助: 0 (幼儿园)
金币: 3.5
帖子: 1
在线: 1.5小时
虫号: 3000540
注册: 2014-02-27
专业: 生物医用高分子材料

[求助] 如何用给定的带有参数的函数拟合已知数据？已有1人参与

已知实验数据和函数原型如下：
data = {{8.5, 0.828}, {12.5, 0.78}, {16.5, 0.741}, {20.4, 0.714}, {23.1, 0.66}, {29.4, 0.558}, {39.4, 0.491}, {49.4, 0.454}};
model = Integrate[(1 - x/r)^2*Exp[-(Log[r/b*Sqrt[1 + (a*a)/(b*b)]])^2/(2 Log[1 + (a*a)/(b*b)])]/(r*Sqrt[2*Pi*Log[1 + (a*a)/(b*b)]]), {r,0, Infinity}];
这是在Mathematica中输入的格式，函数原型可以参见上传的附件。其中a, b是参数，x是自变量。用Findfit求a和b时，Mathematica给出的结果不对：
{a -> 160.339, b -> 186.173}
合理值应该在a=299, b=247左右，在用所给值验证时，所得到的函数值为虚数：
f[a_, b_, x_, r_] :=Integrate[(1 - x/r)^2*Exp[-(Log[r/b*Sqrt[1 + (a*a)/(b*b)]])^2/(2 Log[1 + (a*a)/(b*b)])]/(r*Sqrt[2*Pi*Log[1 + (a*a)/(b*b)]]), {r,0, Infinity}];
f[160.339, 186.173, 8.5, r]
输出为：0.851972 - 4.76418*10^-15 I
不知道哪里可以改进，求高人指点！

另外我尝试用1st opt，但是也不成功：
Constant h=10000;
Parameters a,b;
Variables x,y;
Function y=int(sqr(1-x/r)*exp(-sqr(ln(r*sqrt(1+(a*a)/(b*b))/b))/(2*ln(1+(a*a)/(b*b))))/(r*sqrt(2*Pi*ln(1+(a*a)/(b*b)))),r=0,h);
Data;
8.5 0.828
12.5 0.78
16.5 0.741
20.4 0.714
23.1 0.66
29.4 0.558
39.4 0.491
49.4 0.454
提示函数表达式错误！求指点！！

函数.png

回复此楼

» 猜你喜欢

国社科又开始会评了，不知道这次命运如何已经有4人回复
售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急已经有8人回复
售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急已经有7人回复
售SCI文章，我:8O.5.5.1O.54,科目全,可十急已经有8人回复
售SCI一区T0P文章，我:8O.55.1.O.54,科目全,可伽急已经有8人回复
售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急已经有8人回复
售SCI一区文章，我:8.O.551.O.5.4,科目全,可伽急已经有5人回复
找导师已经有6人回复
售一区SCI文章T0P，我:8O.551.O54,科目全,可十急已经有7人回复
售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急已经有9人回复

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

origin8.0 函数关系中参数拟合求助？已经有4人回复
急！在用origin进行非线性拟合怎样创建新的拟合函数？已经有14人回复
origin 里 “三次样条函数拟合” 与“三次多项式拟合”是一回事吗？如图已经有3人回复
想用GULP拟合纯钼的morse势函数参数已经有13人回复
求助该图形用什么函数曲线拟合，已经有3人回复
Aguado-Paniagua函数拟合已经有6人回复
1stopt 拟合隐函数的四参数已经有7人回复
如何用已知函数去拟合实验数据，得出参数已经有10人回复
关于matlab拟合已经有4人回复
在用origin进行非线性拟合怎样去选择拟合函数？已经有6人回复
通过实验数据拟合，求解公式中的参数已经有14人回复
求助一个逻辑斯蒂函数的拟合已经有8人回复
自定义函数拟合请假大侠帮忙！非常非常感激不尽！！！已经有9人回复
matlab高手帮忙做个自定义复杂函数拟合，含有复数，已经有25人回复
已知函数形式，如何拟合其中的多个系数。已经有3人回复
matlab的残差大小跟什么因素有关？拟合的函数选择？还是给定值的范围？已经有4人回复
跪求1stOpt2.5以上版本，或者请高人帮忙拟合包含复数的函数的参数已经有8人回复
origin进行非线性拟合如何选择函数或者自定义合适的函数已经有19人回复
请问一个有理函数的曲线拟合过程如何做，给出数据，请求拟合，已经有14人回复
【公式拟合问题】如何拟合下面这样的一个公式，多谢大牛。已经有6人回复
【求助】急求关于matlab做非线性拟合的函数已经有8人回复
【求助】matlab函数拟合，拟合结果并保存，方便调用已经有4人回复

1楼 2014-10-27 02:47:26

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

xzczd

木虫 (小有名气)

应助: 72 (初中生)
金币: 1941.4
红花: 11
帖子: 290
在线: 144.2小时
虫号: 3305534
注册: 2014-07-03

【答案】应助回帖

楼主你确定真值在{a->299, b->247}附近？

下图是{a -> 160.339, b -> 186.173}时的拟合效果图：

这拟合效果不算差。

下图是你给的参数对应的图：

囧。

然后，关于直接反代导致虚数的问题：这只是一点数值误差罢了，如你所见，虚部的值非常微小。想要避免这个问题，你可以把全体参数都改成准确值（也就是别用小数），也可以把Integrate改成NIntegrate然后再加个高一点的WorkingPrecision。

赞一下