VPM (Все сообщения пользователя)

function Cached.SWISS_ARMY_KNIFE_INDICATOR()--SWISS ARMY KNIFE INDICATOR (John Ehlers)

--[[ выполняет все обычные функции обычных индикаторов, такие как сглаживание и генерация импульса.
 Также делает некоторые необычные вещи, такие как остановка полосы и фильтрация отклонения полосы. 
 
 Как только вы запрограммируете этот индикатор, вы сможете использовать практически любой метод технического анализа. 
 Этот уникальный общий индикатор является результатом общих концепций цифровой обработки сигналов (DSP) для дискретных сетей сигналов,
 которые появляются в различных формах в техническом анализе.

 Передаточная функция дискретной линейной системы - это отношение выхода системы, деленное на вход. 
 --]]
    --local rounding=rounding
    --local Value=Value;

    local pi=math.pi;
   local exp=math.exp;
   local cos=math.cos;
   local sin=math.sin;
   local sqrt=math.sqrt;
   local floor=math.floor;
        local Filt,Price={},{};
   --local fPeriodDC=Cached.PeriodDC();
   
   return function(I,FSettings,ds )
      local I = I or 1;
      local ds = ds or nil;
      local FSettings = FSettings or {};
      local v_t = FSettings.v_t or 'C';
      local P = FSettings.period or 48;
      local N = FSettings.N or 8;
      local n = FSettings.n or 0;
      local wid = FSettings.wid or "BP";
      local delta1 = FSettings.delta1 or 0.1;

      local avtoP = FSettings.avtoP or 0;
      local rev = FSettings.rev or 1;
      local wg = FSettings.bw_hp or 1; --полоса пропускания
      local bw_ss = FSettings.bw_ss or 1;
      local metod = FSettings.metod or 2; -- адаптации period данных
      --Settings.pole_hp = 0
      --Settings.pole_ss = 2
      local x = FSettings.x or nil;
      --------------------------------------
      local price = x or Value(I, v_t, ds) or 0;
      --------------------------------------
      if I == 1 then 
         Price,Filt={},{};
         Price[I],Filt[I]=p0,p0;
         --fPeriodDC=Cached.PeriodDC()
      end
      --------------------------------------
      local p1=Price[I-1] or price;
      local p2=Price[I-2] or p1;
      local pn=Price[I-n] or p2;
      Price[I] = price;

        --Vars:
        local c0=1;
        local c1=0;
        local b0=1;
        local b1=0;
        local b2=0;
        local a1=0;
        local a2=0;
        local alpha=0;
        local beta1=0;
        local gamma1=0;

      ---------------------------------
      if avtoP == 1  then 
         local DC = fPeriodDC(I, {['v_t']=v_t,['period']=P,['N']=N,['wg']=wg,['bw']=bw_ss,['metod']=metod}, ds)
         --if DC < 6 then DC=6 end
         --if DC > P then DC=P end
         P = round(floor(rev*DC),0)
      end
      ---------------------------------

      local filt=0;
      local filt1=Filt[I-1] or 0;
      local filt2=Filt[I-2] or filt1;
      local fi=2*pi/P;

  if wid == "EMA" then --INDICATOR
    if I <= n then filt,filt1 = price,price; end
      alpha = (cos(fi) + sin(fi) - 1) / cos(fi);
      b0 = alpha;
      a1 = 1 - alpha;
  
  elseif wid == "SMA" then --INDICATOR
     if I <= n then filt,filt1 = price,price; end
        c1 = 1 / n;
        b0 = 1 / n;
        a1 = 1;
      
  elseif wid == "Smooth" then --INDICATOR
        if I <= n then filt,filt1 = price,price; end
      c0 = 1 / 4;
        b1 = 2;
        b2 = 1; 
      
  elseif wid == "Gauss" then --INDICATOR 
        -- фильтр Гаусса является фильтром нижних частот, 
      -- позволяя проходить более низким частотам и ослаблять более высокие частоты.
      -- Использование гауссовых фильтров - это движение к двойной цели: 
      -- уменьшить отставание и уменьшить отставание высокочастотных компонентов относительно отставания низкочастотных компонентов. 
    if I <= n then filt,filt1 = price,price; end
        beta1 = 2.415*(1 - cos(fi));
        alpha = -beta1 + sqrt(beta1*beta1 + 2*beta1);
        c0 = alpha*alpha;
        a1 = 2*(1 - alpha);
        a2 = -(1 - alpha)*(1 - alpha);
  
  elseif wid == "Butter" then --INDICATOR 
        -- Сглаживание, реализованное с помощью этого фильтра, является чрезвычайно скромным, отклоняя только компонент цикла 2 бара.
      if I <= n then filt,filt1 = price,price; end
        beta1 = 2.415*(1 - cos(fi));
        alpha = -beta1 + sqrt(beta1*beta1 + 2*beta1);
        c0 = alpha*alpha / 4;
        b1 = 2;
        b2 = 1;
        a1 = 2*(1 - alpha);
        a2 = -(1 - alpha)*(1 - alpha);
  
  elseif wid == "HP" then 
     if I <= n then filt,filt1 = 0,0; end
       alpha = (cos(fi) + sin(fi) - 1) / cos(fi);
       c0 = 1 - alpha / 2;
       b1 = -1;
       a1 = 1 - alpha;
  
  elseif wid == "2PHP" then 
        --является фильтром высоких частот, позволяя проходить более высоким частотам и ослаблять более низкие частоты.
     if I <= n then filt,filt1 = 0,0; end
      beta1 = 2.415*(1 - cos(fi));
        alpha = -beta1 + sqrt(beta1*beta1 + 2*beta1);
        c0 = (1 - alpha / 2)*(1 - alpha / 2);
      b1 = -2;
        b2 = 1;
        a1 = 2*(1 - alpha);
        a2 = -(1 - alpha)*(1 - alpha);
  
  elseif wid == "BP" then 
     if I <= n then filt,filt1 = price,price; end
      beta1 = cos(fi);
      gamma1 = 1 / cos(4*pi*delta1 / P);
      alpha = gamma1 - sqrt(gamma1*gamma1 - 1);
      c0 = (1 - alpha) / 2;
      b2 = -1;
      a1 = beta1*(1 + alpha);
      a2 = -alpha;
  
  elseif wid == "BS" then --INDICATOR
     if I <= n then filt,filt1 = price,price; end
        beta1 = cos(fi);
        gamma1 = 1 / cos(4*pi*delta1 / P);
        alpha = gamma1 - sqrt(gamma1*gamma1 - 1);
      c0 = (1 + alpha) / 2;
        b1 = -2*beta1;
        b2 = 1;
        a1 = beta1*(1 + alpha);
        a2 = -alpha;
  end;

    if I > n then 
    filt = c0*(b0*price + b1*p1 + b2*p2) + a1*filt1 + a2*filt1 - c1*pn;
    end;
   Filt[I] = filt;

   return filt1,filt
   end
end