通信原理第九次实验报.doc

通信原理第九章实验报告 9.1 %n重复码性能 n=5; m = 0:-0.5:-3; pe = 10.^m; %信道 d = (sign(randn(1,100000))+1)/2; %编码输入 s = [d;d;d;d;d]; %重复5次 s = reshape(s,1,5*length(d)); %将s变成一个序列 %信道误码的仿真 for k=1:length(pe) err = rand(1,length(d)*5); err = errpe(k); %信道误码的随机图样 r = rem(s+err,2); r = reshape(r,5,length(d)); dd = sum(r) 2; %大数判决 error(k) = sum( abs(dd-d) ) / length(d); end loglog(pe,error); xlabel(信道误码率);ylabel(译码后误码率); grid on 9.2 %(7,4)分组码 clear all; close all; H = [1 1 1 0 1 0 0; 0 1 1 1 0 1 0; 1 1 0 1 0 0 1]; G = gen2par(H); %调用Matlab函数求与H对应的生成矩阵G Msg = [0 0 0 0;0 0 0 1;0 0 1 0;0 0 1 1;0 1 0 0;0 1 0 1;0 1 1 0;0 1 1 1; 1 0 0 0;1 0 0 1;1 0 1 0;1 0 1 1;1 1 0 0;1 1 0 1;1 1 1 0;1 1 1 1]; C = rem(Msg*G,2) 9.3a %(7,4)分组码 clear all; close all; H = [1 1 1 0 1 0 0; 0 1 1 1 0 1 0; 1 1 0 1 0 0 1]; G = gen2par(H); %调用Matlab函数求与H对应的生成矩阵G Msg = [0 0 0 0;0 0 0 1;0 0 1 0;0 0 1 1;0 1 0 0;0 1 0 1;0 1 1 0;0 1 1 1; 1 0 0 0;1 0 0 1;1 0 1 0;1 0 1 1;1 1 0 0;1 1 0 1;1 1 1 0;1 1 1 1]; C = rem(Msg*G,2); E = dec2bin([0:2^7-1],7)-48; %调用Matlab函数将整数转换成比特表示 S = rem(H*E,2); %不同的E对应的校验子S %对校验子归类，每个S对应不同的16个E s = S(1,:)*4+S(2,:)*2+S(3,:); for k=1:8 e(k,:)=find(s==k-1); end e0 = dec2bin( e(1,:)-1,7 )-48; %对应S=000的E，编码码字 e1 = dec2bin( e(2,:)-1,7 )-48; %对应S＝001的E e2 = dec2bin( e(3,:)-1,7 )-48; %对应S＝010的E e3 = dec2bin( e(4,:)-1,7 )-48; %对应S＝011的E e4 = dec2bin( e(5,:)-1,7 )-48; %对应S＝100的E e5 = dec2bin( e(6,:)-1,7 )-48; %对应S＝101的E e6 = dec2bin( e(7,:)-1,7 )-48; %对应S＝110的E e7 = dec2bin( e(8,:)-1,7 )-48; %对应S＝111的E 9.3b %循环码 clear all;close all; n = 15; k = 4; p = cyclpoly(n,k,all) %调用Matlab函数得到所有的生成多项式,p是升幂的格式 %编码 [H,G] = cyclgen( n,p(1,:) ); %产生g1(x)对应的监督矩阵和生成矩阵 Msg = [0 0 0 0;0 0 0 1;0 0 1 0;0 0 1 1;0 1 0 0;0 1 0 1;0 1 1 0;0 1 1 1;... 1 0 0 0;1 0 0 1;1 0 1 0;1 0 1 1;1 1 0 0;1 1 0 1;1 1 1 0;1 1 1 1]; C = rem(Msg*G,2