Форум программистов [Powered by Invision Power Board]

Версия для печати темы
Нажмите сюда для просмотра этой темы в оригинальном формате

Форум программистов > Java: Общие вопросы > Обработка массива дробных чисел

Автор: JustAStudent 24.10.2006, 02:52

Здравствуйте.

У меня возникла проблема при обработке статистической зависимости методом Лагранжа. Обработка выполняется следующим методом:

Код


    public double yt (double xt, double[][] array) {
    // определение yt методом Лагранжа
        double Y = 1.0;
        for (int j = 0; j<array.length; j++) {
            double Chisl = 1.0, Znam = 1.0;
            for (int i=0; i < array.length; i++) {
                if (i != j) { 
                    Chisl *= (double) xt - array[i][0]; 
                    Znam *= (double) array[j][0] - array[i][0];
                } 
            } 
            Y += (Chisl/Znam) * array[j][1];
        }
        return Y;
    }

Проблема в непонятно откуда берущейся погрешности вычислений при обработке маленьких (дробных) чисел. Входные данные:
x 3.9 4.2 4.9 5.1 5.5 5.9 6.3 6.9 7.4
y -0.6 -0.3 0.3 0.7 1.0 1.7 2.1 2.5 3.0
Метод должен вернуть значение yt при заданном xt на промежутке [3.9; 7.4], если заданное xt равно x[i] метод должен вычислить и вернуть yt равное y[i]. Но этого не происходит smile

Вместо нужных y-ков он возвращает другие значения:
xt 3.9 4.2 4.9 5.1 5.5 5.9 6.3 6.9 7.4
yt 0.4 0.7 1.3 1.7 2.0 2.7 3.1 3.5 4.0

В программировании я новичек, не знаком со всеми хитростями, сам разобратся не смог smile

Надеюсь на помощь smile

Автор: Metal_Heart 24.10.2006, 09:04

Цитата(JustAStudent @ 24.10.2006, 02:52

)

если заданное xt равно x[i] метод должен вычислить и вернуть yt равное y[i].

1. я не понял, где это у тебя проверяется, либо не понял условия
2. предлагаю, для отладки не использовать массив, а взять только первые его элементы, т.е.:

Цитата

x 3.9
y -0.6

Автор: Kangaroo 24.10.2006, 10:25

может ты сам алгоритм не правильно запрограммил?

Распиши его по-подробней, а то непонятно... может тогда сможем помочь

Автор: Shurr 24.10.2006, 10:42

Простой анализ результатов показывает, что переменная Y должна инициализироваться нулем а не единицей.

P.S. А ведь я уже почти вспомнил полином Лагранджа...

Автор: JustAStudent 24.10.2006, 15:49

Shurr, прав 100% smile

Код


double Y = 0.0;

До меня самого это дошло, но только утром, на паре в универе smile

Всем спасибо smile