部分源碼：

% 首先加載G , H
clear all
load G.mat;
load H.mat;
max_iter=50;
L_frame=size(G,1);
n_frame=100;
start=0;
step=2;
finish=12;
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是單獨加在同相（正交相）的噪聲功率普密度
? ? ? ?[EbN0,number_frame]
? ? ? ??
? ? ? ? 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 時使用
? ? ? ??
? ? ?% ? if success==0 ? ? ? ? ? ? ? ? ? ? ? ? ? % used by test
? ? ?% ? ? ?err_list(err_list_index,1)=number_frame; ? ? ? ? ?% 第一列代表含有錯誤的frame的號數(shù)
? ? ?% ? ? ?err_list(err_list_index,2)=sum(x_msg~=x_dec_msg); % 第二列代表錯誤的個數(shù)
? ? ?% ? ? ?err_list_index=err_list_index+1;
? ? ?% ? end
? ? ? ??
? ? ? ? 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

semilogy(axis_EbN0,plot_pe,'b*-')

xlabel('Eb/N0 信號功率/噪聲功率(db)')
ylabel('BER')
title('誤碼率')
grid on

誤碼率曲線如下：

?D170

總結(jié)

以上是生活随笔為你收集整理的基于MATLAB的LDPC编译码仿真，调制为64QAM的全部內(nèi)容，希望文章能夠幫你解決所遇到的問題。

如果覺得生活随笔網(wǎng)站內(nèi)容還不錯，歡迎將生活随笔推薦給好友。