TA的每日心情 | 怒 2019-11-19 15:34 |
|---|
签到天数: 1 天 [LV.1]初来乍到
|
EDA365欢迎您登录!
您需要 登录 才可以下载或查看,没有帐号?注册
x
* Y4 v* O$ }, X1 `; I# n* r# Q5 |
msg=gf(randint(1,7),1); %information code
8 a( N$ j) C/ ?; d: Y: l# Dc=bchenc(msg,15,7); %encode
# ^3 e. Z4 E5 A' \e=gf([1 0 0 0 0 0 0 0 0 0 0 1 0 0 0],1);%2 bit errors: d3 v3 h0 z3 i6 o
r1=c+e;% receive code+ V6 E& U) x. \% Z
r=r1;
5 T, e+ K( G$ g9 I4 Bg=gf([1 1 1 0 1 0 0 0 1],1);%generator polynomial
* P1 H7 x7 T+ n4 [' B0 um=3;n=0;k=0;
. H* ~! K& |4 g! d' @8 e7 k% Gwhile m>2 & k<30
8 X. E: M% f* V- t: Tk=k+1;
/ N( m$ Q0 d5 q' Z4 m: X# iIF r(1)==15 g8 f3 y O. q5 R0 t* @ P
for j=1:9
; L. o4 q! ]( _' X, f$ w! nr(j)=r(j)+g(j);" L1 V' n( W$ B7 P- t# Z) `
end
" v: i9 x- M+ N1 Dend
6 m4 [! J4 L$ u: \m=sum(r==1);4 ~' _6 M4 l- ?8 _6 X9 R
while r(1)==0 & m>2 & k<30 %left shift
$ c: V- l4 Y0 Q1 [; B, g8 R( h/ w' f; kn=n+1; " W: y4 c& h% a4 `1 t
for i=1:14' D6 `% C' C* j) g, y9 P
r(i)=r(i+1);
7 r! Q- _1 _1 `1 Uend
( _' i( B$ Z" O# y& {* Q# gr(15)=0;8 d0 \# g( k$ C, h: n, [
end
3 m [1 s+ H0 Aend
# b* e2 J4 q& \" \" H# K* Ws=mod(n,15);
: e ^7 ]. {' n& P6 b! r0 J& P# O Dfor k=1:s %cyclic right shift
& H& \; r4 `: Dt=r(15);
" V- F6 V1 M1 V) w- _for i=0:13' v8 w! h' p, ]9 t- o
r(15-i)=r(14-i);
2 o. M6 y( E& n9 Y$ ?2 K' J2 cend& V) X: C' s8 S3 R' n
r(1)=t; 9 Q7 b- L! d1 ^0 l/ G& W, i0 [ H
end
. H( t+ k; K! Z1 Le1=r;2 i4 Z( {/ T g) h
m1=sum(e1==1);0 M# b6 `& O& Z% l4 U4 r
if m1<3; W! D" b% F2 x2 c
disp('OK e(x)=r(x) mod(g(x))');
* d$ b: _1 r1 ~* r8 ~+ _end6 z! o9 j' [5 H! X( X8 q5 N" V0 H
, `, w3 N7 W. W$ N5 N
: V+ ?% O+ G7 W! O0 z5 o& z; W. C: T
想问一下这段代码是什么意思,最好能具体到每一个循环 是能够直接运行的嘛? 他的输出有哪些?最后disp那个输出的是什么意思? |
|