test_uncode.asv

目录结构 ├─WuYuFei ├─WuYufei_matlab ├─cap_turbo ├─mother └─paper constituent_decoder_SemiTh.m co

%****************************************************************
% 内容概述：未编码传输AWGN信道测试
%          每帧的输入是相同的，（不是每帧都不一样，这样可以大幅降低计算量）
%          在帧比特足够多的情况下，应该可以保证随机性（未理论证实）
% 创 建 人：朱殿荣/QQ:235347/MSN:[e-mail]
% 单    位：南京邮电大学，通信工程系
% 创建时间：2005年9月10日
% 修改时间：
% 参考文献：
% 版权声明：任何人均可复制、传播、修改此文件，同时需保留原始版权信息。
%****************************************************************
clear;
clc;
time_begin=datestr(now);
fading_a=1;             %Fading amplitude
snr=0:0.1:6;                            %snr的采样点
EbNoLinear=10.^(snr.*0.1);
num_block_size=1000000000;                    %测试的块尺寸，指包含尾比特的软输入系统系列长度
err_counter=0;        %初始化错误比特计数器
ber=0;                 %初始化错误比特率

random_in=round(rand(1,num_block_size));  %随机数

for nEN=1:length(snr)
    sigma=1/sqrt(2*EbNoLinear(nEN));
    noice=randn(1,num_block_size);    %噪声
    %awgn_in=random_in+sigma*noice;            %信息噪声叠加
    awgn_in=awgn(random_in,snr(nEN),'measured');            %信息噪声叠加
    hard_out(awgn_in>=0)=1;
    hard_out(awgn_in<0)=-1;              %译码
    
    errs=length(find(hard_out(1:num_block_size)~=random_in));%当前点错误bit数
            if errs>0 
                err_counter(1,nEN)=err_counter(1,nEN)+errs;
                nferr(1,nEN)=nferr(1,nEN)+1;
            end
        end
        ber(1,nEN) = err_counter(1,nEN)/nframe/(num_block_size);%误比特率
        %fer(iter,nEN) = nferr(iter,nEN)/nframe; %误帧率
    else
        ber(1,nEN)=1.0e-7;
    end
    fprintf('snr：%1.2f；误码率：%8.4e；\n',...
        snr(nEN),ber(1,nEN));
    %save cap算法06_WYF噪声_max_log_map.mat snr ber;
end
semilogy(snr,ber(1,:));
xlabel('SNR(dB)');
ylabel('Bit Error Rate');
title('未编码信息通过AWGN信道性能图,1024帧长，WYF噪声加法');
%legend('1次迭代','2次迭代','3次迭代');

time_end=datestr(now);
disp(time_begin);
disp(time_end);