计算程序和运行结果如下:
Private Sub Command1_Click()
y = 2376000
p = 0.1345
i = 0
Print " 现在产值:"y
Do
i = i + 1
y = y * (1 + p)
Print i"年后敬拍仔产值亮汪:"贺扰y
Loop While y <2376000 * 2
Print i"年后产值可以翻一番。"
End Sub
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('在加性高斯白噪声下的误码率')
循环码的编译过程
循环码是线慎团性分组码的一种,所以它具有线性分组码的一般特性,此外还具有循环性。循环码的编码和解码设备都不太复杂,且检(纠)错能力强。它不但可以检测随机的错误,还可以检错突发的错误。(n,k)循环码可以检测长为n-k或更短的任何突发错误,包括首尾相接突发错误。\n循环码是一种无权码,循环码编排的特点是相邻两个数码之间符合卡诺图中的邻接条件,即相邻两个数码之间只有一位码元不同,码元就是组成数码的单元。循袜好环码又称格雷码。
告孝铅 一种循环码的全部码字
欢迎分享,转载请注明来源:内存溢出
评论列表(0条)