GMSK.rar_GMSK解调_GMSK调制方面的_gmsk_gmsk 解调_gmsk解调调制

% File_C6: gmsk.m %**** 对于给定信噪比的通信系统，估计差错概率 **** %数据调制采用GMSK clear all; %**** 参数设置部分 **** sr=256000;br=sr;nd=20000;snrindB=5;IPOINT=8; %****滤波器设计**** irfn=21; B=0.25*sr;B2=0.6*sr; %发送器滤波器设计 [xh]=gaussf(B,irfn,ipiont,sr); %接受滤波器设计 [xh2]=gaussf(B2,irfn,ipoint,sr); % 估计差错概率 noe=0; nod=0; %****GMSK调制**** data1=rand(1,nd)>0.5; data11=2*data1-1; data2=oversamp(data11,IPOINT); data3=conv(data2,xh); %定义附加相位序列，I通道信号和Q通道信号 th=zeros(1,length(data3)+1); ich2=zeros(1,length(data3)+1); qch2=zeros(1,length(data3)+1); %产生附加相位信号 for ii=2:length(data3)+1; th(1,ii)=th(1,ii-1)+pi/2*data3(1,ii-1)./IPOIT; end %产生I通道信号和Q通道信号 ich2=cos(th); qch2=sin(th); %计算信道加性白噪声 spow=sum(ich2.*ich2+qch2.*qch2)/nd; attn=0.5*spow*sr/br*10.^(-snrindB/10); attn=sqrt(attn); %****接收机接收到I通道和Q通道信号**** [ich3,qch3]=comb(ich2,qch2,attn); %接收端滤波 ich4=conv(ich3,xh2); qch4=conv(qch3,xh2); %接收端下采样 synncpoint=irfn*IPOINT-INPOINT/2+1; ich5=ich4(syncpoint:IPOINT:length(ich4)); qch5=qch4(syncpoint:IPOINT:length(qch4)); %*****GMSK解调**** demoddata2(1,1)=1; for n=2:nd+1 demoddata2(1,n)=1000*(ich5(1,n-1)*qch5(1,n)-ich5(1,n)*qch5(1,n-1))>0; end demodata=demoddata2(1,2:nd+1); %****估计差错概率（BER）**** noe=sum(abs(data1-demodata)); %统计差错个数 nod=length(data1); %计算发送数据符号的个数 ber=noe/nod; %计算差错概率 fprintf('差错概率估计值=%e\n',noe/nod); %输出计算结果 %End of script file. %File-C6: gaussf.m function[xh]=gaussf(B,irfn,ipoint,sr) %****高斯滤波器设计**** % irfn:设置滤波器的阶数 % ipoint：符号采样点数 % sr:符号传输速率 %B:滤波器的系数 point=ipoint; tr=sr; n=ipoint*irfn; mid=n/2+1; fo=B/sqrt(2*log(2)); for i=1:n icon=i-mid; ym=icon; xt=1/2*(erf(-sqrt(2/log(2))*pi*B*(ym/ipoit-1/2)/tr)+erf(sqrt(2/log(2))*pi*B*(ym/ipoint+1/2)/tr)); end; %End of script file. %File-C6: comb.m function[iout,qout]=comb(idata,qdata,attn) %****产生加性高斯白噪声**** %idata:输入I通道信号， qdata:输入Q通道信号 %iout:输出I通道信号，qout:输出Q通道信号 %attn:attenuation level caused by Eb/No iout=randn(1,length(idata)).*attn; qout=randn(1,length(qdata)).*attn; iout=iout+idata(1,length(idata)); qout=qout+qdata(1,length(qdata)); %End of script file. %File-C6: oversamp.m function[out]=oversamp(indata,sample) %****过采样（在输入数据序列添加零数据） **** %indata:输入数据序列，sample:过采样的点数 nsym=length(indata); out=zeros(1,nsym*sample); for i=1:nsym out(1+(i-1)*sample)=indata(i); end %End of script file.