常用数字信号的产生(C实现)-均匀分布&正态分布
程序员文章站
2022-03-18 14:48:03
小白博主开始学DSP之路,为了激励自己能坚持下去,写一系列博客来记录下来,也欢迎与大家一起讨论。介于我能力所限,这里学习目标定为,学习一个基础知识,写一个C语言程序实现, 最后会形成一个C的函数库,方便调用。对于学习的知识内容,不会面面覆盖,一般只学习一种实现方法,较少涉及复杂数学理论的证明,但保证 ......
小白博主开始学DSP之路,为了激励自己能坚持下去,写一系列博客来记录下来,也欢迎与大家一起讨论。介于我能力所限,这里学习目标定为,学习一个基础知识,写一个C语言程序实现,
最后会形成一个C的函数库,方便调用。对于学习的知识内容,不会面面覆盖,一般只学习一种实现方法,较少涉及复杂数学理论的证明,但保证有理论可以证明,可以查阅数学相关知识自寻证明。
对于所有的C程序,我会调通,供大家参考讨论,完整程序我放到github上,文章中会给出关键代码,程序只使用C语言,因为涉及数学与实现过程,C虽然繁琐,但是易于理解。
这里从数字信号的产生做为出发点,欢迎大家关注,一起学习!
这是第一篇博客,介绍均匀分布的随机数的产生和正态分布随机数的产生。
一、均匀分布随机数产生方法简介
1 //混合同余法,产生(a,b)区间的均匀分布随机数 2 //a:区间下界 b:区间上界 *seed:随机数种子 3 double uniform(double a,double b,long int* seed) 4 { 5 double t; 6 *seed = 2045 * (*seed) + 1; 7 *seed = *seed - (*seed/1048576)*1048576; 8 t = (*seed)/1048576.0; 9 t = a + (b - a) * t; 10 return t; 11 }
二、正态分布随机数产生方法简介
1 //产生均值mean方差sigma的高斯分布随机数 2 #include"uniform.h" 3 double gauss(double mean,double sigma,long int* seed) 4 { 5 int i; 6 double x,y; 7 8 for(x=0,i=0;i<12;i++) 9 x = x + uniform(0.0,1.0,seed); 10 x = x - 6.0; 11 y = mean + x * sigma; 12 return y; 13 }
完整代码参考:https://github.com/HeYingnan/dsp