资源描述:
《GMSK系统的建模与仿真MATLAB程序.doc》由会员上传分享,免费在线阅读,更多相关内容在教育资源-天天文库。
1、主程序:sr=256000.0;%符号率ml=1;%解调电平数br=sr.*ml;%比特率nd=100;%每次循环中的符号个数ebn0=15;%信噪比IPOINT=8;%过采样点数%*************************滤波器初始化***************************irfn=21;%阶数B=0.25*sr;B2=0.6*sr;[xh]=gaussf(B,irfn,IPOINT,sr,1);%发送器滤波器设计[xh2]=gaussf(B2,irfn,IPOINT,sr,0);%接收器滤波器设计%*******************衰落初始化*********
2、***********tstp=1/sr/IPOINT;%时间分辨率itau=[0];%直达波形dlvl=[0];n0=[6];%瑞利衰落th1=[0.0];itnd0=nd*IPOINT*100;itnd1=[1000];now1=1;fd=320;flat=1;%********************开始计算*************************************nloop=1000;%循环数noe=0;%错误数nod=0;%传送数据数目foriii=1:nloop%***************************产生数据*******************
3、*************data1=rand(1,nd.*ml)>0.5;%rand:builtinfunctionsubplot(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);
4、title('基带信号通过高斯白噪声后')th=zeros(1,length(data3)+1);ich2=zeros(1,length(data3)+1);qch2=zeros(1,length(data3)+1);forii=2:length(data3)+1th(1,ii)=th(1,ii-1)+pi/2*data3(1,ii-1)./IPOINT;endich2=cos(th);qch2=sin(th);fich2=fft(ich2,800);fqch2=fft(qch2,800);subplot(4,2,3),plot(abs(fich2));title('同相支路频谱');su
5、bplot(4,2,4),plot(abs(fqch2));title('正交支路频谱');%**************************衰减计算***********************spow=sum(ich2.*ich2+qch2.*qch2)/nd;attn=0.5*spow*sr/br*10.^(-ebn0/10);attn=sqrt(attn);%**********************衰落信道**********************[ifade,qfade]=sefade(ich2,qch2,itau,dlvl,th1,n0,itnd1,now1,leng
6、th(ich2),tstp,fd,flat);%产生的数据送入衰落模拟器itnd1=itnd1+itnd0;%更新衰落计数器%*********************加性高斯白噪声**********************[ich3,qch3]=comb(ifade,qfade,attn);%加入高斯白噪声[ich4,qch4]=compconv(ich3,qch3,xh2);syncpoint=irfn*IPOINT-IPOINT/2+1;ich5=ich4(syncpoint:IPOINT:length(ich4));qch5=qch4(syncpoint:IPOINT:lengt
7、h(qch4));fich5=fft(ich5,100);fqch5=fft(qch5,100);subplot(4,2,5),plot(abs(fich5));title('加入噪声后的同相支路频谱');subplot(4,2,6),plot(abs(fqch5));title('加入噪声后的正交支路频谱');%****************************GMSK解调********************