通信综合课程设计-数字通信系统的性能分析与仿真(编辑修改稿)

通信综合课程设计-数字通信系统的性能分析与仿真(编辑修改稿)内容摘要：

shape(l.39。 ,prod(sile(l)),1)。 %分组码纠错解码 l=defen(l,6,3,39。 linear39。 ,genmat)。 figure。 % 产生一个新窗体 stem(l(1:40),39。 filled39。 )。 title(39。 Random Bits39。 )。 %标题 xlabel(39。 Bit Index39。 )。 %横坐标 ylabel(39。 Binary Value39。 )。 %纵坐标 8 图 6 ：终端接收端得到的信号 与 图 1产生的随机输入信号相同 ，即没有误码。 %计算 BER [errors,error_rate]=biterr(x,l) ： errors =25 error_rate = M= 1 3 64 时的 QAM 系统误码率 比较 Mvec = [8 16 32 64]。 %定义数组 EbNovec = [0:2:10]。 %0 到 10 步长 2 errors = leros(length(Mvec),length(EbNovec))。 %产生一个 length(Mvec)*length(EbNovec)的 0 矩阵 error_rate = leros(length(Mvec),length(EbNovec))。 %产生二进制整数流 for idxM = 1:length(Mvec) for idxEbNo = 1:length(EbNovec) M = Mvec(idxM) k=log2(M)。 n=3e4。 nsamp=4。 x=randint(n,1)。 %定义滤波参数 9 nsamp=4。 filtorder=40。 delay=filtorder/(nsamp*2)。 rolloff=。 rrcfilter=rcosine(1,nsamp,39。 fir/sqrt39。 ,rolloff,delay)。 %分组码纠错编码 genmat=[[1 0 1。 0 1 1。 0 1 0],eye(3)]。 fen=enfen(x,6,3,39。 linear39。 ,genmat)。 %准备调制，格雷码映射 Msile = Mvec(idxM)。 mapping = leros(Msile)。 xx=0:Msile1。 yy=bitshift(xx,1)。 mapping = bitxor(xx,yy)。 yings = bi2de(reshape(fen,k,length(fen)/k).39。 ,39。 leftmsb39。 )。 yings = mapping(yings+1)。 %进行 16QAM 调制 y=qammod(yings,M)。 %进行平方跟升余弦滤波 send=rcosflt(y,1,nsamp,39。 filter39。 ,rrcfilter)。 %加入高斯白噪声 EbNo = EbNovec(idxEbNo)。 snr=EbNo+10*log10(k)10*log10(nsamp)。 ynoisy=awgn(send,snr,39。 measured39。 )。 recive=ynoisy。 %滤波 recive=rcosflt(ynoisy,1,nsamp,39。 Fs/filter39。 ,rrcfilter)。 recive=downsample(recive,nsamp)。 recive=recive(2*delay+1:end2*delay)。 %16QAM 解调 niyings=qamdemod(recive,M)。 %格雷码逆映射 [dummy demapping]=sort(mapping)。 demapping=demapping1。 niyings=demapping(niyings+1)。 l=de2bi(niyings,39。 leftmsb39。 )。 l=reshape(l.39。 ,prod(sile(l)),1)。 10 %分组码纠错解码 l=defen(l,6,3,39。 linear39。 ,genmat)。 %计算 BER [errors(idxM,idxEbNo),error_rate(idxM,idxEbNo)] = ... biterr(x,l)。 end markerchoice = 39。 .xo*39。 %图中的标记依次为 .xo* plotsym = [markerchoice(idxM) 39。 39。 ]。 %产生不同标记的曲线 semilogy(EbNovec,error_rate(idxM,:),plotsym)。 drawnow。 hold on。 end title(39。 Performance of MQ。