TA的每日心情 | 怒 2019-11-19 15:34 |
|---|
签到天数: 1 天 [LV.1]初来乍到
|
EDA365欢迎您登录!
您需要 登录 才可以下载或查看,没有帐号?注册
x
8 K9 a- s% s3 ^, \' M% lmsg=gf(randint(1,7),1); %information code/ s6 I: i; _6 I1 `, h/ m& J
c=bchenc(msg,15,7); %encode& M/ H) v. u& j
e=gf([1 0 0 0 0 0 0 0 0 0 0 1 0 0 0],1);%2 bit errors
: s/ c5 p% \/ j6 j9 T7 dr1=c+e;% receive code
8 _0 c h! j' Kr=r1;8 i+ J9 B/ F7 ?1 q# z: b
g=gf([1 1 1 0 1 0 0 0 1],1);%generator polynomial
{$ `. @7 X/ n& Y8 t5 Rm=3;n=0;k=0;
. U3 ^' p0 J) p6 `$ wwhile m>2 & k<30) S% C" r$ m& p. {5 Y% J; o
k=k+1;* G% }& H) o# K1 |) B
IF r(1)==1; z' x3 H+ T0 Y( B8 T: F, a3 H, V
for j=1:99 e* ]1 y d6 k; w* X3 a. _9 g
r(j)=r(j)+g(j);3 M1 \) c4 R* Z, V
end
( C- b9 M1 R4 n% U& |: e- Nend
" x2 ~6 H4 ]% Y+ m, @# im=sum(r==1);) v- j+ G( b, A) ~2 `7 t# G
while r(1)==0 & m>2 & k<30 %left shift' I! z; K! F! q- E( R$ U* g
n=n+1;
, Q! W# {# v5 r/ bfor i=1:14
) n p% I1 a2 E0 ^6 er(i)=r(i+1);
; ]- P# B# u2 o; `, }end
+ v7 D4 g5 @7 B2 y/ E, kr(15)=0;5 i2 M0 s- s1 z3 R; I
end
, I- H- R' r( U4 [; H {: mend2 B3 E" Z5 V; G
s=mod(n,15);
6 G. q8 t) ], |for k=1:s %cyclic right shift2 n* e) V0 k7 l- |3 N
t=r(15);2 A9 }5 q( N/ ~' O+ l% r4 ? V0 A
for i=0:13
4 I* Z# e. [! s6 p$ Hr(15-i)=r(14-i);
, m& Y X8 { ` e5 N& Nend
- M$ U: q+ y1 o0 W: rr(1)=t; 0 M, K# z0 p# j+ v
end3 r+ P# N% K R
e1=r;4 [' e8 p. c' ^6 U# f0 i& e# K
m1=sum(e1==1);: ]6 K7 }" d: m2 T& U
if m1<3+ Z, H, W1 N8 I; U) m3 o* S
disp('OK e(x)=r(x) mod(g(x))');
4 r4 c' J' ~2 iend1 p4 B: Y( @# f5 n3 ^
) }; ?/ J" {' l6 ^. }: Q* E" A/ f1 O) Y8 \: {2 \! {8 s
想问一下这段代码是什么意思,最好能具体到每一个循环 是能够直接运行的嘛? 他的输出有哪些?最后disp那个输出的是什么意思? |
|