[随机信号实验1.docxVIP

1、熟悉并练习使用下列Matlab 的函数，给出各个函数的功能说明和内部参数的意义，并给出至少一个使用例子和运行结果：1） rand(m，n) 生成m×n 随机矩阵，其元素在（0，1）内例： y=rand(2,3)y = 0.9649 0.9706 0.48540.1576 0.9572 0.80032）randn(m，n) 产生随机数数组或矩阵，其元素服从均值为0，方差为1的正态分布例：Y = randn(3,3)Y = -0.2620 -0.8314 -0.5336 -1.7502 -0.9792 -2.0026 -0.2857 -1.1564 0.96423） normrnd()（1）y = normrnd(m,n，a，b)产生均值为m，标准差为n的正态随机过程，a和b是y的维数。例：y= normrnd(2,1,3,4)y = 3.4367 0.7922 3.3790 1.7275 0.0391 4.9080 0.9418 3.0984 1.8023 2.8252 1.5314 1.72214） y= mean(A) A的均值。例：A = [2 2 3; 3 4 6; 4 5 8; 3 9 7];M=mean(A)M = 3 5 65）var()求方差例：X=[1:1:5;1:2:10];V=var(X,1)V = 0 0.2500 1.0000 2.2500 4.00006）xcorr(x，y)计算x,y的互相关，当x=y时，计算的则是自相关。例：x=normrnd(3,2,1,2); y=normrnd(3,2,1,2);z=xcorr(x,y)z = 2.2262 7.1017 3.35917）periodogram(x) 计算x的功率谱密度例： X=[-30:2:30];Y=periodogram(X);plot(Y)8）fft(x，n)离散傅里叶变换用快速傅里叶算法得到的离散傅里叶变换，返回矩阵每一列的傅里叶变换返回n点的离散傅里叶变换，如果X的长度小于n，X的末尾填零。如果X的长度大于n，则X被截断。当X是一个矩阵时，列的长度也服从同样的操作。X=[0:0.2:1];Y = fft(X,5)Y = 2.0000 -0.5000 + 0.6882i -0.5000 + 0.1625i -0.5000 - 0.1625i -0.5000 - 0.6882i9）normpdf(x，m，a)求正态分布概率密度函数值，m为均值，a为方差。x=[-8:0.1:8];y=normpdf(x,1,2);plot(x,y)（9图）（10图）10）normcdf(x，m，a)求正态分布概率分布函数值，参数为m,ax = (-6:0.1:10); y = normcdf(x,0,1);plot(x,y)11）unifpdf(x，a，b)参数为a,b的均匀分布函数值x = [0:0.1:2];y = unifpdf(x,1,2)y = 0 0 0 0 0 0 0 0 0 0 1 1 1 1 1 1 1 1 1 1 112）unifcdf(x，a，b)求连续均匀分布的概率分布函数值,参数为a,b的均匀分布累计分布函数值 x=[0;0.1;1] ; y=unifcdf(x,-1,1);plot(x,y)（12图）（13图）13）raylpdf(x，a)求瑞利概率密度分布函数值,参数为a x = [-0.5:0.1:8];p = raylpdf(x,2);plot(x,p)14）raylcdf(x，a)求瑞利分布的概率分布函数值,参数为a. x = [0:0.1:8];p = raylcdf(x,2);plot(x,p)（14图）（15图）15）exppdf(x，m)求指数分布的概率密度函数值x=[-1:0.1:8]; y = exppdf(x,2);plot(x,y)16）expcdf(x，m)求指数分布的概率分布函数值，m为均值x = [0:0.2:6];y = expcdf(x,1);plot(x,y)（16图）（18图）17）chol()对称正定矩阵的Cholesky分解（1）R=chol(X) 产生一个上三角阵R，使RR=X。若X为非对称正定，则输出一个出错