问题：如何从数据里估算普瓦松分布的均值？

tanis

晨枫 发表于 2019-2-4 13:33
唉，MATLAB里有histfit命令，干的正是我要的，可惜没法“偷”过来用啊

7 k9 E6 t3 K- c! |Matlab 直接有poissfit啊。。。

+ i; D5 m4 j- I% g' N/ x我感觉人类脑补的水平估计和取中值也差不了多少了既然连忘了看也只是可能出问题，我脚的中值差不多 ）（缺乏工匠精神的我。。。）

松叶牡丹

晨枫 发表于 2019-2-4 13:33
唉，MATLAB里有histfit命令，干的正是我要的，可惜没法“偷”过来用啊

完全不会FORTRAN，这样啊，要不就得自己写算法了？

晨枫

松叶牡丹 发表于 2019-2-3 23:21
6 I+ d* v6 [( ]- b降低精度，把一系列正态分布离散化为有限几个histogram，每次测试点先第一步进行范围判断，然后归到不同 ...

) v  A( u. f/ a: T; x唉，MATLAB里有histfit命令，干的正是我要的，可惜没法“偷”过来用啊

晨枫

tanis 发表于 2019-2-3 23:27
哦哦。。。 那就取最高和第二个的中间位置 有人工干预的时候就用人工，人忘了看就用这个。。。 ...

1 [; L% `8 _  h$ h3 x这也是个办法，总是有点“心有不甘”……

晨枫

松叶牡丹 发表于 2019-2-3 23:21
0 I( T5 F8 {: s. Y0 n降低精度，把一系列正态分布离散化为有限几个histogram，每次测试点先第一步进行范围判断，然后归到不同 ...

1 m4 [2 L0 l7 H& \+ T
  \4 W4 e" L( R. E3 k% K一把抓住：我的问题正是从histogram变换为某种分布，最想poison分布或者log normal distribution，不行的话正态分布也行。有什么现成的计算方法吗?

tanis

晨枫 发表于 2019-2-4 13:23
" ?) T# Z* ^, n2 z  Z. I; U塔板是等间距分布的，这样我还不如直接取温度最高的那一块了

3 H7 b* t1 n- ?1 b0 j# R哦哦。。。 那就取最高和第二个的中间位置 有人工干预的时候就用人工，人忘了看就用这个。。。

晨枫

tanis 发表于 2019-2-3 23:12
( H8 v0 ]7 t- S- V) I! \找温度第二，第三高的两个点，拟合个直线，取中间的那个位置

. j! }$ D' @4 O- D) q( V) h塔板是等间距分布的，这样我还不如直接取温度最高的那一块了

松叶牡丹

3 K/ S6 G  c. i# x3 o' _4 y; A

晨枫 发表于 2019-2-4 13:18
正是：看数据点然后脑补。人工就是这么干的。但操作工要照管的事情太多，一会儿不看，可能就错过了，就可 ...

& u# F' M2 Y6 \8 G0 @5 H8 r
降低精度，把一系列正态分布离散化为有限几个histogram，每次测试点先第一步进行范围判断，然后归到不同的histogram里面。再反向把histogram变换为对应的大致正态分布，这样可行吗？

晨枫

松叶牡丹 发表于 2019-2-3 23:12
您是说每块塔板都是15秒测量一次，所以每15秒要同时给出20个塔板的20个正态分布式子？ ...

& F2 F* b0 c* m1 c每15秒要对20块塔板的温度拟合一条样子像正态分布的曲线。15秒只是计算间隔，每次实际计算其实只有毫秒级的窗口，因为同一个CPU还要分时管好多别的事呢。

晨枫

tanis 发表于 2019-2-3 23:09
5 y* A7 Y1 g% D) j1 N目视是指看数据点脑补么？ 所以一共20个数据点（塔板），然后最高温度可能在点之间？ 15s一次，速度很慢 ...

正是：看数据点然后脑补。人工就是这么干的。但操作工要照管的事情太多，一会儿不看，可能就错过了，就可能坏事。

6 R: A! v9 t; m' {* t" J化工上DCS的采样分基本回路和先进回路，前者一般每秒一次，甚至1/4秒一次，后者每15、30、60秒一次就够了。过程时间常数动辄几十分钟，太快的采样没必要。另一方面，CPU的速度（还是Morotola 68040）和网络带宽（4 megabits 烹per second）有限，不宜太快，还要扫描几百、几千个回路，不宜太快。安全连锁回路更快，在毫秒级，但那不用作连续控制，两回事。

工控级的CPU和网络都用特别成熟（翻译：特别古老）的技术，和实验室、离线是两回事。
" q2 e1 X+ ~' O/ f
3 _5 I/ F. l1 @/ R' j. t9 v在线计算是因为我要把这个峰值点用到控制回路里。拍照、图像识别、离线拟合，这绕的弯子太大，中间节点太多，可靠性根本靠不住，哪里打一个嗝，就可能造成process upset，就是几十万、上百万的损失，那就划不来了。

松叶牡丹

晨枫 发表于 2019-2-4 13:06
切，我有那东西，还费这个事吗？Matlab、Maple，什么不能上啊！问题是我需要在线计算，而且只有FORTRAN一 ...

您是说每块塔板都是15秒测量一次，所以每15秒要同时给出20个塔板的20个正态分布式子？

tanis

晨枫 发表于 2019-2-4 13:06
7 U4 C3 c* H) y+ J4 U4 \  h切，我有那东西，还费这个事吗？Matlab、Maple，什么不能上啊！问题是我需要在线计算，而且只有FORTRAN一 ...

: Q0 d, K2 T- g0 E  z  V% d找温度第二，第三高的两个点，拟合个直线，取中间的那个位置

tanis

晨枫 发表于 2019-2-4 13:03
咳咳，这个其实不是统计问题，是从有限的温度测量估计温度分布曲线的问题。吸收塔一共20块塔板，每块塔板 ...

% {3 J$ M, @' f/ ?! V目视是指看数据点脑补么？ 所以一共20个数据点（塔板），然后最高温度可能在点之间？ 15s一次，速度很慢啊，为啥不能离线算？ 系统太古老了？
6 _4 e* Z/ B/ P& r5 t+ L1 Y
- ^5 ^2 `7 T& G) K5 D如果是系统太老无法整合还有个鬼畜方法。。。 既然肉眼能看到每个点在哪里，可以用电脑拍照把图截下来，然后从图像里把点的值拿到，接着在电脑里拟合。

晨枫

松叶牡丹 发表于 2019-2-3 23:00
做成在线计算器，那就不能查表了，反而太麻烦。坐等观摩高手们过招。

估计吃完年夜饭，就该忙着数红包了。

晨枫

tanis 发表于 2019-2-3 23:00
1 r4 A' |8 X* [3 |原来是泊松分布。。。 惭愧，连最小二乘法都属于复杂算法了。。。还真不知道有什么其他的简便方法。这一 ...

5 o- w! H7 z" ^7 h; p切，我有那东西，还费这个事吗？Matlab、Maple，什么不能上啊！问题是我需要在线计算，而且只有FORTRAN一级的语言从头现编，还不能太耗费CPU时间，耽误了实时控制就偷鸡不成蚀把米了。

晨枫

煮酒正熟 发表于 2019-2-3 22:58
0 u1 t1 X  l$ v  v统计学我只明白最最基本的，泊松分布也理解不深。问几个傻问题吧：一共多少块塔板？为什么每块塔板只能取到 ...

- h' l1 e9 x- [1 l: e& {. P咳咳，这个其实不是统计问题，是从有限的温度测量估计温度分布曲线的问题。吸收塔一共20块塔板，每块塔板只有一个温度测量点，连续测量（每15秒采样一次）。

; Y7 k' l1 ^$ E, w! ^4 ?这温度分布随工艺条件而变，所以这个峰值和凸起一直在动来动去。目视很容易看出峰值在那里，但要用数值办法实时确定，反而犯难了。

松叶牡丹

8 d2 V! `$ l9 o# M
  t. t  Z/ A  W做成在线计算器，那就不能查表了，反而太麻烦。坐等观摩高手们过招。

tanis

晨枫 发表于 2019-2-4 12:46
是我描述得不好。再来一遍。

! K6 A; n/ f% M' R' A+ R4 z# y我有一条样子像泊松分布的温度分布曲线，但只有几个稀疏的点，想用类似泊松 ...

+ o7 j3 z. H% b原来是泊松分布。。。 惭愧，连最小二乘法都属于复杂算法了。。。还真不知道有什么其他的简便方法。这一块早已经还给老师。。。 现在计算能力强大，集成性高，一般直接一个pasfit之类的函数就解决问题了。。。哈哈哈

晨枫

松叶牡丹 发表于 2019-2-3 22:52
2 k8 m; b( }6 r) O) }' W; [如果精度不要求很高的话，不知道查表法的精度能不能达到要求。那个倒是最快。 ...

查表再内插，这要做成在线计算还是有点麻烦了？

煮酒正熟

统计学我只明白最最基本的，泊松分布也理解不深。问几个傻问题吧：一共多少块塔板？为什么每块塔板只能取到一个data point? 我猜是取样也有成本问题？那样儿的话，您知不知道每块板上的随机取样，与这块板的温度峰值是不是最接近、比其他板的随机取样更接近本板的峰值？