%++++++++++++++++++++++++++++++++++++++++++++++++++++++
%	Filename:	gmsk.m
%
%	Description:	Generate baseband GMSK signal
%
% Author: Jeff Laster, MPRG, Virginia Tech 
%++++++++++++++++++++++++++++++++++++++++++++++++++++++

t0 = clock;

%++++++++++++++++++++++++++++++++
% 	Constant Definition
%++++++++++++++++++++++++++++++++

DRate = 1;			% data rate
M = 18;				% no. of sample per bit

%++++++++++++++++++++++++++++++++
%	Variable Definition
%++++++++++++++++++++++++++++++++

N = 12753;			% no. of bits for simulation
%EbNodB = 10;			% Eb/No (dB scale)
BT = 0.3;			% Bandwidth*Period (cannot change )

%++++++++++++++++++++++++++++++++++
% 	Program Control Variables
%++++++++++++++++++++++++++++++++++

fig_no = 1;			% initialize fig no.

%+++++++++++++++++++++++++++++++++++
%	Derived variables
%+++++++++++++++++++++++++++++++++++

T = 1/DRate;				% data period
Ts = T/M;				% sampling period ( discrete-time )
%EbNo = 10^(EbNodB/10);			% Eb/No (linear scale )

%+++++++++++++++++++++++++++++++++++++
%	Generates Input Random Data
%+++++++++++++++++++++++++++++++++++++
[data, RandData] = pn_gen(N,M);

%dodata

%+++++++++++++++++++++++++++++++++++++
%	Generate GMSK Pulse Response
%+++++++++++++++++++++++++++++++++++++
[g,fig_no] = pulse_re(BT,T,Ts,M,fig_no);

%+++++++++++++++++++++++++++++++++++++++++
%	Generate Gaussian Filtered Pulses
%+++++++++++++++++++++++++++++++++++++++++
[b,delay_GF,fig_no] = gauss_pu(g,data,1,fig_no);

%+++++++++++++++++++++++++++++++++++++++++
%	Generate GMSK signal
%+++++++++++++++++++++++++++++++++++++++++
[s,fig_no] = gmsk_sig(b,Ts,M,fig_no);

clear T DRate BT

%data=zeros(size(data));
%for i=1:N
%	for k=1:M
%		data(i*M+M/2-1+k)=Rdiff_mse(i);
%	end
%end

Time = etime(clock,t0)