lianghua.m

基于matlab的gmsk信号的调制

for jj2=[1:1024]
 if     fcos(jj2)<1/1023,fcinn(jj2)=0;
 elseif fcos(jj2)>=2045/1023,fcinn(jj2)=1023;
 else
    for jj1=[1:1022]
    if fcos(jj2)>=(2*jj1-1)/1023 & fcos(jj2)<(2*jj1+1)/1023,fcinn(jj2)=jj1;end
    end
 end
end
%matlab中矩阵是按列存取,第一列从1到n,然后第二列从n+1到2n;
%bitget(8,4)=1,bitget(8,1)=0
for jj4=[1:1024]
 for jj3=[1:8]
     fcbin1((jj4-1)*8+jj3)=bitget(fcinn(jj4),jj3);
 end
 for jj5=[9:10]
    fcbin2((jj4-1)*8+jj5-8)=bitget(fcinn(jj4),jj5);
 end
end
%以上为cos的bin文件的整型存储
%以下则是sin的
for jj6=[1:1024]
 if     fsin(jj6)<1/1023,fsinn(jj6)=0;
 elseif fsin(jj6)>=2045/1023,fsinn(jj6)=1023;
 else
    for jj7=[1:1022]
    if fsin(jj6)>=(2*jj7-1)/1023 & fsin(jj6)<(2*jj7+1)/1023,fsinn(jj6)=jj7;end
    end
 end
end %量化成整型的

for jj8=[1:1024]
 for jj9=[1:8]
     fsbin1((jj8-1)*8+jj9)=bitget(fsinn(jj8),jj9);
 end
 for jj10=[9:10]
    fcbin2((jj8-1)*8+jj10-2)=bitget(fsinn(jj8),jj10);
 end
end
%以下是将二进制数据转化成十进制的数据存储
for jj=[1:1024]
    fcbin(jj)=fcbin1((jj-1)*8+1)+2*fcbin1((jj-1)*8+2)+4*fcbin1((jj-1)*8+3)+8*fcbin1((jj-1)*8+4)+16*fcbin1((jj-1)*8+5)+32*fcbin1((jj-1)*8+6)+64*fcbin1((jj-1)*8+7)+128*fcbin1((jj-1)*8+8);
    %fbinm(jj)=fcbin2((jj-1)*8+1)+2*fcbin2((jj-1)*8+2)+4*fcbin2((jj-1)*8+3)+8*fcbin2((jj-1)*8+4)+16*fcbin2((jj-1)*8+5)+32*fcbin2((jj-1)*8+6)+64*fcbin2((jj-1)*8+7)+128*fcbin2((jj-1)*8+8);
    fbinm(jj)=fcbin2((jj-1)*8+1)+2*fcbin2((jj-1)*8+2)+64*fcbin2((jj-1)*8+7)+128*fcbin2((jj-1)*8+8);
    fsbin(jj)=fsbin1((jj-1)*8+1)+2*fsbin1((jj-1)*8+2)+4*fsbin1((jj-1)*8+3)+8*fsbin1((jj-1)*8+4)+16*fsbin1((jj-1)*8+5)+32*fsbin1((jj-1)*8+6)+64*fsbin1((jj-1)*8+7)+128*fsbin1((jj-1)*8+8);
end

fid1=fopen('rom1.bin','w');
fwrite(fid1,fcbin,'uint64');
fclose(fid1);
fid2=fopen('rom2.bin','w');
fwrite(fid2,fbinm,'uint64');
fclose(fid2);
fid3=fopen('rom3.bin','w');
fwrite(fid3,fsbin,'uint64');
fclose(fid3);