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

今天这篇文章给大家介绍下我编写的产生随机数的函数。

函数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的编程，感兴趣的话可以点击下面的链接或者左下角的【阅读原文】：


免责声明：如果侵犯了您的权益，请联系站长，我们会及时删除侵权内容，谢谢合作！