用matlab模拟FEC和交织两种方式.docx

用matlab模拟FEC和交织两种方式FEC（前向纠错方式）：交织 s1=1,2,3,4,5,6,7,8,9,10,11,12,13,14,15,16,17,18,19,20,21,22,23,24,25,26,27,28,29,30,31,32,33,34,35,36,37,38,39,40,41,42,43,44,45,46,47,48,49,50,51,52,53,54,55,56; x1=(reshape(s1,8,7); x1(2,2)=0;x1(2,3)=0; x1(2,4)=0;x1(2,5)=0; s2=reshape(x1,1,56); x2=reshape(s2,7,8); x2(2,2)=10;x2(2,3)=11; s3=reshape(x2,1,56); a=s1,s2,s3; plot(s1,s2);x1 = 1 2 3 4 5 6 7 8 9 0 0 0 0 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 4849 50 51 52 53 54 55 56L=1000; M=4;%每个符号的比特数N=2M-1;%编码后码字长度K=N-4;%信息长度MSG=randint(L,1);%随机产生L比特信号TP=gftuple(-1:N-1,M);%产生加罗华域元素PG=rspoly(N,K);%产生生成式CODE,ADDED=rsenco(MSG,TP,K);%编码NOI=rand(length(CODE)/M,1) msg=randint(k*500,1,2,st1); code = encode(msg,n,k,hamming/binary); errors = zeros(size(code); errors(n-2:n+3) = 1 1 1 1 1 1; inter = randintrlv(code,st2); inter_err = bitxor(inter,errors); deinter = randdeintrlv(inter_err,st2); decoded = decode(deinter,n,k,hamming/binary); disp(Number of errors and error rate, with interleaving:);Number of errors and error rate, with interleaving: code_err = bitxor(code,errors); decoded = decode(code_err,n,k,hamming/binary); disp(Number of errors and error rate, without interleaving:);Number of errors and error rate, without interleaving: number_without,rate_without = biterr(msg,decoded);L=1000;M=4;%每个符号的比特数N=2M-1;%编码后码字长度K=N-4;%信息长度 MSG=randint(L,1);%随机产生L比特信号TP=gftuple(-1:N-1,M);%产生加罗华域元素PG=rsgenpoly(N,K);%产生生成式CODE,ADDED=rsenco(MSG,TP,K);%编码R=0:0.01:1;%错误率for j=1:length(R) %循环修改错误率 NOI=rand(length(CODE)/M,1)R(j);%加入R%的噪声NOI=(NOI*ones(1,M);%产生突发错误NOI=NOI(:);CODE_NOI=rem(CODE+NOI,2);%噪声加入信号DEC,ERR,CCODE,ERR_C=rsdeco(CODE_NOI,TP,K);%译码MSG=MSG;zeros(ADDED,1);%调整长度error=0;%错误码计数初值Q