wuyima.m

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

% 卷积码的性能，（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;
       %BPSK经过AWGN信道
       r=(2*x-1)+sigma(k)*randn(1,length(x));
        %r经过光纤信道
       r(2:5000)=0.3*r(1:4999)+0.7*r(2:5000);
       r(find(r < 0)) = 0;
       r(find(r > 0)) = 1
       %误码计数
       Pe = length(find(x - r))/x_num;
       Plot_Pe = [Plot_Pe Pe];
end
Plot_Pe
%画图
semilogy(plot_EbNo,Plot_Pe,'m*:')
xlabel('Eb/No (dB)')
ylabel('BER')