基于64QAM调制解调的LDPC编译码算法误码率matlab仿真
2023-09-14 09:14:38 时间
目录
1.算法描述
LDPC译码从译码算法的实现角度可以将译码类型分为硬判决译码和软判决译码两种类型。其中,硬判决译码方式实现过程较为简单,其通过一个预先设置的阈值对译码信息进行判决,如果大于阈值则判决输出“1”,否则判决输出“0”。软件判决译码方式则根据码元错误概率最小的方式进行译码。因此,软件判决译码方式具有更优的译码性能。
BP译码算法是一种基于概率域的译码算法。其译码步骤如下所示:
虽然基于对数域的BP译码算法在一定程度上降低了BP译码算法的复杂度,但是由于其涉及到正切等运算。随着码长的增加,LLR-BP译码算法的复杂依旧较大。由此产生了MS最小和译码算法。 最小和译码算法的步骤如下:
2.仿真效果预览
matlab2022a仿真结果如下:
3.MATLAB核心程序
load G.mat;% 首先加载G , H
load H.mat;
max_iter=25;
L_frame=size(G,1);
n_frame=200;
start=0;
step=2;
finish=10;
r=size(G,1)/size(G,2);
M=6;
Es=42; % 一个64QAM符号能量 也是信号功率
Eb=Es/M;
plot_pe=[];
Q=1;
% err_list_index=1; % used by test
axis_EbN0=start:step:finish;
%******************************************************
for EbN0=start:step:finish;
Liner_EbN0=10^(EbN0/10);
pe_number=0;
variance=0.5*(Eb/Liner_EbN0)/r; % 噪声方差, 来源于a^2/(2*delta^2)=Liner_EbN0,a^2信息幅值的平方即比特能量, 2*delta^2=N0. N0/r得到编码后的噪声功率谱密度
for number_frame=1:1:n_frame % variance是单独加在同相(正交相)的噪声功率普密度
x_msg = (sign(randn(1,size(G,1)))+1)/2; % random bits
x_code_msg= mod(x_msg*G,2); % coding
in=x_code_msg;
x_tx_msg= QAM64m(in); % 16QAM modulation
real_tx=real(x_tx_msg); imag_tx=imag(x_tx_msg);
waveform= x_tx_msg+sqrt(variance)*( randn(size(real_tx)) + randn(size(imag_tx))*i ) ; % AWGN transmission
[f0,f1]=QAM64d(waveform,2*variance);
LLR=log(f0./(f1+eps));
[z_hat, success, k] = log_ldpc_decode(H,LLR,max_iter);; % likelihoods
x_dec_msg = z_hat(1:size(G,1)); % 码率1/2 时使用
pe_number=pe_number+sum(x_msg~=x_dec_msg);
current_time=fix(clock);
fprintf('i am working %g,total(%g) %g年 %g月 %g日 %g时 %g分 %g秒\n\n',Q,(((finish-start)/step)+1)*n_frame,current_time(1),current_time(2),current_time(3),current_time(4),current_time(5),current_time(6))
Q=Q+1;
fprintf('\n\n')
end
pe=pe_number/(L_frame*n_frame);
plot_pe=[plot_pe,pe];
end
figure;
semilogy(axis_EbN0,plot_pe,'b*-')
xlabel('Eb/N0 信号功率/噪声功率(db)')
ylabel('BER')
title('误码率')
grid on
A_070
4.完整MATLAB
V
相关文章
- MATLAB中求矩阵的逆矩阵方法(2种)「建议收藏」
- matlab图像处理初步,MATLAB数字图像处理初步「建议收藏」
- 手眼标定 matlab源码,经典手眼标定算法matlab
- matlab创建控件失败win7,Win7系统无法安装ActiveX控件的解决方法
- 【嵌入式案例分享】使用Matlab生成可供TMS320C6748开发板使用的算法
- matlab griddata插值太慢,griddata在二维插值中的应用问题,有待解决~
- matlab怎么对语音信号处理,语音信号处理MATLAB程序
- Matlab粒子群算法(PSO)优化程序——经典实例
- matlab中如何求插值点,MATLAB插值「建议收藏」
- 最短路径Dijkstra算法原理及Matlab实现「建议收藏」
- matlab多重比较lsd法,多重比较法-LSD I 附赠统计学最全思维导图~[通俗易懂]
- 深入浅出PID控制算法(一)————连续控制系统的PID算法及MATLAB仿真[通俗易懂]
- MATLAB自带的dwt2和wavedec2函数实现基于小波变换的自适应阈值图像边缘检测
- 标准粒子群算法(PSO)及其Matlab程序和常见改进算法_粒子群算法应用实例
- matlab保存所有图,Matlab中图片保存的5种方法
- Mathworks Matlab R2022b 简体中文 (含激活补丁)
- matlab的fprintf语句_matlab绝对值函数
- matlab循环读取txt文件
- Matlab赋值_matlab二维数组赋值
- matlab命令,应该很全了!「建议收藏」
- matlab 马赫带效应,matlab图像处理基础实例
- Matlab中的偏最小二乘法(PLS)回归模型,离群点检测和变量选择|附代码数据
- MATLAB数据挖掘用改进的K-Means(K-均值)聚类算法分析高校学生的期末考试成绩数据
- 分享一个很不错的车辆垂向动力学matlab工具箱
- 基于Matlab的有限元网格自动生成算法 | Q4、Q8、Abaqus单元网格
- 人工智能算法:Matlab遗传算法工具箱使用方法
- 人工智能算法:基于Matlab遗传算法的实现示例
- 【MATLAB】三维绘图 ( 三维绘图步骤 )
- 快速安装Linux上的MATLAB(linux安装matlab)