pn_gen.m

MATLAB 程序设计教程 对初学者很容易上手 内有实例 讲解

function [data, RandData] = pn_gen(N,M)

%++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
%	Filename:	pn_gen.m 
%
%	Description:
%
%	1) This routine generates N random polar impulse train w/M samples/bit
%	2) Generate N random samples
%
%	input :		N  - no. of bits
%			M  - no. of sample per bit
%	output:		data(N*M, 1)  - random data w/ M sample/bit
%			RandData  - random data
%	local var:	xtemp
%
% Author : Jeff Laster, MPRG, Virginia Tech 
%+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++

data  = zeros((N+3)*M,1);
rand('seed',sum(100*clock));		% resets the seed value of the random number generator
xtemp = rand(N,1) - 0.5;		% vector with +/-
RandData  = sign(xtemp);		% take the sign

for i = 1:N
	data( (i-1)*M+1,1 ) = RandData(i);		% insert impulse at beginning
end