在好例子网,分享、交流、成长!
您当前所在位置:首页MATLAB 开发实例MATLAB语言基础 → 线性编组码matlab实现

线性编组码matlab实现

MATLAB语言基础

下载此实例
  • 开发语言:MATLAB
  • 实例大小:2.68KB
  • 下载次数:2
  • 浏览次数:65
  • 发布时间:2020-12-18
  • 实例类别:MATLAB语言基础
  • 发 布 人:xyaiyyqx
  • 文件格式:.m
  • 所需积分:2
 相关标签: 组码

实例介绍

clc;clear all;
%(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关系');


标签: 组码

实例下载地址

线性编组码matlab实现

不能下载?内容有错? 点击这里报错 + 投诉 + 提问

好例子网口号:伸出你的我的手 — 分享

网友评论

发表评论

(您的评论需要经过审核才能显示)

查看所有0条评论>>

小贴士

感谢您为本站写下的评论,您的评论对其它用户来说具有重要的参考价值,所以请认真填写。

  • 类似“顶”、“沙发”之类没有营养的文字,对勤劳贡献的楼主来说是令人沮丧的反馈信息。
  • 相信您也不想看到一排文字/表情墙,所以请不要反馈意义不大的重复字符,也请尽量不要纯表情的回复。
  • 提问之前请再仔细看一遍楼主的说明,或许是您遗漏了。
  • 请勿到处挖坑绊人、招贴广告。既占空间让人厌烦,又没人会搭理,于人于己都无利。

关于好例子网

本站旨在为广大IT学习爱好者提供一个非营利性互相学习交流分享平台。本站所有资源都可以被免费获取学习研究。本站资源来自网友分享,对搜索内容的合法性不具有预见性、识别性、控制性,仅供学习研究,请务必在下载后24小时内给予删除,不得用于其他任何用途,否则后果自负。基于互联网的特殊性,平台无法对用户传输的作品、信息、内容的权属或合法性、安全性、合规性、真实性、科学性、完整权、有效性等进行实质审查;无论平台是否已进行审查,用户均应自行承担因其传输的作品、信息、内容而可能或已经产生的侵权或权属纠纷等法律责任。本站所有资源不代表本站的观点或立场,基于网友分享,根据中国法律《信息网络传播权保护条例》第二十二与二十三条之规定,若资源存在侵权或相关问题请联系本站客服人员,点此联系我们。关于更多版权及免责申明参见 版权及免责申明

;
报警