test_OverLoad_menu.m
上传用户:hnyfjx
上传日期:2013-06-30
资源大小:2149k
文件大小:4k
- %****************************************************************
- % 内容概述:TURBO译码计算速对测试(有菜单)
- % 设定计算时长,检测单位时间内的译码bit数
- % 创 建 人:朱殿荣/QQ:235347/MSN:njzdr@msn.com
- % 单 位:南京邮电大学,通信工程系
- % 创建时间:2005年10月29日
- % 修改时间:
- % 参考文献:
- % 版权声明:任何人均可复制、传播、修改此文件,同时需保留原始版权信息。
- %****************************************************************
- 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 =1;
- 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); %信息噪声叠加
- switch algorithm
- case 1
- [hard_out,soft_out]=decoder_logmap(soft_in,alphain,iter); %log-map译码
- case 2
- [hard_out,soft_out]=decoder(soft_in,alphain,iter); %max-log-map译码
- case 3
- [hard_out,soft_out]=decoder_SemiTh(soft_in,alphain,iter); %门限max-log-map译码
- end
- 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');