fabao2.m

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

close all
T=0.000003692;
 t=[-2.5*T:0.0001*T:2.5*T];
 fun=inline('normcdf(6.1323e5*t+6.1323e5*1.846e-6)*1.3543e5-normcdf(6.1323e5*t-6.1323e5*1.846e-6)*1.3543e5','t')
 su=quadl(fun,-2.5*T,2.5*T);
  q=zeros(5,8);
for i=0:7, q(1,i+1)=quadl(fun,-2.5*T,1.5*T+i*T/8);end,
for i=0:7, q(2,i+1)=quadl(fun,-2.5*T,0.5*T+i*T/8);end,
for i=0:7, q(3,i+1)=quadl(fun,-2.5*T,-0.5*T+i*T/8);end,
for i=0:7, q(4,i+1)=quadl(fun,-2.5*T,-1.5*T+i*T/8);end,
for i=0:7, q(5,i+1)=quadl(fun,-2.5*T,-2.5*T+i*T/8);end,
w=[-1,-1,-1,-1,-1;-1,-1,-1,-1,1;-1,-1,-1,1,-1;-1,-1,-1,1,1;-1,-1,1,-1,-1;-1,-1,1,-1,1;-1,-1,1,1,-1;-1,-1,1,1,1;-1,1,-1,-1,-1;-1,1,-1,-1,1;-1,1,-1,1,-1;-1,1,-1,1,1;-1,1,1,-1,-1;-1,1,1,-1,1;-1,1,1,1,-1;-1,1,1,1,1;1,-1,-1,-1,-1;1,-1,-1,-1,1;1,-1,-1,1,-1;1,-1,-1,1,1;1,-1,1,-1,-1;1,-1,1,-1,1;1,-1,1,1,-1;1,-1,1,1,1;1,1,-1,-1,-1;1,1,-1,-1,1;1,1,-1,1,-1;1,1,-1,1,1;1,1,1,-1,-1;1,1,1,-1,1;1,1,1,1,-1;1,1,1,1,1];
m=w*q;
 r=m*pi;
 r1=r;
 r2=r+pi/2;
  r3=r+pi;
   r4=r+pi*3/2;
   n=zeros(128,8);
   for i=1:32,
       n(4*(i-1)+1,(1:8))=r1(i,(1:8));
       n(4*(i-1)+2,(1:8))=r2(i,(1:8));
       n(4*(i-1)+3,(1:8))=r3(i,(1:8));
       n(4*(i-1)+4,(1:8))=r4(i,(1:8));
   end
    b1=cos(n)
    b2=sin(n)
 M=2^20;
Na=10;
L=0;
time=M/Na;
time=fix(time);
a=round(rand(1,M));
d=[0:Na*8-1];
s=[0,a(1),a(2),a(3),a(4)];
x=5;
 tt=[0:Na*8-1];
hold on
  for i=1:500
      for j=0:9
        s=[s(2:5),a(x)];
        if x==5,
        L=0;
        else a(x-5)=2*a(x-5)-1;       
        L=L+a(x-5); 
        end,
        x=x+1;        
        L=rem(L,4);
        if L<0,L=L+4;end,
        su=s(1)*64+s(2)*32+s(3)*16+s(4)*8+s(5)*4+L+1;
        bb=b1(su,(1:8));
        d(1+j*8:8+j*8)=bb;
    end
    plot(tt,d);
end