gmsk_mod.m

GSM全速率语音突发的的发送端总程序

%
% gmsk_mod:   This function accepts a GSM burst bit sequence and
%             performs a GMSK modulation of the sequence. The
%             modulation is according to the GSM 05.05 recommendations
%
% SYNTAX:     [i,q] = gmsk_mod(burst,Tb,osr,BT)
%
% INPUT:      burst   A differential encoded bit sequence (-1,+1)
%             Tb      Bit duration (GSM: Tb = 3.692e-6 Sec.)
%             osr     Simulation oversample ratio. osr determines the
%                     number of simulation steps per information bit
%             BT      The bandwidth/bit duration product (GSM: BT = 0.3)
%
% OUTPUT:     i,q     In-phase (i) and quadrature-phase (q) baseband
%                     representation of the GMSK modulated input burst 
%                     sequence
%
% SUB_FUNC:   ph_g.m  This sub-function is required in generating the
%                     frequency and phase pulse functions.
%
% WARNINGS:   Sub-function ph_g.m assumes a 3xTb frequency pulse 
%             truncation time
%
% TEST(S):    Function tested using the following relations
%             
%             i)
%
%             I^2 + Q^2 = Cos(a)^2 + Sin(a)^2 = 1
%
%             ii)
%
%             When the input consists of all 1's the resulting baseband
% 	      outputs the function should return a sinusoidal signal of
%             frequency rb/4, i.e. a signal having a periode time of 
%             approximately 4*Tb = 4*3.692e-6 s = 1.48e-5 s for GSM
%
% AUTHOR:   Jan H. Mikkelsen / Arne Norre Ekstr鴐
% EMAIL:    [e-mail] / [e-mail]
%
% $Id: gmsk_mod.m,v 1.5 1998/02/12 10:50:10 aneks Exp $

% ACCUIRE GMSK FREQUENCY PULSE AND PHASE FUNCTION
%
BURST = [0 0 0 0 1 0 1 0 1 0 0 0 1 0 1 1 1 1 1 0 1 0 1 0 0 0 1 0 1 1 1 1 1 0 1 0 1 0 0 0 1 0 1 1 1 1 0 0 1 0 1 0 0 0 1 0 1 1 1 0 1 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1 0 0 1 0 1 0 0 0 0 0 1 0 1 1 1 0 1 0 1 0 0 0 0 0 1 1 1 0 1 0 1 0 1 0 1 0 0 0 1 0 1 0 1 0 1 0 1 0 1 0 0 0 1 0 1 0 1 0 0 0];
DIFF_ENC_DATA = diff_enc(BURST)
[g,q] = ph_g

% PREPARE VECTOR FOR DATA PROCESSING
%
BURST = DIFF_ENC_DATA;
bits = length(BURST);
OSR = 4;
f_res = zeros(1,(bits+2)*OSR);

% GENERATE RESULTING FREQUENCY PULSE SEQUENCE
%
for n = 1:bits,
  f_res((n-1)*OSR+1:(n+2)*OSR) = f_res((n-1)*OSR+1:(n+2)*OSR) + BURST(n).*g;
end

% CALCULATE RESULTING PHASE FUNCTION
%
theta = cumsum(f_res);
theta = pi * theta
% PREPARE DATA FOR OUTPUT
%
I = cos(theta);
Q = sin(theta);

Qbig = [0 -804 -3212 -7962 -15446 -22594 -27466 -28898 -27466 -22005 -12539 0 30273 23170 12539 0 30273 23170 12539 0 30273 23170 12539 0 30273 24279 17869 15446 27466 22005 12539 0 -30273 -24279 -17869 -15446 -27466 -22005 -12539 0 30273 23170 12539 0 30273 23170 12539 0 30273 23170 12539 0 30273 23170 12539 0 30273 23170 12539 0 30273 23170 12539 0 30273 23170 12539 0 30273 23170 12539 0 30273 24279 17869 15446 27466 22005 12539 0 -30273 -24279 -17869 -15446 -27466 -22005 -12539 0 30273 23170 12539 0 30273 23170 12539 0 30273 23170 12539 0 30273 23170 12539 0 30273 24279 17869 15446 27466 22005 12539 0 -30273 -24279 -17869 -15446 -27466 -22005 -12539 0 30273 23170 12539 0 30273 23170 12539 0 30273 24279 17869 15446 27466 22005 12539 0 -30273 -24279 -17869 -15446 -27466 -22005 -12539 0 30273 23170 12539 0 30273 23170 12539 0 30273 23170 12539 0 30273 23170 12539 0 30273 24279 17869 15446 27466 22005 12539 0 -30273 -24279 -17869 -15446 -27466 -22005 -12539 0 30273 24279 17869 15446 27466 22005 12539 0 -30273 -23170 -12539 0 -30273 -23170 -12539 0 -30273 -23170 -12539 0 -30273 -23170 -12539 0 -30273 -23170 -12539 0 -30273 -23170 -12539 0 -30273 -24279 -17869 -15446 -27466 -22005 -12539 0 30273 23170 12539 0 30273 23170 12539 0 30273 23170 12539 0 30273 23170 12539 0 30273 23170 12539 0 30273 23170 12539 0 30273 24279 17869 15446 27466 22005 12539 0 -30273 -24279 -17869 -15446 -27466 -22005 -12539 0 30273 23170 12539 0 30273 23170 12539 0 30273 23170 12539 0 30273 23170 12539 0 30273 23170 12539 0 30273 23170 12539 0 30273 23170 12539 0 30273 23170 12539 0 30273 23170 12539 0 30273 23170 12539 0 30273 23170 12539 0 30273 23170 12539 0 30273 23170 12539 0 30273 23170 12539 0 30273 23170 12539 0 30273 23170 12539 0 30273 23170 12539 0 30273 23170 12539 0 30273 23170 12539 0 30273 23170 12539 0 30273 23170 12539 0 30273 23170 12539 0 30273 23170 12539 0 30273 24279 17869 15446 27466 22005 12539 0 -30273 -24279 -17869 -15446 -27466 -22005 -12539 0 30273 23170 12539 0 30273 24279 17869 15446 27466 22005 12539 0 -30273 -24279 -17869 -15446 -27466 -22005 -12539 0 30273 23170 12539 0 30273 23170 12539 0 30273 24279 17869 15446 27466 22005 12539 0 -30273 -24279 -17869 -15446 -27466 -22005 -12539 0 30273 23170 12539 0 30273 23170 12539 0 30273 23170 12539 0 30273 23170 12539 0 30273 24279 17869 15446 27466 22005 12539 0 -30273 -24279 -17869 -15446 -27466 -22005 -12539 0 30273 23170 12539 0 30273 23170 12539 0 30273 23170 12539 0 30273 23170 12539 0 30273 23170 12539 0 30273 23170 12539 0 30273 23170 12539 0 30273 23170 12539 0 30273 24279 17869 15446 27466 22005 12539 0 -30273 -23170 -12539 0 -30273 -23170 -12539 0 -30273 -23170 -12539 0 -30273 -23170 -12539 0 -30273 -24279 -17869 -15446 -27466 -22005 -12539 0 30273 24279 17869 15446 27466 22005 12539 0 -30273 -24279 -17869 -15446 -27466 -22005 -12539 0 30273 23170 12539 0 30273 23170 12539 0 30273 24279 17869 15446 27466 22005 12539 0 -30273 -24279 -17869 -15446 -27466 -22005 -12539 0 30273 23170 12539 0 30273 23170 12539 0 30273 23170 12539 0 30273 23170 12539 0 30273 24279 17869 15446 27466 22005 12539 0 -30273 -23170 -12539 0 -30273 -23170 -12539 0 -30273 -23170 -12539 0 -30273 -23170 -12539 0 -30273 -24279 -17869 -15446 -27466 -22005 -12539 0 30273 23731 15446 7962 15446 23170 28898 31785]

Qcom = Qbig / 32768