GMSK系统地建模与仿真MATLAB程序.pdfVIP

主程序： 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); % 接收器滤波器设计 %******************* 衰落初始化 ******************** 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; % 传送数据数目 for iii=1:nloop %*************************** 产生数据 ******************************** data1=rand(1,nd.*ml)0.5; % rand: built in function 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); 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)./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)); 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,length(ich 2),tstp,fd,flat); % 产生的数据送入衰落模拟器 itnd1 = itnd1+ itnd0; % 更新衰落计数器 %*********************