近期在SDR平台OAI中做数据信道的编译码,即5G中的相关标准,将LTE中的数据信道编码turbo码改成ldpc码。基于在做完编译码之后需要测试ldpc编译码的snr-bler性能,在此,将测得的snr bler数据画出snr-bler图是工作的一小部分,为了便于工作,在matlab平台上写了一个m文件来实现此功能。
m文件中的代码如下:
%作图一
a=importdata('mcs00.txt') %调用文件中的数据并赋值给数组a
A=a(:,1) %取第一列作为横坐标
B=a(:,2) %取第二列作为纵坐标
hold on
plot(A,B,'-*','Color',[0.0 0.0 0.0]) %画图,-*代表图像以*作为点,后面的三个数值是三色控制该线段的颜色,取值范 围是0到1
%作图二
a=importdata('mcs01.txt')
A=a(:,1)
B=a(:,2)
hold on
plot(A,B,'-*','Color',[0.0 0.0 0.0])
%作图三
a=importdata('mcs02.txt')
A=a(:,1)
B=a(:,2)
hold on
plot(A,B,'-*','Color',[0.0 0.0 0.0])
%作图四
a=importdata('mcs03.txt')
A=a(:,1)
B=a(:,2)
hold on
plot(A,B,'-*','Color',[0.0 0.0 0.0])
%作图五
a=importdata('mcs04.txt')
A=a(:,1)
B=a(:,2)
hold on
plot(A,B,'-*','Color',[0.0 0.0 0.0])
%作图六
a=importdata('mcs05.txt')
A=a(:,1)
B=a(:,2)
hold on
plot(A,B,'-*','Color',[0.0 0.0 0.0])
%作图七
a=importdata('mcs06.txt')
A=a(:,1)
B=a(:,2)
hold on
plot(A,B,'-*','Color',[0.0 0.0 0.0])
%作图八
a=importdata('mcs07.txt')
A=a(:,1)
B=a(:,2)
hold on
plot(A,B,'-*','Color',[0.0 0.0 0.0])
%作图九
a=importdata('mcs08.txt')
A=a(:,1)
B=a(:,2)
hold on
plot(A,B,'-*','Color',[0.0 0.0 0.0])
%作图十
a=importdata('mcs09.txt')
A=a(:,1)
B=a(:,2)
hold on
plot(A,B,'-*','Color',[0.0 0.0 0.0])
%作图十一
a=importdata('mcs10.txt')
A=a(:,1)
B=a(:,2)
hold on
plot(A,B,'-*','Color',[0.0 0.0 0.0])
%作图十二
a=importdata('mcs11.txt')
A=a(:,1)
B=a(:,2)
hold on
plot(A,B,'-*','Color',[0.0 0.0 0.0])
%作图十三
a=importdata('mcs12.txt')
A=a(:,1)
B=a(:,2)
hold on
plot(A,B,'-*','Color',[0.0 0.0 0.0])
%作图十四
a=importdata('mcs13.txt')
A=a(:,1)
B=a(:,2)
hold on
plot(A,B,'-*','Color',[0.0 0.0 0.0])
%加网格
grid on
%将纵坐标设置为log
set(gca,'yscale','log')
%设置x轴 y轴 标题
xlabel ('SNR/dB')
ylabel ('BLER')
title ('SNR-BLER LDPC MCS-0-13 Rc=1/3 AWGN HARQ=1')
%设置x轴的范围
xlim([-6 6])
测试效果:
此次测试共画了14个图像,数据文件如下:
每个文件中的数据格式如下:
其中第一列数据为snr,第二列数据为bler
之后运行m文件,则可以生成如下图所示的效果:
东西很简单,不过至少方便了以后的工作,大家喜欢的可以拷贝代码进行测试。欢迎邮件交流:[email protected]