Сбросьте версию AlgLib TestInterface.

Arthur Albano 2021.10.07 17:28 #41

Привет!

Я совсем новичок в ALGLIB. Итак, я пытался выполнить следующую задачу - вычислить k² с помощью скрипта. Я могу легко сделать это в Excel, но не представляю, как использовать нелинейную модель уравнений с одним боксовым ограничением.

void OnStart()
  {
   const double x[22] = {28.98, 30.48, 30.98, 31.48, 31.98, 32.48, 32.98, 33.48, 33.98, 34.48, 34.98, 35.48, 35.98, 36.48, 36.98, 37.48, 37.98, 38.48, 38.98, 39.98, 40.98, 41.98};
   const double y[22] = {8.26, 7.26, 5.71, 6.15, 5.90, 5.15, 5.03, 4.30, 4.15, 3.30, 2.90, 2.71, 1.96, 1.57, 1.27, 0.90, 0.70, 0.47, 0.33, 0.20, 0.09, 0.06};
   const double w[22] = {2, 8, 2, 2, 2, 1, 11, 12, 24, 12, 6, 12, 21, 7, 74, 51, 159, 385, 171, 109, 18, 18}; //весы
   const double x0 = 37.74; //теорема о центральной стоимости

//модельное уравнение
//y = (1/2)*((x^2-2x0x+x0^2+4k^2)^(1/2) + x+x0) - x
//k^2 был вычислен как 0.582018800686398, после использования Solver в Excel, где общая ошибка была минимизирована суммой w(y-y0)^2
// ограничение на коробку k^2>0;
  };

Заранее спасибо за прочтение!

Спрашивайте! Запуск скрипта через советник Scalp_net

Arthur Albano 2021.10.09 22:36 #42

Arthur Albano #:

Привет!

Я совсем новичок в ALGLIB. Итак, я пытался выполнить следующую задачу - вычислить k² с помощью скрипта. Я могу легко сделать это в Excel, но не представляю, как использовать нелинейную модель уравнений с одним боксовым ограничением.

Заранее спасибо за чтение!

Пока что выглядит уродливо...

//+------------------------------------------------------------------+
//|AlgLibExample.mq5 |
//|Артур Альбано |
//| https://www.mql5.com/ru/users/arthuralbano |
//+------------------------------------------------------------------+
#property copyright "Arthur Albano"
#property link      "https://www.mql5.com/ru/users/arthuralbano"
#property version   "1.00"
//+------------------------------------------------------------------+
//| Функция запуска программы сценария|
//+------------------------------------------------------------------+

#include <Math\Alglib\alglib.mqh>
#include <Math\Alglib\interpolation.mqh>
#include <Arrays\ArrayDouble.mqh>

//+------------------------------------------------------------------+
//||
//+------------------------------------------------------------------+
void OnStart()
  {
   const double x0 = 37.74; //теорема о центральной стоимости

//модельное уравнение
//f(x,k²,x0) = (1/2)*((x^2-2x0x+x0^2+4k^2)^(1/2) + x+x0) - x
//k^2 был рассчитан как 0.582018800686398, после использования Solver в Excel
// box contraint k^2>0;

// https://www.alglib.net/translator/man/manual.cpp.html#example_lsfit_d_nlfb

   const double x[22] = {28.98, 30.48, 30.98, 31.48, 31.98, 32.48, 32.98, 33.48, 33.98, 34.48, 34.98, 35.48, 35.98, 36.48, 36.98, 37.48, 37.98, 38.48, 38.98, 39.98, 40.98, 41.98};
   double w[22] = {2, 8, 2, 2, 2, 1, 11, 12, 24, 12, 6, 12, 21, 7, 74, 51, 159, 385, 171, 109, 18, 18}; //весы
   double y[22] = {8.26, 7.26, 5.71, 6.15, 5.90, 5.15, 5.03, 4.30, 4.15, 3.30, 2.90, 2.71, 1.96, 1.57, 1.27, 0.90, 0.70, 0.47, 0.33, 0.20, 0.09, 0.06};
   
   const int n = ArraySize(x);

   double c[1] = {x0 / 100.};
   CLSFitState state;

   CMatrixDouble xx(n, 1);
   for(int i = 0; i < n; i++)
     {
      xx[i].Set(0,x[i]);
     };

   CLSFit::LSFitCreateWFGH(xx, y, w, c, n,1,1,state);

   const double epsf = 0.000001;
   const double epsx = 0.000001;
   const int maxits = 0;
   CLSFit::LSFitSetCond(state,epsf,epsx,maxits);
   
   double bndl[1] = {0.};
   double bndu[1] = {x0 / 100.};
   CLSFit::LSFitSetBC(state,bndl,bndu);
   
   CNDimensional_PFunc *Pfunc = new CNDimensional_PFunc;

   //CAlglib::LSFitFit(state,Pfunc,,,,rep,state);
   
   int info;
   CLSFitReport rep;
   CLSFit::LSFitResults(state,info,c,rep);
  };
//+------------------------------------------------------------------+
//||
//+------------------------------------------------------------------+
double func(const double _x, const double _k, const double _x0)
  {
   return((1. / 2.) * (pow(pow(_x, 2.) - 2 * _x0 * _x + pow(_x0, 2.) + 4 * pow(_k, 2.), 1. / 2.) + _x + _x0) - _x);
  };
//+------------------------------------------------------------------+

Arthur Albano 2022.01.11 02:23 #43

Готово.

В основном копируем и вставляем из AlgLib TestInterface.mqh.

//+------------------------------------------------------------------+
//|AlgLibTest.mq5 |
//|Copyright 2022, MetaQuotes Ltd. |
//| https://www.mql5.com |
//+------------------------------------------------------------------+
#property copyright "Copyright 2022, MetaQuotes Ltd."
#property link      "https://www.mql5.com"
#property version   "1.00"
//+------------------------------------------------------------------+
//| Функция запуска программы сценария|
//+------------------------------------------------------------------+

#include <Math\Alglib\alglib.mqh>

//+------------------------------------------------------------------+
//||
//+------------------------------------------------------------------+
void OnStart()
  {
   Print(TEST_LSFit_D_NLF());
  };

//+------------------------------------------------------------------+
//| Производный класс от CNDimensional_PFunc|
//+------------------------------------------------------------------+
class CNDimensional_CX_1_Func : public CNDimensional_PFunc
  {
public:
   //--- конструктор, деструктор
                     CNDimensional_CX_1_Func(void);
                    ~CNDimensional_CX_1_Func(void);
   //--- метод
   virtual void      PFunc(double &c[], double &x[], double &func, CObject &obj);
  };
//+------------------------------------------------------------------+
//| Конструктор без параметров|
//+------------------------------------------------------------------+
CNDimensional_CX_1_Func::CNDimensional_CX_1_Func(void)
  {

  }
//+------------------------------------------------------------------+
//| Деструктор|
//+------------------------------------------------------------------+
CNDimensional_CX_1_Func::~CNDimensional_CX_1_Func(void)
  {

  }
//+------------------------------------------------------------------+
//| Этот обратный вызов вычисляет f(c,x), где x - |
//| положение на оси X, а c - регулируемый параметр |
//+------------------------------------------------------------------+
void CNDimensional_CX_1_Func::PFunc(double &c[], double &x[], double &func, CObject &obj)
  {
   double x0 = 37.74;
   func = (1.0/2.0)*(sqrt(x[0]*x[0] - 2.0*x[0]*x0 + x0*x0 + 4*c[0]) + x[0] + x0) - x[0];
  }

//+------------------------------------------------------------------+
//| Нелинейная подгонка с использованием только значения функции |
//+------------------------------------------------------------------+
double TEST_LSFit_D_NLF()
  {
//--- создайте переменные
   CMatrixDouble           x;
   double                  y[];
   double                  c[];
   double                  w[];
   CObject                 obj;
   CNDimensional_CX_1_Func fcx1func;
   CNDimensional_Rep       frep;
   x.Resize(22,1);
//--- инициализация
   x[0].Set(0,28.98);
   x[1].Set(0,30.48);
   x[2].Set(0,30.98);
   x[3].Set(0,31.48);
   x[4].Set(0,31.98);
   x[5].Set(0,32.48);
   x[6].Set(0,32.98);
   x[7].Set(0,33.48);
   x[8].Set(0,33.98);
   x[9].Set(0,34.48);
   x[10].Set(0,34.98);
   x[11].Set(0,35.48);
   x[12].Set(0,35.98);
   x[13].Set(0,36.48);
   x[14].Set(0,36.98);
   x[15].Set(0,37.48);
   x[16].Set(0,37.98);
   x[17].Set(0,38.48);
   x[18].Set(0,38.98);
   x[19].Set(0,39.98);
   x[20].Set(0,40.98);
   x[21].Set(0,41.98);
//--- распределение
   ArrayResize(y,22);
//--- инициализация
   y[0]=8.26;
   y[1]=7.26;
   y[2]=5.71;
   y[3]=6.15;
   y[4]=5.90;
   y[5]=5.15;
   y[6]=5.03;
   y[7]=4.30;
   y[8]=4.15;
   y[9]=3.30;
   y[10]=2.90;
   y[11]=2.71;
   y[12]=1.96;
   y[13]=1.57;
   y[14]=1.27;
   y[15]=0.90;
   y[16]=0.70;
   y[17]=0.47;
   y[18]=0.33;
   y[19]=0.20;
   y[20]=0.09;
   y[21]=0.06;
//--- распределение
   ArrayResize(c,1);
//--- инициализация
   c[0] =  0.5;
//--- создайте переменные
   int maxits=0;
   int info;
   CLSFitStateShell state;
   double diffstep=0.0001;
//--- Подгонка по весам
   CAlglib::LSFitCreateF(x,y,c,diffstep,state);
//--- вызов функции
   double epsf = 0.0000001;
   double epsx = 0.0000001;
   CAlglib::LSFitSetCond(state,epsf,epsx,maxits);
//--- вызов функции
   CAlglib::LSFitFit(state,fcx1func,frep,0,obj);
//--- вызов функции
   CLSFitReportShell rep;
   CAlglib::LSFitResults(state,info,c,rep);
//--- Подгонка с весами
//--- (вы можете изменить вес и посмотреть, как изменится результат)
   ArrayResize(w,22);
//--- инициализация
   w[0]=2;
   w[1]=8;
   w[2]=2;
   w[3]=2;
   w[4]=2;
   w[5]=1;
   w[6]=11;
   w[7]=12;
   w[8]=24;
   w[9]=12;
   w[10]=6;
   w[11]=12;
   w[12]=21;
   w[13]=7;
   w[14]=74;
   w[15]=51;
   w[16]=159;
   w[17]=385;
   w[18]=171;
   w[19]=109;
   w[20]=18;
   w[21]=18;
//--- Очистить весы
// ArrayFill(w,0,ArraySize(w),1);
//--- вызов функции
   CAlglib::LSFitCreateWF(x,y,w,c,diffstep,state);
//--- вызов функции
   CAlglib::LSFitSetCond(state,epsf,epsx,maxits);
//--- вызов функции
   CAlglib::LSFitFit(state,fcx1func,frep,0,obj);
//--- вызов функции
   CAlglib::LSFitResults(state,info,c,rep);
   return(c[0]);
  }
//+------------------------------------------------------------------+

Yevgeniy Koshtenko 2022.06.03 23:00 #44

Сбросьте пожалуйста свежую версию, не компилируется!(

Renat Fatkhullin 2022.06.05 10:25 #45

Koshtenko #:
Сбросьте пожалуйста свежую версию, не компилируется!(

Используйте штатную версию Alglib, которая давно уже идет в поставке терминала. Класс комплексных чисел в нем адаптирован:

Daying Cao 2022.10.14 16:51 #46

qingyouwei #:
Как перемножить две матрицы с помощью lib alglib?

Как перемножить две матрицы с помощью lib alglib?

Можете ли вы скомпилировать ALGLIB и передать его?

При компиляции возникла следующая ошибка: 'complex' - unexpected token complex.mqh 37 8

'complex' - unexpected token complex.mqh 37 8

'{' - ожидается имя complex.mqh 38 3

'complex' - неожиданная лексема, возможно, отсутствует тип? complex.mqh 79 1

'im' - неожиданная лексема, возможно, тип отсутствует? complex.mqh 79 31

'im' - ожидается тип класса complex.mqh 79 31

'}' - ожидается точка с запятой complex.mqh 82 3

'}' - выражения не допускаются в глобальной области complex.mqh 82 3

'complex' - неожиданная лексема, возможно, отсутствует тип? complex.mqh 86 1

'im' - неожиданная лексема, возможно, тип отсутствует? complex.mqh 86 41

'im' - переменная уже определена complex.mqh 86 41

см. предыдущее объявление переменной 'im' complex.mqh 79 31

'im' - ожидается тип класса complex.mqh 86

Библиотеки: MT4Orders Скользящая средняя по Халлу Обсуждение статьи "Набор инструментов

Ivan Shishikin 2023.04.03 15:32 #47

подскажите, а когда планируется добавить функции по спектральному анализу
SSA ?
в ALGLIB они есть, а в МТ, до сих пор нету.
https://www.alglib.net/translator/man/manual.cpp.html#example_ssa_d_basic

Stanislav Korotky 2023.04.03 20:59 #48

Ivan Shishikin #:

подскажите, а когда планируется добавить функции по спектральному анализу
SSA ?
в ALGLIB они есть, а в МТ, до сих пор нету.
https://www.alglib.net/translator/man/manual.cpp.html#example_ssa_d_basic

Может быть пока стандартную библиотеку не пополнят (а это может быть нескоро), устроят имеющиеся исходники (1, 2)?

Ivan Shishikin 2023.04.04 10:03 #49

к сожалению там нет кода для предикатов.
есть только сама модель . но этого мало.

Denis Kirichenko 2023.04.04 10:22 #50

Ivan Shishikin #:

подскажите, а когда планируется добавить функции по спектральному анализу
SSA ?
в ALGLIB они есть, а в МТ, до сих пор нету.
https://www.alglib.net/translator/man/manual.cpp.html#example_ssa_d_basic

Имхо, нет гарантии, что появится вообще... может, чем ждать, лучше напрямую использовать в С-коде?

Библиотеки: ALGLIB - библиотека численного анализа - страница 5