实例介绍
%(n,k)线性分组码
k=3;n=7 %矩阵的行数
m=zeros(2^k,k);
for i = 1:2^k %产生输入矩阵
m(i,:) = bitget(i-1,k:-1:1); %返回i-1的三位二进制值
end
disp(m);
%% 把生成矩阵变为系统矩阵
G=[0,1,0,1,0,1,0
1,1,1,0,0,0,0
1,1,0,0,1,1,0];
temp=G;
G(1,:)=xor(temp(1,:),temp(3,:));
G(2,:)=temp(1,:);
G(3,:)=xor(temp(2,:),temp(3,:));
%%
c=mod(m*G,2);
a=find(c);
%% 把码字变为BPSK
a1=zeros(2^k,n)
for i=1:length(a)
a1(a(i))=-1;
end
b1=find(~c)
for i=1:length(b1)
a1(b1(i))=1
end
%%
H=gen2par(G);%生成矩阵变为校验矩阵
p=G(:,4:7);
i=eye(4)
H=[p',i];%奇偶校验矩阵
for i = 1:2^n %产生输入矩阵
E(i,:) = bitget(i-1,n:-1:1); %返回i-1的三位二进制值
end
s=mod(E*H',2);
s=s 0;
e(1:16,:)=E(1:16,:);
for i=1:2^(n-k)
for j=1:n
if length(find(e(i,:)==1))>length(find(E(i 16*(j-1),:)==1))
e(i,:)=E(i 16*(j-1),:);
end
end
end
%% 标准阵列译码表
C=cell(8,1);
for i=1:8
C{i}=zeros(16,7);
end
for i=1:n
for j=1:16
C{i}(j,:)=xor(e(j,:),c(i,:));
end
end
%% 均值为0,方差为s2的高斯白噪声信道N
flag=0;
for sigma2=1:0.1:100
N=cell(1,10);
r=cell(1,10);
transfernum=100;
for order=1:transfernum
N{order}=normrnd(0,sqrt(sigma2),size(a1,1),size(a1,2));%产生噪声
r{order}=a1 N{order};%加入噪声的序列
%% 对接受序列r解调
for i=1:8
for j=1:7
if r{order}(i,j)<0
b{order}(i,j)=1;
else
b{order}(i,j)=0;
end
end
end
yima=b;
for i=1:8
for j=2:16
if b{order}(i,:)==C{i}(j,:)
yima{order}(i,:)=C{i}(1,:);
end
end
end
%% 计算误码率
mayuantotal=56;
mayuanerror=zeros(1, transfernum);
for i=1:8
for j=1:7
if yima{order}(i,j)~=c(i,j)
mayuanerror(order)=mayuanerror(order) 1;
end
end
end
pe(order)= mayuanerror(order)/ mayuantotal;
end
flag=flag 1;
avpe(flag)=sum(pe)/ transfernum;
end
semilogy(1:0.1:100,avpe,'r');
grid on; % 坐标轴开启
axis([0 100 0 1]); % 限制作图范围
xlabel('方差'); % 横坐标
ylabel('BER'); % 纵坐标
title('噪声方差与BER关系');
标签: 组码
小贴士
感谢您为本站写下的评论,您的评论对其它用户来说具有重要的参考价值,所以请认真填写。
- 类似“顶”、“沙发”之类没有营养的文字,对勤劳贡献的楼主来说是令人沮丧的反馈信息。
- 相信您也不想看到一排文字/表情墙,所以请不要反馈意义不大的重复字符,也请尽量不要纯表情的回复。
- 提问之前请再仔细看一遍楼主的说明,或许是您遗漏了。
- 请勿到处挖坑绊人、招贴广告。既占空间让人厌烦,又没人会搭理,于人于己都无利。
关于好例子网
本站旨在为广大IT学习爱好者提供一个非营利性互相学习交流分享平台。本站所有资源都可以被免费获取学习研究。本站资源来自网友分享,对搜索内容的合法性不具有预见性、识别性、控制性,仅供学习研究,请务必在下载后24小时内给予删除,不得用于其他任何用途,否则后果自负。基于互联网的特殊性,平台无法对用户传输的作品、信息、内容的权属或合法性、安全性、合规性、真实性、科学性、完整权、有效性等进行实质审查;无论平台是否已进行审查,用户均应自行承担因其传输的作品、信息、内容而可能或已经产生的侵权或权属纠纷等法律责任。本站所有资源不代表本站的观点或立场,基于网友分享,根据中国法律《信息网络传播权保护条例》第二十二与二十三条之规定,若资源存在侵权或相关问题请联系本站客服人员,点此联系我们。关于更多版权及免责申明参见 版权及免责申明
网友评论
我要评论