24小时热门版块排行榜

返回列表

当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

艹日大日

木虫 (小有名气)

应助: 0 (幼儿园)
金币: 2410.8
帖子: 112
在线: 67.9小时
虫号: 1108138

[交流] 用origin或1stopt拟合一条隐函数方程曲线

隐函数方程如下：sinh（ay）=bx^c,其中，a，b，c为待拟合参数，x为自变量，y为因变量。拟合数据如下：
x y
1.62255E10 34.359
1.31022E11 45.789
1.62255E11 51.233
1.31022E12 60.937
1.47945E12 67.739
1.62255E12 70.455
1.31022E13 80.543
1.47945E13 82.43
1.62255E13 84.038
2.55026E13 92.766
1.31022E14 101.09
1.47945E14 104.95
2.55026E14 107.12
1.47945E15 116.61
2.55026E15 128.71
2.55026E16 154.61
现在要确定a，b，c三个参数的值，求高手解答

回复此楼

» 猜你喜欢

【复旦大学】二维材料方向招收2026年博士研究生1名已经有0人回复
北京纳米能源与系统研究所王中林院士/曹南颖研究员课题组2026级硕/博/博后招生已经有10人回复
物理学I论文润色/翻译怎么收费? 已经有74人回复
荷兰Utrecht University超快太赫兹光谱王海教授课题招收2026 CSC博士生已经有23人回复
反铁磁体中的磁性切换：两种不同的机制已成功可视化已经有0人回复
求标准粉末衍射卡号 ICDD 01-076-1802 已经有0人回复
新西兰Robinson研究所招收全奖PhD 已经有0人回复
石墨烯转移--二氧化硅衬底石墨烯已经有0人回复
笼目材料中量子自旋液体基态的证据已经有0人回复

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

求教origin拟合简单的渗流方程（又称逾渗方程），做复合材料的应该会遇到吧已经有20人回复
origin8.0拟合langmuir和freundlich原始方程时，未出现曲线而是一条横线已经有2人回复
origin 曲线拟合求助，高手请进已经有11人回复
求ORIGIN拟合包含隐函数的定积分方法已经有6人回复
origin多组数据如何用一条曲线拟合已经有11人回复
【其他】origin8.1 如何找到拟合曲线中任一点对应的X值或Y值已经有5人回复
origin多组数据用一条曲线拟合已经有20人回复
怎么用Origin拟合出过原点的曲线已经有4人回复
用1stopt进行曲线拟合每次给出结果不一样已经有6人回复
origin，sigmaplot等软件自定义曲线拟合（主要是一些动力学方程的模拟）已经有2人回复
求助Origin8.0拟合方程已经有6人回复
求助：用Origin做微分方程的曲线拟合已经有7人回复
【求助】origin拟合外推方程已经有8人回复
origin7.5建立函数拟合为什么会这样？已经有11人回复

» 抢金币啦！回帖就可以得到:

查看全部散金贴

1楼 2011-07-24 11:32:47

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

lipenggg

专家顾问 (正式写手)

★
小木虫: 金币+0.5, 给个红包，谢谢回帖

迭代数: 45
计算用时(时:分:秒:微秒): 00:00:00:762
优化算法: 标准简面体爬山法 + 通用全局优化法
计算结束原因: 达到收敛判断标准
均方差(RMSE): 0
残差平方和(SSE): 0
相关系数(R): 1
相关系数之平方(R^2): 1
决定系数(DC): 1
卡方系数(Chi-Square): 0

参数最佳估算
---------- -------------
a       -7.81770145551721E-18
b       2.16617212990508
c       -2.04267149405294

赞一下

回复此楼

5楼2015-08-05 14:02:29

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

查看全部 5 个回答

锐利的碎片

木虫 (正式写手)

应助: 136 (高中生)
金币: 3637.1
帖子: 988
在线: 1224.9小时
虫号: 961933

★
小木虫(金币+0.5):给个红包，谢谢回帖
艹日大日(金币+1): 2011-07-26 10:28:04

-1.54522485e+000 -1.28408841e-109 1.25973273e+001

CODE:

#!/usr/bin/env python

import scipy.optimize as so

import numpy as np

import pylab as pl

x=np.loadtxt('te.txt')[:,0]

y=np.loadtxt('te.txt')[:,1]

def fun(p):

    return np.arcsinh(p[1]*x**p[2])/p[0]-y

pout=so.leastsq(fun,[1.0,1.0,1.0],maxfev=100000)[0]

print 'a,b,c is ',str(pout)

pl.plot(x,y)

x0=np.linspace(x.min(),x.max(),100)

pl.plot(x0,np.arcsinh(pout[1]*x0**pout[2])/pout[0])

pl.savefig('te.png')