chouyanglianghua.asv

基于matlab的gmsk信号的调制

%function f=f(t)
T=1/270833;
N=8;
dt=T/N;
t=[4*T:0.125*T:4.875*T];                                        
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;   %5*32的矩阵，代表b的可能取值
   -1,-1, 1,-1,-1;   -1,-1, 1,-1, 1;   -1,-1, 1, 1,-1;   -1,-1, 1, 1, 1;   %-1代表b取之为0
   -1, 1,-1,-1,-1;   -1, 1,-1,-1, 1;   -1, 1,-1, 1,-1;   -1, 1,-1, 1, 1;
   -1, 1, 1,-1,-1;   -1, 1, 1,-1, 1;   -1, 1, 1, 1,-1;   -1, 1, 1, 1, 1;
    1,-1,-1,-1,-1;    1,-1,-1,-1, 1;    1,-1,-1, 1,-1;    1,-1,-1, 1, 1;
    1,-1, 1,-1,-1;    1,-1, 1,-1, 1;    1,-1, 1, 1,-1;    1,-1, 1, 1, 1;
    1, 1,-1,-1,-1;    1, 1,-1,-1, 1;    1, 1,-1, 1,-1;    1, 1,-1, 1, 1;
    1, 1, 1,-1,-1;    1, 1, 1,-1, 1;    1, 1, 1, 1,-1;    1, 1, 1, 1, 1;    ];
f=zeros(1,1024); 
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);    %对fi积分，取5个区间内的值
       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];%1023个区间
cutecos=quantiz(c,partition);         %量化余弦值
cutesin=quantiz(s,partition);         %量化正弦值

%------------------------------------高2低8位分开------------------------------------------------

    low_s=mod(cutesin,256);                       %低8位正弦值表
    high_s=floor(cutesin/256);                    %高2位正弦值表
    low_c=mod(cutecos,256);                       %低8位余弦值表
    high_c=floor(cutecos/256);                    %高2位余弦值表
    high=high_s*64+high_c;                        %较高位上芯片值，前两位代表高2位正弦值。后两位代表高2位余弦值
%－－－－－－－－－－－－－－－－－－－－－－－－－－－－－－－－－－－－－－－－－－－－－－－－－－－    
low_s_rom=zeros(8192,1);                          %芯片上的存储空间为8K
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');                %生成bin文件
    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);