fabao3.m

一个有关GMSK基带调制和实现的MATLAB程序

a1=b1';
a2=b2';
z1=reshape(a1,1,128*8);
z1=z1+1;
z1=z1*512;
z2=reshape(a2,1,128*8);
z2=z2+1;
z2=z2*512;
x1=zeros(1,128*8);
x2=zeros(1,128*8);
z1=z1/1024*1023;
z1=round(z1);
z2=z2/1024*1023;
z2=round(z2);
for i=1:1024,
    if z1(i)>=768,z1(i)=z1(i)-768;x1(i)=768;
    elseif z1(i)>=512,z1(i)=z1(i)-512;x1(i)=512;
        elseif z1(i)>=256,z1(i)=z1(i)-256;x1(i)=256;
             else z1(i)=z1(i);x1(i)=x1(i);
             end
         end
         for i=1:1024,
    if z2(i)>=768,z2(i)=z2(i)-768;x2(i)=768;
    elseif z2(i)>=512,z2(i)=z2(i)-512;x2(i)=512;
        elseif z2(i)>=256,z2(i)=z2(i)-256;x2(i)=256;
             else z2(i)=z2(i);x2(i)=x2(i);
             end
         end
 x1=x1/256;
 x2=x2/4;
 x=x1+x2;
 fid1=fopen('z1.bin','w')
fwrite(fid1,z1,'ubit64')
fclose(fid1)
fid2=fopen('z2.bin','w')
fwrite(fid2,z2,'ubit64')
fclose(fid2)
 fid3=fopen('x.bin','w')
fwrite(fid3,x,'ubit64')
fclose(fid3)