传真(Fax)编程

开发平台：
Matlab

test_OverLoad_menu.m：源码内容
							%****************************************************************
% 内容概述：TURBO译码计算速对测试（菜单式）
%          设定计算时长，检测单位时间内的译码bit数
%          该程序用于定性分析各种译码算法在不同条件下的计算速度
%          其计算比特率是对算法复杂性的直观反映。
% 创 建 人：朱殿荣/QQ:235347/MSN:njzdr@msn.com
% 单    位：南京邮电大学，通信工程系
% 创建时间：2005年10月29日
% 修改时间：
% 参考文献：《数字通信－－基础与应用》
%          《High performace parallelised 3GPP Turbo Decoder》
%          《改进的Turbo码算法及其FPGA实现过程的研究》,天津大学，张宁，赵雅兴
%       	K.K.Loo, T.Alukaidey, S.A.Jimaa “High Performance Parallelized
%           3GPP Turbo Decoder”, Personal Mobile Communications
%       	Conference 2003. 5th European (Conf. Publ. No. 492)
%       	3GPP TS 25.212 V6.6.0 (2005-09)
%       	3GPP TS 25.222 V6.2.0 (2004-12) 
%       	刘东华。Turbo码原理与应用技术。电子工业出版社，2004.1
% 版权声明：任何人均可复制、传播、修改此文件，同时需保留原始版权信息。
%****************************************************************
clc;
clear all;
timerLimit = input('测试时长【60秒】:');
if isempty(timerLimit)
   timerLimit = 60;
end
length_interleave = input('交织长度=帧长－尾比特长度【1024】:');
if isempty(length_interleave)
   length_interleave = 1024;
end
iter = input('迭代次数【10】:');
if isempty(iter)
   iter =10;
end
algorithm = input('译码算法【1:LOG-MAP,2:MAX-LOG-MAP(缺省),3:SEMITH-LOG-MAP)】:');
if isempty(algorithm)
   algorithm =2;
end
save_mat = input('是否保存测试结果到MAT文件 【1-保存,0-不保存（缺省）】:');
if isempty(save_mat)
   save_mat=0;
end
if save_mat==1
    matFileName = input('MAT文件名 【''运算时长临时测试数据.mat''】:');
    if isempty(matFileName)
        matFileName='运算时长临时测试数据.mat';
    end
end
fprintf('----------------------测试参数----------------------n'); 
fprintf(' 交织长度=%4dbit；迭代次数=%2dn',length_interleave,iter);
fprintf(' 运算时长=%4d秒n',timerLimit);
switch algorithm
    case 1
        fprintf(' 译码算法：LOG-MAPn');
    case 2
        fprintf(' 译码算法：MAX-LOG-MAPn');
    case 3
        fprintf(' 译码算法：门限MAX-LOG-MAPn');
end
if save_mat==1
    fprintf(' 保存测试结果到 ： %4sn',matFileName);
else
    fprintf(' 不保存测试结果到文件n');
end    
fprintf('----------------------------------------------------n'); 
rate=1/3;           %码率
m=3;                    %尾比特数
fading_a=1;             %Fading amplitude
EbNo=1.0;                            %EbNo的采样点
EbNoLinear=10.^(EbNo.*0.1);
num_block_size=length_interleave+m;     %测试的块尺寸，指包含尾比特的软输入系统系列长度
random_in=round(rand(1,length_interleave));  %随机数
[turbod_out,alphain]=turbo(random_in);      %编码
L_c=4*fading_a*EbNoLinear*rate;
sigma=1/sqrt(2*rate*EbNoLinear);
nframe = 0;    % clear counter of transmitted frames
time_begin=clock;
while etime(clock,time_begin)<timerLimit        %nferr:当前迭代次数、EbNo点的错误帧数
    nframe = nframe + 1; 
    noice=randn(4,num_block_size);    %噪声
    soft_in=L_c*(turbod_out+sigma*noice);            %信息噪声叠加
    [hard_out,soft_out]=decoder_all_algorithm(soft_in,alphain,iter,algorithm)%译码
    %注意：提高速度请给上句加上分号－－看不到译码结果，但可以更快译码。
end
if save_mat==1
    save (matFileName,'length_interleave','iter','nframe','timerLimit');
end
fprintf('----------------------测试结果----------------------n');  
fprintf(' 交织长度=%4dbit； 迭代次数=%2dn',length_interleave,iter);
fprintf(' 运算帧数=%2d；运算bit数=%6dn',nframe,nframe*length_interleave);
fprintf(' 运算时长=%4d秒；运算bit速率=%6.2e bit/sn',timerLimit,nframe*length_interleave/timerLimit);
switch algorithm
    case 1
        fprintf(' 译码算法：LOG-MAPn');
    case 2
        fprintf(' 译码算法：MAX-LOG-MAPn');
    case 3
        fprintf(' 译码算法：门限MAX-LOG-MAPn');
end
if save_mat==1
    fprintf(' 保存测试结果到 ： %4sn',matFileName);
else
    fprintf(' 未保存测试结果到文件n');
end    
fprintf('----------------------------------------------------n'); 

						