Статьи

Рациональность и MS Excel (и другие калькуляторы)

Этим утром Матье получил хороший опыт в своем курсе  вычислительных методов в актуарной науке . Но давайте начнем с некоторых математических формальных определений.

Во-первых, вспомните, что  http://latex.codecogs.com/gif.latex?y^x это — как-то — стандартное выражение. Никто не должен удивляться, увидев такое выражение. Как правило (как описано в  http://en.wikipedia.org/…  ), эта функция определяется только тогда, когда  http://latex.codecogs.com/gif.latex?y\in\mathbb{R}_+. Идея заключается в том , что  определение  о  http://latex.codecogs.com/gif.latex?y^x том , что

http://latex.codecogs.com/gif.latex?y^x%20=%20\exp\left(x\log[y>\right)

И это  определение . Такая функция существует только если  http://latex.codecogs.com/gif.latex?y\in\mathbb{R}_+ (возможно, исключая  http://latex.codecogs.com/gif.latex?0). Это было бы стандартным определением в реальном анализе.

Теперь эта функция «сила» появляется и в комплексном анализе, когда речь идет о единичных корнях. От случая, если  http://latex.codecogs.com/gif.latex?z=y^{\frac{1}{k}}e^{i%20\frac{2n\pi}{k}}, где  http://latex.codecogs.com/gif.latex?y\in\mathbb{R}_+ и  http://latex.codecogs.com/gif.latex?k\in\mathbb{N}_\star, для некоторых  http://latex.codecogs.com/gif.latex?n\in\mathbb{N}, тогда  http://latex.codecogs.com/gif.latex?z^k=y. Таким образом, в комплексном анализе может быть сложнее определить правильно,  http://latex.codecogs.com/gif.latex?y^x поскольку он может быть не уникальным. Но мы можем связать (иногда, когда  http://latex.codecogs.com/gif.latex?x инверсия целого числа или, может быть, рациональное число?) С корнями полиномиальных функций. Пока что ничего нового…

Вернемся к проблеме Матье. На самом деле, в своем курсе, он хотел вычислить  http://latex.codecogs.com/gif.latex?(-8)^{\frac{1}{3}}. С французской версией Excel введите

ты получаешь  http://latex.codecogs.com/gif.latex?-2. Если вы посмотрите на окно «Помощь», у вас есть еще несколько деталей

Похоже , эта  шляпа  функция может быть использована для определения объектов , таких как  http://latex.codecogs.com/gif.latex?y^x. Но с

ты получаешь

(имеется в виду, что это проблема …). Также можно использовать   функцию power  ( puissance in French) Excel,

Здесь вы также получите

Странная часть здесь является то, что в окне «помощи», вы можете прочитать , что эта  сила  функция может быть использована с любым числом в  http://latex.codecogs.com/gif.latex?\mathbb{R}.

Еще один момент … как насчет   ? Каким-то образом это просто квадрат предыдущего (с дробью) … Здесь, набрав

ты получаешь

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

 

ты получаешь

 Но если вы печатаете (что должно быть близко, из свойства непрерывности степенной функции)

 

ты получаешь

и аналогично

В  Wolfram Mathworld введите

Mathematica признает, что мы пытаемся иметь дело с единичными корнями: результат здесь

с — как и ожидалось — числовым приближением

С Matlab Матье получил то же самое, что и Mathematica (его десятичное приближение). И в заключение, с R, Матье получил

> (-8)^(1/3)
[1] NaN
> (-8)^(.333333333333333)
[1] NaN

Поэтому для R вы не можете использовать эту функцию с отрицательными числами.

Теперь, как мы можем интерпретировать эти результаты?

1) Я понимаю, что с MS Excel,  http://latex.codecogs.com/gif.latex?x^{ab}\neq%20\left(x^a\right)^bпоскольку

http://latex.codecogs.com/gif.latex?(-8)^{\frac{2}{3}}\neq%20\left((-8)^{\frac{1}{3} } \ справа) ^ 2

что проблематично. Например, в страховании с ежемесячными скидками у нас есть такие функции, как  http://latex.codecogs.com/gif.latex?u^{\frac{k}{12}}. Что если

http://latex.codecogs.com/gif.latex?u^{\frac{k}{12}}\neq%20\left(u^{\frac{1}{12}}\right)^k

2) Проблема возникает — вероятно (MS Excel не является открытым программным обеспечением, поэтому его может быть трудно проверить) — из-за того факта, что  http://latex.codecogs.com/gif.latex?y^{\frac{1}{n}} интерпретируется как обратная (возможно) биективная функция. Чтобы быть более конкретным,  http://latex.codecogs.com/gif.latex?x=y^{\frac{1}{n}}значит это  http://latex.codecogs.com/gif.latex?x^n=y. Когда  http://latex.codecogs.com/gif.latex?n нечетное целое число, то (в реальном анализе) существует однозначное обратное, и, таким образом,  http://latex.codecogs.com/gif.latex?y^{\frac{1}{n}} оно однозначно определено, поскольку  http://latex.codecogs.com/gif.latex?x\mapsto%20x^n является биективной  http://latex.codecogs.com/gif.latex?\mathbb{R}\rightarrow\mathbb{R} функцией. Это то , что делают MS Excel (и Google):  http://latex.codecogs.com/gif.latex?x\mapsto%20x^3 биективная  http://latex.codecogs.com/gif.latex?\mathbb{R}\rightarrow\mathbb{R} функция, так  http://latex.codecogs.com/gif.latex?(-8)^{\frac{1}{3}} означает , что мы должны найти  уникальное значение (реального)   таким образом, что  . Таким образом, почему-то имеет смысл вернуться  .http://latex.codecogs.com/gif.latex?xhttp://latex.codecogs.com/gif.latex?x^3=-8http://latex.codecogs.com/gif.latex?-2

3) Все еще есть проблемы с Google и Mathematica. Это нормально, чтобы вернуть единичные корни в  http://latex.codecogs.com/gif.latex?\mathbb{C}. Но как получается, что есть только одна ценность? Я имею в виду, да  http://latex.codecogs.com/gif.latex?1+\sqrt{3}%20\%20i , это возможный ответ, так как

http://latex.codecogs.com/gif.latex?(1+\sqrt{3}%20\%20i)^3=-8

но можно также наблюдать это, и аналогично,  http://latex.codecogs.com/gif.latex?(-2)^3=-8 и

http://latex.codecogs.com/gif.latex?(1-\sqrt{3}%20\%20i)^3=-8

Можно проверить с

С R, так как здесь мы имеем дело не со степенной функцией, а с корнями, если мы хотим найти  http://latex.codecogs.com/gif.latex?x такую  http://latex.codecogs.com/gif.latex?x^3=-8, что функция

> polyroot(c(8,0,0,1))
[1]  1+1.732051i -2+0.000000i  1-1.732051i

Что отличается … Странно, не правда ли?