untitled610.m

来自「这部分是用MATLAB产生随机过程的一系列源代码.含产生高斯分布随机变量,多变量」· M 代码 · 共 22 行

%%带通随机过程抽样的产生
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);
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是假设的（参考P137说明）
M=50;
% 确定带通过程的自相关系数和频谱
bpp_autocorr=Rx_est(band_pass_process,M); 
bpp_spectrum=fftshift(abs(fft(bpp_autocorr)));
% 绘图
f=[0:1:50];
plot(f,bpp_spectrum) % 如果用stem来画图，得到的是离散的图形