function [signal_out,I_out,Q_out,phase,gt,qt] = mod_gmsk(data,data_len,sample_number,Rb,alpha)
%GMSK基帶調制
%**************************************************************************
% data              調制信號
% data_len          碼元個數
% sample_number     每個碼元采樣點數
% Rb                碼元速率
% signal_out        基帶調制輸出
% I_out             I路輸出
% Q_out             Q路輸出
% alpha             alpha = Bb*Tb = Bb/Rb
%**************************************************************************

% data_len = 10;                %碼元個數
% sample_number = 8;            %采樣點數
% Rb = 16000;                   %碼元速率
% alpha = 0.25;                 %BbTb值
% data1 = randint(1,data_len);  
% data = 2*data1-1;             %傳輸的序列

Tb = 1/Rb;                      %碼元時間
Bb = alpha/Tb;                  %調制高斯濾波器3dB帶寬
fs = Rb*sample_number;          %采樣速率

%--------------------------------------------------------------------------
%高斯濾波器
[gt,qt] = gauss_filter(Bb,Tb,fs,sample_number);
%**************************************************************************

%--------------------------------------------------------------------------
%相位
phase = gmsk_phase(data,data_len,sample_number,qt);
%**************************************************************************

%--------------------------------------------------------------------------
%調制信號產生
I_out = cos(phase);
Q_out = sin(phase);
signal_out = I_out + j*Q_out;
%**************************************************************************