matlab中rand函数是产生0到1的随机分布,matlab中randn函数是产生标准正态分布,randint是产生整数随机数,默认为0和1 >> rand(3) ans = 0.8147 0.9134 0.2785 0.9058 0.6324 0.5469 0.1270 0.0975 0.9575 >> randn(3) ans = -0.4326 0.2877 1.1892 -1.6656 -1.1465 -0.0376 0.1253 1.1909 0.3273 >> randint(3) ans = 1 1 0 0 0 0 1 1 1 也可直接用randn函数产生高斯分布序列,例如: 程序代码:
y=randn(1,2500);
y=y-mean(y);
y=y/std(y);
a=0.0128;
b=sqrt(0.9596);
y=a+b*y;
1) rand产生的是[0,1]上的均匀分布的随机序列 2) randn产生均值为0,方差为1的高斯随机序列,也就是白噪声序列; 也就是说,可以直接使用上面两个函数对原始信号添加噪声(例如y=x+rand(length(x),1)或者y=x+randn(length(x),1)) ============================================================= 产生指定方差和均值的随机数: 设某个随机变量x均值为mu,方差为var^2,若要产生同样分布的随机变量y,但使新的随机变量参数改变,均值为mu_1,方差为var_1^2,可以用如下公式进行变换:
- y=var_1/var*(x-mu)+mu_1,其中x为随机变量,其余为常数(原分布参数)。具体到正态分布,若要产生均值为u,方差为o^2的M*N的随机数矩阵,可以用y=o*randn(M,N)+u得到。
- 对于均匀分布,若要产生[a,b]区间的均匀分布的M*N的随机数矩阵,则可以用y=rand(M,N)*(b-a)+a得到。
- y = wgn(m,n,p) 产生一个m行n列的高斯白噪声的矩阵,p以dBW为单位指定输出噪声的强度。
- y = wgn(m,n,p,imp) 以欧姆(Ohm)为单位指定负载阻抗。
- y = wgn(m,n,p,imp,state) 重置RANDN的状态。
- y = wgn(…,POWERTYPE) 指定p的单位。POWERTYPE可以是'dBW', 'dBm'或'linear'。线性强度(linearpower)以瓦特(Watt)为单位。
- y = wgn(…,OUTPUTTYPE) 指定输出类型。OUTPUTTYPE可以是'real'或'complex'。
- y = awgn(x,SNR) 在信号x中加入高斯白噪声。信噪比SNR以dB为单位。x的强度假定为0dBW。如果x是复数,就加入复噪声。
- y = awgn(x,SNR,SIGPOWER) 如果SIGPOWER是数值,则其代表以dBW为单位的信号强度;如果SIGPOWER为'measured',则函数将在加入噪声之前测定信号强度。
- y = awgn(x,SNR,SIGPOWER,STATE) 重置RANDN的状态。
- y = awgn(…,POWERTYPE) 指定SNR和SIGPOWER的单位。POWERTYPE可以是'dB'或'linear'。如果POWERTYPE是'dB',那么SNR以dB为单位,而SIGPOWER以dBW为单位。如果POWERTYPE是'linear',那么SNR作为比值来度量,而SIGPOWER以瓦特为单位
转载于:https://my.oschina.net/itfanr/blog/195662