一个turbo码的仿真程序

				function [output, state] = encode_bit(g, input, state)								% the rate is 1/n				% k is the constraint length				% m is the amount of memory				[n,k] = size(g);				m = k-1;								% determine the next output bit				for i=1:n				   output(i) = g(i,1)*input;				   for j = 2:k				      output(i) = xor(output(i),g(i,j)*state(j-1));				   end;				end								state = [input, state(1:m-1)];