Matlab设计卷积编码和Viterbi已实现_functions.docVIP

%%---------------------------------------testViterbi.m--------------------------------------------------------------%% %%% Annotated by Quan Xu on Nov. 22, 2011 clear all; clc; close; delete *.mat; EsNodB = 0:1:10; EsNo = 10.^(EsNodB/10); N_samples = 100; N_bits = 512; for k = 1 : length(EsNodB) for m = 1 : N_samples % noise power WGNPower = 1/EsNo(k); % for uncoded system WGNPower1 = 2/EsNo(k); % for 1/2 convolutional coded system % random message msg = randint(1,N_bits); % unencode data b = msg; % Encode data bits trel = poly2trellis(3,[7 5]); % Define trellis b1 = convenc(msg,trel); % Encode % modulate x = 1 - 2*b; x1 = 1 - 2*b1; % AWGN channel h = ones(size(x)); h1 = ones(size(x1)); % Gaussian noise generation w = sqrt(WGNPower/2)*(randn(size(x)) + sqrt(-1)*randn(size(x))); w1 = sqrt(WGNPower1/2)*(randn(size(x1)) + sqrt(-1)*randn(size(x1))); % recieved signals y = h.*x + w; y1 = h1.*x1 + w1; % demodulate uncoded data outUncod = uncode(y,h); % Soft-decision viterbi decoding outVit1 = simeudistVit(trel, y1); errrate(m) = mean(abs(outUncod-msg)); errrate1(m) = mean(abs(outVit1-msg)); end % for m ber(k) = mean(errrate) ber1(k) = mean(errrate1) end % for k semilogy(EsNodB,ber,k,EsNodB,ber1,-r); grid on axis([1 10 10^-5 1]) xlabel(10log(Eb/No)) ylabel(BER) title(AWGN channel viterbi decoder); h = legend(uncoded, [2,1,2]convolutional encoding); set(h,Interpreter,none) save conViterbi.mat; % EOF %%---------------------------------------- testViterbi.m ------------------------------------------------------%% %%----------------------------------------simEudist.m ------------------------------------------------------%% %% Calulating Euclidean distance %% (x-y).^2 -- xy ; sin