![]() |
| ||
|
|
Новости |
Поиск |
FAQ |
Статьи |
Книги |
Новости VCL
| | Форумы ЧАТ |
Орешник |
|
| |||
| |||
« предыдущая ветвь | форум | следующая ветвь »
Большие числа с плавающей точкой [D6, WinXP]
POMAH © (2002-02-27 13:04:00)Требуется библиотека, содержащая реализацию и процедуры работы с сабжем Задача - решение систем уравнений больших порядков -> стандартных типов уже не хватает. Рассмотрю любые предложения :-)
Заранее спасиба.
pasha_golub (2002-02-27 14:22:00) [1]Домножь числа так, чтобы не было дробных и используй longint, после этого переведи обратно.
Виктор Щербаков (2002-02-27 14:20:00) [2]На самом деле метод Гаусса тоже неидеален. У полученного с его помощью решения может быть плохая точность если матрица плохо обусловлена. Т.е. незначительные изменения элементов матрицы (неизбежные в машинных вачислениях) приводят к неприемлемым отклонениям решения. Лучше использовать метод регуляризации Тихонова или SVD.
О проблеме точности возникающей при решении систем линейных уравнений и методах её решения подробно написано в
Форсайт Дж., Малькольм М., Моулер Е. Машинные методы математических вычислений. М.: Мир, 1980.
POMAH (2002-02-27 14:12:00) [3]Да, я как раз об этом думаю. Но всётаки вопрс по прежнему открыт. В данный момент нужен также интерфейс маткада(хоть через него сделать чтоли).
Виктор Щербаков (2002-02-27 14:05:00) [4]Для решения систем линейных уравнений не обязательно вычислять определитель матрицы системы. Попробуй порешать методом Гаусса с выбором ведущего элемента.
POMAH (2002-02-27 13:51:00) [5]:-) Да. Причём всё работает нормально, пока размер мантиссы числа <= ABS(порядка экспонеты). Проверято на маткаде и екселе. Маткад считает правильно, а у моей процедуры и у екселя результат одинаковый(неверный). Юзал extended - с даблом ещё хуже :-(
Виктор Щербаков (2002-02-27 13:40:00) [6]Хорошо, а как ты обратную матрицу находишь?
Вычисляешь ли для этого определитель?
POMAH (2002-02-27 13:28:00) [7]Нужна большая точность(порядка 10 знаков после запятой), которая пропадает, например, в процессе нахождения обратной матрицы
Alx2 (2002-02-27 13:22:00) [8]Так, наверное, нормировать систему можно, чтобы числа большими не были?
POMAH (2002-02-27 13:19:00) [9]Матричный.
Виктор Щербаков (2002-02-27 13:12:00) [10]А какие методы используешь для решения систем?