function H = RaisedCosSpectrum(L,fsT,beta)
%
% Time-stamp: <96/02/22 10:35:51 pparis>
%
% Computes the spectrum of an ideal raised cosine pulse. beta is the roll-off 
% factor. (see Stuber, page 164) L indicates the maximum frequency (in
% multiples of 1/T) to be evaluated. fsT indicates the number of samples
% per 1/T.
%
% Syntax: H=Raised(L, fsT,beta)

% $Id: RaisedCosSpectrum.m,v 1.1 2003/09/11 19:08:59 pparis Exp $

f = -L:1/fsT:L;

H= 10^(-14)*ones(size(f));

i1 = find(abs(f) <= (1-beta)/2);
H(i1) = 1;

i2 = find(abs(f) < (1+beta)/2 & abs(f) > (1-beta)/2);
H(i2) = (1-sin(pi/beta*(abs(f(i2))-0.5)))/2;
