发信人: jzxu (自然), 信区: Physics
标 题: 问个多参数fitting的问题 (转载)
发信站: BBS 未名空间站 (Wed May 14 22:19:44 2008), 站内

【 以下文字转载自 EE 讨论区 】
发信人: jzxu (自然), 信区: EE
标 题: 问个多参数fitting的问题
发信站: BBS 未名空间站 (Wed May 14 22:18:54 2008), 站内

发信人: jzxu (自然), 信区: Computation
标 题: 问个多参数fitting的问题
发信站: BBS 未名空间站 (Wed May 14 22:17:06 2008), 转信

公式很简单

dM1/dt = -qM1-c1M1+p1c2M2

dM2/dt = -qM2-c2M2+c1M1+c3M3

dM3/dt = -qM3-c3M3+p2c2M2

initial condition: c1M1+c2M2=c3M3, p1+p2=1

S(q)=M1+M2+M3

现在已有数据系列数据S(q)，请问如何做data fitting 得到参数c1,c2,c3和p1,p2? 本
来打算用lsqcurvfit和ode45，但是总是搞不定，大家有什么好办法吗？



--

※ 来源:·BBS 未名空间站 海外: mitbbs.com 中国: mitbbs.cn·[FROM: 66.79.]

发信人: Tecnai (绿洲), 信区: Physics
标 题: Re: 问个多参数fitting的问题 (转载)
发信站: BBS 未名空间站 (Thu May 15 21:18:12 2008), 转信

5个参数？
我记得费曼似乎可以用他们来拟合一头大象

【 在 jzxu (自然) 的大作中提到: 】
: 【 以下文字转载自 EE 讨论区 】
: 发信人: jzxu (自然), 信区: EE
: 标 题: 问个多参数fitting的问题
: 发信站: BBS 未名空间站 (Wed May 14 22:18:54 2008), 站内
: 发信人: jzxu (自然), 信区: Computation
: 标 题: 问个多参数fitting的问题
: 发信站: BBS 未名空间站 (Wed May 14 22:17:06 2008), 转信
: 公式很简单
: dM1/dt = -qM1-c1M1+p1c2M2
: dM2/dt = -qM2-c2M2+c1M1+c3M3
: ...................

--
于是,我似乎看到多年后,你站在一扇蓝色大门前……
下午三点的阳光,你仍然有几颗青春痘,你笑着,我跑向你问你好不好,你点点头
三年,五年,甚至更久以后,我们会变成什么样的大人呢?
虽然,我闭上眼睛也看不到自己,但是,我却可以看见你.



※ 来源:·BBS 未名空间站 海外: mitbbs.com 中国: mitbbs.cn·[FROM: 128.180.]

发信人: dingle (快准狠), 信区: Physics
标 题: Re: 问个多参数fitting的问题 (转载)
发信站: BBS 未名空间站 (Fri May 16 11:51:32 2008)

我也想知道。我以前做的时候就是估计其中一两个参数的合理值，然后
基本上就是试出来的。不是很清楚有没有比较系统的方法。


【 在 jzxu (自然) 的大作中提到: 】
: 发信人: jzxu (自然), 信区: EE
: 标 题: 问个多参数fitting的问题
: 发信站: BBS 未名空间站 (Wed May 14 22:18:54 2008), 站内
: 发信人: jzxu (自然), 信区: Computation
: 标 题: 问个多参数fitting的问题
: 发信站: BBS 未名空间站 (Wed May 14 22:17:06 2008), 转信
: 公式很简单
: dM1/dt = -qM1-c1M1+p1c2M2
: dM2/dt = -qM2-c2M2+c1M1+c3M3
: dM3/dt = -qM3-c3M3+p2c2M2
: ...................




--

※ 来源:·WWW 未名空间站 海外: mitbbs.com 中国: mitbbs.cn·[FROM: 64.241.]

发信人: rossby (五十岚已夜), 信区: Physics
标 题: Re: 问个多参数fitting的问题 (转载)
发信站: BBS 未名空间站 (Fri May 16 14:19:55 2008), 转信

An old saying in physics is something like this

"Giving me two parameters, I can fit the world."


【 在 Tecnai (绿洲) 的大作中提到: 】
: 5个参数？
: 我记得费曼似乎可以用他们来拟合一头大象



--

※ 来源:·BBS 未名空间站 海外: mitbbs.com 中国: mitbbs.cn·[FROM: 141.213.]

发信人: Tecnai (绿洲), 信区: Physics
标 题: Re: 问个多参数fitting的问题 (转载)
发信站: BBS 未名空间站 (Fri May 16 14:42:36 2008), 转信

真的两个就可以了么
帮我fit一个吧

【 在 rossby (五十岚已夜) 的大作中提到: 】
: An old saying in physics is something like this
: "Giving me two parameters, I can fit the world."


--
于是,我似乎看到多年后,你站在一扇蓝色大门前……
下午三点的阳光,你仍然有几颗青春痘,你笑着,我跑向你问你好不好,你点点头
三年,五年,甚至更久以后,我们会变成什么样的大人呢?
虽然,我闭上眼睛也看不到自己,但是,我却可以看见你.



※ 来源:·BBS 未名空间站 海外: mitbbs.com 中国: mitbbs.cn·[FROM: 128.180.]

发信人: xiongxiong (熊熊), 信区: Physics
标 题: Re: 问个多参数fitting的问题 (转载)
发信站: BBS 未名空间站 (Fri May 16 16:04:54 2008)

用mathematica,

先
DSolve[{m1'[t] == - q m1[t] - c1 m1[t] + p1 c2 m2[t],
m2'[t] == -q m2[t] - c2 m2[t] + c1 m1[t] + c3 m3[t],
m3'[t] == -q m3[t] - c3 m3[t] + p2 m2[t]}, {m1[t], m2[t], m3[t]},
t]

然后用FindFit函数.

--

※ 来源:·WWW 未名空间站 海外: mitbbs.com 中国: mitbbs.cn·[FROM: 132.177.]

发信人: mifepristone (弃婴), 信区: Physics
标 题: Re: 问个多参数fitting的问题 (转载)
发信站: BBS 未名空间站 (Fri May 16 23:13:18 2008)

我做过，比你还复杂一点，大概是100个关联微分方程，我当时搞了无数方法
都不能有效收敛，最后用的是shooting method, 微分方程用runga-kutta积分，
在我那台1997年生产的超级垃圾计算机上大概5分钟就收敛了


【 在 jzxu (自然) 的大作中提到: 】
: 发信人: jzxu (自然), 信区: EE
: 标 题: 问个多参数fitting的问题
: 发信站: BBS 未名空间站 (Wed May 14 22:18:54 2008), 站内
: 发信人: jzxu (自然), 信区: Computation
: 标 题: 问个多参数fitting的问题
: 发信站: BBS 未名空间站 (Wed May 14 22:17:06 2008), 转信
: 公式很简单
: dM1/dt = -qM1-c1M1+p1c2M2
: dM2/dt = -qM2-c2M2+c1M1+c3M3
: dM3/dt = -qM3-c3M3+p2c2M2
: ...................



--

※ 来源:·WWW 未名空间站 海外: mitbbs.com 中国: mitbbs.cn·[FROM: 129.105.]

发信人: shadowman (hiahia), 信区: Physics
标 题: Re: 问个多参数fitting的问题 (转载)
发信站: BBS 未名空间站 (Fri May 16 23:47:20 2008), 转信

少了点，7个差不多
你记得也许是
“给我一根支点，我能撬动地球。”吧

【 在 rossby (五十岚已夜) 的大作中提到: 】
: An old saying in physics is something like this
: "Giving me two parameters, I can fit the world."



--

※ 来源:·BBS 未名空间站 海外: mitbbs.com 中国: mitbbs.cn·[FROM: 140.180.]

发信人: rossby (五十岚已夜), 信区: Physics
标 题: Re: 问个多参数fitting的问题 (转载)
发信站: BBS 未名空间站 (Fri May 16 23:57:05 2008), 转信


No, this was made by a famous physicist, once a bull in my field
cited it to mock a 2-parameter fitting talk.

【 在 shadowman (hiahia) 的大作中提到: 】
: 少了点，7个差不多
: 你记得也许是
: “给我一根支点，我能撬动地球。”吧



--

※ 来源:·BBS 未名空间站 海外: mitbbs.com 中国: mitbbs.cn·[FROM: 68.40.]

发信人: lynhard (寥若晨星), 信区: Physics
标 题: Re: 问个多参数fitting的问题 (转载)
发信站: BBS 未名空间站 (Sat May 17 00:26:27 2008), 转信

I think it was Landau who once said that 5 parameters could
fit an elephant. Something like that. Correct me if I am wrong.

【 在 rossby (五十岚已夜) 的大作中提到: 】
: No, this was made by a famous physicist, once a bull in my field
: cited it to mock a 2-parameter fitting talk.



--

※ 修改:·lynhard 于 May 17 00:27:28 修改本文·[FROM: 18.97.]
※ 来源:·BBS 未名空间站 海外: mitbbs.com 中国: mitbbs.cn·[FROM: 18.97.]

发信人: rossby (五十岚已夜), 信区: Physics
标 题: Re: 问个多参数fitting的问题 (转载)
发信站: BBS 未名空间站 (Sat May 17 01:20:27 2008), 转信


呵呵，google了一下，应该是Eugene Wigner说的，
不过到底说的是四个还是两个还真google不出来准确解，一个版本是
Give me four parameters and I can fit an elephant.
Give me five and I can make it wag its tail.

另外一个版本就是把四和五替换为二和三。


【 在 lynhard (寥若晨星) 的大作中提到: 】
: I think it was Landau who once said that 5 parameters could
: fit an elephant. Something like that. Correct me if I am wrong.



--

※ 来源:·BBS 未名空间站 海外: mitbbs.com 中国: mitbbs.cn·[FROM: 68.40.]

发信人: Microsystem (clam), 信区: Physics
标 题: Re: 问个多参数fitting的问题 (转载)
发信站: BBS 未名空间站 (Sat May 17 02:56:04 2008), 转信


呵呵，这些做理论的！实验中大部分的数据都是多参数拟合，甚至这个标准模型也是建
立在很多物理常量的误差范围内的。这个拟合大象的说法比较偏激。

【 在 rossby (五十岚已夜) 的大作中提到: 】
: 呵呵，google了一下，应该是Eugene Wigner说的，
: 不过到底说的是四个还是两个还真google不出来准确解，一个版本是
: Give me four parameters and I can fit an elephant.
: Give me five and I can make it wag its tail.
: 另外一个版本就是把四和五替换为二和三。


--

曾梦想仗剑走天涯
看一看世界的繁华



※ 来源:·BBS 未名空间站 海外: mitbbs.com 中国: mitbbs.cn·[FROM: 69.235.]

发信人: xiumu (朽木长风), 信区: Physics
标 题: Re: 问个多参数fitting的问题 (转载)
发信站: BBS 未名空间站 (Sat May 17 21:21:37 2008)

多参数也要看啥方程，要不我给你十个参数你给我拟合个大象来？
--

※ 来源:·WWW 未名空间站 海外: mitbbs.com 中国: mitbbs.cn·[FROM: 70.56.]