function y = mod_gmsk_030610 (x,n,M)

global k Fc Fd Fs Fs_n;

% x 為長度為1024*5=5120點的數字信號，每個碼元里有Fs_n/Fd＝120個數，共42+2/3個碼元。

% 首先將基帶數字信號 x 變成雙極性信號（非零信號）

x=-1+2*x;   %  “0”——“－1”  “1”——“1”

% 再將雙極性基帶信號通過高斯低通濾波器。


B=0.3*Fd;
a=sqrt(log(2)/2)/B;
t=n; 
h=sqrt(pi)/a*exp(-(pi/a*t).^2);

y0=fftfilt(h,x);
y0=y0/max(y0);

% GMSK just like FSK2 Modulation Procedure, 
% M == 2
Tone = Fd/2 ;
y = cos( 2*pi* (Fc-Tone/2 + Tone*y0').*n );   % -Fd/4, +Fd/4, h=0.5; GMSK