eye.m

基于matlab的gmsk信号的调制

times=256;%总的屏数
na=10;%每屏的码元数
L2=8;%每码元的抽样点数
N=times*na*L2;%总的显示的点数
%ts=1;
len=na*L2;%每屏的显示点数
da=zeros(1,N);%每个点对应的纵坐标组成的矩阵
ats=zeros(1,times*na+5);
ats2=zeros(1,times*na+6);
ats1=round(rand(1,times*na+6));%伪随机序列,多出6个,dan极性
for ii3=[2:times*na+6]
    ats2(ii3)=xor(ats1(ii3-1),ats1(ii3));
end
for ii4=[1:times*na+5]
    ats(ii4)=ats2(ii4+1);
end
atd=ats;
for ii1=[1:times*na+5]
    if ats(ii1)==0,atd(ii1)=-1;end %ats->单极性
end
for ii=[4:times*na+3] %有用的序列值,120个
at1=atd(1:ii-3);
ls=sum(at1);%保证初始相位为0
l=mod(ls,4); %决定象限
n=ats(ii-2)*16+ats(ii-1)*8+ats(ii)*4+ats(ii+1)*2+ats(ii+2)*1;
 for ii2=[1:8]
da((ii-4)*8+ii2)=fcos(n*32+l*8+ii2);
 end
end 
%tt=linspace(0,na*ts,len);
tt=[1:na*L2];
hold on
for ii5=[1:len:N-len]
    %figure(1);
    plot(tt,da(ii5-1+(1:len)));
    xlabel('Time');ylabel('sin(fai(k))');
    title('The eyechart of GMSK');
end