利用混合乘同余法产生(0,1)均匀分布随机数.doc

概念 混合同余法： 用一常数与已产生的伪随机数之积的同余式做递推公式,产生均匀伪随机数的方法。 即： 第一步：递推式 其中： M=2^k,k2整数；A=1(mod4),即 为正整数,初值为非负整数。 第二步： 可以证明， 是周期为2k的伪随机数。 混合同余法（0，1）均匀分布的随机数 L=60,取k=2,c=1,n=5,l=60,=1 公式原理： 见上第一第二步。 程序： clear clc k=8; c=1; n=5; x=1; M=2^k; A=2^n+1; store=[ ]; for i=1:60 ee=x/M; store(i)=ee; x=mod(A*x+c,M); end plot(1:60,store) title(混合同余法产生随机数) store= 0.00390625000000000 0.132812500000000 0.386718750000000 0.765625000000000 0.269531250000000 0.898437500000000 0.652343750000000 0.531250000000000 0.535156250000000 0.664062500000000 0.917968750000000 0.296875000000000 0.800781250000000 0.429687500000000 0.183593750000000 0.0625000000000000 0.0664062500000000 0.195312500000000 0.449218750000000 0.828125000000000 0.332031250000000 0.960937500000000 0.714843750000000 0.593750000000000 0.597656250000000 0.726562500000000 0.980468750000000 0.359375000000000 0.863281250000000 0.492187500000000 0.246093750000000 0.125000000000000 0.128906250000000 0.257812500000000 0.511718750000000 0.890625000000000 0.394531250000000 0.0234375000000000 0.777343750000000 0.656250000000000 0.660156250000000 0.789062500000000 0.0429687500000000 0.421875000000000 0.925781250000000 0.554687500000000 0.308593750000000 0.187500000000000 0.191406250000000 0.320312500000000 0.574218750000000 0.953125000000000 0.457031250000000 0.0859375000000000 0.839843750000000 0.718750000000000 0.722656250000000 0.851562500000000 0.105468750000000 0.484375000000000 第二题：LS 递推，L=15，输入4位M序列，v（K)为均值c的白噪声序列， z(k)-2.1z(k-1)+1.2z(k-2)=u(k-1)+0.4u(k-2)+v(k) 幅值：0.03 程序： %FLch3RLSeg3 clear%清理工作间变量 L=15;% M序列的周期 y1=1;y2=1;y3=1;y4=0;%四个移位积存器的输出初始值 for i=1:L;%开始循环，长度为L x1=xor(y3,y4);%第一个移位积存器的输入是第3个与第4个移位积存器的输出的“或” x2=y1;%第二个移位积存器的输入是第3个移位积存器的输出 x3=y2;%第三个移位积存器的输入是第2个移位积存器的输出 x4=y3;%第四个移位积存器的输入是第3个移位积存器的输出 y(i)=y4;%取出第四个移位积存器幅值为0和1的输出信号， if y(i)0.5,u(i)=-0.03;%如果M序列的值为1时,辨识的输入信号取“-0.03” else u(i)=0.03;%当M序列的值为0时,辨识的输入信号取“0.03” end%小循环结束 y1=x1;y2=x2;y3=x3;y4=x4;%为下一次的输入信号做准