подскажите как связать это с решением СЛАУ методом Жордана-Гаусса?
тоисть заставить искать базисные и опорные решения






Это сообщение отредактировал(а) Веталька - 24.1.2010, 21:56

У меня та же тема...
Хотелось бы знать, как реализовать на с++ нахождение базисных решений СЛАУ и как вообще дать программе знать, когда какое решение (общее, базисное, опорное) искать...

Ответьте на вопрос: что такое "общее решение" для вырожденной непротиворечивой СЛАУ?
Просто я даже не знаю, что Вам писать. Просто если вы не понимаете математику, речь о С++ водить бесполезно.

ну, если слау вырожденая, то нужно проверить систему на совместность. если система совместна, то может быть тривиальное решение или базисное. Вроде так... Если нет, для начала буду рад ссылке на толковое разъяснение даной темы для бестолочи (меня)...

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

Нет, решение СЛАУ одно. Если система вырождена и совместна, то решение есть вектор X некоего подпространства. Небоходимо найти 
1) размерность подпространства решений 
2) выразить решение (вектор X) в виде л.к.:
X = X0 + a1X1 + a2X2 + .. + anXn, где X1 - Xn - базис подпространства, X0 - постоянный вектор, an - произвольные коэффициентики

Деление решений на базисные, небазисные и прочую фигню - исключительно искусственно.

З.Ы. Я не вижу в приведенном коде даже проверки на совместность.

Добавлено через 1 минуту и 22 секунды
Толковой ссылки не знаю. Учебники по ЛинАлу, наверно. Пишите конкретно, чего не ясно, здесь ответят.

Для точности: дано задание на лабораторную: 1)а) решить СЛУ методом Жордана-Гаусса; б)найти все базисные решения; 2) найти все опорные (неотрицательные базисные) решения СЛУ.
Приведенный вверху код решает систему методом Жордана-Гаусса, если количество уравнений равно количеству неизвестных (или ранг основной матрицы равен количеству неизвестных-так наверное точнее сказано).
Если количество уравнений меньше количества неизвестных, то именно тогда, я так понимаю, нужно искать базисные решения.
По поводу проверки на совместность: нужно сравнить ранг основной матрицы и расширенной. Но как вписать поиск ранга в имеющийся код?

Что значит "как вписать"? У тебя есть матрицы-массивы. Есть гугл и запрос "нахождение ранга матрицы". Что конкретно не понятно?

Ну вот, ранг теперь вычисляется. 

Далее: если ранг равен количеству переменных, то можна выводить результат решения методом ж-г--последний столбец;
если ранг меньше количества переменных, то нужно искать общее решение (в виде выражения базисных неизвестных через свободные) и базисное решение. Это я взял с небольших указаний к лабораторной... 
Вы пишете, что:
 

это будет С(n^m), где n-количество неизвестных, а m-количество базисных неизвестных.
базисные неизвестные-это те, которые после прохода массива методом ж-г встречаются только в одном уравнении (тоесть 1 при (i=j) и 0 в других ячейках столбика j).



1)Допустим, система АХ=В. постоянный вектор-это В?
2)Если что не так, прошу исправить меня.
3)Теперь самый важный на сегодня вопрос: мне теперь разбивать вектор, который есть, и потом каким-то образом выражать общее решение (прошу подсказать как), или можно обойтись без этого?

Ну вот смотри:
x+4y+8z = 5

Соответственно ранг = 1
Кол-во переменных = 3

Базис подпространства решений = 3-1 = 2 базисные функции
x = 5 - 4y -8z

Соответственно решение:
(x,y,z) = (5-4y-8z, y, z) =
= (5, 0, 0) + (-4, 1, 0)y + (-8, 0, 1)z

Это чтобы был наглядный пример, на который ориентироваться при работе с общим видом.

спасибо), пробую дальше...

Итак, программа решает системы с рангом равным количеству неизвестных, а также находит общее решение для системы с рангом меньшим за количество неизвестных. Остается найти базисные... Также отреставрировать сортировку...

Вопросы:
1) может кто знает, почему стопорится сортировка?
2) как исходя из общих, найти базисные решения?

базисные ты уже нашел. Базисные - это общие при значениях свободных переменных = 0, вроде бы.

если взять такую матрицу:
 3х5
    1 -1 3 1 4
    2 3 -1 -1 6
    1 -6 10 4 6

то базисных решений здесь будет 6.
Находятся они путем замены поочередно неизвестных нулями. При этом количество заменяемых на ноль (вроде бы) равно рангу матрицы. сейчас мучаю это решение на с++.

Добавлено через 4 минуты и 19 секунд
для предыдущей матрицы базисные решения: 
Х1=(9, 2, 0, 0)
Х2=(0, 11, 9, 0)
Х3=(11, 0, -2, 0)
Х4=(10, 0, 0, 5)
Х5=(0, 20, 0, -45)
Х6=(0, 0, 20, 55)

Добавлено через 6 минут и 7 секунд
А как на счет сортировки? Почему вектор ругается?

Подожди.
(С) Базисным решением системой уравнений называется такое решение при котором свободные переменные равны нулю.

Возможно... Но у меня в методичке в указаниях лабы решен примерчик (тот, который выше, с шестью базисными). Судя по этому решению, нужно заменять все пары по очереди (тоесть, комбинация с m по n). "От руки" такие решения я научился находить (так же как и в методичке). Задача стоит научить программу это делать.