example7.m

或许这个代码会对大家有用

clear all;
clc;
%定义进制 ,码长 ,和信息长度
M=16; N =15; K =11 ;
%定义载波频率 ,数字传输频率和模拟频率
Fc =10000; Fd =1000 ; Fs =100000 ;
%产生一个 1 1 0 0 X1包含待传输信息的矩阵
% randint产生随机整数矩阵
msg=randint(K*100,1);
%使用 BCH编码 .码长为 N
%信息长度为 K
code=encode(msg, N, K, ' bch' );
%使用 1 6- QASK数字调制方法
%进制为 M=16
modu=dmod(code, Fc, Fd, Fs, 'qask',M);
%假定加性白高斯噪声的 STD为 0 . 1
std_value=0.1 ;
modu_noise =modu+randn(length(modu), 1 )*std_value;
%解调过程
demo =ddemod(modu_noise, Fc, Fd, Fs, 'qask' , M);
% BCH解码
%码长为 N,信息长度为 K
msg_r =decode(demo, N, K, ' bch' );
rate=biterr(msg, msg_r, M)