Ehlers - SineTrend

["Henry Amand" <h.amand@xxxxxxxxx>]

Hello,

was reading the book by Ehlers and he had a SineTrend system in there.
Produced some nice results and i typed it into PowerEditor.

However, it does not seem to work. Ehlers made some changes for
it on his website, wich i have included, but still it does not work.
No signals are triggered.


anyone any idea where this code might go wrong ? I unfortunatly
can not seem to find it...

greetings



-----------------------------------------------------------------

{***************************************

 Written by: Ehlers - Rocket science for traders
             typed by Henry Amand

 Description: SineTrend system, page 120

****************************************}

Inputs: Price((H+L)/2);

Vars: Smooth (0),
      Detrender (O) , 
      I1 (0), 
      Q1 (O), 
      jI (O), 
      JQ (0), 
      I2 (0), 
      Q2 (0), 
      Re (0), 
      Im (0), 
      Period (0), 
      SmoothPeriod (O), 
      SmoothPrice (O), 
      DCPeriod (O), 
      RealPart (0), 
      Imagpart (O), 
      count (0), 
      DCPhase (O), 
      DCSine (O), 
      LeadSine (O), 
      Itrend (O), 
      Trendline (O), 
      Trend (0), 
      DaysinTrend (O);



If CurrentBar  >  5  then begin


Smooth =   (4*Price  +  3*Price[1]   + 2*Price[2]   + Price[3])   /  10; 
Detrender =   (.0962 * Smooth + .5769 * Smooth[2] - .5769 * 
Smooth[4] - .0962*Smooth[6]) * (.075 * Period[1] + .54);


{Compute  InPhase and Quadrature components.} 

Q1  =   (.0962 * Detrender + .5769 * Detrender[2] - .5769 
* Detrender[4] - .0962 * Detrender[6]) * (.075 * Period[1] + .54); 
I1  = Detrender[3];


{Advance the phase of  I1 and Q1 by 90 degrees} 

jI = (.0962 * I1 + .5769 * I1[2] - .5769 * I1[4] - 
.0962 * I1[6]) * (.075 * Period[1] + .54); 
JQ = (.0962 * Q1 + .5769 * Q1[2] - .5769 * Q1[4] - 
.0962 * Q1[6]) * (.075 * Period[1] + .54);


{Phasor addition for 3 bar averaging)} 

I2 = I1 - JQ; 
Q2 = Q1 + jI;


{Smooth the  I and Q components before applying 
the discriminator} 

I2 = .2 * I2 + .8 * I2[1];
Q2 = .2 * Q2 + .8 * Q2[1];


{Homodyne Discriminator}

Re = I2 * I2[1] + Q2 * Q2[1];
Im = I2 * Q2[1] - Q2 * I2[1];
Re = .2 * Re + .8 * Re[1];
Im = .2 * Im + .8 * Im[1];

If Im <> 0 and Re <> 0 then Period = 360 / ArcTangent (Im/Re); 
If Period > 1.5 * Period[1] then Period = 1.5 * Period[1];
If Period < .67 * Period[1] then Period = .67 * Period[1];
If Period < 6 then Period = 6;
If Period > 50 then Period = 50;
Period = .2 * Period + .8 * Period[1];
SmoothPeriod = .33 * Period + .67 * SmoothPeriod[1];


{Compute Dominant Cycle  Phase}

SmoothPrice = (4 * price + 3 * Price[1] + 2 * Price[2] + 
Price[3]) / 10;
DCPeriod = IntPortion(SmoothPeriod + .5); 
RealPart = 0; 
ImagPart = 0; 


For count = 0 To DCPeriod - 1 begin 

RealPart = RealPart + sine (360 * count / DCPeriod) * 
(SmoothPrice[count]);
ImagPart = imagPart + CoSine (360 * count / DCPeriod) * 
(SmoothPrice[count]); 

End; 

If AbsValue(ImagPart) > 0 then DCPhase = 
Arctangent(RealPart / ImagPart); 
If AbsValue(ImagPart) <= .001 then DCPhase = 
DCPhase + 90 * Sign(RealPart);
DCPhase = DCPhase + 90;


{Compensate for one bar lag of the Weighted Moving Average}

DCPhase = DCPhase + 360 / SmoothPeriod;

If ImagPart < 0 then DCPhase = DCPhase + 180;
If DCPhase > 315 then DCPhase = DCPhase - 360;


{Compute the Sine and LeadSine Indicators}

DCSine = Sine(DCPhase); 
LeadSine = Sine(DCPhase + 45);


{Compute Trendline as simple average over the measured 
dominant cycle period} 

ITrend = 0; 

For count = 0 to DCPeriod - 1 begin
ITrend = ITrend + Price[count]; 
End;

If DCPeriod > 0 then ITrend = ITrend / DCPeriod; 
Trendline = (4 * ITrend + 3 * ITrend[1] + 2 * ITrend[2] 
+ ITrend[3]) / 10; 
If CurrentBar < 12 then Trendline = Price;


{Assume Trend Mode} 

Trend = 1;


{Measure days in trend from last crossing of the 
Sinewave Indicator lines}

If Sine(DCPhase) Crosses Over Sine(DCPhase + 45) or 
Sine(DCPhase) Crosses Under Sine(DCPHase + 45) Then begin
DaysInTrend = 0; 
Trend = 0;
End;

DaysInTrend = DaysInTrend + 1;
If DaysInTrend < .5 * SmoothPeriod then Trend = 0;


{Cycle Mode If delta phase is +/- 50% of dominant cycle 
change of phase}

If SmoothPeriod <> 0 and (DCPhase - DCPhase[1] > .67 * 360 / 
SmoothPeriod and DCPhase - DCPhase[1] < 1.5 * 360 / 
SmoothPeriod) then Trend = 0;


{Declare a Trend Mode if the SmoothPrice is more than 1.5% 
from the Trendline} 

If AbsValue ((SmoothPrice - Trendline) / Trendline) >= 
.015 then Trend = 1;

If Trend = 1 then begin

If Trend [1] = 0 then begin

If MarketPosition = -1 and SmoothPrice >= Trendline then buy;
If MarketPosition = 1 and SmoothPrice < Trendline then sell;

End;

If SmoothPrice Crosses Over Trendline then buy;
If SmoothPrice Crosses Under Trendline then sell;
End;

If Trend = 0 then begin

If LeadSine Crosses Over DCSine then buy;
If LeadSine Crosses Under DCSine then sell;

End;
End;