15.2.3. Регрессия специального вида

Кроме рассмотренных, в Mathcad встроено еще несколько видов трехпара-метрической регрессии. Их реализация несколько отличается от приведенных выше вариантов регрессии тем, что для них, помимо массива данных, требуется задать некоторые начальные значения коэффициентов а,ь,с. Используйте соответствующий вид регрессии, если хорошо представляете себе, какой зависимостью описывается Ваш массив данных. Когда тип регрессии плохо отражает последовательность данных, то ее результат часто бывает неудовлетворительным и даже сильно различающимся в зависимости от выбора начальных значений. Каждая из функций выдает вектор уточненных параметров а,b,с.

  • expfit(x,y,g) —регрессия экспонентой f(x)=aebx+c;
  • lgsfit(x,y,g) —регрессия логистической функцией f (x)=a/ (1+bесх);
  • sinf it (x,y,g) —регрессия синусоидой f(x) =asin (х+b)+с;
  • pwfit(x,y,g) — регрессия степенной функцией f(x)=axb+c;
  • iogfit(x,y,g) — рефессия логарифмической функцией f(x) =aln(х+b)+с;
  • lnfit(x,y) — регрессия двухпараметрической логарифмической функцией f(x)=aln(x)+b;
    • х — вектор действительных данных аргумента;
    • у — вектор действительных значений того же размера;
    • g — вектор из трех элементов, задающий начальные значения а,b,с.

Правильность выбора начальных значений можно оценить по результату регрессии — если функция, выданная Mat head, хорошо приближает зависимость у (х), значит они были подобраны удачно.

Пример расчета одного из видов трехпараметрической регрессии (экспоненциальной) приведен в листинге 15.13 и на рис. 15.17. В предпоследней строке листинга выведены в виде вектора вычисленные коэффициенты а,ь,с, а в последней строке через эти коэффициенты определена искомая функция f (х).

Листинг 15.13. Экспоненциальная регрессия

Многие задачи регрессии данных различными двухпараметрическими зависимостями у (х) можно свести к более надежной, с вычислительной точки зрения, линейной регрессии. Делается это с помощью соответствующей замены переменных.

Рис. 15.17. Экспоненциальная регрессия (листинг 15.13)