用c++或VB编写CRC码【（7,3）循环码】生成程序，生成码为10111

主要有拍启函数void code(char m,char* in,char* out)实现编码

m为生成码，这里10111对应为23

in为输入3位二进制码字符串的指针

out为输出7位二进制码字符串的指针

#include<stdio.h>

#define M char(23) //对应生成码10111

void code(char ,char*,char*)

void main()

{

char incode[8],outcode[8]

outcode[7]='颂贺态\0'

printf("输入三位由0，1组成的二进制码：")

scanf("%s",&incode)

code(M ,incode,outcode)

printf("输出码为：%s",outcode)

}

void code(char m,char* in,char* out)

{

int i

char j,temp=0

for(i=0i<3i++)

{

temp=temp<<1

if(in[i]=='1')

temp=temp^m

}

for(i=0,j=64i<7i++,j>野源>=1)

{

if(temp&j)

out[i]='1'

else

out[i]='0'

}

}

1、循环码编者磨码与解码Matlab源程序（实验以（7，4）循环码进行分析）

m

=

3

n

=

2^m-1

%定义码长

k

=

n-m

%信息位长

msg

=

randint(k*4,1,2)

%随机提取信号，引起一致地分布的任意整数矩阵

subplot(2,2,1)

stem(msg)

title('编码器输入信号')

p=cyclpoly(n,k)

%循环码生成多项式，n=7，k=4

code

=

encode(msg,n,k,'cyclic',p)

%编码函数，对信号进行差错编码

subplot(2,2,2)

stem(code)

title('编码器输出信号')

recode=decode(code,n,k,'cyclic',p)

%对信号进行译码，对接收到的码字进行译码，恢复

出原始的信息，译码参数和方式必须和编码时采用的严格相同

subplot(2,2,3)

stem(recode)

title('译码器输出信号')

t=-1:0.01:1

x=recode

%将recode赋值给x，并进行长度与fft设定

N=length(x)

fx=fft(x)

df=100/N

n=0:N/2

f=n*df

subplot(2,2,4)

plot(f,abs(fx(n+1))*2/念胡N)

grid

title('频谱图')

2、误码率与信噪比之间的关系程序（以（首高斗3，2）循环码进行测试）

m

=

2

n

=

2^m-1

%定义码长

k

=

n-m

%信息位长

Fs=40

%系统采样频率

Fd=1

%码速率

N=Fs/Fd

M=2

%进制数

for

SNRpBit=1:100%信噪比

SNR=SNRpBit/log2(M)

%制造100个信息组，每组k位

msg

=

randint(100,k,[0,1])

code

=

encode(msg,n,k,'cyclic/binary')

%加入噪声

%在已调信号中加入高斯白噪声

noisycode=awgn(code,SNR-10*log10(0.5)-10*log10(N),'measured',[],'dB')

%将浮点数转化为二进制，波形整形过程

for

i=1:100

for

a=1:k+1

if

noisycode(i,a)<0.5

noisycode(i,a)

=

0

else

noisycode(i,a)

=

1

end

end

end

%译码

newmsg

=

decode(noisycode,n,k,'cyclic')

%计算误码率

[number,ratio]=biterr(newmsg,msg)

result(SNRpBit)=ratio

disp(['The

bit

error

rate

is',num2str(ratio)])

end

%不同信噪比下循环码经过加性高斯白噪声信道的误码率

figure(1)

stem(result)

title('循环码在不同信噪比下的误码率')

legend('误码率','*')

xlabel('信噪比')

ylabel('在加性高斯白噪声下的误码率')

---