分组码.zip

close all clear all clc %% 主函数 N = 1000; %每个信噪比条件下要发送多少个样本点 snrdB_min = -5; snrdB_max = 5;% SNR范围是0-5dB snrdB = snrdB_min:1:snrdB_max; mentercarluo = 100;% 蒙特卡罗实验次数 ber_nocode = zeros(mentercarluo,length(snrdB)); %临时存储矩阵 ber_code = zeros(mentercarluo,length(snrdB)); %临时存储矩阵 for m=1:mentercarluo %开始蒙特卡洛实验 sym_initial = round(rand(1,N));% 生成原始序列 % 获得汉明码编码C_hamming G = [1 0 0 0 1 1 0;... 0 1 0 0 1 0 1;... 0 0 1 0 0 1 1;... 0 0 0 1 1 1 1]; % 生成矩阵G [k,n] = size(G); C_ham = hamming(sym_initial,G); [groups,~] = size(C_ham); C_ham_stream = reshape(C_ham.',[1,length(C_ham(:))]);% 转换成码流 % 两种序列都通过BPSK-AWGN信道 [~,BER_dir] = BPSK(sym_initial,snrdB);% 未编码码流 [RX_ham,~] = BPSK(C_ham_stream,snrdB);% 汉明编码码流 % 汉明码解码并计算SNR errors_ham = zeros(1,length(snrdB)); % 预分配错误内存 for i=1:length(snrdB) RX_ham1 = reshape(RX_ham(i,:),[n,groups]).'; C_result1 = decode(RX_ham1,G); errors_ham(i) = sum(sum(mod(C_result1+C_ham(:,1:k),2))); end BER_ham = errors_ham/N; % 得到BER ber_nocode(m,:) = BER_dir; %每次蒙特卡罗实验 存入临时矩阵 ber_code(m,:) = BER_ham; %每次蒙特卡罗实验 存入临时矩阵 end ber_nocode_final = mean(ber_nocode,1); %对蒙特卡罗实验结果，按列求均值 ber_code_final = mean(ber_code,1); %对蒙特卡罗实验结果，按列求均值 % 计算两种方案的BER-SNR曲线 figure semilogy(snrdB,ber_nocode_final,'r-*',snrdB,ber_code_final,'b-s');grid on; set(gca,'xtick',-5:1:5); xlabel('信噪比 SNR / dB');ylabel('误码率 BER'); legend('未编码','汉明编码');