-=ЮрА=-

Форум программистов и сисадминов Киберфорум

Войти

Регистрация
Восстановить пароль

Блоги

Сообщество

Поиск

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

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

Показов 2259 Комментарии 5

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

Итак задача : нам следует возвести в степень два вещественных числа, а именно необходимо определить значение с
$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 -=ЮрА=- вне форума

-=ЮрА=-

Регистрация 08.08.2009
Адрес Украина, Донецк
Сообщений 10,325
Записей в блоге 24

Новые блоги и статьи Все статьи Все блоги /
Часы электронные Uhbif79 12.08.2026 Выкладываю программу часов. Программа позволяет: 1. Использовать системное время и дату, 2. Есть возможность вводить время и дату вручную. 3. Реализованы 2 будильника: начало и конец рабочего дня. . . .	Часы с будильником на основе класса QLCDNumber Uhbif79 12.08.2026 Всем добрый день, выкладываю программу часов с будильником на основе класса QLCDNumber. Здесь я пробовал самостоятельно создавал классы, впервые столкнулся с видимостью переменной одного класса из. . .	Установка MinGW GCC 16.2 и CMake 8Observer8 10.08.2026 VK Видео: https:/ / vkvideo. ru/ video-240781534_456239017 YouTube: eY5-5PyI9NM Текстовая версия	Неделя из жизни имитационной модели склада: мои кривые руки растут, откуда надо anaschu 10.08.2026 Неделя из жизни имитационной модели склада: как я почти написал неправильную логику и что с этим делать Работаю сейчас над учебно-рабочим проектом: строю в AnyLogic имитационную модель процессов. . .
Калькулятор для расчета родства russiannick 07.08.2026 1. Задача: Создать калькулятор для расчета родства. Родственных связей существует 8 ступеней, такие как: p - отец P - мать q - муж Q - жена b - брат B - сестра s - сын S - дочь	Мир по моей воле kumehtar 07.08.2026 Когда-то кажется, что всё просто. Ты весь такой светлый. Причиняешь добро. Борешься за справедливость в этом тёмном мире. Потом начинаешь замечать одну неприятную вещь. Почти каждый хороший. . .	Кредитный калькулятор Maks 05.08.2026 Решение задачи по прикладной информатике средствами 1С. Задача: Напишите приложение-калькулятор, которое помогает рассчитывать параметры кредита для аннуитетного и дифференцированного видов. . .	У нас сейчас поговорку "Опять 25" нужно переделать на "Опять +35". kumehtar 04.08.2026 С ностальгией вспоминаю времена моего детства, когда у нас и правда +25 - была максимальная температура летом. Раньше +25 °C реально казались вершиной жары, когда можно было весь день пропадать на. . .

Наверх

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

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