一、实验目的:
(1) 熟悉 Matlab 工作环境及工具箱;
(2) 掌握线性分组码的编码、译码原理以及纠错原理。
二、实验仪器:
Matlab2015a
三、实验原理:
四、实验步骤:
1.先利用Matlab编写函数;
2.再根据函数将实验要求生成实验结果;
3.根据实验内容填写实验报告;
五、实验内容及数据整理:
用Matlab软件编程实现线性分组码的信道编码和译码
六、实验结果及讨论
clc;clear;close all
%编码
G=input('请输入生成矩阵G,例如:G=[1 0 1 1 1;0 1 1 0 1]\n G=');
[k,n]=size(G);
r=n-k;
m=input('请输入需传送信息m,如m=[0 0 0 1 1 0 1 1]\n m=');
l=length(m);
if(mod(l,k))
disp('输入的信息有误');
else
ge=l/k;
%将输入序列转化成矩阵m
temp1=[];
for i=1:ge
temp1(i,:)=m(k*(i-1)+1:i*k);
end
m=temp1;
%求校验矩阵H
c=mod(m*G,2);
P=G(:,k+1:n);
H=[P',eye(r)];
disp('生成矩阵');G
disp('校验矩阵');H
disp('编码矩阵');c
end
disp('敲回车键继续'); pause
%解码
y=input('输入接收序列y,如:y=[0 0 0 0 0 0 1 1 0 1 1 0 1 1 1 1 0 0 1 0]\n y=');
temp2=[];
for i=1:ge
temp2(i,:)=y(1,n*(i-1)+1:i*n);
end
y=temp2
s=mod(y*H',2);
e=s*pinv(H');
for i=1:ge
for j=1:n
if(e(i,j)>0.5-eps)
e(i,j)=1;
else
e(i,j)=0;
end
end
end
cc=mod(y+e,2);
%cc=xor(y,e)
sc=cc(:,1:k);
disp('差错图样'); e
disp('估计值'); cc
disp('译码序列'); sc