function [gt,qt]=gauss_filter(Bb,Tb,fs,sample_number)
%高斯濾波器
%**************************************************************************
%Bb             濾波器的3dB帶寬
%Tb             碼元時間
%fs             采樣速率
%sample_number  采樣個數(shù)
%**************************************************************************

% Rb = 16000;
% Tb = 1/16000;
% sample_number = 8;
% fs = sample_number*Rb;
% Bb = 0.25/Tb;
%--------------------------------------------------------------------------
%g(t)
irfn = 3;                        %截取irfn個Tb區(qū)間
n = irfn * sample_number;
mid = (n./2);
t = -mid/fs:1/fs:mid/fs;
gt = 1/2*(erf(-sqrt(2/log(2))*pi*Bb*(t-1/2*Tb))+erf(sqrt(2/log(2))*pi*Bb*(t+1/2*Tb)));
%**************************************************************************

%--------------------------------------------------------------------------
%q(t)
qt=zeros(1,irfn*sample_number);
for i=1:irfn*sample_number
    for j=1:i
       qt(i)=qt(i)+gt(j)/sample_number/2;
    end;
end;
%**************************************************************************