楼主思考的很深入啊，来支持

mark

楼主有没有广泛的测试过? 看看在生成大量随机数后会不会出现特殊情况(如http://www.zhihu.com/question/20222653)

楼主根本就没理解这句“痴语”。

一、说计算机不能生成，指的是计算机通过计算不能生成。 这一点在目前是肯定的。

二、楼主你的驱动读外部数据，这明显就不是计算生成的了。  那是计算机外部的设备生成的了。

这个谁不知道啊，纯计算和计算机无关，就像1+1=2一样，你还想算出3吗？
“专业路过”说的情况是有可能的，处理好就没有问题了。生成数据能通过各种测试。

附件内容
1 以秒为单位用时间戳函数连续测量Beep（0,0）的数据
2 以周期数为单位连续测量Beep（0,0）的数据

上传的附件：

• QQ图片20140720155952.jpg （10.25kb，3次下载）
• QQ图片20140720160056.jpg （8.48kb，1次下载）

#include <stdio.h>
#include <stdlib.h>

#define php_uint32 unsigned int//!!
#define php_int32 int//!!
#define PHPAPI
#define MT_N (624)
#define M             (397)                /* a period parameter */
#define N             MT_N                 /* length of state vector */
#define BG(x) (x)
#define loBit(u)      ((u) & 0x00000001U)  /* mask all but lowest    bit of u */
#define loBits(u)     ((u) & 0x7FFFFFFFU)  /* mask     the highest   bit of u */
#define hiBit(u)      ((u) & 0x80000000U)  /* mask all but highest   bit of u */
#define mixBits(u, v) (hiBit(u)|loBits(v)) /* move hi bit of u to hi bit of v */
#define twist(m,u,v)  (m ^ (mixBits(u,v)>>1) ^ ((php_uint32)(-(php_int32)(loBit(u))) & 0x9908b0dfU))
#define TSRMLS_D
php_uint32        gstate[MT_N+1] = {0};  /* state vector + 1 extra to not violate ANSI C */
php_uint32        *next = NULL;       /* next random value is computed from here */
int                        left = 0;        /* can *next++ this many times before reloading */

unsigned int asm_get_rndseed()
{
        asm("rdtsc");
}

static inline void php_mt_initialize(php_uint32 seed, php_uint32 *state)
{
        /* Initialize generator state with seed
           See Knuth TAOCP Vol 2, 3rd Ed, p.106 for multiplier.
           In previous versions, most significant bits (MSBs) of the seed affect
           only MSBs of the state array.  Modified 9 Jan 2002 by Makoto Matsumoto. */

        register php_uint32 *s = state;
        register php_uint32 *r = state;
        register int i = 1;

        *s++ = seed & 0xffffffffU;
        for( ; i < N; ++i ) {
                *s++ = ( 1812433253U * ( *r ^ (*r >> 30) ) + i ) & 0xffffffffU;
                r++;
        }
}

static inline void php_mt_reload()                                                         
{      
        /* Generate N new values in state
           Made clearer and faster by Matthew Bellew (matthew.bellew@home.com) */                  
        
        register php_uint32 *state = BG(gstate);                                                     
        register php_uint32 *p = state;                                                            
        register int i;                                                                             
        
        for (i = N - M; i--; ++p)
                *p = twist(p[M], p[0], p[1]);                                                      
        for (i = M; --i; ++p)
                *p = twist(p[M-N], p[0], p[1]);                                                     
        *p = twist(p[M-N], p[0], state[0]);                                                         
        BG(left) = N;
        BG(next) = state;                                                                           
}

PHPAPI php_uint32 php_mt_rand()
{
        /* Pull a 32-bit integer from the generator state
           Every other access function simply transforms the numbers extracted here */

        register php_uint32 s1;

        if (BG(left) == 0) {
                php_mt_reload();
        }
        --BG(left);

        s1 = *BG(next)++;
        s1 ^= (s1 >> 11);
        s1 ^= (s1 <<  7) & 0x9d2c5680U;
        s1 ^= (s1 << 15) & 0xefc60000U;
        return ( s1 ^ (s1 >> 18) );
}

int main()
{
        unsigned int array[1000] = {0};
        php_mt_initialize(asm_get_rndseed(), gstate);
        for (int i = 0; i < 10000000; i++)
        {
        php_mt_initialize(asm_get_rndseed(), gstate);
                unsigned int rr = php_mt_rand();
                unsigned int rr_ = rr % 1000;
                //printf("%u\n", rr_);
                array[rr_]++;
        }

        for (int i = 0; i < 1000; i++)
        {
                printf("%d:%d\n", i, array[i]);
        }
}
用这个吧，直接从php扒出来的，Mersenne Twister算法

标记一下

谢谢8楼介绍的伪随机数算法。

就你这样2个数据根本说明不了问题，还妄想真随机，你觉得你这思路别人想不到么？为什么没有普及开呢？

我手拿一颗球，往地上丢弹起来，每次的时间都不一样，但也不能称为随机。

　　是不是真随机数，要拿数据来说啦，你觉得不是你可以证明它不是，例如你可估算它后面的值，如果总是说对了，那它肯定是伪随机数了。
　　楼上扔球的时间肯定是随机数了，但要将这些随机数直接使用是不行的，需要用数学方法将随机性提取出来才好作其它使用。

同意你这个说法  而且我觉得lz的这个办法  虽然看出去好像是真随机   但是从科学的角度很不严谨   我不是数学专业  但是我起码了解高斯分布曲线   这东西可以反映数据分布的集中性   我觉得理想状态的随机是不可能的    就算是数学上的真随机   【我觉得】也只是对数据本身的特性规定一个指标   当这个特性达到或高于这个指标时候  我们就可以认为数据是真随机   而且  我以前看过相关的帖子   有点点印象记得数据的随机【程度】  的确是可以用某种散点图来体现的    具体可以去查阅相关资料
     第二点就是，我个人觉得lz的办法出来的数据可能会有如下情况：由于随机数的生成取决于执行时间，首先，计算机的机器时间周期不是连续且大致稳定的，函数的执行时间也肯定是机器时间周期的整数倍，毕竟函数的执行本质上就是机器指令的执行。所以我觉得这个办法出来的随机数， 从宏观来看，会在机器时间周期的整数倍那里相对集中。
    当然就像我第一点所说的，如果说从严谨的科学角度看，虽然数据是有部分集中，但是如果数据的某种特性达到某个规定真随机的指标，那么我们也是可以说办法可行，的确是真随机

非常赞同！

　　上面是利用声卡驱动生成随机数的例子，方法很简单随便编个小程序就可以验证它，随机数根据考查是均匀分布的，是不是真随机数需要有人去证实，我还没有观察到反例。如果用来生成少量随机数还行，数量大了效率就太低了。
　　上面是利用计算机硬件来获得随机数。用纯粹计算能得到真随机数吗？我想也是可以的但程序不会太简单，可以利用模拟技术来实现仅举一例，例如模拟256个小球在矩形空间做二维刚体碰撞运动，如果初速度和运动方向都是不同的，它们将演绎出形形色色的运动状态，当运动达到相对稳定后可以从小球的坐标提取随机量...

　　8楼acqqer的方法也可以生成真随机数，他的方法是采用系统周期计数值作为种子，启动一个周期达到10的6000次方的随机函数取值作为输出。再优化一下是可以实用的。

楼主的思路反了。不要被现象蒙蔽本质

　　不要总是指出不对，请用道理或数据打倒楼主的谬论。楼主认为真随机数无非是一些没有规律的乱码，即使理论上能找出规律，实际上很困难或找不出规律的都可以一样放心使用。
　　有人曾对我说投掷骰子也不能生成随机数，理由是你如果把这个骰子的所有动力学参数都搞清楚，那么骰子的运动就是必然的事情毫无随机性可言，但实际上骰子的运动状态由多种因素决定，你的测量未必能在不干扰骰子运动的条件下完成，所以根本不可能实现，即使测量实现了你也无法控制那些因素，例如投掷者的力道方向等，所有那些说道毫无意义，只是搞笑而已。

随机数貌似可以结合混沌就可以啦 ux(1-x)，当u=4时，可以达到好的随机性，然后就是组装组装其他公式～

不错，不错

　随机函数再好也是源于同样的规则，生成的数据无法和真随机数比，尽管周期大的了不得，如果大量引用据说能产生一些不相干的东西使你的运算目标受到影响。

貌似windows是多任务的，任何一个函数的执行时间可能恒定吗？

http://zh.wikipedia.org/wiki/%E9%9A%8F%E6%9C%BA%E6%95%B0

　　做过一些测试，一般函数算式等执行时间都是严格一致的。但例外的也有就像上面说的。据说有的计算机已经有专门生成随机数的芯片，但还没见到。实际上计算机的指令寄存器的数据实际是很随机的，但是无法将其读出，否则将是不错的随机数源头。

不太懂这些东西，坐看打什么讨论