Форум программистов, компьютерный форум, киберфорум
Наши страницы

С++ для начинающих

Войти
Регистрация
Восстановить пароль
 
 
Рейтинг: Рейтинг темы: голосов - 15, средняя оценка - 4.80
Pein95
3 / 3 / 0
Регистрация: 02.01.2013
Сообщений: 116
#1

Свести матрицу к треугольному виду - C++

02.01.2013, 20:44. Просмотров 2077. Ответов 22
Метки нет (Все метки)

C++
1
2
3
4
5
for (int i = 1; i < N; ++i)
        for (int k = i; k < N; ++k)
            for (int j = N-1; j >= 0; --j)
                if (mas[k][i-1] !=0 &&  mas[i-1][i-1] !=0  && mas[i-1][j] != 0)
                mas[k][j] -= ceil((mas[k][i-1] / mas[i-1][i-1] * mas[i-1][j]));
Правильный ли ето алгоритм?
0
Надоела реклама? Зарегистрируйтесь и она исчезнет полностью.
Similar
Эксперт
41792 / 34177 / 6122
Регистрация: 12.04.2006
Сообщений: 57,940
02.01.2013, 20:44
Здравствуйте! Я подобрал для вас темы с ответами на вопрос Свести матрицу к треугольному виду (C++):

Привести матрицу к треугольному виду - C++
Коэффициенты системы линейных уравнений заданы в виде прямоугольной матрицы. С помощью допустимых преобразований привести систему к...

Задача: привести матрицу к треугольному виду - C++
Блин, я никак не разберусь, может тут кто знает? Задание: Дана квадратная матрица размера m*n. С помощью допустимых преобразований...

Как привести матрицу к треугольному виду - C++
Нужно привести матрицу к треугольному виду, можно написать прогу , а то какие нашел не работают ни черта

Как привести матрицу к треугольному виду по методу гауса? - C++
Может у кого есть прога? Нужна срочно!

Привести матрицу к треугольному виду Гауссом (не могу найти работающий код) - C++
Нужен код для приведения матрицы к треугольному виду (матрица статическая, целочисленная). Пробовал гуглить, в том числе и тут....

Как привести матрицу к верхнему треугольному виду? С++ (Имеется код нахождения нижней треуг. матрицы) - C++
Нахождение нижней треуг. матрицы: for (i=0, k=nn-1; i&lt;nn&amp;&amp;k&gt;=0; i++, k--) for (j=0; j&lt;i;...

22
-=ЮрА=-
Заблокирован
Автор FAQ
02.01.2013, 21:19 #2
Pein95, посмотри сюда - там готовый код
Двумерные массивы (СЛАУ). Привести систему к треугольному виду
0
Pein95
3 / 3 / 0
Регистрация: 02.01.2013
Сообщений: 116
03.01.2013, 20:41  [ТС] #3
Как то там все сложно очень)) Мне нужно только свести матрицу к треугольному виду)

Добавлено через 23 часа 14 минут
C++
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
void triangle_view(double a[][N],int n)
{
    double v;
    for (int k=0,i,j,im; k<n-1; k++)
    {
 
        
        im = k;
        for (i=k+1; i<n;i++)
        {
            if (fabs(a[im][k]) < fabs(a[i][k]))
                im = i;
        }
        
        if (im != k)
        {
            for (j=0; j<n; j++)
            {
                v = a[im][j]; 
                a[im][j] = a[k][j];
                a[k][j] = v;
            }
 
        }
 
    
        for (i=k+1; i<n; i++)
        {
            
            v = 1.0*a[i][k]/a[k][k];
            a[i][k] = 0;
            if (v!=0)
                for(j=k+1; j<n; j++)
                    a[i][j] = a[i][j] - v*a[k][j];
 
        }
    }
}
вот я немножко переделал алгорит из той статьи. Но он не работает коректно(

мне нужно свести матирцу к треугольному виду для того, чтоб потом найти её ранг.
0
-=ЮрА=-
Заблокирован
Автор FAQ
03.01.2013, 20:49 #4
Цитата Сообщение от Pein95 Посмотреть сообщение
Как то там все сложно очень)) Мне нужно только свести матрицу к треугольному виду)
- Там метод Гаусса с перестановкой строк(о которой все почему то забывают), ничего там сложного нет (абсолютно), идёт прямой ход (преобразование матрицы к треугольному виду) и после обратный ход, который позволяет найти корни. Переделывать мой код прямого хода не надо было - он отточен почти до безупречности,
бери в работу
Ещё раз код преобразования к треугольнику
Кликните здесь для просмотра всего текста
C++
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
void PryamoiHod(int n, double **a, double *b)
{
        double v;
        for(int k = 0,i,j,im; k < n - 1; k++)
        {
                im = k;
                for(i = k + 1; i < n; i++)
                {
                        if(fabs(a[im][k]) < fabs(a[i][k]))
                        {
                                im = i;
                        }
                }
                if(im != k)
                {
                        for(j = 0; j < n; j++)
                        {
                                v                = a[im][j];
                                a[im][j] = a[k][j];
                                a[k][j]  = v;
                        }
                        v     = b[im];
                        b[im] = b[k];
                        b[k]  = v;
                }
                for(i = k + 1; i < n; i++)
                {
                        v               = 1.0*a[i][k]/a[k][k];
                        a[i][k] = 0;
                        b[i]    = b[i] - v*b[k];
                        if(v != 0)
                        for(j = k + 1; j < n; j++)
                        {
                                a[i][j] = a[i][j] - v*a[k][j];
                        }
                }
        }
}

править его не надо - каждая строка имеет скрытый смысл, даже если он не лежит на поверхности...
1
Pein95
3 / 3 / 0
Регистрация: 02.01.2013
Сообщений: 116
03.01.2013, 20:52  [ТС] #5
но у меня нет вектора b... я его и забрал. и то что у меня осталось работает не правильно.
0
-=ЮрА=-
Заблокирован
Автор FAQ
03.01.2013, 20:57 #6
Цитата Сообщение от Pein95 Посмотреть сообщение
но у меня нет вектора b... я его и забрал. и то что у меня осталось работает не правильно.
- исключи из параметров функции и кода всё что относится к b а вот пересатновку строк, как ты сделал,
убирать не следует!
0
Pein95
3 / 3 / 0
Регистрация: 02.01.2013
Сообщений: 116
03.01.2013, 21:00  [ТС] #7
C++
1
2
3
4
5
6
7
8
9
10
 if (im != k)
        {
            for (j=0; j<n; j++)
            {
                v = a[im][j]; 
                a[im][j] = a[k][j];
                a[k][j] = v;
            }
 
        }
а ето разве не перестановка?
0
-=ЮрА=-
Заблокирован
Автор FAQ
03.01.2013, 21:10 #8

Не по теме:

Pein95, да, прости я писал зайдя с телефона, там вроди как был пробел на месте перестановки строк. Хорошо сейчас посомтрим что у тебя не срослось с кодом...



Добавлено через 1 минуту
Pein95, покажи весь свой код, думаю дело в том что возможно у тебя несоответствие размерности
Цитата Сообщение от Pein95 Посмотреть сообщение
void triangle_view(double a[][N],int n)
,
0
Pein95
3 / 3 / 0
Регистрация: 02.01.2013
Сообщений: 116
03.01.2013, 21:14  [ТС] #9
Допустим у меня матрица 5х5 заполнена елементами от 1..25. Ранг такой матрицы должен быть 2. Так как количество ненулевых строк и будет рангом. Но если использовать етот алгоритм у меня не будет ни одной нулевой строки.

C++
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
void triangle_view(double a[][N],int n)
{
    double v;
    for (int k=0,i,j,im; k<n-1; k++)
    {
 
        
        im = k;
        for (i=k+1; i<n;i++)
        {
            if (fabs(a[im][k]) < fabs(a[i][k]))
                im = i;
        }
        
        if (im != k)
        {
            for (j=0; j<n; j++)
            {
                v = a[im][j]; 
                a[im][j] = a[k][j];
                a[k][j] = v;
            }
 
        }
 
    
        for (i=k+1; i<n; i++)
        {
            
            v = 1.0*a[i][k]/a[k][k];
            a[i][k] = 0;
            if (v!=0)
                for(j=k+1; j<n; j++)
                    a[i][j] = a[i][j] - v*a[k][j];
 
        }
    }
}
функция

C++
1
2
3
const int N = 5;
double a[N][N];
triangle_view(a, 5);
0
-=ЮрА=-
Заблокирован
Автор FAQ
03.01.2013, 21:20 #10
Вобщем вод упрощённый код прямого хода со статическим массивом передаваемым в функцию
C++
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
#include <cmath>
#include <iostream>
using namespace std;
 
void PryamoiHodSimple(int n, double **a);
 
int main()
{
    double a[][4] = 
    { 
        {1.2,0,3.4,5.6},
        {-15,16.6,0,25},
        {0,-3,8,19},
        {3,2,125.5,16}
    };
    int i, j;
    int n = sizeof(a[0])/sizeof(a[0][0]);
    cout<<"\tINPUT:"<<endl;
    for(i = 0; i < n; i++)
    {
        for(j = 0; j < n; j++)
            cout<<a[i][j]<<"\t";
        cout<<endl;
    }
    double *ptr[4] = {0};//Массив указателей - позволит использовать синтаксис ** 
    //без привязки к раземеру матрицы
    for(i = 0; i < n; i++)
        ptr[i] = a[i];
    PryamoiHodSimple(n, ptr);
    cout<<"\tOUTPUT:"<<endl;
    for(i = 0; i < n; i++)
    {
        for(j = 0; j < n; j++)
            cout<<a[i][j]<<"\t";
        cout<<endl;
    }
    return 0;
}
 
void PryamoiHodSimple(int n, double **a)
{
    double v;
        for(int k = 0,i,j,im; k < n - 1; k++)
        {
                im = k;
                for(i = k + 1; i < n; i++)
                {
                        if(fabs(a[im][k]) < fabs(a[i][k]))
                        {
                                im = i;
                        }
                }
                if(im != k)
                {
                        for(j = 0; j < n; j++)
                        {
                                v                = a[im][j];
                                a[im][j] = a[k][j];
                                a[k][j]  = v;
                        }
                }
                for(i = k + 1; i < n; i++)
                {
                        v       = 1.0*a[i][k]/a[k][k];
                        a[i][k] = 0;
                        for(j = k + 1; j < n && v != 0; j++)
                        {
                                a[i][j] = a[i][j] - v*a[k][j];
                        }
                }
        }
}
0
Миниатюры
Свести матрицу к треугольному виду  
-=ЮрА=-
Заблокирован
Автор FAQ
03.01.2013, 21:22 #11
Цитата Сообщение от Pein95 Посмотреть сообщение
Допустим у меня матрица 5х5 заполнена елементами от 1..25. Ранг такой матрицы должен быть 2.
вот собственно почитай о ранге http://ru.wikipedia.org/wiki/Ранг_матрицы
0
Pein95
3 / 3 / 0
Регистрация: 02.01.2013
Сообщений: 116
03.01.2013, 21:27  [ТС] #12
Если несложно попробуйте свести матрицу 5х5 заполненую от 1..25 к треугольному виду. А то у меня опять чтото не то...
должно быть примерно так.
1 0 -1 -2 -3
0 1 2 3 4
0 0 0 0 0
0 0 0 0 0
0 0 0 0 0
0
-=ЮрА=-
Заблокирован
Автор FAQ
03.01.2013, 21:29 #13
Вот кстати на счёт 1-25
C++
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
#include <cmath>
#include <iostream>
using namespace std;
 
void PryamoiHodSimple(int n, double **a);
 
int main()
{
    double a[][5] = 
    { 
        {1,2,3,4,5},
        {6,7,8,9,10},
        {11,12,13,14,15},
        {16,17,18,19,20},
        {21,22,23,24,25}
    };
    int i, j;
    int n = sizeof(a[0])/sizeof(a[0][0]);
    cout<<"\tINPUT:"<<endl;
    for(i = 0; i < n; i++)
    {
        for(j = 0; j < n; j++)
            cout<<a[i][j]<<"\t";
        cout<<endl;
    }
    double *ptr[5] = {0};//Массив указателей - позволит использовать синтаксис ** 
    //без привязки к раземеру матрицы
    for(i = 0; i < n; i++)
        ptr[i] = a[i];
    PryamoiHodSimple(n, ptr);
    cout<<"\tOUTPUT:"<<endl;
    for(i = 0; i < n; i++)
    {
        for(j = 0; j < n; j++)
            cout<<a[i][j]<<"\t";
        cout<<endl;
    }
    return 0;
}
 
void PryamoiHodSimple(int n, double **a)
{
    double v;
        for(int k = 0,i,j,im; k < n - 1; k++)
        {
                im = k;
                for(i = k + 1; i < n; i++)
                {
                        if(fabs(a[im][k]) < fabs(a[i][k]))
                        {
                                im = i;
                        }
                }
                if(im != k)
                {
                        for(j = 0; j < n; j++)
                        {
                                v                = a[im][j];
                                a[im][j] = a[k][j];
                                a[k][j]  = v;
                        }
                }
                for(i = k + 1; i < n; i++)
                {
                        v       = 1.0*a[i][k]/a[k][k];
                        a[i][k] = 0;
                        for(j = k + 1; j < n && v != 0; j++)
                        {
                                a[i][j] = a[i][j] - v*a[k][j];
                        }
                }
        }
}
Всё что подчёркнуто на скриншоте розовым следует воспринимать как ноль (там не ноль ввиду особенностей машинного представления чисел с плавающей точкой)
Да теперь согласен - ранг у полученной матрицы равен 2-м
0
Миниатюры
Свести матрицу к треугольному виду  
-=ЮрА=-
03.01.2013, 21:30
  #14

Не по теме:

Цитата Сообщение от Pein95 Посмотреть сообщение
должно быть примерно так.
1 0 -1 -2 -3
0 1 2 3 4
0 0 0 0 0
0 0 0 0 0
0 0 0 0 0
- с перестановкой строк будет как нас скриншоте

0
Pein95
3 / 3 / 0
Регистрация: 02.01.2013
Сообщений: 116
03.01.2013, 21:32  [ТС] #15
ну а как заставить програму те подчеркнутые цыфры щитать как ноль?)
0
03.01.2013, 21:32
MoreAnswers
Эксперт
37091 / 29110 / 5898
Регистрация: 17.06.2006
Сообщений: 43,301
03.01.2013, 21:32
Привет! Вот еще темы с ответами:

Привести вектор к треугольному виду - C++
Я решаю систему лин. уравнений одной малоизвестной модификацией метода Гаусса. В чем суть: =&gt; есть вектор с таким содержимым: 1...

приведение матрицы к треугольному виду - C++
вводится размерность матрицы,нужно привести ее к треугольному виду. Помогите,пожалуйста. прочитала все темы из киберфорума,но только там...

Двумерные массивы (СЛАУ). Привести систему к треугольному виду - C++
Коэффициенты системы линейных уравнений заданы в виде прямоугольной матрицы.С помощью допустимых преобразований привести систему к...

С помощью допустимых изменений привести систему к треугольному виду... - C++
Надо выполнить вот это задание &quot;: Коефіцієнти системи лінійних рівнянь задані у вигляді прямокутної матриці. З допомогою допустимих...


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

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

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