【编码译码】基于matlab实现汉明码仿真（含编码译码）

最新推荐文章于 2023-08-06 21:59:23 发布

matlab科研助手

最新推荐文章于 2023-08-06 21:59:23 发布

阅读量887

点赞数 2

分类专栏：信号处理文章标签： matlab 开发语言

本文链接：https://blog.csdn.net/matlab_dingdang/article/details/131873252

版权

信号处理专栏收录该内容

638 篇文章 125 订阅

订阅专栏

该文介绍了使用Matlab进行汉明码的仿真，包括编码和译码过程，对比了软判决和硬判决在AWGN信道下的性能，展示了误码率曲线，并提供了相关代码示例，涉及QAM调制和解调。此外，还提到了作者提供的其他Matlab仿真服务，如优化算法、机器学习、图像处理和无人机路径规划等领域。

摘要由CSDN通过智能技术生成

✅作者简介：热爱科研的Matlab仿真开发者，修心和技术同步精进，matlab项目合作可私信。

🍎个人主页：Matlab科研工作室

🍊个人信条：格物致知。

更多Matlab仿真内容点击👇

智能优化算法神经网络预测雷达通信无线传感器电力系统

信号处理图像处理路径规划元胞自动机无人机

⛄ 内容介绍

利用matlab软件进行汉明码仿真，包含了编码、译码。同时译码分别采用了软判决+最小欧式距离译码、硬判决+最小汉明距离译码方法。

⛄ 部分代码

BitSeq=randi([0,1],[5120000,2]); %生成比特序列SymbolSeq=bi2de(BitSeq,'left-msb');  %比特序列转化为符号序列ModulatedSeq=qammod(SymbolSeq,4);  %调制后的信号序列%QPSK实际误码率：error=zeros([1,11]);for ebno=0:1:10    AWGNChannel = comm.AWGNChannel('NoiseMethod','Signal to noise ratio (Eb/No)','EbNo',ebno,'BitsPerSymbol',2,'SignalPower',mean(abs(ModulatedSeq).^2));    outsignal=AWGNChannel(ModulatedSeq);    recv=qamdemod(outsignal,4);    rec=de2bi(recv,'left-msb');    [~,error(ebno+1)]=biterr(rec,BitSeq);endsemilogy(0:1:10,error,'o-');grid on axis([0,10,1e-5,1])hold onxlabel('Eb/No');ylabel('BER');title('AWGN');%%BitSeq=randi([0,1],[512000,4]); %生成比特序列BitSeq_1 = reshape(BitSeq',[2048000,1]);BitSeq_Encode = encode(BitSeq_1,7,4,'hamming/binary');ModulatedSeq=qammod(BitSeq_Encode,4,'InputType','bit');  %调制后的信号序列error=zeros([1,11]);for ebno=0:1:10    AWGNChannel = comm.AWGNChannel('NoiseMethod','Signal to noise ratio (Eb/No)','EbNo',ebno,'BitsPerSymbol',2,'SignalPower',mean(abs(ModulatedSeq).^2)*7/4);    outsignal=AWGNChannel(ModulatedSeq);    recv = qamdemod(outsignal,4,'OutputType','bit');    rec = decode(recv,7,4,'hamming/binary');    [~,error(ebno+1)]=biterr(rec,BitSeq_1);endsemilogy(0:1:10,error,'o-');legend('Uncoded','Soft decision','Hard decision')%%function [res] = SoftD(Seq)% 软判决函数% INPUT：接收到的编码后的符号% OUTPUT：软判决后的比特序列    global ModulatedY;    global base;    dist = sum(dis_complex(ModulatedY,Seq),2);    [~,index] = min(dist);    res = base(index,:);endfunction [res] = dis_complex(a,b)% 复数欧式距离函数    res = sqrt((imag(a)-imag(b)).^2+(real(a)-real(b)).^2);end