[西门子] 西门子SCL编程实例——产生随机数的算法

[复制链接]
查看92391 | 回复0 | 2024-5-15 21:26:35 | 显示全部楼层 |阅读模式
今天这篇文章给大家介绍下我编写的产生随机数的函数。



函数FC5006_RandomReal可根据需要产生一个指定范围内的随机数,它采用线性同余算法产生随机数,其公式如下:

X[n+1]=(a*X[n] +c) % m

这是一个递归公式,利用X[n]产生X[n+1]。其中:

    a:系数,取值范围:0<a<m;

    c:增量,取值范围:0<=c<m;

    m:模,取值范围:m>0;

说明:
①m应尽可能的大,通常至少大于2^30(2的30次方);
②如果m选取为2的幂(即m=2^n),则a的值通常应满足a mod 8=5;
③当m和a的选取比较合理时,对于c的约束并不强,但要保证c与m互质,例如c可以选择1或11;
这个递归公式中,X[n]用来产生X[n+1],因此,X[n]被称为种子。也就是说,前一个随机数是后一个随机数的种子。在实际使用时,第一个种子一般是作为参数赋值给函数的。这种递归公式计算对于C/C++等语言很适合,但是对于PLC的编程语言有所困难。因此,在实际PLC编程时,可以使用系统时间戳作为种子,即:

randomNumber=(a*seed +c) % m;

下面介绍下我编写的函数FC5006_RandomReal,它有两个输入参数:

    minValue:需要产生随机数范围的最小值;

    maxValue:需要产生随机数范围的最大值;

三个常数:

    M_MOD:公式中的m;

    A_FACTOR:公式中的a;

    C_INCREMENT:公式中的c;


所有参数如下图所示:



代码如下:



我编写了一个函数测试随机数算法,产生一个1~1000之间的实数,代码如下:



多次测试结果如下:



如果你对这个函数有任何问题,欢迎留言讨论。

我的书《西门子S7-1200/1500 PLC SCL语言编程》比较详细的介绍了SCL的编程,感兴趣的话可以点击下面的链接或者左下角的【阅读原文】:

本帖子中包含更多资源

您需要 登录 才可以下载或查看,没有账号?注册哦

x
您需要登录后才可以回帖 登录 | 注册哦

本版积分规则