center of gravity macd

Viewing 3 posts - 1 through 3 (of 3 total)

Author

Posts
11/16/2023 at 9:14 AM #223783 quote

ARLEQUIN49
Participant

Veteran

Bonjour à tous,

J’aurais besoin d’aide pour convertir cet indicateur de Métatrader 4 qui est le centre de gravité avec le MACD qui est assez bien efficient.

Merci d’avance

//+——————————————————————+
//| COGMACD.mq4 |
//| Original Code from NG3110@latchess.com |
//| Linuxser 2007 for TSD http://www.forex-tsd.com/ |
//| Mod by Brooky @ Brooky-Indicators.com |
//+——————————————————————+
#property copyright “ANG3110@latchess.com”
//———ang_pr (Din)——————–
#property indicator_separate_window
#property indicator_buffers 8
#property indicator_color1 RoyalBlue
#property indicator_color2 Red
#property indicator_color3 Green
#property indicator_color4 Orange
#property indicator_color5 LimeGreen
#property indicator_color6 Olive
#property indicator_color7 DarkBlue
#property indicator_color8 Tomato

#property indicator_level1 0
#property indicator_levelstyle 0
#property indicator_levelcolor DimGray

#property indicator_width1 2
#property indicator_width2 2
#property indicator_width3 2
#property indicator_width4 1
#property indicator_width5 1
#property indicator_width6 2
#property indicator_width7 2
#property indicator_width8 2

#property indicator_style1 0
#property indicator_style2 0
#property indicator_style3 0
#property indicator_style4 2
#property indicator_style5 2
#property indicator_style6 0
#property indicator_style7 0
#property indicator_style8 0
//———————————–
extern int bars_back = 240;
extern int fma = 12;
extern int sma = 26;
extern int sigma = 9;
extern double fib.band = 0.618;

extern int m = 3;
extern int i = 0;
extern double kstd = 3.618;

//———————–
double fx[], sqh[], sql[], stdh[], stdl[], stochdata[],stochdata2[], stochsdata[];
double ai[10,10], b[10], x[10], sx[20];
double sum;
int ip, p, n, f;
double qq, mm, tt;
int ii, jj, kk, ll, nn;
double sq, std;
//*******************************************
int init()
{
IndicatorShortName(“COGMACD: Mod by Brooky-Indicators.com”);

SetIndexBuffer(0, fx);SetIndexStyle(0, DRAW_LINE);
SetIndexBuffer(1, sqh);
SetIndexBuffer(2, sql);
SetIndexBuffer(3, stdh);
SetIndexBuffer(4, stdl);
SetIndexBuffer(5, stochdata);SetIndexStyle(5, DRAW_HISTOGRAM);
SetIndexBuffer(6, stochsdata);SetIndexStyle(6, DRAW_LINE);
SetIndexBuffer(7, stochdata2);SetIndexStyle(7, DRAW_HISTOGRAM);
p = MathRound(bars_back);

nn = m + 1;

return(0);
}
//———————————————————-
int deinit()
{

}
//**********************************************************************************************
int start()
{
int mi;
//——————————————————————————————-

p = bars_back;
sx[1] = p + 1;
SetIndexDrawBegin(0, Bars – p – 1);
SetIndexDrawBegin(1, Bars – p – 1);
SetIndexDrawBegin(2, Bars – p – 1);
SetIndexDrawBegin(3, Bars – p – 1);
SetIndexDrawBegin(4, Bars – p – 1);
//———————-sx——————————————————————-

int rlimit;
int rcounted_bars=IndicatorCounted();
//—- check for possible errors
if(rcounted_bars<0) return(-1);
//—- the last counted bar will be recounted
if(rcounted_bars>0) rcounted_bars–;
rlimit=Bars-rcounted_bars;
//—- main loop
for(int ri=0; ri<rlimit; ri++)
{
if(iMACD(NULL,0,fma,sma,sigma,PRICE_CLOSE,MODE_MAIN,ri)>0)stochdata[ri]= iMACD(NULL,0,fma,sma,sigma,PRICE_CLOSE,MODE_MAIN,ri);
else stochdata2[ri]= iMACD(NULL,0,fma,sma,sigma,PRICE_CLOSE,MODE_MAIN,ri);
stochdata[ri]= iMACD(NULL,0,fma,sma,sigma,PRICE_CLOSE,MODE_MAIN,ri);//iStochastic(NULL,0,fma,sma,sigma,MODE_SMA,0,MODE_MAIN,ri);
stochsdata[ri]= iMACD(NULL,0,fma,sma,sigma,PRICE_CLOSE,MODE_SIGNAL,ri);//iStochastic(NULL,0,fma,sma,sigma,MODE_SMA,0,MODE_SIGNAL,ri);
}

for(mi = 1; mi <= nn * 2 – 2; mi++)
{
sum = 0;
for(n = i; n <= i + p; n++)
{
sum += MathPow(n, mi);
}
sx[mi + 1] = sum;
}
//———————-syx———–
for(mi = 1; mi <= nn; mi++)
{

sum = 0.00000;
for(n = i; n <= i + p; n++)
{

if(mi == 1)
sum += ((iMACD(NULL,0,fma,sma,sigma,PRICE_CLOSE,MODE_MAIN,n)+iMACD(NULL,0,fma,sma,sigma,PRICE_CLOSE,MODE_SIGNAL,n))+0.0000001)/2;//rsi_period iRSI(NULL,0,rsi_period,prICE_CLOSE,n)
else
sum += (((iMACD(NULL,0,fma,sma,sigma,PRICE_CLOSE,MODE_MAIN,n)+iMACD(NULL,0,fma,sma,sigma,PRICE_CLOSE,MODE_SIGNAL,n))+0.0000001)/2) * MathPow(n, mi – 1);
}
b[mi] = sum;
}
//===============Matrix=======================================================================================================
for(jj = 1; jj <= nn; jj++)
{
for(ii = 1; ii <= nn; ii++)
{
kk = ii + jj – 1;
ai[ii, jj] = sx[kk];
}
}
//===============Gauss========================================================================================================
for(kk = 1; kk <= nn – 1; kk++)
{
ll = 0; mm = 0;
for(ii = kk; ii <= nn; ii++)
{
if(MathAbs(ai[ii, kk]) > mm)
{
mm = MathAbs(ai[ii, kk]);
ll = ii;
}
}
if(ll == 0)
return(0);

if(ll != kk)
{
for(jj = 1; jj <= nn; jj++)
{
tt = ai[kk, jj];
ai[kk, jj] = ai[ll, jj];
ai[ll, jj] = tt;
}
tt = b[kk]; b[kk] = b[ll]; b[ll] = tt;
}
for(ii = kk + 1; ii <= nn; ii++)
{
qq = ai[ii, kk] / ai[kk, kk];
for(jj = 1; jj <= nn; jj++)
{
if(jj == kk)
ai[ii, jj] = 0;
else
ai[ii, jj] = ai[ii, jj] – qq * ai[kk, jj];
}
b[ii] = b[ii] – qq * b[kk];
}
}
x[nn] = b[nn] / ai[nn, nn];
for(ii = nn – 1; ii >= 1; ii–)
{
tt = 0;
for(jj = 1; jj <= nn – ii; jj++)
{
tt = tt + ai[ii, ii + jj] * x[ii + jj];
x[ii] = (1 / ai[ii, ii]) * (b[ii] – tt);
}
}
//===========================================================================================================================
for(n = i; n <= i + p; n++)
{
sum = 0;
for(kk = 1; kk <= m; kk++)
{
sum += x[kk + 1] * MathPow(n, kk);
}
fx[n] = x[1] + sum;
}
//———————————–Std———————————————————————————–
sq = 0.0;
for(n = i; n <= i + p; n++)
{
sq += MathPow((((iMACD(NULL,0,fma,sma,sigma,PRICE_CLOSE,MODE_MAIN,n)+iMACD(NULL,0,fma,sma,sigma,PRICE_CLOSE,MODE_SIGNAL,n))+0.0000001)/2) – fx[n], 2);
}
sq = MathSqrt(sq / (p + 1)) * kstd;
std = iStdDevOnArray(stochdata,0,p,0,MODE_SMA,i) * kstd;
for(n = i; n <= i + p; n++)
{
sqh[n] = fx[n] + sq;
sql[n] = fx[n] – sq;
stdh[n] = fx[n] + (fib.band*std);
stdl[n] = fx[n] – (fib.band*std);
}
//——————————————————————————-
//ObjectMove(“sstart” + sName, 0, Time[p], fx[p]);
//—————————————————————————————————————————-
return(0);
}
//==========================================================================================================================

Graph-MT4.png

07/23/2024 at 12:29 PM #235662 quote
Iván González
Moderator

Legend
Bonjour Voici ce que j’ai réalisé :

Center of Gravity of the MACD

Ca a l’air pas mal mais c’est très lourd…
```
//----------------------------------------------------------------//
// PRC COGMACD
//version = 0
//23.07.2024
//Iván González @ www.prorealcode.com
//Sharing ProRealTime knowledge
//----------------------------------------------------------------//
//-----Inputs-----------------------------------------------------//
m = 3
p = 240

fma = 12
sma = 26
sigma = 9

fibband = 0.618
kstd = 3.618
//----------------------------------------------------------------//
//-----MACD calculation-------------------------------------------//
// Calculate the MACD and Signal line
myMacd = macd[fma, sma, sigma](close)
signal = MACDSignal[fma, sma, sigma](close)
MacdMiddle = (myMacd + signal) / 2
//----------------------------------------------------------------//
nn = m + 1

// Calculate the sx array which contains the sum of powers of n
for mi = 1 to nn * 2 - 2 do
sum = 0
for n = 0 to p do
sum = sum + pow(n, mi)
next
$sx[mi + 1] = sum
next
//----------------------------------------------------------------//
// Calculate the b array using MacdMiddle values
for mi = 1 to nn do
sum = 0
for n = 0 to p do
if mi = 1 then
sum = sum + MacdMiddle[n]
else
sum = sum + MacdMiddle[n] * pow(n, mi - 1)
endif
next
$b[mi] = sum
next
//----------------------------------------------------------------//
// Populate the ai matrix using the sx values
for jj = 1 to nn do
for ii = 1 to nn do
kk = ii + jj - 1
$ai[(ii - 1) * nn + jj] = $sx[kk]
next
next
//----------------------------------------------------------------//
// Apply Gauss method to solve the system of equations
for kk = 1 to nn - 1 do
ll = 0
mm = 0
qq = 0
tt = 0
for ii = kk to nn do
if abs($ai[(ii - 1) * nn + kk]) > mm then
mm = abs($ai[(ii - 1) * nn + kk])
ll = ii
endif
next
if ll <> kk then
for jj = 1 to nn do
tt = $ai[(kk - 1) * nn + jj]
$ai[(kk - 1) * nn + jj] = $ai[(ll - 1) * nn + jj]
$ai[(ll - 1) * nn + jj] = tt
next
tt = $b[kk]
$b[kk] = $b[ll]
$b[ll] = tt
endif
for ii = kk + 1 to nn do
qq = $ai[(ii - 1) * nn + kk] / $ai[(kk - 1) * nn + kk]
for jj = 1 to nn do
if jj = kk then
$ai[(ii - 1) * nn + jj] = 0
else
$ai[(ii - 1) * nn + jj] = $ai[(ii - 1) * nn + jj] - qq * $ai[(kk - 1) * nn + jj]
endif
next
$b[ii] = $b[ii] - qq * $b[kk]
next
next
//----------------------------------------------------------------//
// Calculate the coefficients of the polynomial using back substitution
$x[nn] = $b[nn] / $ai[(nn - 1) * nn + nn]

for ii = nn - 1 downto 1 do
tt = 0
for jj = 1 to nn - ii do
tt = tt + $ai[(ii - 1) * nn + (ii + jj)] * $x[ii + jj]
next
$x[ii] = (1 / $ai[(ii - 1) * nn + ii]) * ($b[ii] - tt)
next
//----------------------------------------------------------------//
// Calculate the Center of Gravity (COG) and the standard deviation bands
sq = 0.0

for n = 0 to p do
sum = 0
for kk = 1 to m do
sum = sum + $x[kk + 1] * pow(n, kk)
next

// Center of Gravity calculation
$fx[n] = $x[1] + sum

// Accumulate the squared differences for standard deviation
sq = sq + pow(MacdMiddle[n] - $fx[n], 2)

// Calculate standard deviation
st = sqrt(sq / (p + 1)) * kstd

// Calculate the bands based on standard deviation
$sqh[n] = $fx[n] + st
$sql[n] = $fx[n] - st
$stdh[n] = $fx[n] + (fibband * st)
$stdl[n] = $fx[n] - (fibband * st)
next
//----------------------------------------------------------------//
return $fx[p] as "Center of Gravity" coloured("blue") style(line,2), $sqh[p] as "Upper Band" coloured("red") style(line,2), $sql[p] as "Lower Band" coloured("green") style(line,2), $stdh[p] as "Upper Fibonacci Band" coloured("orange") style(dottedline,1), $stdl[p] as "Lower Fibonacci Band" coloured("orange") style(dottedline,1), myMacd coloured("yellow") style(line,2), 0 as "zero"
```
07/23/2024 at 4:32 PM #235676 quote

ARLEQUIN49
Participant

Veteran

Merci beaucoup
Author

Posts