%****************************************************************
% 內(nèi)容概述：AWGN信道測試
%          每個(gè)SNR采樣點(diǎn)的測試幀數(shù)是隨snr指數(shù)上升的，
%          這樣不僅減少了計(jì)算量而且可以保證精度。
%          每幀的輸入是相同的，（不是每幀都不一樣，這樣可以大幅降低計(jì)算量）
%          在幀比特足夠多的情況下，應(yīng)該可以保證隨機(jī)性（未理論證實(shí)）
% 創(chuàng) 建 人：朱殿榮/QQ:235347/MSN:njzdr@msn.com
% 單    位：南京郵電大學(xué)，通信工程系
% 創(chuàng)建時(shí)間：2005年7月24日
% 修改時(shí)間：
% 參考文獻(xiàn)：
% 版權(quán)聲明：任何人均可復(fù)制、傳播、修改此文件，同時(shí)需保留原始版權(quán)信息。
%****************************************************************
clear;
clc;
rate=1/3;           %碼率
wear=1;             %損耗系數(shù)
time_begin=datestr(now);
snr=0:0.5:3;                              %snr的采樣點(diǎn)
EbNo=10.^(snr.*0.1);

num_frame=round(10.^(snr+1));         %測試的幀數(shù)；
num_block_size=1024;                    %測試的塊尺寸，指包含尾比特的軟輸入系統(tǒng)系列長度
err_counter(3,1:length(snr))=0;         %初始化錯(cuò)誤比特計(jì)數(shù)器
guage=0;                                %初始化計(jì)算量
counter_guage=0;                        %初始化計(jì)算量計(jì)數(shù)器
for num_it=1:1:3
    for ii=1:length(snr)
        guage=guage+num_frame(ii);      %統(tǒng)計(jì)計(jì)算量
    end
end
random_in=round(rand(1,num_block_size-3));  %隨機(jī)數(shù)
[turbod_out,alphain]=turbo(random_in);      %編碼

for num_it=1:1:3
    for ii=1:length(snr)

        L_c=4*wear*EbNo(ii)*rate;
        sigma=1/sqrt(2*rate*EbNo(ii));
               
        for frame=1:num_frame(ii)
        
            noice=randn(3,num_block_size);    %噪聲
            soft_in=turbod_out+sigma*noice;            %信息噪聲疊加
            soft_in=L_c*soft_in;
            %soft_in=awgn(turbod_out,snr(ii),'measured');    %通過AWGN信道，加入噪聲
            counter_guage=counter_guage+1;
            fprintf('總幀數(shù)：%6.0f；完成幀數(shù)：%6.0f；當(dāng)前snr點(diǎn)總幀數(shù)：%4.0f；當(dāng)前snr點(diǎn)完成幀數(shù)：%4.0f；計(jì)算進(jìn)度：%2.1f%% \n',...
                guage,counter_guage,num_frame(ii),frame,counter_guage*100/guage);            
            [hard_out,soft_out]=deturbo_cap(soft_in,alphain,num_it);
            err_counter(num_it,ii)=err_counter(num_it,ii)+...
            length(find(hard_out(1:num_block_size-3)~=random_in));
        end
        ber(num_it,ii)=err_counter(num_it,ii)/((num_block_size-3)*num_frame(ii));
    end
end
semilogy(snr,ber(1,:),snr,ber(2,:),snr,ber(3,:));
xlabel('SNR(dB)');
ylabel('BER');
title('3GPP標(biāo)準(zhǔn) max-log-map譯碼算法 譯碼性能圖,1024交織長度，楚宛琴噪聲加法');
legend('1次迭代','2次迭代','3次迭代');
save cap算法06_楚宛琴噪聲加法.mat snr ber soft_out;
time_end=datestr(now);
disp(time_begin);
disp(time_end);
%plot(snr,ber);
%counter_ber
%soft_in(1,1:10);
%hard_out(1:10);
%soft_out(1:10);

%--------------------------
%soft_in2(1,:)=soft_in(1,alphain);
%soft_in2(2,:)=soft_in(3,:);
%soft_out=decoder_3GPP_MAX_new(soft_in2);
%soft_out(alphain)=soft_out;
%----------------------------
%soft_out=decoder_3GPP_MAX_new(soft_in);
%for ii=1:num_block_size
%    if soft_out(ii)>0
%        hard_out(ii)=1;
%    else
%        hard_out(ii)=-1;
%    end
%end
%----------------------------