matlab 窄带高斯过程仿真,通信原理基于matlab的计算机仿真_源代码

《通信原理：基于matlab的计算机仿真》书上的实例源代码。。

[ff,glf]=T2F( tt, glt );

figure(2)

plot(ff,abs(glf));

xlabel('频率(Hz)'); ylabel('窄带高斯过程样本的幅频特性')

figure(3)

subplot(411);

plot(tt,real(gt));

title('窄带高斯过程样本')

subplot(412)

plot(tt,real(glt).*cos(2*pi*fc*tt)-imag(glt).*sin(2*pi*fc*tt))

title('由等效基带重构的窄带高斯过程样本')

subplot(413)

plot(tt,real(glt));

title('窄带高斯过程样本的同相分量')

subplot(414)

plot(tt,imag(glt));

xlabel('时间t(秒)'); title('窄带高斯过程样本的正交分量')

%求窄带高斯信号功率;注：由于样本的功率近似等于随机过程的功率，因此可能出现一些偏差

P_gt=sum(real(gt).^2)/T;

P_glt_real = sum(real(glt).^2)/T;

P_glt_imag = sum(imag(glt).^2)/T;

%验证窄带高斯过程的同相分量、正交分量的正交性

a = real(glt)*(imag(glt))'/T;

用到的子函数

function [t,st]=bpf(f,sf,B1,B2)

%This function filter an input at frequency domain by an ideal bandpass filter %Inputs:

% f: frequency samples

% sf: input data spectrum samples

% B1: bandpass's lower frequency

% B2: bandpass's higher frequency

%Outputs:

% t: frequency samples

% st: output data's time samples

df = f(2)-f(1);

T = 1/df;

hf = zeros(1,length(f));

bf = [floor( B1/df ): floor( B2/df )] ;