Модераторы: LSD, AntonSaburov
  

Поиск:

Ответ в темуСоздание новой темы Создание опроса
> Обработка массива дробных чисел, Недопустимая погрешность 
V
    Опции темы
JustAStudent
Дата 24.10.2006, 02:52 (ссылка) | (нет голосов) Загрузка ... Загрузка ... Быстрая цитата Цитата


Новичок



Профиль
Группа: Участник
Сообщений: 11
Регистрация: 24.10.2006

Репутация: нет
Всего: 1



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

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

    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 
PM MAIL   Вверх
Metal_Heart
Дата 24.10.2006, 09:04 (ссылка) | (нет голосов) Загрузка ... Загрузка ... Быстрая цитата Цитата


а почему бы и нет?
**


Профиль
Группа: Участник
Сообщений: 728
Регистрация: 31.3.2005
Где: Москва

Репутация: 4
Всего: 12



Цитата(JustAStudent @  24.10.2006,  02:52 Найти цитируемый пост)
если заданное xt равно x[i] метод должен вычислить и вернуть yt равное y[i].


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

Цитата

x 3.9
y -0.6




--------------------
 не стыдно учиться, а стыдно не учиться 
PM ICQ   Вверх
Kangaroo
Дата 24.10.2006, 10:25 (ссылка) | (нет голосов) Загрузка ... Загрузка ... Быстрая цитата Цитата


AA - Aussie Animal
****


Профиль
Группа: Участник Клуба
Сообщений: 2042
Регистрация: 7.10.2006
Где: US

Репутация: 21
Всего: 104



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

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


--------------------
Lost....
PM MAIL MSN   Вверх
Shurr
Дата 24.10.2006, 10:42 (ссылка) | (нет голосов) Загрузка ... Загрузка ... Быстрая цитата Цитата


Шустрый
*


Профиль
Группа: Участник
Сообщений: 112
Регистрация: 24.10.2006

Репутация: 5
Всего: 7



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

P.S. А ведь я уже почти вспомнил полином Лагранджа...
PM MAIL   Вверх
JustAStudent
Дата 24.10.2006, 15:49 (ссылка) | (нет голосов) Загрузка ... Загрузка ... Быстрая цитата Цитата


Новичок



Профиль
Группа: Участник
Сообщений: 11
Регистрация: 24.10.2006

Репутация: нет
Всего: 1



Shurr, прав 100% smile 
Код

double Y = 0.0;

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

Это сообщение отредактировал(а) JustAStudent - 24.10.2006, 16:16
PM MAIL   Вверх
  
Ответ в темуСоздание новой темы Создание опроса
Правила форума "Java"
LSD   AntonSaburov
powerOn   tux
javastic
  • Прежде, чем задать вопрос, прочтите это!
  • Книги по Java собираются здесь.
  • Документация и ресурсы по Java находятся здесь.
  • Используйте теги [code=java][/code] для подсветки кода. Используйтe чекбокс "транслит", если у Вас нет русских шрифтов.
  • Помечайте свой вопрос как решённый, если на него получен ответ. Ссылка "Пометить как решённый" находится над первым постом.
  • Действия модераторов можно обсудить здесь.
  • FAQ раздела лежит здесь.

Если Вам помогли, и атмосфера форума Вам понравилась, то заходите к нам чаще! С уважением, LSD, AntonSaburov, powerOn, tux, javastic.

 
0 Пользователей читают эту тему (0 Гостей и 0 Скрытых Пользователей)
0 Пользователей:
« Предыдущая тема | Java: Общие вопросы | Следующая тема »


 




[ Время генерации скрипта: 0.0439 ]   [ Использовано запросов: 22 ]   [ GZIP включён ]


Реклама на сайте     Информационное спонсорство

 
По вопросам размещения рекламы пишите на vladimir(sobaka)vingrad.ru
Отказ от ответственности     Powered by Invision Power Board(R) 1.3 © 2003  IPS, Inc.