基于卷积码的差错控制系统仿真

%实现(2,1,9)的编码、软判决滑动窗维特比译码，其生成多项式为G0=561(八进制),G1=753(八进 %制)，调制方式为BPSK，信道为AWGN，译码深度=10 clear all; close all SNRdB=-10:0.5:3; %设置信噪比范围，0-3dB declen=40;%译码深度 SNRnum=length(SNRdB); %信噪比数目 iter=10;%每个信噪比下的迭代次数 for i=1:SNRnum %循环内计算每个信噪比下的误码率 for j=1:iter %每一个信噪比下迭代计算误码率10次，再求平均误码率 trel=poly2trellis(9,[561 753]); %卷积码(2,1,9)网格图，约束长度为9 siglen=1000000; %设置信号长度 msg=randint(siglen,1); %生成0，1序列，长度同信号长度 encode=convenc(msg,trel,0); %从0状态开始做卷积编码 I=0.5*ones(siglen*2,1); y=encode-I; bpsk=sign(y);%BPSK调制 channelout=awgn(bpsk,SNRdB(i)); %添加高斯白噪声 debpsk=channelout*0.5+0.5; %解调 parti=0:.15:.9;%设置量化等级划分 codebk=0:7;%设输出等级 [x,qcode]=quantiz(debpsk,parti,codebk); %量化，准备维特比软判决 decode= vitdec(qcode',trel,declen,'cont','soft',3);%维特比译码，量化级数=2^3 [errorbit,errorrate(j)] = biterr(decode(declen+1:end),msg(1:end-declen)); %计算本次BER end BER(i)=sum(errorrate)/iter;%求平均BER end semilogy(SNRdB,BER);%画不同信噪比下的误码率 xlabel('信噪比SNR(dB)'); ylabel('误码率'); grid on;