TA的每日心情 | 怒 2019-11-19 15:34 |
|---|
签到天数: 1 天 [LV.1]初来乍到
|
EDA365欢迎您登录!
您需要 登录 才可以下载或查看,没有帐号?注册
x
" t. \- k( P9 K- B2 h
msg=gf(randint(1,7),1); %information code1 S' p5 U3 o; f1 f) s7 l
c=bchenc(msg,15,7); %encode
/ l9 G# x+ Y# w0 Oe=gf([1 0 0 0 0 0 0 0 0 0 0 1 0 0 0],1);%2 bit errors
; c2 ~0 i" V" p9 r* X7 E- x# @/ Lr1=c+e;% receive code# Z4 v+ Q: l& Z
r=r1;
6 O! p6 s' P0 x- ?g=gf([1 1 1 0 1 0 0 0 1],1);%generator polynomial
5 O1 k: r! T5 s+ q0 Um=3;n=0;k=0;
0 K( N% U, k1 \, F8 |6 H2 U+ o5 ~while m>2 & k<300 j. D3 _5 X0 c4 o
k=k+1;0 ?) ]! v, C0 @; h' n
IF r(1)==13 D* K: D4 u5 {4 G' _
for j=1:9. }1 A h- z7 @- b1 ?1 W3 q% V# P
r(j)=r(j)+g(j);. U* R; o6 _1 B1 B, X6 w. u* h
end
/ n# |, h- }, P7 d( B0 I- vend0 O# C4 Q& D/ v! a
m=sum(r==1);2 S& G% U. l6 y
while r(1)==0 & m>2 & k<30 %left shift
# K* B4 k w" d* {$ W8 In=n+1;
, _2 L, |9 C3 @/ ]+ @( n" yfor i=1:14
* ] r3 Z4 n3 U& t/ _/ cr(i)=r(i+1);4 v2 V$ V5 N0 z# d/ x7 G v
end
F, l0 N0 ]2 ?$ {( kr(15)=0;5 y* P7 w% C8 o/ P8 Q/ R
end $ j; u$ p0 ^ L0 ?# Y. U; B7 v
end% O; T$ v0 A( p1 Z) I% V5 Q5 d
s=mod(n,15);
' m3 D+ ~* t0 V! Z4 C$ ^4 Ofor k=1:s %cyclic right shift
8 |0 j! p# r5 w B! Ht=r(15);
: o) l# R+ z, v, K) k7 ffor i=0:13
5 a4 x2 b' u- h: T* Mr(15-i)=r(14-i);
. E. J! `) _2 `2 B% Hend
9 _. p6 O" V+ Q) {/ q0 Xr(1)=t;
* Z* N+ B- N2 W+ {" v* tend# h$ O8 g c' \& j6 Y
e1=r;+ [5 m( e) _; v- H0 H
m1=sum(e1==1);
t0 }: U& V# E) z Gif m1<37 k3 H; V- c9 s) P
disp('OK e(x)=r(x) mod(g(x))');
4 \& V& s& V( H* L1 jend/ }. {4 |1 }! j1 u$ x4 M9 G ~
4 X+ B9 N1 r) {4 E. s& b' A: l1 V, l$ S3 C- F6 ?
想问一下这段代码是什么意思,最好能具体到每一个循环 是能够直接运行的嘛? 他的输出有哪些?最后disp那个输出的是什么意思? |
|