anafsk.m

? bpsk信号的产生可用与仿真的输入信号

function [y,iflaw]=anafsk(N,Ncomp,Nbf);
%ANAFSK Frequency Shift Keying (FSK) signal.
%	[Y,IFLAW]=ANAFSK(N,NCOMP,NBF) simulates a phase coherent
%	Frequency Shift Keying (FSK) signal. This signal is a succession
%	of complex sinusoids of NCOMP points each and with a normalized
%	frequency uniformly chosen between NBF distinct values between
%	0.0 and 0.5.  
% 	Such signal is only 'quasi'-analytic.
%	
%	N     : number of points
%	NCOMP : number of points of each component (default: N/5)
%	NBF   : number of distinct frequencies     (default: 4  )
% 	Y     : signal
% 	IFLAW : instantaneous frequency law  (optional).
%
%	Example :
%	 [signal,ifl]=anafsk(512,64,5); clf; figure(gcf);
%  	 subplot(211); plot(real(signal)); subplot(212); plot(ifl);
%
%	See also ANABPSK, ANAQPSK, ANAASK.

%	O. Lemoine - June 1995, F. Auger - August 1995.
%	Copyright (c) 1996 by CNRS (France).
%
%	------------------- CONFIDENTIAL PROGRAM -------------------- 
%	This program can not be used without the authorization of its
%	author(s). For any comment or bug report, please send e-mail to 
%	[e-mail] 

if (nargin == 0),
 error('The number of parameters must be at least 1.');
elseif (nargin == 1),
 Ncomp=round(N/5); Nbf=4;
elseif (nargin == 2),
 Nbf=4;
end;

if (N <= 0),
 error ('The signal length N must be strictly positive' );
end;

MatlabVersion=version; MatlabVersion=str2num(MatlabVersion(1));
if (MatlabVersion==4), rand('uniform'); end;

m=ceil(N/Ncomp);
freqs=0.25+0.25*(floor(Nbf*rand(m,1))/Nbf-(Nbf-1)/(2*Nbf));
iflaw=kron(freqs,ones(Ncomp,1)); iflaw=iflaw(1:N,1);
y=exp(j*2*pi*cumsum(iflaw));