function y=g(t)
%global a T B
%B為3db帶寬
%GMSK中B*T=0.3
%T為碼元周期
%a是決定filter性能的參數(shù)
T=1/270833;%s
B=0.3/T;%t=[-2.5/270833:0.0000001:2.5/270833];
a=sqrt(log(2)/2)/B;
y=1/4/T*erfc(pi/a*(t-T/2))-1/4/T*erfc(pi/a*(t+T/2));
%example is as follows:
%t=[-0.00000923:0.00000001:0.00000923];
%g=1/4/0.000003692*erfc(pi/0.000007245*(t-0.00000134))-1/4/0.000003692*erfc(pi/0.000007245*(t+0.00000134));
%plot(t,y)
%quad(@g,-2.5*T,2.5*T)
end