资源描述：
《GMSK系统的建模与仿真MATLAB程序.docx》由会员上传分享，免费在线阅读，更多相关内容在行业资料-天天文库。
1、主程序：sr=256000.0;%符号率ml=l;%解调电平数br=sr.*ml;%比特率nd=100;%每次循环中的符号个数cbn()=15;%信噪比1POINT二&%过采样点数irfn=21;%阶数滤波器初始化************************水**E二0・25*sr;B2=0.6*sr;[xh]=gaussf(B,irfn,lP()lNT,sr,l);%发送器■滤波器设计[xh2]=gaussf(B2,irFn,IP()INT,sr,0);%接收器滤波器设计%*戎************
2、****取衰落初女台彳匕********************lstp=1/sr/lPOINT;%时间分辨率itau=[0];%直达波形dlvl=[0];n0=[6];%瑞利衰落th1=[().()];imdO=nd*IP()INTnOO;itndl=[1000];nowl=1;fd二320;精品精品flat=1;%********************开始计算**************************宜**窘*******nk)op=1000;%循环数noc=0;%错误数nod=0;%传送数
3、据数冃foriii=l:nloop<^)*戎****************安********产生^^,**********官******************戎**daval=rand(l,nd*ml)>0.5;%rand:builtinfunctionsubplot(4,2,1),stairs(datal);产生的原始信号)*********************官**********axis([0,50,-2,2]);%*事******官*********官********GMSKi同制Idata
4、l1=2*datal-1;data2=ovcrsamp(datal1Jcnglh(dalal1),1POINT);subplot(4,2,2),plot(data3);dava3=conv(dala2,xh);曲X基带信号通过高斯白噪声后')th=zcros(lJcngth(data3)+1);ich2=zcros(1,length(data3)+1);qch2=zcros(l,length(data3)+l);forii=2:lcngth(data3)+lth(l4i)=th(l,ii-l)+pi/2*d
5、ata3(l,ii-l)./IPOINT;end精品ich2=cos(th);qch2=sin(th);fich2=fft(ich2,800);fqch2=fft(qch2,800);subplot(4,2,3),plot(abs(fich2));titlcC同相支路频谱)；subplot(4,2,4),plot(abs(fqch2));讹正交支路频谱)；%**戎***************官******客衰减j十算***********官***********Spow=sum(ich2.*ich2+qc
6、h2.*qch2)/nd;attn=0.5*spc)w*sr/br*lO.A(-cbn()/10);(%****************avtn=sqrt(attn);衰落信道*******************：fo卅ffadc,qfadc]=scfadc(ich2,qch2,itau,dlvl,thl,n0,itndl^lowl,lcngth(ich2),lstp,fd,flat);%产生的数据送入衰落棋拟器imdl=itndl+imdO;%更新衰落计数器忆*********************加
7、枠盲出斤白瓏吉****官******************[ich3,qch3]=comb^fadc^fadc^atm);%加入高斯白噪声[ich4,qch4]=compconv(ich3,qch3,xh2);syncpoint=irfh*IPOINT-IPOlNT/2+1;ich5=ich4(syncpoinv：IP()lNT:lcnglh(ich4));精品qch5=qch4(syncpoint:IPOINT:lcngth(qch4));fich5=fft(ich5JOO);fqch5=fft(qch
8、5,100);subplot(4,2,5),plot(abs(fich5));tklcC加入噪声后的同相支路频谱')；subplot(4^,6),plot(abs(fqch5));tklcC加入噪声后的正交支路频谱')；%****************************gmsk解调*****************************dcmoddala2(l,1)=-1；fork=3:2:nd*ml+
