chouyanglianghua.m

基于matlab的gmsk信号的调制

%function f=f(t)
T=1/270833;
N=8;
dt=T/N;
t=[4*T:0.125*T:4.875*T];                                        
f=zeros(1,1024); 
close all;
%hold on
%for i=1:50                                              %画图50次
%b=sign(randn(1,2000));
b=[-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;];
for  L=0:3
   for bb=1:32
       for n=1:8
       k=4;                                             
       f1=pi*b(bb,k-3)*quad(@g,-2.5*T, 1.5*T+t(n)-k*T);    
       f2=pi*b(bb,k-2)*quad(@g,-2.5*T, 0.5*T+t(n)-k*T);
       f3=pi*b(bb,k-1)*quad(@g,-2.5*T,-0.5*T+t(n)-k*T);
       f4=pi*b(bb,k  )*quad(@g,-2.5*T,-1.5*T+t(n)-k*T);
       f5=pi*b(bb,k+1)*quad(@g,-2.5*T,-2.5*T+t(n)-k*T);
       f(n+(bb-1)*32+L*8)=f1+f2+f3+f4+f5+L*pi/2;
       end
   end
end

s=sin(f);
c=cos(f);

partition=[-1+2/1024:2/1024:1-2/1024];%量化
cutecos=quantiz(c,partition);
reshape(cutecos,16,64);
cutesin=quantiz(s,partition);
reshape(cutesin,16,64);
%------------------------------------高2低8位分开------------------------------------------------

    low_s=mod(cutesin,256);
    high_s=floor(cutesin/256);
    low_c=mod(cutecos,256);
    high_c=floor(cutecos/256);
    high=high_s*64+high_c;
%－－－－－－－－－－－－－－－－－－－－－－－－－－－－－－－－－－－－－－－－－－－－－－－－－－－    
low_s_rom=zeros(8192,1);
low_c_rom=zeros(8192,1);
team_rom =zeros(8192,1);
for n1=1:1024
     low_s_rom(n1*8-7,1)=low_s(n1);
     low_c_rom(n1*8-7,1)=low_c(n1);
     team_rom(n1*8-7,1)=high(n1);
  
end
 
    fid = fopen('low_s_rom.bin', 'wb');
    fwrite(fid,low_s_rom);    
    
    fid = fopen('low_c_rom.bin', 'Wb');
    fwrite(fid,low_c_rom);  
    
    fid = fopen('team_rom.bin', 'Wb');
    fwrite(fid,team_rom);