爱吱声

标题: 问题：如何从数据里估算普瓦松分布的均值？ [打印本页]
作者: 晨枫    时间: 2019-2-4 11:38
标题: 问题：如何从数据里估算普瓦松分布的均值？
本帖最后由 晨枫 于 2019-2-3 22:42 编辑
$ h7 C; z5 `/ f" k0 I
6 w, C' h2 M( p8 Z9 v: X, I9 l一般估算均值是简单的算术平均，但我的问题比较特别，其实都不是统计问题。) i8 l$ u2 H* e7 R9 O
( G: A1 b1 w$ j3 w" q* X
是这样的，我有一个吸收塔，从塔顶到塔底有一个温度分布，形状大体像普瓦松分布：3 V, S  j5 J& o& p
6 J% u! {$ ~5 F% Y; `# c  j! j' G

; J- W% i, j7 l8 v1 V
3 O) e0 a! e' W" L: W我要估算峰值出现在哪块塔板，所以想到用统计的办法。实际上，统计里这是概率密度，曲线高低代表在这一点的采样“数量”，但我在每一个点只有一个数据点。说到底，我就是要拟合一个“钟形曲线”，然后找峰值所在点。相当于上图中红线（或者lambda=2）里横轴1-2之间的位置。+ |3 ]7 ]  l- h. l0 @9 K; `; N" Q

% B2 {) r) _6 J6 V3 U' I: L正态分布有现成的计算办法，但这样的“扭曲正态分布”或者普瓦松分布有什么简单办法吗？我需要能在DCS上实现，所以不能用太复杂的离线算法。
' {6 @4 a0 Y, f* \" m0 y+ {$ S$ R/ R$ b0 \) j6 C( W6 S+ M, S* [
爱坛里博士多如狗，教授满地走，想象起来，或许有谁在工作中碰到过类似的问题？
作者: 晨枫    时间: 2019-2-4 12:30
没人理我？都在忙着吃年夜饭？6 Y2 z9 p  C' d' v) Q
1 ?+ Q. z8 w& z; V
@煮酒正熟 @holycow @tanis @关中农民 @老马丁 @Dracula
作者: holycow    时间: 2019-2-4 12:43
泊松分布？我还是没理解你釆样拿到了哪些数据点
作者: 松叶牡丹    时间: 2019-2-4 12:45
没怎么看明白题目，您是说已知y值，要求对应正态分布？
作者: 晨枫    时间: 2019-2-4 12:46
holycow 发表于 2019-2-3 22:43
  D- ?& m/ X1 N3 K泊松分布？我还是没理解你釆样拿到了哪些数据点
* a# S/ t; S& H9 ^( `
是我描述得不好。再来一遍。
( K$ f8 p2 J' r% b0 A3 t3 e5 P2 M4 G* O, O
我有一条样子像泊松分布的温度分布曲线，但只有几个稀疏的点，想用类似泊松分布的曲线内插，然后计算峰值所在的横坐标点。
: r5 H( N+ u1 C5 c  A
: l) t0 k' ~2 Q+ E' I这要用于在线计算，不能搞非线性最小二乘之类的复杂算法。
4 @/ n) T5 Z1 h! M  z( J
7 @! |) L. l& U. V2 b4 m" @) U怎么弄？
作者: 晨枫    时间: 2019-2-4 12:47
本帖最后由 晨枫 于 2019-2-3 22:48 编辑
) e' r& h$ o& a
松叶牡丹 发表于 2019-2-3 22:453 p( ~# O. A" F( h& _; d& G/ @* }3 V
没怎么看明白题目，您是说已知y值，要求对应正态分布？

9 P1 i, s. B+ N6 @4 \2 y& x* L; U  t  d- {; L
抱歉抱歉，差不多就是这个意思，不过不是传统的对称正态分布，而是歪向一边的skewed normal distribution，请见5楼，不知道是否描述清楚一点？
作者: 晨枫    时间: 2019-2-4 12:50
松叶牡丹 发表于 2019-2-3 22:45
" \2 D/ _$ q# T5 T& N( |" F% [没怎么看明白题目，您是说已知y值，要求对应正态分布？

! M2 o+ c$ V6 l如果太麻烦，最后可能也就用对称的正态分布近似了。反正都是近似……
作者: 松叶牡丹    时间: 2019-2-4 12:52
晨枫 发表于 2019-2-4 12:50; r% Y2 R  r: |% b5 b1 y
如果太麻烦，最后可能也就用对称的正态分布近似了。反正都是近似…… ...
+ J7 Y2 T) }! g  V  p, G
如果精度不要求很高的话，不知道查表法的精度能不能达到要求。那个倒是最快。
作者: 煮酒正熟    时间: 2019-2-4 12:58
统计学我只明白最最基本的，泊松分布也理解不深。问几个傻问题吧：一共多少块塔板？为什么每块塔板只能取到一个data point? 我猜是取样也有成本问题？那样儿的话，您知不知道每块板上的随机取样，与这块板的温度峰值是不是最接近、比其他板的随机取样更接近本板的峰值？
作者: 晨枫    时间: 2019-2-4 12:59
松叶牡丹 发表于 2019-2-3 22:52! {4 A% V  F6 L' W" b
如果精度不要求很高的话，不知道查表法的精度能不能达到要求。那个倒是最快。 ...
2 @% ]* ?7 \) q0 n; f
7 [) k' V6 @6 V/ l3 T! {3 t% q
查表再内插，这要做成在线计算还是有点麻烦了？
作者: tanis    时间: 2019-2-4 13:00
晨枫 发表于 2019-2-4 12:46
3 n/ ~; k" h- t. o2 W. d& H是我描述得不好。再来一遍。1 b" _, Y5 S! U4 @/ C: ~$ I

' U( b: g2 ^  T# a我有一条样子像泊松分布的温度分布曲线，但只有几个稀疏的点，想用类似泊松 ...
' p6 U. G! y' t5 B& _# ^/ c" [
原来是泊松分布。。。 惭愧，连最小二乘法都属于复杂算法了。。。还真不知道有什么其他的简便方法。这一块早已经还给老师。。。 现在计算能力强大，集成性高，一般直接一个pasfit之类的函数就解决问题了。。。哈哈哈
作者: 松叶牡丹    时间: 2019-2-4 13:00
本帖最后由 松叶牡丹 于 2019-2-4 13:06 编辑 * g! P+ K/ C$ v3 y8 s. t4 O0 I

+ n1 F9 |8 c) O: `* L0 A0 v做成在线计算器，那就不能查表了，反而太麻烦。坐等观摩高手们过招。
作者: 晨枫    时间: 2019-2-4 13:03
煮酒正熟 发表于 2019-2-3 22:587 w# {, `! A, H
统计学我只明白最最基本的，泊松分布也理解不深。问几个傻问题吧：一共多少块塔板？为什么每块塔板只能取到 ...

+ z# y, m1 S1 K咳咳，这个其实不是统计问题，是从有限的温度测量估计温度分布曲线的问题。吸收塔一共20块塔板，每块塔板只有一个温度测量点，连续测量（每15秒采样一次）。
, W0 x7 Q3 j( x5 F: n
+ L& d( Q9 Y& c6 l1 [1 D! [这温度分布随工艺条件而变，所以这个峰值和凸起一直在动来动去。目视很容易看出峰值在那里，但要用数值办法实时确定，反而犯难了。
作者: 晨枫    时间: 2019-2-4 13:06
tanis 发表于 2019-2-3 23:00  E* b/ V/ a) E* R1 `4 s; q7 @' N
原来是泊松分布。。。 惭愧，连最小二乘法都属于复杂算法了。。。还真不知道有什么其他的简便方法。这一 ...
7 A% r' C4 i6 G& v6 t
切，我有那东西，还费这个事吗？Matlab、Maple，什么不能上啊！问题是我需要在线计算，而且只有FORTRAN一级的语言从头现编，还不能太耗费CPU时间，耽误了实时控制就偷鸡不成蚀把米了。
作者: 晨枫    时间: 2019-2-4 13:06
松叶牡丹 发表于 2019-2-3 23:006 E  \5 \+ W4 P& X( N% |9 w; Z
做成在线计算器，那就不能查表了，反而太麻烦。坐等观摩高手们过招。
4 B4 D0 b1 X$ X6 ^5 l8 w6 o
估计吃完年夜饭，就该忙着数红包了。
作者: tanis    时间: 2019-2-4 13:09
晨枫 发表于 2019-2-4 13:030 |# X" {  ^6 J: k8 a
咳咳，这个其实不是统计问题，是从有限的温度测量估计温度分布曲线的问题。吸收塔一共20块塔板，每块塔板 ...
; U3 D% M6 V( f/ w* h
目视是指看数据点脑补么？ 所以一共20个数据点（塔板），然后最高温度可能在点之间？ 15s一次，速度很慢啊，为啥不能离线算？ 系统太古老了？
5 \. N4 h/ w# C- a# f6 B
4 G% N+ P! `$ T7 {: u如果是系统太老无法整合还有个鬼畜方法。。。 既然肉眼能看到每个点在哪里，可以用电脑拍照把图截下来，然后从图像里把点的值拿到，接着在电脑里拟合。
作者: tanis    时间: 2019-2-4 13:12
晨枫 发表于 2019-2-4 13:06
8 Z6 l% m, b- N% u& y/ S; a切，我有那东西，还费这个事吗？Matlab、Maple，什么不能上啊！问题是我需要在线计算，而且只有FORTRAN一 ...
, k$ d- k) A: X9 D3 B5 N
找温度第二，第三高的两个点，拟合个直线，取中间的那个位置
作者: 松叶牡丹    时间: 2019-2-4 13:12
晨枫 发表于 2019-2-4 13:06$ `- ^+ p' t# `9 O9 d, O
切，我有那东西，还费这个事吗？Matlab、Maple，什么不能上啊！问题是我需要在线计算，而且只有FORTRAN一 ...
2 V, ~9 p- R8 W$ c, \6 B
您是说每块塔板都是15秒测量一次，所以每15秒要同时给出20个塔板的20个正态分布式子？
作者: 晨枫    时间: 2019-2-4 13:18
tanis 发表于 2019-2-3 23:098 P. K2 A1 ?4 M& U
目视是指看数据点脑补么？ 所以一共20个数据点（塔板），然后最高温度可能在点之间？ 15s一次，速度很慢 ...
3 r( ~% r0 w  o5 M% ^/ m1 t$ @/ E" D
正是：看数据点然后脑补。人工就是这么干的。但操作工要照管的事情太多，一会儿不看，可能就错过了，就可能坏事。! |9 Z3 k$ P: L

6 G6 {- N- |/ g+ G2 O化工上DCS的采样分基本回路和先进回路，前者一般每秒一次，甚至1/4秒一次，后者每15、30、60秒一次就够了。过程时间常数动辄几十分钟，太快的采样没必要。另一方面，CPU的速度（还是Morotola 68040）和网络带宽（4 megabits 烹per second）有限，不宜太快，还要扫描几百、几千个回路，不宜太快。安全连锁回路更快，在毫秒级，但那不用作连续控制，两回事。
8 Q) ]  y2 B4 H0 I; t+ F0 Y% {8 D/ L& l
工控级的CPU和网络都用特别成熟（翻译：特别古老）的技术，和实验室、离线是两回事。% `! `2 O3 A' u$ e6 Q

+ f" C; `/ |& q, S在线计算是因为我要把这个峰值点用到控制回路里。拍照、图像识别、离线拟合，这绕的弯子太大，中间节点太多，可靠性根本靠不住，哪里打一个嗝，就可能造成process upset，就是几十万、上百万的损失，那就划不来了。
作者: 晨枫    时间: 2019-2-4 13:21
松叶牡丹 发表于 2019-2-3 23:12  ~7 ]' T' ~* c; \: r
您是说每块塔板都是15秒测量一次，所以每15秒要同时给出20个塔板的20个正态分布式子？ ...
# M1 T6 w4 f& G) r
每15秒要对20块塔板的温度拟合一条样子像正态分布的曲线。15秒只是计算间隔，每次实际计算其实只有毫秒级的窗口，因为同一个CPU还要分时管好多别的事呢。
作者: 松叶牡丹    时间: 2019-2-4 13:21
本帖最后由 松叶牡丹 于 2019-2-4 13:23 编辑 0 H% k- T# {1 k! N2 `% X2 H
晨枫 发表于 2019-2-4 13:185 |2 H7 A+ s8 k9 ~6 }$ K
正是：看数据点然后脑补。人工就是这么干的。但操作工要照管的事情太多，一会儿不看，可能就错过了，就可 ...

+ P# [* R1 J5 V' b5 s( B, v& ?; O3 [" l. S
降低精度，把一系列正态分布离散化为有限几个histogram，每次测试点先第一步进行范围判断，然后归到不同的histogram里面。再反向把histogram变换为对应的大致正态分布，这样可行吗？
作者: 晨枫    时间: 2019-2-4 13:23
tanis 发表于 2019-2-3 23:12
, o  T. B& b7 ~找温度第二，第三高的两个点，拟合个直线，取中间的那个位置

- j: H! \5 T: W$ ^& R' n5 b塔板是等间距分布的，这样我还不如直接取温度最高的那一块了
作者: tanis    时间: 2019-2-4 13:27
晨枫 发表于 2019-2-4 13:23/ \5 h) B( D# Y5 u, g
塔板是等间距分布的，这样我还不如直接取温度最高的那一块了
) B! a+ g$ O) S0 I
哦哦。。。 那就取最高和第二个的中间位置 有人工干预的时候就用人工，人忘了看就用这个。。。
作者: 晨枫    时间: 2019-2-4 13:28
松叶牡丹 发表于 2019-2-3 23:21
; z* \% x- [3 @' H- w% C& S! w降低精度，把一系列正态分布离散化为有限几个histogram，每次测试点先第一步进行范围判断，然后归到不同 ...
1 Z" K$ i/ T+ R. ?/ `( b7 ~* b
5 {, y) k& F8 s% z( Q6 X& F3 f
一把抓住：我的问题正是从histogram变换为某种分布，最想poison分布或者log normal distribution，不行的话正态分布也行。有什么现成的计算方法吗?
作者: 晨枫    时间: 2019-2-4 13:29
tanis 发表于 2019-2-3 23:27
6 C2 v1 x# e7 F" f" h哦哦。。。 那就取最高和第二个的中间位置 有人工干预的时候就用人工，人忘了看就用这个。。。 ...

" X" P/ h" u. {& T5 n  F这也是个办法，总是有点“心有不甘”……
作者: 晨枫    时间: 2019-2-4 13:33
松叶牡丹 发表于 2019-2-3 23:21
# x" w  A8 R8 ]$ W降低精度，把一系列正态分布离散化为有限几个histogram，每次测试点先第一步进行范围判断，然后归到不同 ...

4 V: @: j! O& y- h: O+ F6 e, J, G9 r5 z唉，MATLAB里有histfit命令，干的正是我要的，可惜没法“偷”过来用啊
作者: 松叶牡丹    时间: 2019-2-4 13:36
晨枫 发表于 2019-2-4 13:33+ b# K, p4 K/ f
唉，MATLAB里有histfit命令，干的正是我要的，可惜没法“偷”过来用啊

: h, N% |2 `- m9 O2 T" c' s2 `& {* ^1 k' ]完全不会FORTRAN，这样啊，要不就得自己写算法了？
作者: tanis    时间: 2019-2-4 13:37
晨枫 发表于 2019-2-4 13:33
( \$ ]6 y% k6 D4 }唉，MATLAB里有histfit命令，干的正是我要的，可惜没法“偷”过来用啊

/ E  g" ?2 N: p7 A/ |8 ^- jMatlab 直接有poissfit啊。。。
3 ^' x+ w" L7 Z& I1 w2 a4 o: F! O( B( a1 N9 G
我感觉人类脑补的水平估计和取中值也差不了多少了既然连忘了看也只是可能出问题，我脚的中值差不多 ）（缺乏工匠精神的我。。。）9 Y( r. ^, P! S
作者: 晨枫    时间: 2019-2-4 13:40
tanis 发表于 2019-2-3 23:37
, m! @  w  q1 S3 \$ i! l$ t  s6 A# ]9 DMatlab 直接有poissfit啊。。。# p' L/ ?  b% r4 d* N
' a5 x5 U+ ]. W# F1 I& H% W
我感觉人类脑补的水平估计和取中值也差不了多少了既然连忘了看也 ...
. r# t8 C- P# N# Z6 Z& A
问题是我没法把MATLAB装进DCS里用啊。嗯，如果不行的话，可能也只能内插取中值了。
作者: 松叶牡丹    时间: 2019-2-4 13:41
您这个采样数据如果有之前很长时间的大批量数据的话，要不选一部分建模，建立比如10个distribution式子，用剩下数据validate，看看是不是符合六西格玛质量控制要求？如果不行就提高精度，设20，30或者更多式子？我这是笨办法。
作者: 晨枫    时间: 2019-2-4 13:43
松叶牡丹 发表于 2019-2-3 23:41
7 T, k1 S, y/ K您这个采样数据如果有之前很长时间的大批量数据的话，要不选一部分建模，建立比如10个distribution式子，用 ...

& T" _$ `) r3 `+ }$ P这个好像不行？我需要的是实时估计当前的峰值点，历史数据能帮什么忙呢？
作者: gnomegordon    时间: 2019-2-4 13:55
如果问题是拟合若干点的光滑曲线 找到最大值，kernel density 就好，matlab里的histfit 也是基于这个算法。fortran里面library浩若烟海，查查肯定能找到。速度绝对够快，不需要POISSON 假设。8 O' p7 F3 D# C5 e. ?( s& J3 o
7 U7 D0 X8 a$ p: N
试着搜一下，这个fortran code 试试？： https://ideas.repec.org/c/boc/bocode/g00004.html
作者: 松叶牡丹    时间: 2019-2-4 14:07
晨枫 发表于 2019-2-4 13:434 y  a2 ~" E- P! \' j4 \+ ^
这个好像不行？我需要的是实时估计当前的峰值点，历史数据能帮什么忙呢？ ...
5 x! ~  m. h+ `( N" |: ?
历史数据建模自己写算法做验证用。不过如果可以像楼上那位说的可以找到现成library调用，kernel density estimation和 guassian mixture model应该都可以。
作者: 晨枫    时间: 2019-2-4 14:16
gnomegordon 发表于 2019-2-3 23:55- K* K* U) f! @. v
如果问题是拟合若干点的光滑曲线 找到最大值，kernel density 就好，matlab里的histfit 也是基于这个算法。 ...
. c7 _/ t3 d4 p' J3 O8 H
十分感谢！打开看到的是SRC文件，还在找FORTRAN……
作者: 晨枫    时间: 2019-2-4 14:19
松叶牡丹 发表于 2019-2-4 00:07: G. H# y" Z, _4 {- m
历史数据建模自己写算法做验证用。不过如果可以像楼上那位说的可以找到现成library调用，kernel density  ...
6 t% R% I6 l4 i
多谢！这个kernel density以前没有玩过，正在研究中……
作者: Dracula    时间: 2019-2-4 14:24
本帖最后由 Dracula 于 2019-2-4 14:28 编辑
4 R; Z) C0 q, j% p
4 \8 x0 v6 Z4 M2 S! u2 r; f说实话，你这个题目我没看懂。
) e- D& g& F, a9 C: G" Z: I0 G5 H) u; u
但如果你你单纯就是想实现
: x8 a2 B% L6 N1 \/ y3 @, V
6 C, _; \4 C" b
我就是要拟合一个“钟形曲线”，然后找峰值所在点。相当于上图中红线（或者lambda=2）里横轴1-2之间的位置。
( A; ?* p, x3 u5 Q; o9 i  B  q
2 E! k" p3 p5 s3 I) s( N3 N3 l
其实应该挺容易。假设这个分布是log normal，如果你知道这个分布的参数， μ和 σ，你要想计算峰值（英语是mode），有现成公式，套公式就可以了（公式你可以上网查，比如Wikipedia上就有）。估计这个分布怎么做，我没看懂你的题目，不好说。但既然正态分布你知道能怎么做，log-normal没什么本质区别，就是公式不一样。简单的使用method of moments应该很容易就能推导出公式。实现就是简单的算术计算，不需要什么数值方法。$ ~; L, ]4 d& ~/ Y# Z1 n: T

( n9 C+ M5 O8 G6 p  G
作者: 晨枫    时间: 2019-2-4 14:34
Dracula 发表于 2019-2-4 00:24
0 q1 m/ y5 C% E0 T7 ^) d说实话，你这个题目我没看懂。5 O4 f4 v: W+ u8 O$ {) J: b

0 k; Z; s/ l' Q5 ]但如果你你单纯就是想实现
8 d6 g! o5 x3 t; O4 t
有了 μ和 σ想计算峰值就容易了，我的问题是如何从histogram计算log normal的 μ和 σ。看来这也是个办法，把数据log一下，再按normal做。我也想到过这个办法，就是不大确定。明天试试，上班的活儿用下班的时间试，可不划不来嘛。
作者: 松叶牡丹    时间: 2019-2-4 14:35
Dracula 发表于 2019-2-4 14:24+ @1 w) ?9 F. k/ p
说实话，你这个题目我没看懂。
* f( A, s. m4 b/ V: E( U4 H8 Y4 Y9 O6 a. f% ~9 P
但如果你你单纯就是想实现

4 Y: g8 r; J  o& F看题目意思应该是不知道参数的。
作者: 松叶牡丹    时间: 2019-2-4 14:36
晨枫 发表于 2019-2-4 14:19
) u4 l+ y& q5 e% r1 \3 i多谢！这个kernel density以前没有玩过，正在研究中……
$ s. g9 J% a& Q2 H" ^+ w
晨大辛苦，您太客气了。祝新年快乐！
作者: gnomegordon    时间: 2019-2-4 14:39
晨枫 发表于 2019-2-4 14:16
5 ^8 J: _2 x. c: _& w十分感谢！打开看到的是SRC文件，还在找FORTRAN……

" D9 C9 y9 B. d0 U5 v* q! Tapologize. 网上搜code太麻烦，还得验证。最好有本书可以翻翻 或者搜library
作者: Dracula    时间: 2019-2-4 14:40
本帖最后由 Dracula 于 2019-2-4 14:45 编辑
, f6 I; C, r0 ?4 m  R+ s8 y
晨枫 发表于 2019-2-4 14:343 `' S! |- i' g: X/ R- E
有了 μ和 σ想计算峰值就容易了，我的问题是如何从histogram计算log normal的 μ和 σ。看来这也是个办 ...

+ {3 s' k8 ^7 _& h9 b9 a$ u9 [
, D+ K% z$ n$ T) n+ d( G! J9 q怎么计算分布参数的问题，你的题目我没看明白，不好说，但是正态分布你会做，log-normal 没有任何本质区别，一样的办法，就是数学公式不一样就是的了。应该不难。/ M1 \# ]7 }# r, I8 f

7 `+ Z5 l. \* \& `% V6 Y- R* V5 c; u, j1 `. O' |7 p. r
（标准的统计学问题，估计log-normal分布参数也是有公式的，你到网上去查个公式就是的了。）3 H+ q) J% z/ g7 P- T$ V2 \
0 i. y+ r% m  b8 b% W& G
作者: 数值分析    时间: 2019-2-4 14:56
本帖最后由 数值分析 于 2019-2-4 15:05 编辑 ) b% w" k& @. I$ e* G

) _0 e# h9 q0 U% ?, ]( O6 M你应该不用拟合分布函数吧？你只想知道峰值的位置，然后你又知道（或者说你假设）是泊松分布，所以峰值的位置一定是 x=lambda,（这里lambda不一定是整数），那么剩下的就是从样本里推断lambda了，这是个典型的估计啊. 对于泊松分布，lambda正好是期望，所以一般来用样本均值估计期望。, ~0 i8 X# G8 H9 J
你给每一个板子从最左边顺序编个号，i=0，1，2，3。。。，然后设每块板子i的对应温度样本值xi，，然后计算sum（i*xi）/n [即累加所有的（板号乘以对应温度）然后除以板数】 （因为你的分布曲线可能和泊松分布差一个常数，所以最后结果得scale一下）不就可以了么？当然，这得假设你的histogram真的得长得像泊松分布分布。
5 |" }1 e5 k/ H2 x/ ?
作者: 视觉错误    时间: 2019-2-4 17:40
石化行业的DCS应该Honeywell多吧，这年头不支持OPC的很少了。
9 j4 J* n  C) B: [% j7 T- v9 U数据送到电脑上算，算了以后在送回去。
作者: 视觉错误    时间: 2019-2-4 18:02
晨枫 发表于 2019-2-4 13:33
4 }# i! T2 [1 l0 ^% i) _" j唉，MATLAB里有histfit命令，干的正是我要的，可惜没法“偷”过来用啊
, S4 ^1 s4 z' z. q# J: G' e0 \
我记得MATLAB支持OPC
作者: 数值分析    时间: 2019-2-4 18:34
数值分析 发表于 2019-2-4 14:56
6 j" d& n" F, ^" A$ ?; u9 o! H5 `你应该不用拟合分布函数吧？你只想知道峰值的位置，然后你又知道（或者说你假设）是泊松分布，所以峰值的位 ...
* I1 {1 _2 J; C& O* z+ C7 E; ^; p
多解释一句scale那块儿。因为泊松分布曲线下面的面积是1，而你的histogram显然不是，所以你的histogram和泊松分布差一个常数。你求出来的lambda的估计要用你histogram的面积归一一下。
作者: 小刀    时间: 2019-2-4 20:39
可以试试GMM Guassian Mixed Model去拟合统计分布8 U6 |9 N; N1 L! u" w& x
作者: 关中农民    时间: 2019-2-4 21:47
晨枫 发表于 2019-2-4 12:30. P0 J" D  l1 k7 O
没人理我？都在忙着吃年夜饭？1 \' V9 j/ J/ F! ^

7 s' U  |3 Y  ?* p3 q- k% `- D, C@煮酒正熟 @holycow @tanis @关中农民 @老马丁 @Dracula  ...

8 n6 i2 K9 K0 H9 s3 V& ~9 ]晨大，这得数学博士才中啊，额完全外行了，看见这个只能联想到面条
作者: Dracula    时间: 2019-2-4 22:07
我又看了一下你这个题，终于看明白了。你的问题是一条曲线类似于统计学上Poisson或者log-normal的density function的形状，在这条曲线上你已知坐标是x=1，x=2，... x=20，这20个点的值，你想要知道的是曲线的最大值是在那个点上。不知道我这个理解对不对。6 T& s' m9 a: T/ i) t  W6 Z

$ C1 V3 P7 H" O, f" e" h如果我的理解是对的话，这不是个统计学问题。你画的那个也不是histogram，因为histogram的纵坐标是在每个值观测到的sample size，而你的图的纵坐标是温度，不是一回事。因此统计学的书你不用查，查了也没用。解决这个问题最显而易见的办法就是最小二乘法，但应该是没有分析解，你不能用。我好奇的是如果假设假设曲线的形状类似于正态分布的density function，你们是怎么解的，使用最小二乘法应该是一样没有分析解。如果解正态分布有特别的巧妙的办法的话，或许稍微修改一下就可以用到log-normal的情况。" k5 d* s! S0 ?1 k  T2 |
0 Q( U' |- I+ x* ^8 G
9 j3 q+ H/ b" X  S3 ~

2 f$ u- O* \) K3 _7 N7 D  P( Z( f( L6 p) s4 d7 t* u+ x
作者: 晨枫    时间: 2019-2-4 22:17
数值分析 发表于 2019-2-4 00:566 @: N4 K/ b1 r( }
你应该不用拟合分布函数吧？你只想知道峰值的位置，然后你又知道（或者说你假设）是泊松分布，所以峰值的位 ...
& C8 S9 I. V* H
这个办法好！回头试一下！我是打算用这个办法当正态分布处理的，没想到也可以相当直接地套到泊松分布。可能这就解决我的问题了！多谢！
作者: 晨枫    时间: 2019-2-4 22:20
Dracula 发表于 2019-2-4 08:07% A3 b3 |! ]2 O6 T
我又看了一下你这个题，终于看明白了。你的问题是一条曲线类似于统计学上Poisson或者log-normal的density f ...

% [7 @6 l! s' \2 X8 i3 p' z对，就是这个意思。我也提到了，不是统计问题，只是“形似”，想看看统计里有没有现成的办法。楼上42楼就是我一开始想到的办法，但只想到那能用于正态分布，正想改造为对数正态，没想到可以直接套泊松。这就好了。
作者: 晨枫    时间: 2019-2-4 22:21
小刀 发表于 2019-2-4 06:395 _/ H0 z4 @$ m% v3 l2 D+ H
可以试试GMM Guassian Mixed Model去拟合统计分布

( Y7 d# Y( ]& D" ?5 {# _这个还是太复杂了。用在控制回路里，必须KISS。但还是要谢一个！
作者: 晨枫    时间: 2019-2-4 22:23
视觉错误 发表于 2019-2-4 03:40# B! k! v) O/ s" F) ]
石化行业的DCS应该Honeywell多吧，这年头不支持OPC的很少了。) V% X" d# N* v* ?% K0 q7 b3 V  k
数据送到电脑上算，算了以后在送回去。 ...
4 b7 e; J, i) o  X7 g8 C3 P4 s) c+ K
我们有OPC，问题是可靠性。用以下层基本的回路控制一般不用OPC，当机或者“交通堵塞”的后果太大。这是惯例。只有上层的APC可以用OPC，当了就自动shed到基本控制。
作者: 晨枫    时间: 2019-2-4 22:25
视觉错误 发表于 2019-2-4 04:022 U4 h+ p, Q8 j/ J/ l
我记得MATLAB支持OPC
. [4 I" ?1 T' ^7 _8 ?/ S
是的，我以前还试过用MATLAB C通过OPC与DCS相连，在技术上这是做得到的，但可靠性达不到要求。OPC是不作为可靠的控制信息通道使用的，只能传送点监视数据或者一般数据采集。
作者: 晨枫    时间: 2019-2-4 22:26
gnomegordon 发表于 2019-2-4 00:396 ], W' U- S$ J4 {8 i5 c
apologize. 网上搜code太麻烦，还得验证。最好有本书可以翻翻 或者搜library
! G1 U* d( f; y( n; o
再次感谢。楼下45楼有好办法，我先试试那个办法，比kernel density简单多了。
作者: 晨枫    时间: 2019-2-4 22:26
松叶牡丹 发表于 2019-2-4 00:360 M2 p6 b+ i9 }8 m% M9 {) A
晨大辛苦，您太客气了。祝新年快乐！

/ y  v7 b( r$ ^9 K  q松叶MM新年快乐！
作者: 雨楼    时间: 2019-2-4 22:31
晨枫 发表于 2019-2-3 23:46
3 ^$ Z8 b+ s7 Z) ]是我描述得不好。再来一遍。
7 h) i( W: |- E& u8 n8 ~* L: h8 c! p
- b7 I- X; D( N5 D3 R- F( q5 q我有一条样子像泊松分布的温度分布曲线，但只有几个稀疏的点，想用类似泊松 ...
& P- u0 ^; w$ J2 Q3 T1 @7 i! O
就是正态分布然后在x轴上平移么？ 类似Y=(X-a)^2.
作者: Dracula    时间: 2019-2-4 22:33
晨枫 发表于 2019-2-4 22:20
; J+ w' g: F2 F, |3 {+ V对，就是这个意思。我也提到了，不是统计问题，只是“形似”，想看看统计里有没有现成的办法。楼上42楼就 ...
2 ~& `( H% K/ X6 z
42楼那个办法不对。那是把这当成个统计学的问题来处理，但这不是个统计学问题。你的纵坐标是温度，不是sample size，不能这么用。最明显的是，那个办法解出来的量纲是温度，而你想要的应该是具体是那块板，因此那个解和你想要的没什么关系。
4 x6 q$ ]  v0 t! P7 O7 j# {$ p
0 a+ `: Z* y& D
作者: 晨枫    时间: 2019-2-4 22:35
数值分析 发表于 2019-2-4 04:34
1 I) Y2 t# o$ g0 f多解释一句scale那块儿。因为泊松分布曲线下面的面积是1，而你的histogram显然不是，所以你的histogram和 ...
6 d' n3 _7 K/ M0 v1 t, x9 D! M( N3 z
多谢！记住了！
  Q9 D- C3 X0 M+ M8 q
8 B  P) j; s/ P' T. Y; t其实你说的办法我已经试过。我把正态分布一边的尾巴砍掉，至少外观上接近泊松或者对数正态。只要有峰在，估计出来均值就还不错，越对称越准确，就有点窃喜。但对道理不摸底，不敢放手用。除非在数学上站得住脚，否则在线的时候没人看着，给我乱估一个就完蛋了。现在看来，道理就是你说的，这个办法不只适用于正态分布。曲线只有只有半边的话，就有点悬，这个可以理解。一般到不了这个情况，程序里简单判别一下也不难，另作处理。
作者: 雨楼    时间: 2019-2-4 22:37
晨枫 发表于 2019-2-4 00:03
9 m9 t2 U0 }5 o& T9 c2 |7 s$ s9 n咳咳，这个其实不是统计问题，是从有限的温度测量估计温度分布曲线的问题。吸收塔一共20块塔板，每块塔板 ...

/ Z) n' s, R) G) |, Y2 }5 n0 ~1. 20个数据点在分布上有没有规律。比如两头低中间高。" e+ y3 f( j0 v3 S4 i
2。规律稳定么？% g, f/ T7 R+ w: q' O
3。可不可以简化成20个点里找最大值。" a4 J/ }0 O+ Q: p
4。峰值如果不在采样点（塔板），而在塔板之间，只能按相邻塔板的问题计算温度曲线斜率，然后插值，而其要比较峰值塔板两侧的斜率，取较大的。
作者: 晨枫    时间: 2019-2-4 22:38
雨楼 发表于 2019-2-4 08:31+ o# i6 e3 y3 r/ O  V7 F5 Z2 a
就是正态分布然后在x轴上平移么？ 类似Y=(X-a)^2.
# U% [" M7 G2 V8 g1 X, [3 b
差不多。我开始也想过用抛物线然后平移，但平移量本身也要最小二乘出来。可能还是可以线性化然后用简单的最小二乘。我来试试看。
作者: 晨枫    时间: 2019-2-4 22:38
Dracula 发表于 2019-2-4 08:33& f3 S9 C( h' O' k
42楼那个办法不对。那是把这当成个统计学的问题来处理，但这不是个统计学问题。你的纵坐标是温度，不是sa ...
" N- S, R( V0 W9 r7 J
嗯，我再想想。谢谢。
作者: 晨枫    时间: 2019-2-4 22:42
雨楼 发表于 2019-2-4 08:37" A* M' C7 e% v4 t( ]# h% ]8 w
1. 20个数据点在分布上有没有规律。比如两头低中间高。
7 ^! Q/ p* [7 c  n$ K2。规律稳定么？8 q; L$ ^/ A& U0 c( t% Q! w
3。可不可以简化成20个点里找最 ...

' t- ]+ V3 W1 D2 h1 E* u两头低，中间高，但峰值形状随工艺条件而变，可以从“一头歪”的泊松变到对称的正态，然后继续变到往另一头歪的泊松。选最大值然后插值也是可以的，我就是在想是不是有可以一次性计算出来的，而不需要这样搜索。+ G4 F  F! J& y/ \; @& v

' v7 |4 _3 C. J# K) I' f如果42楼的办法最终不行，我可能回过来用你的办法。谢谢！
作者: 视觉错误    时间: 2019-2-4 22:42
晨枫 发表于 2019-2-4 22:25
8 p& t, v! D* \( z; @  ~' ~是的，我以前还试过用MATLAB C通过OPC与DCS相连，在技术上这是做得到的，但可靠性达不到要求。OPC是不作 ...
# T9 j& s6 ^4 u; ]
这个峰值位置还需要用于控制吗？7 F7 @4 @' T) ~8 Z7 o( ?
我理解这个峰值位置计算出来也就是用于参考吧。
作者: 晨枫    时间: 2019-2-4 22:44
视觉错误 发表于 2019-2-4 08:42
% C. m* D0 w: g+ X# p5 D这个峰值位置还需要用于控制吗？
9 {# v, F3 X( o5 b我理解这个峰值位置计算出来也就是用于参考吧。 ...

, Z1 q& U7 c7 F0 o! i" [& f7 w% L我就是想用于控制。这对应于吸收塔里放热反应的热点位置，决定了吸收效率。跑得太偏了，要么浪费能源，要么吸收不达标。传统的单点温度控制效果有限，很容易被上下移动的峰值位置所“误导”。
作者: 视觉错误    时间: 2019-2-4 22:45
同意伯爵的看法，本质上是个曲线拟合。7 q* s( i7 u" v; q# D. t) P% F
这个曲线有点像某种分布曲线，晨大强调这个造成误解了。
作者: 晨枫    时间: 2019-2-4 22:46
视觉错误 发表于 2019-2-4 08:45
' I; k6 w1 X. W' \, y' x同意伯爵的看法，本质上是个曲线拟合。
  w5 Z$ E# Q) y+ V7 H& b9 x; ?2 U这个曲线有点像某种分布曲线，晨大强调这个造成误解了。 ...

4 r6 b" A* t' E+ P3 U1 p0 n8 M8 ]抱歉！我还特意提一句，这不是统计问题，但还是误导了。多谢各位指点。我这会儿有很多办法可以试试了！
作者: Dracula    时间: 2019-2-4 22:53
晨枫 发表于 2019-2-4 22:38% h4 W+ ?$ t8 z  q1 a0 ?
嗯，我再想想。谢谢。
9 K+ x9 |& g0 S1 K* P2 s2 ], d$ D
我曾经想过一个和42楼类似的办法。区别是分母不是板数，而是各个板加在一起的温度的和。如果这条曲线真的是个density function，也就是说，曲线下的面积等于1的话，这个办法是可以的。这个办法就是把温度类比成统计学里的sample size来处理了。但是曲线下的面积等于1这个假设肯定不对。我曾经想过把它scale成1就是的了。但问题是你的纵坐标0度的设置应该纯粹是arbitrary的，但是这个0度的设置会直接影响scaling以及最后结果，因此这个办法也不行。( z  k+ h! F# {. G0 r

) `- _4 T/ |* \. j' D如果想不出数学上的分析解的话，我的建议是你想一想人的直觉是怎么来处理这个问题的。你前面说，人眼一下就能看出来。把人直觉的逻辑想清楚的话，写个if then else的程序实现应该不算很难。
+ i& C& l$ T+ |( N6 u8 S* T. Q' C) e! V' U
作者: 雨楼    时间: 2019-2-4 23:49
１，　计算最大值，% x" `* T6 ^8 D
２，　计算左右的斜率
3 C. Q$ [6 p" m$ v, B３，　计算左右直线的交点。
0 J- x. Z$ `/ O9 X+ B9 ^1 M４，　由交点根据斜率反推峰值的大概值。这个是quick&dirty的解决办法。想要准确，还得曲线拟合。２阶应该就可以了。过高反而会引入误差。但是拟合的资源消耗可能系统付不起。/ `- s9 Y, ~' n9 `# w

# k2 Y; R$ S: _" m. K4 N& m* Thttp://www.aswetalk.net/bbs/home ... um&picid=102923, Y3 L" `2 X% k, S3 A) [
作者: 晨枫    时间: 2019-2-5 00:01
Dracula 发表于 2019-2-4 08:53
" [& Y# ^2 p9 T1 u2 K我曾经想过一个和42楼类似的办法。区别是分母不是板数，而是各个板加在一起的温度的和。如果这条曲线真的 ...

' x  v! ^1 G( i4 D/ _多谢。接下来我会按42楼办法用实际数据多试试，包括scaling，然后和目视结果比较。如果可靠的话，就能用了。否则还要另想办法，如楼下68楼的。) l# E# \  I5 C7 I7 Y& l

, N/ T2 h& t4 @9 q7 B模拟人类思维的办法想过，不大好弄。太复杂的IF...THEN容易把自己绕进去。工程上还是要KISS。
" ^: _8 e5 Q3 T) B+ p) r6 ]0 b/ o1 _5 E5 J
多谢伯爵帮我想这个问题！
作者: 晨枫    时间: 2019-2-5 00:02
雨楼 发表于 2019-2-4 09:49! w5 i  D2 j) C: Z$ D
１，　计算最大值，, [: a* N  Q) x) o" ~9 C6 m6 _9 A% _! Q
２，　计算左右的斜率) A  j0 z5 O0 J- z8 z
３，　计算左右直线的交点。
1 ?" }8 p$ q3 s6 g
这个办法也好！我会试试看。多谢了。
作者: 晨枫    时间: 2019-2-5 00:07
关中农民 发表于 2019-2-4 07:47" C' u% H  L1 c
晨大，这得数学博士才中啊，额完全外行了，看见这个只能联想到面条 ...
  p2 l  @5 w, s/ |. h
哎，不是想着你们数据处理的问题多嘛，可能有经验。反正还是要谢一个！
作者: Dracula    时间: 2019-2-5 00:37
晨枫 发表于 2019-2-5 00:01
1 A3 ]0 J) S1 Y& M多谢。接下来我会按42楼办法用实际数据多试试，包括scaling，然后和目视结果比较。如果可靠的话，就能用 ...
$ W  f, N; e+ G0 ?
我在你这儿再贴几张Melissa Benoist的照片/ T6 G6 k( J. M! _/ A/ U
1 Y" D2 b3 K3 s8 Y6 s% a5 J% `- m
4 Q1 \/ w3 c; y  H/ ~( v& X' m
6 G9 M8 Y: |7 D+ F
( D+ M* _, U$ {9 {; n5 i, R7 c) ^
; h/ C% K" R  W  Z- c# L, I

5 M  D  N4 s- p' P) I& p# T1 h7 F

0 C0 N! X- t* W# ?1 R" w
8 X1 l! ]7 [, S# g1 W# y1 T! Y& s" g1 G, O5 S1 [$ {5 l* e
; m0 u1 }* f! y: k2 K; H
/ a, ^. |8 W7 k, {

: P" H9 K8 X+ @8 @1 z  b1 L 1 V" i. v* y' o) N2 M  j6 i# p" W* ^( _
  G( V: t5 U4 \4 t
3 p8 E" O" [9 g1 K: {% f
4 B1 O- n, X* G- {  M% e
6 P( a) ]2 f) F6 p9 Q) ]
. D' a" g- E' ^
; B" ~: C/ g" Z3 c( ?# q% B

( c1 \7 t. H- W* f) V6 C1 |. m
' X  X! [. P5 ]  p* M& f
7 ?3 ^" i. B* H! v5 b" r* x
, P) ?- |. [9 H, X5 |( E% S
1 r, w: }2 A, [' d5 [. a
8 e9 x' o  t9 o% ?( v6 d" T# p/ K2 O) r4 H6 D& u$ T0 M
5 D) g5 V" ^& P! P( L

3 v5 J8 t' B1 f% ]" u
4 [2 U! L. [: @7 _; {7 X
: b9 y3 N2 c( O4 U ; ~7 d3 {9 B  m7 A

& u3 H, Q6 K3 c: M' m& X) \* n% s# t6 h9 V% T$ H3 G; R! ?: s, _/ g

) I* T1 o' ?: o
作者: tanis    时间: 2019-2-5 00:52
晨枫 发表于 2019-2-4 22:25
. y- {% w8 x- b1 x- T: v2 d7 {  k是的，我以前还试过用MATLAB C通过OPC与DCS相连，在技术上这是做得到的，但可靠性达不到要求。OPC是不作 ...
1 k7 J0 \: N6 u8 s/ @. W
大工程的可靠性果然是我们小lab里完全不能比的。
作者: 数值分析    时间: 2019-2-5 00:52
本帖最后由 数值分析 于 2019-2-5 00:58 编辑
" w' f7 @, N: H1 Q: J7 X# R
Dracula 发表于 2019-2-4 22:534 ^/ |6 F5 _0 w; P" w8 {9 |
我曾经想过一个和42楼类似的办法。区别是分母不是板数，而是各个板加在一起的温度的和。如果这条曲线真的 ...

0 t$ Z& z. _2 n
: [( T$ K# _: d  C. r3 L, u这个和统计其实关系不大.你可以把他想象成求重心问题.已知一条曲线和x轴围成一个形状,如果这个形状是均匀厚度的匀质材料构成的一块板子,那么这个形状重心的x座标是多少?这个x座标(如果存在的话)就是这个分布的数学期望.这其实就是一个加权平均问题.当然,一个任意形状的重心和最高点的x座标当然不一定一致.不过数学上可以证明,高斯曲线和博松分布曲线围成的图形重心的x座标和最高点的x座标正好一样.和统计没关系.你再想想?
作者: 数值分析    时间: 2019-2-5 00:54
本帖最后由 数值分析 于 2019-2-5 01:04 编辑 . @6 K+ D' y. a
晨枫 发表于 2019-2-4 22:38- M, N+ R+ \) _8 |  J4 x1 E+ x
嗯，我再想想。谢谢。
7 B* E: L5 P. E5 b$ R3 A4 I
3 {* o# k& g2 v. p6 s$ H7 P. }
请见74楼回复.谢谢.  o( q7 D$ E* q. r9 a# ~5 F
& N9 C! Q- \1 s2 p0 U4 U
任意偏态分布最高点的位置就不能简单的用样本均值来估计了.不过也有办法,如果已知分布函数可以用矩估计或者最大似然估计.
作者: Dracula    时间: 2019-2-5 01:02
数值分析 发表于 2019-2-5 00:52
" @; V) `, x0 p4 B0 Y/ ?  g这个和统计其实关系不大.你可以把他想象成求重心问题.已知一条曲线和x轴围成一个形状,如果这个形状如果是 ...
% ^; S$ P  C& E# V9 L0 @$ i
这个和零点的选择是有关的。如果把温度类比为具体某块板的sample size的话，统计学的那个解就是以sample size（也就是温度）来加权求平均值。移动纵坐标的0点，相当于所有的sample size都加或减了一个相同的数。但是数学上，分子分母同乘同除同一个数可以，同加同减同一个数，值要变的。这个问题上，0点的选择应该完全是arbitrary的，因此这个办法应该是有问题的。& d& Q: e" w$ n, o5 Y- t4 s

! t  m( c" X6 B% [; u4 a: L% p" s
作者: 晨枫    时间: 2019-2-5 01:03
tanis 发表于 2019-2-4 10:52
& l- I/ i1 A  R" Q( {3 `大工程的可靠性果然是我们小lab里完全不能比的。
* \" ~) ?; q6 d1 S
是啊，实验室规模的可行性和工业规模的可靠性是两个很不同的概念，实施起来的考虑完全不一样。实验室是探路、开路的，工程是修路、维持交通和拉动经济的。
作者: 晨枫    时间: 2019-2-5 01:06
Dracula 发表于 2019-2-4 10:37
6 C+ U) B3 q, H我在你这儿再贴几张Melissa Benoist的照片

! o8 @7 c# ^) u1 o呵呵，好久没见Melissa了。她笑起来还是很charming的！少了点妖气，这是她的长处，还是缺点。Too well rounded, not enough sparkle or something to make you on edge。
作者: 数值分析    时间: 2019-2-5 01:07
Dracula 发表于 2019-2-5 01:02" v( F+ Z1 c/ S  U
这个和零点的选择是有关的。如果把温度类比为具体某块板的sample size的话，统计学的那个解就是以sample  ...
; u. m2 z4 @2 N' M, {1 t/ ^
你可以试试,平移没有问题的.你把他想象成求重心问题,曲线平移x,重心也平移x..
作者: Dracula    时间: 2019-2-5 01:09
数值分析 发表于 2019-2-5 01:07
/ I( p& g- f1 @# U4 X- U: m: C! c你可以试试,平移没有问题的.你把他想象成求重心问题,曲线平移x,重心也平移x.. ...

3 |& h2 @) s  n/ i我在纸上推过公式，结果是随0点的选择变化的。$ l3 \- a8 Q, N$ X! j
) X# Y: r* b: T" ?& Z/ O- v* ^
作者: 数值分析    时间: 2019-2-5 01:17
数值分析 发表于 2019-2-5 01:07
5 e( i" P$ K( G- u) V6 c) M你可以试试,平移没有问题的.你把他想象成求重心问题,曲线平移x,重心也平移x.. ...

( B% L& x) m8 u, F# l: Qintegral f(x)* x*dx=lambda右平移a个单位,则新重心位置integral f(x-a)*x*dx. 设t=x-a, integral f(t)* (t+a)*d(t+a)=integral f(t)* (t+a)* dt=integral f(t)* t* dt+integral f(t)* a* dt=lambda+a (因为integral f(t)*a *dt=a*integral f(t)* dt=a,而 integral f(t)* t* dt=lambda)$ Q# n3 N# j& Q* D  t3 B
形状右平移a个单位,重心也右平移a个单位
作者: 数值分析    时间: 2019-2-5 01:20
数值分析 发表于 2019-2-5 01:17
( N- ?! d$ l; r3 @integral f(x)* x*dx=lambda右平移a个单位,则新重心位置integral f(x-a)*x*dx. 设t=x-a, integral f(t)*  ...

6 z; L8 j" G6 e7 M: Z当然,前提是integral f(x)* dx=1,所以我跟晨风说要归一,否则确实不灵.
作者: Dracula    时间: 2019-2-5 01:21
晨枫 发表于 2019-2-5 01:06' M' _- E! x# _( \$ D
呵呵，好久没见Melissa了。她笑起来还是很charming的！少了点妖气，这是她的长处，还是缺点。Too well ro ...

' A& }3 j4 ]% r. M! M看来你不去我的那个Superhero电视剧美女贴，那儿我最近一个月基本上每隔几天就会贴几张的。昨天庆祝CW宣布Supergirl会有第五季，我刚贴了16张Melissa Benoist的照片。而且那里除了Melissa Benoist以外，别的美女我也经常会贴几张，象最近有Katie McGrath，Emma Watson，Virginia Gardner和Elizabeth Laith。欢迎常去那儿，观看加分。: t% ^; i, E& k5 I1 l

# U' g3 ]$ r  F0 l% I
作者: Dracula    时间: 2019-2-5 01:37
数值分析 发表于 2019-2-5 01:20  _2 Z9 Z9 _3 A
当然,前提是integral f(x)* dx=1,所以我跟晨风说要归一,否则确实不灵.
! `% i- [5 X) f  R& @7 ^4 S
曲线下面的面积等于1，这个条件肯定不满足。因为这本来就不是个概率论的问题。+ V# t5 v( |6 E- X
3 b5 m+ y3 {! V  h. v
那个公式是sum(xi * yi) / sum (yi), 如果纵坐标的零点移动，就是说yi' = yi + t, 你再算 sum(xi * yi) / sum (yi)不等于sum(xi * yi') / sum (yi')
( T! z7 h3 X2 T  z; P2 }! [( U3 _/ H1 `* y' h$ L0 m6 Q% O$ [0 N
作者: 数值分析    时间: 2019-2-5 01:41
Dracula 发表于 2019-2-5 01:37: @/ S7 M' l0 c$ B" D; h
曲线下面的面积等于1，这个条件肯定不满足。因为这本来就不是个概率论的问题。
0 P, ]* D8 Z( k8 c) S2 K
9 s" A$ R6 ?7 I- N那个公式是sum(xi * yi)  ...
2 F# ^- C* }% V  ~% `
所以我和晨风说要归一么.用histogram 面积归一以后,没问题.这实际是个加权平均问题,加权平均要求所以权重加起来和是1.即integral f(x)* dx=1,现在权重是温度,加起来肯定不是1.但只要除以总面积,(这里就是总温度),就还是满足这个关系的.不影响结果.晨风只关心最高点出现的位置,而不关心最高点是多少,这是关键.
作者: 晨枫    时间: 2019-2-5 01:46
Dracula 发表于 2019-2-4 11:218 A& s8 Z& l5 q) `' R
看来你不去我的那个Superhero电视剧美女贴，那儿我最近一个月基本上每隔几天就会贴几张的。昨天庆祝CW宣 ...
1 K$ K1 ]4 q" d( v* M
这等好地方怎么错过了？赶紧去看！
作者: 晨枫    时间: 2019-2-5 01:49
数值分析 发表于 2019-2-4 11:41- d0 u! ^% j- e
所以我和晨风说要归一么.用histogram 面积归一以后,没问题.这实际是个加权平均问题,加权平均要求所以权重 ...

  f+ n2 v$ W1 w# g3 Y话说，如果选“爱坛最学术贴”，这个贴有没有希望当选？我肯定投一票！0 [- K( D4 X8 R: I; D9 b% H

! Z) V6 A9 v7 W  Q多谢各位老大帮忙、指点。正在用Excel抓历史数据验算，看看这办法灵不灵验！
作者: 数值分析    时间: 2019-2-5 01:51
晨枫 发表于 2019-2-5 01:49
( \0 I5 N7 I( \9 `' ]) J2 y( b- C7 p* Z话说，如果选“爱坛最学术贴”，这个贴有没有希望当选？我肯定投一票！
8 [6 O4 t9 I& q$ j* b$ K6 @1 }- V& F* B. [5 s' f6 J0 x4 b# ]; |  X2 |+ x
多谢各位老大帮忙、指点。正在用 ...
: @! p% D' Z/ ~+ E: t$ l- c) L
* J1 @8 u* l7 C. q0 n8 b% p
如果不灵就是你那个偏态曲线和博松分布曲线实际上并不像,即重心和最高点不重合.不过有的修.如果到那一步咱们再谈怎么修.
作者: 数值分析    时间: 2019-2-5 01:54
晨枫 发表于 2019-2-5 01:49
9 R/ \, @8 A5 t) N, ^话说，如果选“爱坛最学术贴”，这个贴有没有希望当选？我肯定投一票！
* i6 F' e) t3 F) ~) H* s* C6 ^7 v
/ ?# g* i. O6 O4 i, ~9 b$ K多谢各位老大帮忙、指点。正在用 ...
4 t+ t6 G( C6 l" q
不过不管灵不灵,晨大可以帮我验证这样一个事儿,即把整个曲线平移n个单位,用同样的算法算完,结果应该是老结果平移n个单位.这个不管是不是博松,只要面积归一一定都灵.
作者: Dracula    时间: 2019-2-5 02:00
本帖最后由 Dracula 于 2019-2-5 03:16 编辑 7 V9 H0 K% t/ o+ J, H
数值分析 发表于 2019-2-5 01:41) ]; k1 V/ W" H  |# F
所以我和晨风说要归一么.用histogram 面积归一以后,没问题.这实际是个加权平均问题,加权平均要求所以权重 ...
! b$ _8 j: S+ [5 ^: J7 z

8 B& p8 n  I* ?  Z1 s假设一个最简单的情况吧。只有两个点，y1和y2，y1<y2，如果你把零点设在y1, 那么y2的权重是1，y1的权重是0，只有第二个点的值决定结果。但是如果零点设在接近于负无穷，那么不管y1, y2的值是多少，都接近相当于两个点的权重都是0.5。零点的选择肯定是对结果有影响的。但是因为零点的选择是arbitrary的，这种情况不应该出现，因此我认为这个算法有问题。
; O2 f; s- r1 a: P; ?/ i1 B1 @8 h/ z" r8 e# ^6 r
0 r+ Q# `5 Z+ L6 [! n& ^7 g* V
作者: 晨枫    时间: 2019-2-5 02:01
数值分析 发表于 2019-2-4 11:543 O5 W/ o3 e, A4 V
不过不管灵不灵,晨大可以帮我验证这样一个事儿,即把整个曲线平移n个单位,用同样的算法算完,结果应该是老 ...
! o3 Z- R9 V' p. k7 X9 [) f
我用“掐尾”正态分布已经试过了，不归一都精度不错。我再归一试试看！
作者: 晨枫    时间: 2019-2-5 02:02
数值分析 发表于 2019-2-4 11:51! k+ H# l$ J$ c% u/ |! l
如果不灵就是你那个偏态曲线和博松分布曲线实际上并不像,即重心和最高点不重合.不过有的修.如果到那一步 ...
* i! I1 ~( g$ Q
多谢！will report back!
作者: holycow    时间: 2019-2-5 02:15
数值分析 发表于 2019-2-4 09:418 Q# y( B. [" @5 _
所以我和晨风说要归一么.用histogram 面积归一以后,没问题.这实际是个加权平均问题,加权平均要求所以权重 ...

$ \* M' U8 _+ Z: k3 E1 |伯爵的意思是说，总温度凭什么以零摄氏度做原点？如果零度不是原点，则和原点的相对温度差之和完全是主观确定的，就不能拿来当scaling的分母
作者: 数值分析    时间: 2019-2-5 02:23
本帖最后由 数值分析 于 2019-2-5 02:25 编辑 $ S9 }6 K; |$ k  {$ R, M
holycow 发表于 2019-2-5 02:15
- w8 x0 X0 A8 w0 A9 r  B/ |伯爵的意思是说，总温度凭什么以零摄氏度做原点？如果零度不是原点，则和原点的相对温度差之和完全是主观 ...
5 a  `3 L, }5 v/ A" b6 W, e

& |( K, y9 z% |% t9 L- c8 B这个答案很简单,因为用零度才像泊松分布,如果上下平移的话,重心还是存在的,只是和最高点不再重合.你可以试想一下把泊松分布加上一,然后重新归一,也能得一个新的分布,这个分布也有期望,但期望很可能就不是最高点了. 不过单峰分布,只要不是骨骼轻奇(偏度skewness特别大),基本上最高点和重心差不太远.
作者: 数值分析    时间: 2019-2-5 02:32
holycow 发表于 2019-2-5 02:15
& _4 c, }% k7 y" O' T& d伯爵的意思是说，总温度凭什么以零摄氏度做原点？如果零度不是原点，则和原点的相对温度差之和完全是主观 ...

( c$ f- [2 q, D! m  \顺便说一下,如果是对称的单峰分布的话,就没有这个问题,随便上下平移,只要归一就可以.
作者: holycow    时间: 2019-2-5 02:42
数值分析 发表于 2019-2-4 10:32% `- N" o( r6 d- q# s# }- X# H, X
顺便说一下,如果是对称的单峰分布的话,就没有这个问题,随便上下平移,只要归一就可以. ...

$ L# b; B5 ^' P  P3 O% L& b. |8 K1. 极值出在哪里，只要估计出lambda即可5 y/ g, U  D% h5 Q: x
2. Lambda的估计需要依赖于归一. ]- I* g* o. F' C
3. 归一的分母是可以主观确定的 (导致曲线下面积变动)2 {& v, ?* Z; b4 e
: `1 y; a# V1 d6 e
就算是对称单峰分布，也要先解决这个峰的陡峭程度才知道这个峰在哪里，恰恰是峰的陡峭程度依赖于归一的分母...
作者: tanis    时间: 2019-2-5 03:26
数值分析 发表于 2019-2-5 02:23, Q  x* {% e3 G, H  V+ z7 i( j
这个答案很简单,因为用零度才像泊松分布,如果上下平移的话,重心还是存在的,只是和最高点不再重合.你可以 ...
& p1 }9 i2 V; H/ t; ]  P( F1 J7 ~
冒昧的问一句，你搞过竞赛么~
- Q" H- y7 A) B* p5 w$ c9 C' d! f  A8 x; `$ X# ~2 O
思维方式挺像的~
作者: Dracula    时间: 2019-2-5 03:43
数值分析 发表于 2019-2-5 02:23  j7 u- h% `' `% N) X
这个答案很简单,因为用零度才像泊松分布,如果上下平移的话,重心还是存在的,只是和最高点不再重合.你可以 ...
( q% c* O9 |8 ?+ c3 a" C! A
问题就是这个0度在哪儿你并不知道。至于曲线下的面积必须是1这一点，只要各个点同乘或同除一个数就都可以做到，这个条件并不能提供任何额外的约束来确定零度这个参数。  _9 {# d# |0 e/ e) M9 l8 B9 O

: K( a& h0 @- W9 i1 n8 O" X/ w
作者: 木不铎    时间: 2019-2-5 05:31
不麻烦啊。查一下维基百科上关于“泊松分布”的页面嘛。( Z; D2 [' u$ N1 a1 u+ A3 L# I. c
( B7 @8 ~# w; E" [, W9 j2 v
泊松分布的概率密度函数为
. \9 C3 W. P, Q! _! R
2 d: C2 H8 o/ w( A2 p其中λ是单位时间（或单位面积）内随机事件的平均发生率，k代表发生某类事件的次数。
7 O# L; G, ], w: l+ V: V这里有一个很好的例子如下：
' y( a9 l) H. B& f/ O3 p2 \- u% W: s5 H# L

7 S8 e9 {! G# k* n7 I6 ?
对某公共汽车站的客流做调查，统计了某天上午10:30到11:47来到候车的乘客情况。假定来到候车的乘客各批（每批可以是1人也可以是多人）是互相独立发生的。观察每20秒区间来到候车的乘客批次，共观察77分钟*3=231次，共得到230个观察记录。其中来到0批、1批、2批、3批、4批及4批以上的观察记录分别是100次、81次、34次、9次、6次。使用极大似真估计（MLE），得到 的估计为λ=（81*1+34*2+9*3+6*4）/231=0.8658。

5 q+ V; O$ ^" b8 k/ Q3 e6 e7 d( w7 R" R. o
也就是说20秒之内平均有0.8658批客人。
0 l+ v! b) X8 Q: F, `2 I
2 c# X! ^) Z8 g; E- q0 M9 z0 |这个例子应该和斯基的问题很类似。根据统计数字，用这个MLE方法，就能得到你的均值λ
作者: 晨枫    时间: 2019-2-5 05:49
木不铎 发表于 2019-2-4 15:31/ s; h# \" h. @
不麻烦啊。查一下维基百科上关于“泊松分布”的页面嘛。
- |4 p6 M% y& R  U6 u. g7 W  {" H3 c! y
泊松分布的概率密度函数为
5 m! o) s1 r3 t! l
谢谢。这和42楼“数值分析”的方法是一样的。



欢迎光临 爱吱声 (http://aswetalk.net/bbs/) Powered by Discuz! X3.2