2

Есть некая функция

alt text

где alt text

Возникла необходимость вычислять обратную функцию

alt text

Сейчас думаю над элегантным и не слишком медленным методом вычисления x(y). Какие у меня есть варианты?

5
  • Прошу прощения, поверил обещаниям, что поддерживается ввод разметки TeX. @Barmaley, переделал в глинянные таблички.
    – karmadro4
    19 апр 2012 в 19:03
  • Ничего не понял - формулу то напиши...
    – Barmaley
    19 апр 2012 в 19:13
  • @karmadro4, Это вопрос на форум Математика. 19 апр 2012 в 20:38
  • 2
    Блин, опять настал формулам капут. @Павел Зайцев, лол, нет. Не пугайтесь, этот матан не так страшен :-)
    – karmadro4
    19 апр 2012 в 21:14
  • @karmadro4, TeX-разметка поддерживается только на форуме "Математика".
    – angry
    19 апр 2012 в 23:07

3 ответа 3

5

Стандартное рекуррентное соотношение не подойдет?

// x(n + 1) = 0.5*(x(n) + n / x(n))

scanf("%u", &test_number);

x_new = test_number;
x_old = 0;

while (x_old != x_new) {
    x_old = x_new;
    x_new = (x_new + test_number / x_new) >> 1;
}

printf("%u\n", x_new);
1
  • О, формула Герона! Как раз сейчас экспериментирую с ней. Кстати, сходится гораздо быстрее если покормить не нулем, а грубой оценкой.
    – karmadro4
    19 апр 2012 в 20:29
4

Если правильно понял насчет натуральных чисел, то вот для целых 32-бит.

unsigned int 
sqrt32(unsigned long n)  
{  
    unsigned int c = 0x8000;  
    unsigned int g = 0x8000;

    for(;;) {  
        if(g*g > n)  
            g ^= c;  
        c >>= 1;  
        if(c == 0)  
            return g;  
        g |= c;  
    }  
}
4
  • Интересный подход, хотя я мало чего понял :-) А ограничение натуральными числами позволяет не волноваться о том, существует ли обратная функция или нет.
    – karmadro4
    19 апр 2012 в 21:30
  • Алгоритм не мой, из интернета. Насколько понимаю программная модель кусочка железяки. По крайней мере так сказал человек, который когда-то разрабатывал железо.
    – avp
    19 апр 2012 в 21:35
  • Спасибо, попробую поиграться с ним и обобщить для разных разрядностей.
    – karmadro4
    19 апр 2012 в 21:57
  • Спасибо, скопировал себе.
    – northerner
    20 апр 2012 в 3:31
0

Разложить в ряд Тейлора

Ваш ответ

By clicking “Отправить ответ”, you agree to our terms of service and acknowledge you have read our privacy policy.

Всё ещё ищете ответ? Посмотрите другие вопросы с метками или задайте свой вопрос.