yimahouyantu.m

来自「信号在信道进行传输时」· M 代码 · 共 35 行

% 卷积码的性能，（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经过光纤信道
       m(2:10004)=0.3*r(1:10003)+0.7*r(2:10004);
       t(find(m < 0)) = 0;
       t(find(m > 0)) = 1;
       %译码
       dd=viterbi(G,1,t);
       %误码计数
       Pe = length(find(x - dd))/x_num;
       Plot_Pe = [Plot_Pe Pe];
end
eyediagram(x,2)
eyediagram(r,2);
eyediagram(dd,2);