TA的每日心情 | 怒 2019-11-19 15:34 |
|---|
签到天数: 1 天 [LV.1]初来乍到
|
EDA365欢迎您登录!
您需要 登录 才可以下载或查看,没有帐号?注册
x
. K3 h! i; C) j# j imsg=gf(randint(1,7),1); %information code
9 I9 {% a, g& C% I( C5 c! P, ic=bchenc(msg,15,7); %encode& g5 b( Z8 f- ]5 L
e=gf([1 0 0 0 0 0 0 0 0 0 0 1 0 0 0],1);%2 bit errors
) H/ [0 C+ ]. i6 i5 P7 ]r1=c+e;% receive code
6 N1 @) z; U/ H `9 `' q/ B% t2 ^r=r1;* P, l" q Y6 K* \1 r
g=gf([1 1 1 0 1 0 0 0 1],1);%generator polynomial O: |2 M3 _* s7 [. |
m=3;n=0;k=0;
# E8 y5 X4 j) ?3 p; nwhile m>2 & k<30
* L- m2 \ W0 C3 mk=k+1;# o8 d& M+ c9 J# f
IF r(1)==18 h5 w' O8 i& A A% a
for j=1:9
# @% c( O- }& u; X0 q/ }; R% Ir(j)=r(j)+g(j);
0 g: t& b5 o. m( M+ H/ uend0 }, B1 Q. i5 E8 j
end3 a7 v& j8 Q& B# x- }" q. E
m=sum(r==1);% g+ T0 f! v$ w1 A* A" p7 p* E
while r(1)==0 & m>2 & k<30 %left shift
1 Y$ ~) c( ^8 X% \/ {9 n0 M# H3 pn=n+1;
. A' n2 h+ M R% Z5 vfor i=1:14
- a- g8 Q% s+ E/ I. U9 g- O6 Cr(i)=r(i+1);
$ e j, p- X4 [9 M5 }! x8 Q+ Uend! F, A3 Y6 c% ]! H5 }. A. H
r(15)=0;
4 x' j9 n4 |2 J& Y1 qend
; u5 e1 Z7 F" nend* D5 H3 a: D3 M5 z( }. B, O7 G
s=mod(n,15);
4 h$ B% e1 ? [0 ~/ rfor k=1:s %cyclic right shift
0 ~% k: k0 w% Z9 e5 j& m$ |# {t=r(15);+ @; l# J* q& H# [
for i=0:135 f4 L3 [8 D# ^( K
r(15-i)=r(14-i);8 W# |5 s8 O9 D/ R! j( l% g
end
0 X4 ^. l) }+ W9 J2 d! n B6 `r(1)=t;
% y o2 ]+ N- Q$ I2 E% D+ L9 `end0 O0 J$ f+ t& w. Q
e1=r;9 F5 r" S+ }3 r
m1=sum(e1==1);) z+ h0 p: X# [% g9 b3 d# U4 H
if m1<30 v& p3 q0 X3 }, V' q0 S" G
disp('OK e(x)=r(x) mod(g(x))');; n% n) s3 u/ b: u3 J
end J1 z/ b/ c7 X/ m8 r
$ o8 \& t4 E ?. u
, _$ j7 [5 S1 V! S9 ]
想问一下这段代码是什么意思,最好能具体到每一个循环 是能够直接运行的嘛? 他的输出有哪些?最后disp那个输出的是什么意思? |
|