Hi guys!

I have an M-File which models a unit-energy Raised Cosine Pulse;

function y = rcpuls(a,t)
tau = 1; % Set the symbol time
t = t + 10^(-7); % This avoids dividing by zero
tpi = pi/tau; atpi = tpi*a;
at = 4*a^2/tau^2;
y = sin(tpi*t).*cos(atpi*t)./(tpi*t.*(1-at*t.^2));

A typical example would be typing in the command window;

>> rcpulse = rcpuls(0.3,[-6:0.01:6]);
>> rcpulse = rcpulse/norm(rcpulse);

This returns a 30% Raised Cosine Pulse.

My problem is that I need a unit-energy Square Root Raised Cosine Pulse!!

Does anyone know what I need to do to edit this code?

Thanks

Will

Don't Worry found an equation and came up with this;
Works well;

function y = rtrcpuls(a,t)
tau = 1;
t = t+0.0000001
tpi = pi/tau;
amtpi = tpi*(1-a);
aptpi = tpi*(1+a);
ac = 4*a/tau;
at = 16*a^2/tau^2;
y = (sin(amtpi*t) + (ac*t).*cos(aptpi*t))./(tpi*t.*(1-at*t.^2));
y = y/sqrt(tau);

Regards 

Will