Полезные функции от KimIV - страница 47

Новый комментарий
 

так у меня тоже совпадает если вызывать вот так

  for (int i=0; i<r; i++)  {
    Y[i]=Close[i+1];
    X[i]=i;
  }
    
  Array_LR(X, Y);
  for (i=0; i<r; i++) {
    SetArrow(170, Blue, "arr"+i+r, Time[i+1], Y[i]);
  }

Правда точки в обоих вариантах накладываються не точно. Но это скорее всего особенность SetArrow()

Вот картинка

 
Prival писал (а) >>
Правда точки в обоих вариантах накладываються не точно. Но это скорее всего особенность SetArrow()

Нет, это особенность графического объекта OBJ_ARROW. Он привязывается не центром массы, а серединой верхней границы.

 

Функция ArrayMo().

Возвращает Моду - максимум кривой плотности распределения. Функция принимает следующие необязательные параметры:

  • x - Массив значений числового ряда.
  • d - Точность значений числового ряда, количество знаков после запятой.
//+----------------------------------------------------------------------------+
//|  Автор    : Ким Игорь В. aka KimIV,  http://www.kimiv.ru                   |
//+----------------------------------------------------------------------------+
//|  Версия   : 21.06.2008                                                     |
//|  Описание : Возвращает Моду - максимум кривой плотности распределения.     |
//+----------------------------------------------------------------------------+
//|  Параметры:                                                                |
//|    x - массив значений числового ряда                                      |
//|    d - точность значений числового ряда, количество знаков после запятой   |
//+----------------------------------------------------------------------------+
double ArrayMo(double& x[], int d=4) {
  double e, s=0;
  double m[][2];             // временный массив:
                             //  столбец 1 - количество значений
                             //  столбец 2 - значения
  int    i, k=ArraySize(x);
  int    n;                  // номер строки временного массива m
  int    r;                  // количество строк во временном массиве m

  if (k>0) {
    for (i=0; i<k; i++) {
      e=NormalizeDouble(x[i], d);
      n=ArraySearchDouble(m, e);
      if (n<0) {
        r=ArrayRange(m, 0);
        ArrayResize(m, r+1);
        m[r][0]++;
        m[r][1]=e;
      } else m[n][0]++;
    }
    ArraySort(m, WHOLE_ARRAY, 0, MODE_DESCEND);
    s=m[0][1];
  } else Print("ArrayMo(): Массив пуст!");

  return(s);
}
 

Пример использования функции ArrayMo().

Определение наиболее часто встречающегося ценового уровня High среди последних 1000-чи баров текущего графика:

#define R 1000
void start() {
  double a[R];
  for (int i=0; i<R; i++) a[i]=High[i];
  Message(ArrayMo(a, 4));
}
ЗЫ. Во вложении скрипт для тестирования функции ArrayMo().
Файлы:
test_arraymo.mq4  5 kb
 

Полностью опубликована библиотека функций b-Array, предназначенная для работы с массивами.

 

Есть еще одна, расчет ковариации

//+----------------------------------------------------------------------------+
//|  Автор    : Сергей Привалов aka Prival,  Skype: privalov-sv                |
//+----------------------------------------------------------------------------+
//|  Версия   : 11.09.2008                                                     |
//|  Описание : Рассчет ковариации массива cvar(X,Y)                           |
//+----------------------------------------------------------------------------+
//|  Параметры:                                                                |
//|    X    - массив значений числового ряда, ось X                            |
//|    Y    - массив значений числового ряда, ось Y                            |
//+----------------------------------------------------------------------------+

double cvar(double &X[], double &Y[])
{
      double mo_X = 0, mo_Y = 0, res = 0;
      int    i,N=ArraySize(X);
     
      if(N>1 && ArraySize(Y)==N)  {
        for( i = 0; i < N; i ++ ) {
            mo_X +=X[i]-X[0];
          mo_Y +=Y[i];
      }
      mo_X /=N;
      mo_Y /=N;
      for( i = 0; i < N; i ++ ) res +=(X[i]-mo_X)*(Y[i]-mo_Y);
      res /=N;
    } else Print("cvar(): Недостаточное количество элементов ряда! N=", N, " или не совпадает размерность");
   return(res);

исправлено.

 
Prival писал (а) >>

Есть еще одна, расчет ковариации

добавте её в библиотеку. Хотя опрераций с массивами (матрицами) намного больше. Но думаю постепенно наполним.

Есть пару вопросов:

1. что такое mo_XY?

2. В строчке накопления МО по Х

mo_X +=X[i]-X[0];
Зачем отнимать X[0]?
3. Зачем массив Х должен быть упорядоченным?
 

1. mo_XY можно удалить, проверял различные варианты расчета. Это осталось от плохого варианта.

2. Этот алгоритм что я привел имеет минимальные шансы получить ошибку в расчетах если в качестве X будет поступать Time[]. Перемножение больших чисел постепенно вызывают накомление ошибки и она вылазит. Именно для этой цели, устранение возможного появления этой ошибки, X (дополнительно) смещается в начало координат вычитанием X[0]

3. Перемудрил, может быть не упорядочен, главное что бы значения введенные в X соответсвовали Y

Сейчас поправлю

 
mo_X += X[0]; // Забыли наверное.
Это лишная операция. Можете перепроверить.
 
TheXpert писал (а) >>

Не соглашусь.

Это хорошее правило - недоверие. Проверте в любом математическом пакете. Результаты выложим. Я счас это сделаю в MathCade.

1...404142434445464748495051525354...128
Новый комментарий
Причина обращения: