1.txt

用Matlab产生一满足条件的窄带随机过程

用Matlab产生一满足条件的窄带随机过程；
y1=wgn(1000,1,2);
y2=wgn(1000,1,2);
windowSize = 5;
a1=filter(ones(1,windowSize)/windowSize,1,y1);
a2=filter(ones(1,windowSize)/windowSize,1,y2);
plot(y1);
title('高斯白噪声');
plot(a1);
title('高斯白噪声通过低通滤波器后');
f=1000;
t=0:4/f:4*pi/f;
m=cos(2*pi*f.*t);
m1=sin(2*pi*f.*t);
y=a2*m1-a1*m;
plot(y);
title('窄带随机过程')

编写子程序计算该随机过程的均值函数，自相关函数，功率谱，包络，包络平方及相位的一维概率密度画出相应的图形并给出解释。
mean(y)
ans = 0.0182    0.0182    0.0182    0.0182
var(y)
ans =0.3363    0.3363    0.3363    0.3363
x=xcorr(y);
plot(x)
title('自相关函数')
f=fft(y);
ab=abs(f);
plot(ab)
title('功率谱密度')
bao=abs(y);
plot(bao);
title('包络')
hist(bao);
title('包络的一维概率密度')
bao1=bao.^2;
plot(bao1); 
title('包络的平方')
hist(bao1);
title('包络平方的一维概率密度')
an=angle(y);
plot(an);
title('相位图形')
hist(an);
title('相位的一维概率密度')即得到自相关函数的图形

自行产生任一随机信号数据（可参考实验一）。（1）绘制该段信号的波形；（2）计算并画图求均值、方差及自相关函数；（3）计算并画图求其功率谱。
r=unifrnd(0,2,1000,1);
>> plot(r);title('信号波形')
g=mean(r)
g =
    1.0016
stem(g) ;title('均值')
h=var(r)
h =
    0.3557
>> stem(h) ;title('方差')
t=xcorr(r);
>> plot(t);title('自相关')
>> y=fft(abs(t));title('功率谱')