yimahou.m

信号在信道进行传输时

% 卷积码的性能，（7，5），R＝1/2
clear all;
close all;
R=1/2;
EbN0dB=0:6;
%BPSK调制
EsN0dB=EbN0dB-10*log10(R);
N0=10.^(-EsN0dB/10);%求awgn 信道的单边功率谱密度
sigma=sqrt(N0/2);
plot_EbNo =0:6;
Plot_Pe = []
x_num=5000;
G=[1 1 1;1 0 1];
for k=1:length(N0)
       %信源产生
       x = bingen(x_num);
       x(find(x < 0)) = 0;
       x(find(x > 0)) = 1;
       %卷码编码
       s=cnv_encd(G,1,x);
       %BPSK经过AWGN信道
       r=(2*s-1)+sigma(k)*randn(1,length(s));
       %r经过光纤信道
       r(2:10004)=0.3*r(1:10003)+0.7*r(2:10004);
       r(find(r < 0)) = 0;
       r(find(r > 0)) = 1
       %译码
       dd=viterbi(G,1,r);
       %误码计数
       Pe = length(find(x - dd))/x_num;
       Plot_Pe = [Plot_Pe Pe];
end
Plot_Pe
%画图
semilogy(plot_EbNo,Plot_Pe,'m*:')
xlabel('Eb/No (dB)')
ylabel('BER')