new.txt

GMSK基带系统仿真 本试验是用数字基带处理方法来实现调试期算法的基带硬件实验。 通过用软件matlab进行了GMSK相位路径和眼图的仿真,之后对任何一种五位随机码

clear;
clc;

cosfai=zeros(4,8,32);
for d=1:32
    for L=0:3
        for t=0:7
            m=change(d-1);
            ph=phase(t/8,m)+pi*L/2; 
            cosfai(L+1,t+1,d)=cos(ph);
            sinfai(L+1,t+1,d)=sin(ph);  
        end
    end
end



N=1000+7;                             %伪随机序列码元数
ran=sign(randn(1,N));                 %伪随机序列
for nn=6:1:N-2
 L=mod(sum(ran(1:nn-5)),4)+1;      %求象限值
    n=locate1(ran(nn:-1:nn-4));         %求状态值
    for ii=1:1:8          
       m1(nn-5,ii)=cosfai(L,ii,n);%读余弦表，取8个抽样的离散值 
       m2(nn-5,ii)=sinfai(L,ii,n);
    end
end        
figure(1)
for count=1:1:(N-4)/10         
 m11=m1(1:10,:);%同一屏幕显示10个码元,取10行
    mmm=reshape(m11',1,80) ;
     t=1:1:80;
    plot(t,mmm);%画一个屏幕的眼图
    m1(1:10,:)=[];
    title('余弦的眼图COSfai的眼图')
   xlabel('X轴')
   ylabel('Y轴')
    hold on; %禁止清屏
end
figure(2)
for ct=1:1:(N-4)/10
    m22=m2(1:10,:);
    sss=reshape(m22',1,80);
     t=1:1:80;
    plot(t,sss);%画一个屏幕的眼图
    m2(1:10,:)=[];
    hold on; %禁止清屏
    title('正弦的眼图Sinfai的眼图')
     xlabel('X轴')
   ylabel('Y轴')
end
%…………………………
s=1;
for n=1:1:32
    for L=1:1:4
        for j=1:1:8
Cos_table(s,:)=change2(cosfai(L,j,n));
Sin_table(s,:)=change2(sinfai(L,j,n));
s=s+1;
        end
    end
end

Tablebincos(:,1:8)=Cos_table(:,3:10);                   %余弦表低8位
Tablebinsin(:,1:8)=Sin_table(:,3:10);                    %正弦表低8位
Tablebincs(:,1:2)=Sin_table(:,1:2);
Tablebincs(:,3:4)=0;
Tablebincs(:,5:6)=0;
Tablebincs(:,7:8)=Cos_table(:,1:2);
TTablebincos=Tablebincos';                           %二进制余弦表和正弦表的生成
TTablebinsin=Tablebinsin';
tempcs1=Tablebincs';
Transbincos=flipud(TTablebincos);
Transbinsin=flipud(TTablebinsin);
Transbincs=flipud(tempcs1);
for i=1:1:1024
    TTransbincos(:,8*(i-1)+1)=Transbincos(:,i);
    TTransbinsin(:,8*(i-1)+1)=Transbinsin(:,i);
    TTransbincs(:,8*(i-1)+1)=Transbincs(:,i);
end
fid=fopen('Tablesin.bin','wb');
fwrite(fid,TTransbinsin,'ubit1');                        %生成二进制正弦表文件（低8位）
fid=fopen('Tablecos.bin','wb');
fwrite(fid,TTransbincos,'ubit1');                       %生成二进制余弦表文件（低8位）
fid=fopen('Tablecs.bin','wb');
fwrite(fid,TTransbincs,'ubit1');                   %生成二进制正弦和余弦表文件（高2位）
fclose('all');