Быстрый алгоритм возведения в степень

Запись от -=ЮрА=- размещена 26.07.2013 в 10:07

Проблема довольно известная и возможно кому то покажется глупым её рассмотрение, НО прежде чем закрыть данную ветку раздела блогов внимательно читаем название и понимаем что будет не так как в учебниках. Заинтрегованы?Тогда читаем дальше!

Итак задача : нам следует возвести в степень два вещественных числа, а именно необходимо определить значение с
$https://www.cyberforum.ru/cgi-bin/latex.cgi?c = {a}^{b}$ и сделать это максимально быстро.
(кто побежал за учебником алгебры в раздел функций комплексного переменного и прочей экзотикой - быстро её закрыли, тут такого не будет, т.к комплексов в явном виде в расчётах не будет)

Давайте посмотрим на равенство внимательно:
Если прологорифмировать обе части равенства получим:
$https://www.cyberforum.ru/cgi-bin/latex.cgi?ln(c) = b*ln(a)$
а теперь снова проведём обратное преобразование:
$https://www.cyberforum.ru/cgi-bin/latex.cgi?c = {e}^{b*ln(a)}$
Кому-то покажется странным замена 1-й матоперации возведения в степень сразу двумя - логорифмированием и потенцированием. Что ж отмечу, что нахождение логорифма и экспоненты через разложение в ряд Тейлора потребует чуть ли не банальных 10-12 операций сложения и является весьма тривиальной задачей, самое главное - что при данном построении решения выигрыш в скорости будет невероятно большим (всем кто сомневается предлагаю проверить)

Вроде хорошо, но вот уже слышу свист и фразу : "а если показатель степени отрицательное?!Логорифма отрицательного числа не бывает!"
Что ж отвечаю $https://www.cyberforum.ru/cgi-bin/latex.cgi?|c| = {e}^{b*ln(|a|)}$
Ага снова недовольство, где наш комплекс зачем нам его модуль?

И вот тут магия - оказалось что
$https://www.cyberforum.ru/cgi-bin/latex.cgi?c = Re({e}^{b*ln(|a|)}) + j*Im({e}^{b*ln(|a|)})$
обладает удивительным свойством, а именно
$https://www.cyberforum.ru/cgi-bin/latex.cgi?atan2(Re({e}^{b*ln(|a|)}), Im({e}^{b*ln(|a|)})) = \frac{\pi }{10}$
Что из этого следует - а следует то что при отрицательном показатели степени нам всего лишь достаточно полученный нами результат $https://www.cyberforum.ru/cgi-bin/latex.cgi?|c| = {e}^{b*ln(|a|)}$ записать вот так
$https://www.cyberforum.ru/cgi-bin/latex.cgi?\hat{c} = |c|*cos(\frac{\pi }{10}) + j*|c|*sin(\frac{\pi }{10})$

Сама проблема рассмотрена в данной теме
Возведение в степень, отрицательные числа

Коды разложений не привожу а даю ссылку на свой FAQ по математике где всё есть (ссылка)

Размещено в Без категории

« Безопасный ввод через функцию-шаблон Главная страница

Надоела реклама? Зарегистрируйтесь и она исчезнет полностью.

Всего комментариев 5

Комментарии

Пи/10 используем лишь когда основание степени a < 0
т.е изначально находим
$https://www.cyberforum.ru/cgi-bin/latex.cgi?|c| = {e}^{b*ln(|a|)}$
Потом смотрим на основание - если оно < 0
то потребуется конечное преобразование с представлением комплекса в алгебраической форме через его модуль и как оказывается магическую фазу Пи/10.
На счёт с - честно скажу не понял, если выражаться простым языком, с - просто значение f(x) = exp(x)

Запись от -=ЮрА=- размещена 30.07.2013 в 23:47 -=ЮрА=- вне форума

Сегодня за долгое время просмотрел свой блог, наткнулся на некогда живую дискуссию без меня

Фаза Пи/10 получается лишь при показателях степени b = 2,1, 4,1, 6,1 и т. д. Возьми в качестве b число, равное 4,2 и "магическая" фаза сразу станет равной Пи/5.

- да прав, но смотреть надо глубже, а не пытаться закидать меня камнями...
Уточнение в алгоритме (справедливо для всех b >= 0)
$https://www.cyberforum.ru/cgi-bin/latex.cgi?atan2(Re({e}^{b*ln(|a|)}), Im({e}^{b*ln(|a|)})) = \pi*(b - floor(b))$

MathCAD файл оновлен по ссылке https://www.cyberforum.ru/cpp-... ost5923921

Не по теме:

тот, кто знает мир науки, поймёт к чему здесь этот смайл :p .

Запись от -=ЮрА=- размещена 19.03.2014 в 22:22 -=ЮрА=- вне форума

Что то не сходится: Re=Im, а угол почему-то . Ваще то синус равен косинусу для совсем другого угла, для . Да и скакого перепугу они уравнялись? Кроме того, каждая из компонент Re и Im не могут быть равны модулю только в том случае, ели вторая равна нолю, так что одновременно они не могут быть равны модулю.

Я написал чёрным по белому, в последнем посте
Уточнение в алгоритме (справедливо для всех b >= 0)
$https://www.cyberforum.ru/cgi-bin/latex.cgi?atan2(Re({e}^{b*ln(|a|)}), Im({e}^{b*ln(|a|)})) = \pi*(b - floor(b))$

Не по теме:

прежде чем писать прошу читать что написано

Запись от -=ЮрА=- размещена 20.03.2014 в 20:55 -=ЮрА=- вне форума

Будь другом, приведи доказательство, что: $https://www.cyberforum.ru/cgi-bin/latex.cgi?atan2(...) = \pi*(b-floor(b));$

Дробная часть показателя описывает фазу : на промежутке 0...1 содержится весь период для тангенса, т.е Пи радиан
Вобщем проводим аналогию (0...Пи)-> (0...1)

Запись от -=ЮрА=- размещена 20.03.2014 в 21:02 -=ЮрА=- вне форума

Зосима, повторюсь магическое не Пи/ 10 а магическое
$https://www.cyberforum.ru/cgi-bin/latex.cgi?atan2(Re({e}^{b*ln(|a|)}), Im({e}^{b*ln(|a|)})) = \pi*(b - floor(b)) = \phi$

Отмечал пару постов выше

Для тех кто только включился поясняю мы можем посчитать экспоненту без формулы Эйлера(http://ru.wikipedia.org/wiki/Формула_Эйлера), верней выйти на её окнчательный резултат путём более простых преобразований (менее ресурсозатратных для ЭВМ разложений)
$https://www.cyberforum.ru/cgi-bin/latex.cgi?|c| = {e}^{b*ln(|a|)}$ записать вот так
$https://www.cyberforum.ru/cgi-bin/latex.cgi?\hat{c} = |c|*cos(\phi) + j*|c|*sin(\phi)$
На вопрос зачем это нужно - отвечаю для численных методов - это быстрей прямого алгоритма в лоб на порядок.

Кликните здесь для просмотра всего текста

Повторяю ссылки на всякий пожарный

https://www.cyberforum.ru/cpp-... ost5923921

Запись от -=ЮрА=- размещена 20.03.2014 в 23:03 -=ЮрА=- вне форума

Новые блоги и статьи Все статьи Все блоги /
Музыка, написанная Искусственным Интеллектом volvo 04.12.2025 Всем привет. Некоторое время назад меня заинтересовало, что уже умеет ИИ в плане написания музыки для песен, и, собственно, исполнения этих самых песен. Стихов у нас много, уже вышли 4 книги, еще 3. . .	От async/await к виртуальным потокам в Python IndentationError 23.11.2025 Армин Ронахер поставил под сомнение async/ await. Создатель Flask заявляет: цветные функции - провал, виртуальные потоки - решение. Не threading-динозавры, а новое поколение лёгких потоков. Откат?. . .	Поиск "дружественных имён" СОМ портов Argus19 22.11.2025 Поиск "дружественных имён" СОМ портов На странице: https:/ / norseev. ru/ 2018/ 01/ 04/ comportlist_windows/ нашёл схожую тему. Там приведён код на С++, который показывает только имена СОМ портов, типа,. . .	Сколько Государство потратило денег на меня, обеспечивая инсулином. Programma_Boinc 20.11.2025 Сколько Государство потратило денег на меня, обеспечивая инсулином. Вот решила сделать интересный приблизительный подсчет, сколько государство потратило на меня денег на покупку инсулинов. . . .	Ломающие изменения в C#.NStar Alpha Etyuhibosecyu 20.11.2025 Уже можно не только тестировать, но и пользоваться C#. NStar - писать оконные приложения, содержащие надписи, кнопки, текстовые поля и даже изображения, например, моя игра "Три в ряд" написана на этом. . .
Мысли в слух kumehtar 18.11.2025 Кстати, совсем недавно имел разговор на тему медитаций с людьми. И обнаружил, что они вообще не понимают что такое медитация и зачем она нужна. Самые базовые вещи. Для них это - когда просто люди. . .	Создание Single Page Application на фреймах krapotkin 16.11.2025 Статья исключительно для начинающих. Подходы оригинальностью не блещут. В век Веб все очень привыкли к дизайну Single-Page-Application . Быстренько разберем подход "на фреймах". Мы делаем одну. . .	Фото: Daniel Greenwood kumehtar 13.11.2025	Расскажи мне о Мире, бродяга kumehtar 12.11.2025 — Расскажи мне о Мире, бродяга, Ты же видел моря и метели. Как сменялись короны и стяги, Как эпохи стрелою летели. - Этот мир — это крылья и горы, Снег и пламя, любовь и тревоги, И бескрайние. . .	PowerShell Snippets iNNOKENTIY21 11.11.2025 Модуль PowerShell 5. 1+ : Snippets. psm1 У меня модуль расположен в пользовательской папке модулей, по умолчанию: \Documents\WindowsPowerShell\Modules\Snippets\ А в самом низу файла-профиля. . .

Сообщение форума

Отменить изменения