资源描述：
《指定功率的白噪声的产生方案.docx》由会员上传分享，免费在线阅读，更多相关内容在应用文档-天天文库。
1、指定功率的白噪声的产生方案12020年4月19日文档仅供参考，不当之处，请联系改正。随机信号课程设计杜华贵起超刘越胡谦徐涛李盼辉一、题目：设计指定功率为20的白噪声产生方案并实现。分析：白噪声过程是服从一定分布的随机过程，它的特点是功率谱密度为常数。平稳白噪声过程均值为零，自相关函数在t=0处为一个冲击，功率谱函数在整个频率范围内是常数。因为其均值为零，因此方差在数值上等于功率。实际上，白噪声都是理想模型，在模拟系统中不可能存在，因此我们利用混合同余法生成随机数的方法产生一个无限近似的白噪声。为了分析简单方便，易于理解，利用混合同余法产生一个
2、服从X~U[-√3,√3]的白噪声。其均值为0，方差和功率为1。二、原理：产生随机数的方法是先用一定的方法产生[0,1]均匀分布的随机数，然后经过一个适当的变换就能够得到符合某一概率模型的随机数。常见的产生[0,1]均匀分布的随机数的方法有乘同余法和混合同余法。运用混合同余法产生在[-√3,√3]均匀分布的随机数,经过调整种子数、乘子和随机数序列的长度，能够优化使其分布更加接近X~U[-√3,√3]。并经过MATLAB编程实现其随机分布，计算并绘制22020年4月19日文档仅供参考，不当之处，请联系改正。出自相关函数，功率谱密度，概率密度曲线
3、的图像。●混合同余法介绍：经过同余运算生成伪随机数的方法称为同余法，常见的同余法包括加同余法、乘同余法、混合同余法、除同余法。其中乘同余法和混合同余法的性能更好，有速度快、内存省、周期长、统计特性好等优点。混合同余法是Lehmer在1951年提出的，其迭代公式为：Xnmod(AXn-1C,M)(0.1)Yn=Xn/M(0.2)公式(0.1)、(0.2)中，mod表示求余函数，A,C,M均为正整数。其中M是模数，A是乘子，C是增量，X0为初始值(0X0M)，即种子数，当C=0时，称此算法为乘同余法；若C0，则称算法为混合同余法，当C取不为零的适
4、当数值时，有一些优点，但优点并不突出，故常取C=0。Xn是在(0,M)内服从均匀分布的随机变量，Yn则是在(0,1)内服从均匀分布的随机变量。式中X0,A,C,M的取值并不是随意的，模M大小是发生器周期长短的主要标志，常见有M为素数，取A为M的原根，则周期T=M-1。试验统计表明，用以下参数进行混合同余法产生的随机序列的统计特性较好：Xn=mod(2045Xn-1+1,220)（1.3）混合同余法产生的随机序列具有以下特点：?Xn重复周期较小，由于Xn取值在(0,M)内，其周期TM，T32020年4月19日文档仅供参考，不当之处，请
5、联系改正。受X0,A,C,M的值的影响，在编程实现时，浮点运算也会对T产生影响?用此方法产生的随机序列，在一个周期内任意两个随机数不可能相等，这往往与实际情况不相符经Hull和Dobell证明，只有X0,A,C,M满足以下一些关系才能实现周期最大化，即T=M，条件如下：?C与M互质（或互素，即它们的最大公约数为1）?设q为某一质数，M分别能被q和4整除，且(A-1)能被q和4整除产生具有最大周期的伪随机序列的混合同余法算法为：Xnmod((4a+1)Xn-1(2c+1),2k)Yn=Xn/2k（1.4）（1.5）由于M=2k,k2时，M只有
6、一个素数因子2，且4也是M的因子，此时A=4a+1，正好满足了T=M的第二个条件；而此时C=2c+1刚好与M互质，即满足T=M的第一个条件。三、验证方案：根据白噪声定义和特性：其均值为0，功率谱密度为常数，方差等于功率，其自相关函数在0处为一个冲击，因此我们需要验证其均值、功率谱密度和自相关函数即可验证其分布。四、流程图：42020年4月19日文档仅供参考，不当之处，请联系改正。图1五、源程序：a=input('请输入所要指定的功率a=:')%混合同余法产生（-3^0.5到3^0.5的均匀分布n=7340;%随机数长度x=zeros(1,n)
7、;%分配存储空间x(1)=224;%设定种子数y=zeros(1,n-1);fori=2:n%循环产生[0,1]均匀分布的随机数x(i)=mod(2045*x(i-1)+1,2^20);y(i-1)=x(i)/2^20;endy=y*(2*3^0.5)-3^0.5;%对[0,1]均匀分布变换，使其均值为零，方差为1%经过变换生成功率为a的白噪声X=a^(1/2)*y;%功率为a的白噪声（因为均值为0）52020年4月19日
