pe.m

MATLBA 信道仿真,完成信道卷积码循环码的编译码

%计算原始误码率和卷积码编译码后的误码率
function [p1,p2]=Pe(snr_in_dB)
%snr_in_dB为输入的信噪比
%p1为原始误码率
%p2为卷积码编译码后的误码率

E=1;
snr=10^(snr_in_dB/10);
sgma=sqrt((E/snr)/2);%根据信噪比计算噪声方差

N=100;%测试的帧数
input=round(rand(1,N*25));%信息序列
%原始误码率，经过信道
u1=2*input-1;
u1=sign(u1+sgma*randn(size(u1))+eps);
u1=(u1+1)/2
%判决，计算误码率p1
p1=(length(find(u1~=input)))/length(input);
%经卷积码编译码后的误码率
G=[1 0 0 1 1;1 1 0 1 1];%编码器的冲激响应
u2=viterbi_encoder1(G,input)%卷积码编码
%经过信道
u2=2*u2-1;
u2=sign(u2+sgma*randn(size(u2))+eps);
u2=(u2+1)/2;
[u,hanming]=VitDec(G, u2, 1)%卷积码译码
u=u';
%判决，计算误码率p2
p2=(length(find(u~=input)))/length(input);