今天为了方便调试上位机程序,主要是VB中曲线更新绘制,本来以为会有现成的函数,毕竟有些工具里自带,比如Matlab。
网上搜了一下,PLC是没有的随机函数的,需要自己编写。于是决定用PLC编程实现,在指定范围内生成一个随机数。
发现一个帖子,是2005年的工控网帖子,讨论的比较思路比较清晰,还有例程。10多年前的就有人思考的问题,我现在才思考,想来自己还差的很远。
原文链接:用S7-200产生一个随机数的方法
转帖如下:
下面看这样一个C程序:
//rand01.c
#include
static unsigned int RAND_SEED;
unsigned int random(void)
{
RAND_SEED=(RAND_SEED*123+59)%65536;
return(RAND_SEED);
}
void random_start(void)
{
int temp[2];
movedata(0x0040,0x006c,FP_SEG(temp),FP_OFF(temp),4);
RAND_SEED=temp[0];
}
main()
{
unsigned int i,n;
random_start();
for(i=0;i<10;i++)
printf("%u\t",random());
printf("\n");
}
时钟中断服务程序中,每当TIMER_LOW转满时,此时,记数器也会转满,记数器的值归零,即TIMER_LOW处的16位二进制归零,而TIMER_HIGH加一。rand01.c中的movedata(0x0040,0x006c,FP_SEG(temp),FP_OFF(temp),4);
正是把TIMER_LOW和TIMER_HIGH两个16位二进制数放进temp数组,再送往RAND_SEED,从而获得了“随机种子”。
现在,可以确定的一点是,随机种子来自系统时钟,确切地说,是来自计算机主板上的定时/计数器在内存中的记数值。这样,我们总结一下前面的分析,并讨论一下这些结论在程序中的应用:
//rand02.cpp
#include
#include
using namespace std;
int main()
{
unsigned int seed=5;
srand(seed);
unsigned int r=rand();
cout< }
//rand03.cpp
#include
#include
using namespace std;
int main()
{
srand((unsigned)time(NULL));
unsigned int r=rand();
cout< return 0;
}
//rand04.cpp
#include
#include
using namespace std;
int main()
{
int rNum,m=20;
char *ch=new char[m];
for ( int i = 0; i //大家看到了,随机种子会随着for循环在程序中设置多次
srand((unsigned)time(NULL));
rNum=1+(int)((rand()/(double)RAND_MAX)*36); //求随机值
switch (rNum){
case 1: ch[i]='a';
break ;
case 2: ch[i]='b';
break ;
case 3: ch[i]='c';
break ;
case 4: ch[i]='d';
break ;
case 5: ch[i]='e';
break ;
case 6: ch[i]='f';
break ;
case 7: ch[i]='g';
break ;
case 8: ch[i]='h';
break ;
case 9: ch[i]='i';
break ;
case 10: ch[i]='j';
break ;
case 11: ch[i]='k';
break ;
case 12: ch[i]='l';
break ;
case 13: ch[i]='m';
break ;
case 14: ch[i]='n';
break ;
case 15: ch[i]='o';
break ;
case 16: ch[i]='p';
break ;
case 17: ch[i]='q';
break ;
case 18: ch[i]='r';
break ;
case 19: ch[i]='s';
break ;
case 20: ch[i]='t';
break ;
case 21: ch[i]='u';
break ;
case 22: ch[i]='v';
break ;
case 23: ch[i]='w';
break ;
case 24: ch[i]='x';
break ;
case 25: ch[i]='y';
break ;
case 26: ch[i]='z';
break ;
case 27:ch[i]='0';
break;
case 28:ch[i]='1';
break;
case 29:ch[i]='2';
break;
case 30:ch[i]='3';
break;
case 31:ch[i]='4';
break;
case 32:ch[i]='5';
break;
case 33:ch[i]='6';
break;
case 34:ch[i]='7';
break;
case 35:ch[i]='8';
break;
case 36:ch[i]='9';
break;
}//end of switch
cout< }//end of for loop
cout< return 0;
}
而运行结果显示的随机字符串的每一个字符都是一样的,也就是说生成的字符序列不随机,所以我们需要把srand((unsigned)time(NULL)); 从for循环中移出放在for语句前面,这样可以生成随机的字符序列,而且每次运行生成的字符序列会不同(呵呵,也有可能相同,不过出现这种情况的几率太小了)。
--------------------------------------------------------------------------------------
litianxi回复:
在很多场合都要用到随机数,比如随机分配职工电影票(我总是最角落,边上一定是公司最漂亮的MM)
随机选定公司出差(公游)的人员(27选2,选来选去总有俺,另一个一定是在12个漂亮妹妹中随机选定)
所以,掌握这个东东还是有点用处的……
[题目:]
制作一个函数,用于输出从1到36(数字)中的任意7个,而且输出的7个数字不能重复(买过体彩的应该清楚)
大家试试啦……
1、对于种子的初始化,如果不是对随机数的应用有非常特殊的要求,不要反复初始化。
对于前面那段 C++ 代码,在循环中重复初始化种子不是很好。我没有调试那段代码,但可以预见,每次循环生成的随数会比较接近,按一定的规律不严格的递增或递减。因为这段代码在计算机中运行速度会非常快,也就是说时间变化很小,每次初始化的种子就很接近,从而使生成的散列有一定的“规律”。
如果不是对随机数的算法和自己应用的要求理解得很清楚,就不要反复初始化种子。前面的代码实际上是"规律“性的初始化,已经有背于”随机“的要求了
2、伪随机数的产生,现在用得较多的是“线性同余法"
就是下面这个式子
R(n+1) = [R(n) * a + b] mod c
为使随机数分布尽量均匀,a、b 均为质数, c 一般取值域内的最大值(mod 是求余数)
从这个式了可以看出,每次产生的随机数都跟上一次产生的数有关系,那么,第一个数是怎么来的呢?这就是线性同余法中必须用的的”种子",也就是说,给定某个种子后,所产生的随机数序列是固定的,在计算机编程中,一般使用系统时间来初始化种子,就是前面代码中的 srand((unsigned)time(NULL)); 这一句了。因为每次运行程序的时间肯定不一样,所以产生散列肯定也不一样,从而达到“随机”的目的。
a,b,c 的取值我用的是 a=3373, b=1, c=32768
下面的两个子程序是我在我的项目(S7-200 226)中产生随机的系统编号用的,因为我的编号中只有4位数采用了随机数,所以下面的程序中用的是整型,最大范围为32767。如果需要更宽范围的随机数,可以采用双字类型,并适当修改程序,代码很简单,就是将上面那个表达式用 S7-200 的指令表示出来就行了。
这两个子程序是从 MicroWIN V4.0 中导出来的,可以将它们用文本编辑器保存为 AWL 文件后直接导入 MicroWIN。
使用时在第一个扫描周期调用 Srand 初始种子,需要随机数的地方调用 Random
Random 有了个最大范围参数,可以限制生成的随机数的最大范围,比如我只需要4位随机数,所以一般这样调用 CALL Random, 10000, vw0,生成的数就在 0-9999 范围内
下面是代码:
SUBROUTINE_BLOCK Srand:SBR17