%第6章例10
N=1000; %抽样的数目
for i=1:2:N;
[X1(i) X1(i+1)]=gngauss;
[X2(i) X2(i+1)]=gngauss;
end; %标准高斯输入噪声过程
A=[1 -0.9]; %低通滤波器参数
B=1;
Xc=filter(B,A,X1); %filter(B,A,X) 一维数字滤波器,X输入,B、A分别为滤波器系数
Xs=filter(B,A,X2);
fc=1000/pi; %载波频率
for i=1:N;
band_pass_process(i)=Xc(i)*cos(2*pi*fc*i)-Xs(i)*sin(2*pi*fc*i);
end; %设T=1;
%确定带通过程的自相关系数和频谱
M=50;
bpp_autocorr=Rx_est(band_pass_process,M);
bbp_spectrum=fftshift(abs(fft(bpp_autocorr)));
%绘图
plot(bbp_spectrum);
xlabel('f');title('功率谱');
figure(2),subplot(1,2,1),plot(X1);title('X1');
subplot(1,2,2),plot(Xc);title('Xc');
figure(3),subplot(1,2,1),plot(X2);title('X2');
subplot(1,2,2),plot(Xs);title('Xs');