samp9_1.m

来自「数字信号处理的Matlab实现的随书光盘」· M 代码 · 共 15 行

%Samp9_1
N=100000;   %数据个数
randn('state',0);   %设置产生随机数的状态
y=randn(1,N);  %产生一个随机序列
disp('平均值:');
yMean=mean(y)   %计算随机序列的均值
disp('平方值:');
y2p=y*y'/N    %计算其均方值,这里利用了矩阵相乘的算法
disp('平方根:');
ysq=sqrt(y2p)     %计算其均方根值
disp('标准差:');
ystd=std(y,1)   %相当于ystd=sqrt(sum((y-yMean).^2)/(N-1))
disp('方差:');
yd=ystd.*ystd