Форум программистов, компьютерный форум, киберфорум
Наши страницы
С++ для начинающих
Войти
Регистрация
Восстановить пароль
 
 
Рейтинг: Рейтинг темы: голосов - 70, средняя оценка - 4.66
Stas0n
3 / 4 / 1
Регистрация: 13.07.2011
Сообщений: 313
#1

Решить систему алгебраических линейных неоднородных уравнени - C++

22.07.2011, 15:58. Просмотров 9371. Ответов 152
Метки нет (Все метки)

У меня есть система линейных уравнений. В ней 4000 уравнений.
Киньте плиз код для её решения. Желательно, чтобы он был максимально быстрым.
0
Надоела реклама? Зарегистрируйтесь и она исчезнет полностью.
Similar
Эксперт
41792 / 34177 / 6122
Регистрация: 12.04.2006
Сообщений: 57,940
22.07.2011, 15:58
Я подобрал для вас темы с готовыми решениями и ответами на вопрос Решить систему алгебраических линейных неоднородных уравнени (C++):

Решить систему линейных алгебраических уравнений
Решить CЛАУ 2A^tX=Y^2,где |2 1 5 6| |3| A= |5 2 2 6|...

Методом Гаусса решить систему n линейных алгебраических уравнений
Помогите пожалуйста написать на с++ методом Гаусса решить систему n линейных...

Решить систему алгебраических уравнений методом Гаусса
Решить систему алгебраических уравнений методом Гаусса: 6X1-X2-X3=11,33,...

Решить систему линейных неравенств
Нужно написать программу, решающую систему неравенств. Программа должна...

Решить систему линейных уравнений вида Ax=b методом Зейделя
2)Решить систему линейных уравнений вида Ax=b методом Зейделя 6.9 0,0319...

Система линейных алгебраических уравнений
Помогите пожалуйта с програмой вопрос жизни и смерти!!! Мне нужно написать...

152
-=ЮрА=-
Заблокирован
Автор FAQ
08.08.2011, 16:55 #121
Вот набросал проект проверки, есть нули, но есть и довольно приличные числа, пока не знаю что явилось причиной не очень точных вычислений, ещё думаю...
1
Вложения
Тип файла: rar check_results.rar (38.3 Кб, 9 просмотров)
Stas0n
3 / 4 / 1
Регистрация: 13.07.2011
Сообщений: 313
08.08.2011, 16:56  [ТС] #122
Что касается быстроты, то меньше 5 секунд все это дело заняло...
Осталось разобраться кто прав))

Добавлено через 18 секунд
Хотя пивас с меня будет))))
0
-=ЮрА=-
Заблокирован
Автор FAQ
08.08.2011, 16:59 #123
На свежую голову вечером подумаю, главное что алгоритм вышел довольно быстрым и уже есть норм результаты, я подозреваю что огрех может быть сдесь
C++
1
2
3
4
5
6
void PryamoiHod(int n, double **a, double *b)
{
....
if(v != 0)
                        for(j = k + 1; j < n; j++)
                                a[i][j] -= v*a[k][j];
Ветьвь
C++
1
if(v != 0)
ввёл для ускорения алгоритма, как по мне логика от этого не страдает а там кто его знает
1
Stas0n
3 / 4 / 1
Регистрация: 13.07.2011
Сообщений: 313
08.08.2011, 17:09  [ТС] #124
Пожалуй надо будет разобраться в твоем алгоритме... пока что не сильно понятно что там и как...
удивительно как быстро работает
на scilab весь процесс (считывание матрицы, столбца, решение, вывод в файл) занял около 200 секунд, а тут 1-2 раза маргнуть успел...

Добавлено через 7 минут
Слушай, а ты не смотрел в inpA присутствуют все диагональные элементы? (даже не важно нулевые или нет)
0
-=ЮрА=-
Заблокирован
Автор FAQ
08.08.2011, 17:14 #125
Цитата Сообщение от hello19 Посмотреть сообщение
Пожалуй надо будет разобраться в твоем алгоритме... пока что не сильно понятно что там и как...
Позже коменты напишу в срр-шнике, быстрота потому что нужно было написать быстрый алгоритм вот и пришлось изголятьсяВобще любая CAD-система работает не столь прямолинейно как наше приложение, юзаются продвинутые алгоритмы решения - они хоть и дают решение в любом случае, но требуют гораздо дольше времени на вычисления...Пришли ответы из scilab, мне охота погрешности сравнить, ведь я сравнивал Е = А*Х - В по сути нашёл абсолютную погрешность, а чтобы относительную получить нужно Е[i] / B[i] как то об этом, только сейчас подумал

Добавлено через 2 минуты
Цитата Сообщение от hello19 Посмотреть сообщение
Слушай, а ты не смотрел в inpA присутствуют все диагональные элементы? (даже не важно нулевые или нет)
Это твоя исходная матрица, как говориться что дал с тем и я решалНолей на диагоналях по моим подсчётам там всего 3 А то что не все первые элементы имеют возрастающий индекс, не говорит о том что этого индекса нет в какой-нибудь строке, ты понял о чём я?
1
Stas0n
3 / 4 / 1
Регистрация: 13.07.2011
Сообщений: 313
08.08.2011, 17:18  [ТС] #126
Про индекс я понял... просто интересно все ли диагональные элементы затронуты в каждой строчке..
Хотя это вообщем то не нужно.
0
grizlik78
Эксперт С++
1982 / 1475 / 191
Регистрация: 29.05.2011
Сообщений: 3,048
08.08.2011, 17:29 #127
Цитата Сообщение от hello19 Посмотреть сообщение
на scilab весь процесс (считывание матрицы, столбца, решение, вывод в файл) занял около 200 секунд, а тут 1-2 раза маргнуть успел...
hello19, мне, конечно, всё равно по-большому счёту, но на всякий случай ещё раз замечу, что scilab решает матрицу из сообщения #30 меньше секунды. Заметно меньше.
0
Stas0n
3 / 4 / 1
Регистрация: 13.07.2011
Сообщений: 313
08.08.2011, 17:33  [ТС] #128
Цитата Сообщение от grizlik78 Посмотреть сообщение
hello19, мне, конечно, всё равно по-большому счёту, но на всякий случай ещё раз замечу, что scilab решает матрицу из сообщения #30 меньше секунды. Заметно меньше.
Молодец, что заметил, но у меня матрица коэффициентов то порядка 3638.
0
grizlik78
Эксперт С++
1982 / 1475 / 191
Регистрация: 29.05.2011
Сообщений: 3,048
08.08.2011, 17:34 #129
hello19, а сообщение я для кого указываю? Какой там порядок матрицы?
0
-=ЮрА=-
Заблокирован
Автор FAQ
08.08.2011, 17:34 #130
Цитата Сообщение от grizlik78 Посмотреть сообщение
hello19, мне, конечно, всё равно по-большому счёту, но на всякий случай ещё раз замечу, что scilab решает матрицу из сообщения #30 меньше секунды. Заметно меньше.
-убери извещения аля printf в срр-нике и получешь время алгоритма не дольше, вчитайся в топики скилаб дал тоже не верный результат...
0
Stas0n
3 / 4 / 1
Регистрация: 13.07.2011
Сообщений: 313
08.08.2011, 17:35  [ТС] #131
to Юра: А вот то, что scilab нарешал. Я кинул на вход матрицу порядка 3638.
Могу в принципе кинуть еще и без нулевых элементов... там не сильно результаты отличаются
0
Вложения
Тип файла: rar 33.rar (15.4 Кб, 7 просмотров)
-=ЮрА=-
Заблокирован
Автор FAQ
08.08.2011, 17:37 #132
Цитата Сообщение от hello19 Посмотреть сообщение
Молодец, что заметил, но у меня матрица коэффициентов то порядка 3638.
- чувак (grizlik78), просто не представляет как геморно работать с таким объёмом данных, уж мы то знаем

Добавлено через 1 минуту
Цитата Сообщение от hello19 Посмотреть сообщение
Я кинул на вход матрицу порядка 3638.
- ок беру в работу
1
grizlik78
Эксперт С++
1982 / 1475 / 191
Регистрация: 29.05.2011
Сообщений: 3,048
08.08.2011, 17:40 #133
Цитата Сообщение от -=ЮрА=- Посмотреть сообщение
убери извещения аля printf в срр-нике и получешь время алгоритма не дольше
Мне не интересна эта программа. Я не пытался сравнивать её с scilab'ом. Успокойся.
Цитата Сообщение от -=ЮрА=- Посмотреть сообщение
скилаб дал тоже не верный результат
Если он решал эту систему 200 секунд, то я, конечно, во что угодно поверю.
Но если правильно использовать инструмент, то эта система scilab'ом решается верно и практически мгновенно.

Добавлено через 1 минуту
Цитата Сообщение от -=ЮрА=- Посмотреть сообщение
- чувак (grizlik78), просто не представляет как геморно работать с таким объёмом данных

Цитата Сообщение от -=ЮрА=- Посмотреть сообщение
уж мы то знаем
Да да. В начале ветки обсуждения я видел.
0
-=ЮрА=-
Заблокирован
Автор FAQ
08.08.2011, 17:42 #134
Х[1] = -0.0000000004
Х[425] = 0.0045016537

по первому выражению (1;1)(425;1) = 0
они должны быть равны.

В моём решении
Х[1] = -46.500
Х[425] = 46.500

Я говорю - матрица А с нулями при решении даст не верные результаты
Чуть позже детально всё сверю
1
Stas0n
3 / 4 / 1
Регистрация: 13.07.2011
Сообщений: 313
08.08.2011, 17:47  [ТС] #135
Цитата Сообщение от -=ЮрА=- Посмотреть сообщение
Х[1] = -0.0000000004
Х[425] = 0.0045016537

по первому выражению (1;1)(425;1) = 0
они должны быть равны.

В моём решении
Х[1] = -46.500
Х[425] = 46.500

Я говорю - матрица А с нулями при решении даст не верные результаты
Чуть позже детально всё сверю
Я согласен что равны. А ты подсчитай руками))) такой косяк)) у меня со scilab тот же косяк...
Там должны не 46 а что то около 200 стоять..

to grizlik78: Слушай, а как ты решаешь систему?
вот я например вот так x = linsolve (a,-b);
0
grizlik78
Эксперт С++
1982 / 1475 / 191
Регистрация: 29.05.2011
Сообщений: 3,048
08.08.2011, 17:51 #136
Цитата Сообщение от hello19 Посмотреть сообщение
to grizlik78: Слушай, а как ты решаешь систему?
Хорошо, сейчас в соответствующем разделе покажу. А пока вот вам правильное решение, полученное в scilab. Когда программа будет давать что-то похожее, можно будет на неё и посмотреть
1
Вложения
Тип файла: zip x.zip (11.8 Кб, 15 просмотров)
Stas0n
3 / 4 / 1
Регистрация: 13.07.2011
Сообщений: 313
08.08.2011, 17:54  [ТС] #137
гризлик, ты результат то видел в каком виде выдаешь)))
забавная змейка))
покажи свой код, как ты решение само делаешь?
0
grizlik78
Эксперт С++
1982 / 1475 / 191
Регистрация: 29.05.2011
Сообщений: 3,048
08.08.2011, 17:57 #138
Цитата Сообщение от hello19 Посмотреть сообщение
гризлик, ты результат то видел в каком виде выдаешь)))
Открой в нормальном редакторе. Это из-за того, что в линуксе переносы строки кодируются не так, как в Windows. Я, конечно, могу перекодировать, но надеюсь ты и сам осилишь.
Цитата Сообщение от hello19 Посмотреть сообщение
покажи свой код, как ты решение само делаешь?
Чуточку терпения, говорю же, в соответствующем разделе сейчас выложу.
1
Stas0n
3 / 4 / 1
Регистрация: 13.07.2011
Сообщений: 313
08.08.2011, 18:13  [ТС] #139
просто скажи чем ты пользовался непосредственно для решения системы? linsolve?

Добавлено через 13 минут
Гризлик, твое решение очень похоже на правильное)
0
grizlik78
Эксперт С++
1982 / 1475 / 191
Регистрация: 29.05.2011
Сообщений: 3,048
08.08.2011, 18:16 #140
Цитата Сообщение от hello19 Посмотреть сообщение
просто скажи чем ты пользовался непосредственно для решения системы? linsolve?
Нет, конечно же. Я сразу сказал, что нужно использовать разреженные матрицы. исползовался lusolve, посмотреть можно здесь.

Добавлено через 1 минуту
Цитата Сообщение от hello19 Посмотреть сообщение
Гризлик, твое решение очень похоже на правильное)
Ну по-крайней мере умножение A на x даёт b с очень небольшими погрешностями. Я смотрел.
1
08.08.2011, 18:16
MoreAnswers
Эксперт
37091 / 29110 / 5898
Регистрация: 17.06.2006
Сообщений: 43,301
08.08.2011, 18:16
Привет! Вот еще темы с решениями:

Система линейных алгебраических уравнений
Решить систему линейных алгебраических уравнений. Метод решения выбрать...

Система линейных алгебраических уравнений
помогите исправить ошибки! #include &quot;stdafx.h&quot; #include &lt;iostream&gt; #include...

Решение системы линейных алгебраических уравнений
Добрый вечер!!! Помогите, пожалуйста, написать программу для решения системы...

Решение СЛАУ - Системы линейных алгебраических уравнений
Задачу почти закончил писать, матрицу приводит к ступенчатому виду, только не...


Искать еще темы с ответами

Или воспользуйтесь поиском по форуму:
140
Ответ Создать тему
Опции темы

КиберФорум - форум программистов, компьютерный форум, программирование
Powered by vBulletin® Version 3.8.9
Copyright ©2000 - 2018, vBulletin Solutions, Inc.
Рейтинг@Mail.ru