GMSK.rar_gmsk

主程序：

sr=256000.0; % 符号率

ml=1; % 解调电平数

br=sr.*ml; % 比特率

nd = 100; % 每次循环中的符号个数

ebn0=15; % 信噪比

IPOINT=8; % 过采样点数

%************************* 滤波器初始化***************************

irfn=21; % 阶数

itau = [0]; %直达波形

dlvl = [0];

n0=[6]; %瑞利衰落

th1=[0.0];

itnd0=nd*IPOINT*100;

itnd1=[1000];

now1=1;

fd=320;

flat =1;

%********************开始计算*************************************

subplot(4,2,1) ,stairs(data1);

title('产生的原始信号')

axis([0,50,-2,2]);

%***************************GMSK 调制********************************

data11=2*data1-1;

data2=oversamp(data11,length(data11),IPOINT);

data3=conv(data2,xh);

subplot(4,2,2),plot(data3);

title('基带信号通过高斯白噪声后')

th=zeros(1,length(data3)+1);

th(1,ii)=th(1,ii-1)+pi/2*data3(1,ii-1)./IPOINT;

end

ich2=cos(th);

qch2=sin(th);

fich2=fft(ich2,800);

fqch2=fft(qch2,800);

subplot(4,2,3),plot(abs(fich2));

title('同相支路频谱');

subplot(4,2,4),plot(abs(fqch2));

%********************** 衰落信道 **********************

[ifade,qfade]=sefade(ich2,qch2,itau,dlvl,th1,n0,itnd1,now1,length(ich

itnd1 = itnd1+ itnd0; %更新衰落计数器

%*********************加性高斯白噪声 **********************

[ich3,qch3]= comb(ifade,qfade,attn);% 加入高斯白噪声

[ich4,qch4] = compconv(ich3,qch3,xh2);

syncpoint =irfn*IPOINT-IPOINT/2+1;

subplot(4,2,6),plot(abs(fqch5));

title('加入噪声后的正交支路频谱');

%****************************GMSK 解调*****************************

demoddata2(1,1)=-1;

for k=3:2:nd*ml+1

demoddata2(1,k)=ich5(1,k)*qch5(1,k-1)*cos(pi*(k))>0;

for n=2:2:nd*ml+1

demoddata2(1,n)=ich5(1,n-1)*qch5(1,n)*cos(pi*(n))>0;