24小时热门版块排行榜

返回列表

q1036210267

木虫 (正式写手)

应助: 26 (小学生)
金币: 2182.2
散金: 806
红花: 4
帖子: 893
在线: 262.1小时
虫号: 1346635
注册: 2011-07-15
性别: GG
专业: 药物化学

引用回帖:

19楼: Originally posted by snoopyzhao at 2011-11-24 09:49:01:
我知道的很少，我自己只用 R 来做，当然，我做的可能比较简单

您是怎么用 R 来做的呢? 能不能教下我啊。

赞一下

回复此楼

生活中坚持快乐的原则不改变，学习工作中踏实上进！`~`

21楼2011-11-24 09:56:14

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

dingd

铁杆木虫 (职业作家)

应助: 1641 (讲师)
金币: 15037.3
散金: 101
红花: 234
帖子: 3410
在线: 1223.5小时
虫号: 291104
注册: 2006-10-28

【答案】应助回帖

q1036210267(金币+2): 2011-12-08 14:41:09

引用回帖:

16楼: Originally posted by q1036210267 at 2011-11-24 08:52:24:
您好，我现在遇到个问题，有两个自变量，两个因变量，两个回归模型共享部分参数。这种非线性方程怎么去解啊？

如下面附图。
有两个方程，y1 和 y2，两个自变量x1和x2。
a, b, c, d 是四个共享的参数。

...

应该是你的版本太老了。

赞一下

回复此楼

22楼2011-11-24 11:50:49

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

q1036210267

木虫 (正式写手)

应助: 26 (小学生)
金币: 2182.2
散金: 806
红花: 4
帖子: 893
在线: 262.1小时
虫号: 1346635
注册: 2011-07-15
性别: GG
专业: 药物化学

内容已删除

回复此楼

生活中坚持快乐的原则不改变，学习工作中踏实上进！`~`

23楼2011-11-24 12:29:47

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

锐利的碎片

木虫 (正式写手)

star watcher

应助: 136 (高中生)
金币: 3637.1
散金: 252
红花: 22
帖子: 988
在线: 1224.9小时
虫号: 961933
注册: 2010-03-05
专业: 凝聚态物性 II ：电子结构

【答案】应助回帖

★ ★
余泽成(金币+2): 谢谢参与应助！ 2011-11-26 23:26:35
q1036210267(金币+2): 2011-12-08 14:42:02

CODE:

import numpy as np

from scipy.optimize import leastsq

dat=np.loadtxt('dat.txt')

x,y=dat[:,0],dat[:,1]

def err(p):

    a,b=p

    return y-a*((b/x)**9-(b/x)**6)

r=leastsq(err,[1,2])

a,b=r[0]

print a,b