main1012.m

放大转发协同通信的有关仿真代码

 %function [source_frame_bit,coded_bit,coded_data] = generate_coded_data( );
source_bit=randint(10,244);
constraint_length = [4];
source_frame_bit =[];
coded_bit =[];
coded_data =[];
codegenerator = [15,17,13,15];
for i=1:10
    CRC_CCITT = CRC_16(source_bit(i,:));%the crc bits
    frame_bit(i,:) = [source_bit(i,:),CRC_CCITT];
    source_frame_bit =[source_frame_bit,frame_bit(i,:)];
    %codegenerator = [15,17,13,15];
    %coded_bit(i,:) = convenc(frame_bit(i,:), poly2trellis(constraint_length, codegenerator));
    %coded_data =[coded_data,coded_bit(i,:)];
end
source_frame_bit;
coded_bit = convenc(source_frame_bit, poly2trellis(constraint_length, codegenerator));%1*10400
for i=1:10
    coded_data = [coded_data;coded_bit((i-1)*1040+1:i*1040)];%10*1040
end
transmit_symbol = coded_data;
transmit_symbol_N1 = zeros(10,780);

transmit_symbol_N1(:,1:3:end) = transmit_symbol(:,1:4:end);
transmit_symbol_N1(:,2:3:end) = transmit_symbol(:,2:4:end);
transmit_symbol_N1(:,3:3:end) = transmit_symbol(:,3:4:end);

N1_symbol=[];
for i=1:10
N1_symbol = [N1_symbol,transmit_symbol_N1(i,:)];
end
symbol_sequence =1 - 2* N1_symbol;

insert_decode_data = zeros(1,10400);
        insert_decode_data(1:4:end) = symbol_sequence(1:3:end);
        insert_decode_data(2:4:end) = symbol_sequence(2:3:end);
        insert_decode_data(3:4:end) = symbol_sequence(3:3:end);
        insert_decode_data(4:4:end) = 0;
        opmode = 'trunc'; 
        tblen = 5;
        decode_symbol_N1= vitdec(insert_decode_data, poly2trellis(constraint_length, codegenerator), tblen, opmode, 'unquant');
        error_count3 = sum(abs(source_frame_bit - decode_symbol_N1))