untitled610.asv

这部分是用MATLAB产生随机过程的一系列源代码.含产生高斯分布随机变量,多变量高斯过程样本的产生,高斯--马尔科夫过程的产生等等..分别运行了这些程序.希望对大家有用

%%带通随机过程抽样的产生
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;
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)