j**u
发帖数: 6059 | 1 【 以下文字转载自 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，但是总是搞不定，大家有什么好办法吗？ |
T****i
发帖数: 2190 | 2 5个参数？
我记得费曼似乎可以用他们来拟合一头大象
【在 j**u 的大作中提到】
: 【 以下文字转载自 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
|
d****e
发帖数: 251 | 3 我也想知道。我以前做的时候就是估计其中一两个参数的合理值，然后
基本上就是试出来的。不是很清楚有没有比较系统的方法。
【在 j**u 的大作中提到】
: 【 以下文字转载自 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
|
r****y
发帖数: 1437 | 4 An old saying in physics is something like this
"Giving me two parameters, I can fit the world."
【在 T****i 的大作中提到】
: 5个参数？
: 我记得费曼似乎可以用他们来拟合一头大象
|
T****i
发帖数: 2190 | 5 真的两个就可以了么
帮我fit一个吧
【在 r****y 的大作中提到】
: An old saying in physics is something like this
: "Giving me two parameters, I can fit the world."
:
|
x********g
发帖数: 595 | 6 用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函数. |
m**********e
发帖数: 12525 | 7 我做过，比你还复杂一点，大概是100个关联微分方程，我当时搞了无数方法
都不能有效收敛，最后用的是shooting method, 微分方程用runga-kutta积分，
在我那台1997年生产的超级垃圾计算机上大概5分钟就收敛了
【在 j**u 的大作中提到】
: 【 以下文字转载自 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
|
s*******n
发帖数: 1474 | 8 少了点，7个差不多
你记得也许是
“给我一根支点，我能撬动地球。”吧
【在 r****y 的大作中提到】
: An old saying in physics is something like this
: "Giving me two parameters, I can fit the world."
:
|
r****y
发帖数: 1437 | 9 No, this was made by a famous physicist, once a bull in my field
cited it to mock a 2-parameter fitting talk.
【在 s*******n 的大作中提到】
: 少了点，7个差不多
: 你记得也许是
: “给我一根支点，我能撬动地球。”吧
|
l*****d
发帖数: 754 | 10 I think it was Landau who once said that 5 parameters could
fit an elephant. Something like that. Correct me if I am wrong.
【在 r****y 的大作中提到】
: No, this was made by a famous physicist, once a bull in my field
: cited it to mock a 2-parameter fitting talk.
|
r****y
发帖数: 1437 | 11 呵呵，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.
另外一个版本就是把四和五替换为二和三。
【在 l*****d 的大作中提到】
: I think it was Landau who once said that 5 parameters could
: fit an elephant. Something like that. Correct me if I am wrong.
|
M*********m
发帖数: 2024 | 12 呵呵，这些做理论的！实验中大部分的数据都是多参数拟合，甚至这个标准模型也是建
立在很多物理常量的误差范围内的。这个拟合大象的说法比较偏激。
【在 r****y 的大作中提到】
: 呵呵，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.
: 另外一个版本就是把四和五替换为二和三。
:
|
x***u
发帖数: 6421 | 13 多参数也要看啥方程，要不我给你十个参数你给我拟合个大象来？ |
