fi.m

基于matlab的gmsk信号的调制

global dt df t f N;
close all
clear
T=1/270833;
%h=0.5;
B=0.3/T;
a=sqrt(log(2)/2)/B;
t=[-2.5*T:5*T/100:2.5*T];
g1=0.5*erfc(pi*(t+T/2)/a);
g2=0.5*erfc(pi*(t-T/2)/a);
g_fun=(g2-g1)/(2*T);
bb=zeros(1,5);
sum=zeros(32,8);
aa=[1,1,1,1,1;
    1,1,1,1,-1;
    1,1,1,-1,-1;
    1,1,1,-1,1;
    1,1,-1,-1,-1;
    1,1,-1,-1,1;
    1,1,-1,1,-1;
    1,1,-1,1,1;
    1,-1,1,1,1;
    1,-1,1,1,-1;
    1,-1,1,-1,1;
    1,-1,1,-1,-1;
    1,-1,-1,1,1;
    1,-1,-1,1,-1;
    1,-1,-1,-1,1;
    1,-1,-1,-1,-1;
    -1,1,1,1,1;
    -1,1,1,1,-1;
    -1,1,1,-1,-1;
    -1,1,1,-1,1;
    -1,1,-1,-1,-1;
    -1,1,-1,-1,1;
    -1,1,-1,1,-1;
    -1,1,-1,1,1;
    -1,-1,1,1,1;
    -1,-1,1,1,-1;
    -1,-1,1,-1,1;
    -1,-1,1,-1,-1;
    -1,-1,-1,1,1;
    -1,-1,-1,1,-1;
    -1,-1,-1,-1,1;
    -1,-1,-1,-1,-1;
];
cc=zeros(1,5);
dd=qua(0.5*T)
for j=0:1:3
  for i=1:1:32
   for n=0:1:7
    bb(1)=qua(1.5*T+T*n/8);
    bb(2)=qua(0.5*T+T*n/8);
    bb(3)=qua(-0.5*T+T*n/8);
    bb(4)=qua(-1.5*T+T*n/8);
    bb(5)=qua(-2.5*T+T*n/8);
   %bb(1)=quad(@g,-2.5*T,1.5*T+T*n/8);bb(2)=quad(@g,-2.5*T,0.5*T+T*n/8);bb(3)=quad(@g,-2.5*T,-0.5*T+T*n/8);bb(4)=quad(@g,-2.5*T,-1.5*T+T*n/8);bb(5)=quad(@g,-2.5*T,-2.5*T+T*n/8);
    cc=bb.*aa(i, : );
    sum(i+j*32,n+1)=((cc(1)+cc(2)+cc(3)+cc(4)+cc(5))+j/2)*pi;
   %sum(n+1)=bb(1)+bb(2)+bb(3)+bb(4)+bb(5)
    end
  end
 end 
sum
end