С Новым годом! Форум программистов, компьютерный форум, киберфорум
С++ для начинающих
Форум программистов и сисадминов Киберфорум > Форум программистов > C/C++ > C++ > С++ для начинающих
Войти
Регистрация
Восстановить пароль
Блоги Сообщество Поиск Заказать работу  
 
0 / 0 / 1
Регистрация: 29.09.2014
Сообщений: 97

Решение системы линейных уравнений

07.01.2015, 00:21. Показов 883. Ответов 0
Метки нет (Все метки)
Студворк — интернет-сервис помощи студентам
Задание:
Решить систему линейных уравнений:
-сформировать матрицу А коэффициентов и матрицу В правых частей (матрицы должны быть размещены в файле in_file.txt);
-вычислить обратную матрицу https://www.cyberforum.ru/cgi-bin/latex.cgi?{A}^{-1};
-вычислить матрицу X - решение системы уравнений;
-напечатать решение в файл o_file.txt


x1-x2+x4=0
x2-x3+x5=0
-x1+x3+x6+0
-5*x1-25*x2-9*x3=-2
25*x2+7*x4-43*x5=17
9*x3+43*x5-31*x6=-20


И вот что у меня вышло (код ниже).
Подскажите, пожалуйста, где у меня там ошибка? Вроде все хорошо, но программа выдает не то, что нужно(если считать даже через калькулятор онлайн или ручкам, то выходит совсем другой результат, так что проблема в программе, но что не так, не пойму), если говорить о том, что она выдает, то вот:
Кликните здесь для просмотра всего текста
3.53562e+041 1.08434e-311 2.69493e-312 3.14281e-317 6.94713e-308
6.94796e-308 1.21696e-008 6.94772e-308 2.1319e-312 1.08434e-311
2.12213e-313 1.4822e-323 1.08434e-311 5.57703e-307 -1.#QNAN
6.36887e-314 5.03947e-322 2.92907e+267 6.94803e-308 7.35627e-314
1.31778e+048 1.08434e-311 4.41376e-312 1.79087e-307 3.05514e+267

1.#QNAN 1.#QNAN 1.#QNAN 0 6.94713e-308
1.#QNAN 1.#QNAN 1.#QNAN 0 1.08434e-311
1.#QNAN 1.#QNAN 1.#QNAN 1.#QNAN 1.#QNAN
6.36887e-314 5.03947e-322 2.92907e+267 6.94803e-308 7.35627e-314
1.31778e+048 1.08434e-311 4.41376e-312 1.79087e-307 3.05514e+267

1.#QNAN 1.#QNAN 1.#QNAN 0 6.94713e-308
1.#QNAN 1.#QNAN 1.#QNAN 0 1.08434e-311
1.#QNAN 1.#QNAN 1.#QNAN 1.#QNAN 1.#QNAN
6.36887e-314 5.03947e-322 2.92907e+267 6.94803e-308 7.35627e-314
1.31778e+048 1.08434e-311 4.41376e-312 1.79087e-307 3.05514e+267

1.#QNAN 1.#QNAN 1.#QNAN 0 6.94713e-308
1.#QNAN 1.#QNAN 1.#QNAN 0 1.08434e-311
1.#QNAN 1.#QNAN 1.#QNAN 1.#QNAN 1.#QNAN
6.36887e-314 5.03947e-322 2.92907e+267 6.94803e-308 7.35627e-314
1.31778e+048 1.08434e-311 4.41376e-312 1.79087e-307 3.05514e+267

1.#QNAN 1.#QNAN 1.#QNAN 0 6.94713e-308
1.#QNAN 1.#QNAN 1.#QNAN 0 1.08434e-311
1.#QNAN 1.#QNAN 1.#QNAN 1.#QNAN 1.#QNAN
6.36887e-314 5.03947e-322 2.92907e+267 6.94803e-308 7.35627e-314
1.31778e+048 1.08434e-311 4.41376e-312 1.79087e-307 3.05514e+267

1.#QNAN 1.#QNAN 1.#QNAN 0 6.94713e-308
1.#QNAN 1.#QNAN 1.#QNAN 0 1.08434e-311
1.#QNAN 1.#QNAN 1.#QNAN 1.#QNAN 1.#QNAN
6.36887e-314 5.03947e-322 2.92907e+267 6.94803e-308 7.35627e-314
1.31778e+048 1.08434e-311 4.41376e-312 1.79087e-307 3.05514e+267

-1 0 1 0 0
1.#QNAN 1.#QNAN 1.#QNAN 0 1.08434e-311
1.#QNAN 1.#QNAN 1.#QNAN 1.#QNAN 1.#QNAN
6.36887e-314 5.03947e-322 2.92907e+267 6.94803e-308 7.35627e-314
1.31778e+048 1.08434e-311 4.41376e-312 1.79087e-307 3.05514e+267

0 0 1 0 0
1.#QNAN 1.#QNAN 1.#QNAN 0 1.08434e-311
1.#QNAN 1.#QNAN 1.#QNAN 1.#QNAN 1.#QNAN
6.36887e-314 5.03947e-322 2.92907e+267 6.94803e-308 7.35627e-314
1.31778e+048 1.08434e-311 4.41376e-312 1.79087e-307 3.05514e+267

0 1 1 0 0
1.#QNAN 1.#QNAN 1.#QNAN 0 1.08434e-311
1.#QNAN 1.#QNAN 1.#QNAN 1.#QNAN 1.#QNAN
6.36887e-314 5.03947e-322 2.92907e+267 6.94803e-308 7.35627e-314
1.31778e+048 1.08434e-311 4.41376e-312 1.79087e-307 3.05514e+267

0 1 -1 0 0
1.#QNAN 1.#QNAN 1.#QNAN 0 1.08434e-311
1.#QNAN 1.#QNAN 1.#QNAN 1.#QNAN 1.#QNAN
6.36887e-314 5.03947e-322 2.92907e+267 6.94803e-308 7.35627e-314
1.31778e+048 1.08434e-311 4.41376e-312 1.79087e-307 3.05514e+267

0 1 -1 0 0
1.#QNAN 1.#QNAN 1.#QNAN 0 1.08434e-311
1.#QNAN 1.#QNAN 1.#QNAN 1.#QNAN 1.#QNAN
6.36887e-314 5.03947e-322 2.92907e+267 6.94803e-308 7.35627e-314
1.31778e+048 1.08434e-311 4.41376e-312 1.79087e-307 3.05514e+267

0 1 -1 0 1
1.#QNAN 1.#QNAN 1.#QNAN 0 1.08434e-311
1.#QNAN 1.#QNAN 1.#QNAN 1.#QNAN 1.#QNAN
6.36887e-314 5.03947e-322 2.92907e+267 6.94803e-308 7.35627e-314
1.31778e+048 1.08434e-311 4.41376e-312 1.79087e-307 3.05514e+267

-1 0 1 0 0
1 -1 0 1 0
1.#QNAN 1.#QNAN 1.#QNAN 1.#QNAN 1.#QNAN
6.36887e-314 5.03947e-322 2.92907e+267 6.94803e-308 7.35627e-314
1.31778e+048 1.08434e-311 4.41376e-312 1.79087e-307 3.05514e+267

0 0 1 0 0
-1 -1 0 1 0
1.#QNAN 1.#QNAN 1.#QNAN 1.#QNAN 1.#QNAN
6.36887e-314 5.03947e-322 2.92907e+267 6.94803e-308 7.35627e-314
1.31778e+048 1.08434e-311 4.41376e-312 1.79087e-307 3.05514e+267

0 1 1 0 0
-1 0 0 1 0
1.#QNAN 1.#QNAN 1.#QNAN 1.#QNAN 1.#QNAN
6.36887e-314 5.03947e-322 2.92907e+267 6.94803e-308 7.35627e-314
1.31778e+048 1.08434e-311 4.41376e-312 1.79087e-307 3.05514e+267

0 1 -1 0 0
-1 0 1 1 0
1.#QNAN 1.#QNAN 1.#QNAN 1.#QNAN 1.#QNAN
6.36887e-314 5.03947e-322 2.92907e+267 6.94803e-308 7.35627e-314
1.31778e+048 1.08434e-311 4.41376e-312 1.79087e-307 3.05514e+267

0 1 -1 0 0
-1 0 1 0 0
1.#QNAN 1.#QNAN 1.#QNAN 1.#QNAN 1.#QNAN
6.36887e-314 5.03947e-322 2.92907e+267 6.94803e-308 7.35627e-314
1.31778e+048 1.08434e-311 4.41376e-312 1.79087e-307 3.05514e+267

0 1 -1 0 1
-1 0 1 0 0
1.#QNAN 1.#QNAN 1.#QNAN 1.#QNAN 1.#QNAN
6.36887e-314 5.03947e-322 2.92907e+267 6.94803e-308 7.35627e-314
1.31778e+048 1.08434e-311 4.41376e-312 1.79087e-307 3.05514e+267

-1 0 1 0 0
1 -1 0 1 0
0 1 0 0 1
6.36887e-314 5.03947e-322 2.92907e+267 6.94803e-308 7.35627e-314
1.31778e+048 1.08434e-311 4.41376e-312 1.79087e-307 3.05514e+267

0 0 1 0 0
-1 -1 0 1 0
-5 1 0 0 1
6.36887e-314 5.03947e-322 2.92907e+267 6.94803e-308 7.35627e-314
1.31778e+048 1.08434e-311 4.41376e-312 1.79087e-307 3.05514e+267

0 1 1 0 0
-1 0 0 1 0
-5 -25 0 0 1
6.36887e-314 5.03947e-322 2.92907e+267 6.94803e-308 7.35627e-314
1.31778e+048 1.08434e-311 4.41376e-312 1.79087e-307 3.05514e+267

0 1 -1 0 0
-1 0 1 1 0
-5 -25 -9 0 1
6.36887e-314 5.03947e-322 2.92907e+267 6.94803e-308 7.35627e-314
1.31778e+048 1.08434e-311 4.41376e-312 1.79087e-307 3.05514e+267

0 1 -1 0 0
-1 0 1 0 0
-5 -25 -9 0 1
6.36887e-314 5.03947e-322 2.92907e+267 6.94803e-308 7.35627e-314
1.31778e+048 1.08434e-311 4.41376e-312 1.79087e-307 3.05514e+267

0 1 -1 0 1
-1 0 1 0 0
-5 -25 -9 0 0
6.36887e-314 5.03947e-322 2.92907e+267 6.94803e-308 7.35627e-314
1.31778e+048 1.08434e-311 4.41376e-312 1.79087e-307 3.05514e+267

-1 0 1 0 0
1 -1 0 1 0
0 1 0 0 1
-25 -9 0 0 0
1.31778e+048 1.08434e-311 4.41376e-312 1.79087e-307 3.05514e+267

0 0 1 0 0
-1 -1 0 1 0
-5 1 0 0 1
0 -9 0 0 0
1.31778e+048 1.08434e-311 4.41376e-312 1.79087e-307 3.05514e+267

0 1 1 0 0
-1 0 0 1 0
-5 -25 0 0 1
0 25 0 0 0
1.31778e+048 1.08434e-311 4.41376e-312 1.79087e-307 3.05514e+267

0 1 -1 0 0
-1 0 1 1 0
-5 -25 -9 0 1
0 25 0 0 0
1.31778e+048 1.08434e-311 4.41376e-312 1.79087e-307 3.05514e+267

0 1 -1 0 0
-1 0 1 0 0
-5 -25 -9 0 1
0 25 0 7 0
1.31778e+048 1.08434e-311 4.41376e-312 1.79087e-307 3.05514e+267

0 1 -1 0 1
-1 0 1 0 0
-5 -25 -9 0 0
0 25 0 7 -43
1.31778e+048 1.08434e-311 4.41376e-312 1.79087e-307 3.05514e+267

-1 0 1 0 0
1 -1 0 1 0
0 1 0 0 1
-25 -9 0 0 0
25 0 7 -43 0

0 0 1 0 0
-1 -1 0 1 0
-5 1 0 0 1
0 -9 0 0 0
0 0 7 -43 0

0 1 1 0 0
-1 0 0 1 0
-5 -25 0 0 1
0 25 0 0 0
0 0 7 -43 0

0 1 -1 0 0
-1 0 1 1 0
-5 -25 -9 0 1
0 25 0 0 0
0 0 9 -43 0

0 1 -1 0 0
-1 0 1 0 0
-5 -25 -9 0 1
0 25 0 7 0
0 0 9 0 0

0 1 -1 0 1
-1 0 1 0 0
-5 -25 -9 0 0
0 25 0 7 -43
0 0 9 0 43

1.#QNAN 1.#QNAN 1.#QNAN 1.#INF 1.#INF 1.#INF
1.#QNAN 1.#QNAN 1.#QNAN -1.#IND -1.#INF 1.#INF
1.#QNAN 1.#QNAN 1.#QNAN -1.#INF -1.#INF 1.#INF
1.#QNAN 1.#QNAN 1.#QNAN -1.#IND 1.#INF 1.#INF
1.#QNAN 1.#QNAN 1.#QNAN 1.#INF -1.#INF -1.#INF
1.#QNAN 1.#QNAN 1.#QNAN -1.#INF 1.#INF -1.#INF
X1=1.#QNAN
X2=1.#QNAN
X3=1.#QNAN
X4=1.#QNAN
X5=1.#QNAN
X6=1.#QNAN


(В самом конце, корни, которые, похоже на то, что делятся на 0)

Исходная матрица(находится в файле in_file) такова:
Кликните здесь для просмотра всего текста
6
1 -1 0 1 0 0 0
0 1 -1 0 1 0 0
-1 0 1 0 0 1 0
-5 -25 -9 0 0 0 -2
0 25 0 7 -43 0 17
0 0 9 0 43 -31 -20


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
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144

#include <iostream>
#include <fstream>
#include <math.h>
#include <stdio.h>
 
using namespace std;
double opr(double* , int,int);
int main()
{ 
    int N,M,i,j,i0,j0;
    double res,c,res1;
    ifstream f("in_file.txt");
    ofstream r("out_f.txt");
    f>>N;
    double mas[N][N], vd[N-1][N-1],d[N][N],x[N];
    double B[N];
    for (i=0;i<N;i++)
    {
        for (j=0;j<N;j++)
        {
            f>>mas[i][j];
        }
        f>>B[i];
        
    }
    for (i0=0;i0<N;i0++)
    {
        for (j0=0;j0<N;j0++)
        {
            
            for (i=0;i<N-1;i++)
            {
                for (j=0;j<N-1;j++)
                {
                    if (i<i0&&j<j0)
                        vd[i][j]=mas[i][j];
                    if (i<i0&&j>=j0)
                        vd[i][j]=mas[i][j+1];
                    if(i<i0&&j<j0)
                        vd[i][j]=mas[i+1][j];
                                
                }
            }
            
            for (i=0;i<N-1;i++)
            {
                for (j=0;j<N-1;j++)
                r<<vd[i][j]<<'\t';
                r<<endl;
            }
            r<<endl;
            
            res1=opr(&vd[0][0], N-1, N-1);
            cout<<"resyltat glavnii= "<<res1*pow((-1),(i0+j0))<<endl;
            d[j0][i0]=(pow((-1),(i0+j0)))*res1;
        }
    }
    
    res=opr(&mas[0][0], N, M);
    cout<<"res = "<<res<<endl;
    
    for (i=0;i<N;i++)
    {
        for (j=0;j<N;j++)
        {
            cout<<mas[i][j]<<'\t';
        }
        cout<<endl;
    }
    for (i=0;i<N;i++)
    {
        x[i]=0;
        for (j=0;j<N;j++)
        {
            r<<d[i][j]/res<<'\t';
            x[i]=x[i]+d[i][j]*B[j]/res;
        }
        r<<endl;
    }
    for (i=0;i<N;i++)
    r<<"X"<<i+1<<"="<<x[i]<<endl;
}
 
        double opr(double* ap,int N,int M)
{ int i,j,d,x=1,p,k;
double opred, koeff,pff;
if(N>1)//цикл считает определитель
{
opred=1;
p=0;
    if( *(ap+(N-1)*M+N-1)==0)
    {
        k=0;
        for(d=N-2; d>=0 && k<1;d--) 
        {
            if( *(ap+d*M+(N-1))!=0)
            {
                p=p+d-(N-1);
                x=1;
                k++;
                for (j=0;j<M;j++)
                {
                    pff= *(ap+(N-1)*M+j);
                    *(ap+(N-1)*M+j)=(*(ap+d*M+j));
                    *(ap+d*M+j)=pff;
                }
            }
            else
            {
                x=0;
            }
            cout<<k<<"\t"<<d<<endl;
        
        }
    }
    
    if(x==0)
    {
        opred=0;
    }
    else
    {
        for(i=0;i<N-1;i++)
        {
            koeff= -( *(ap+i*M+N-1))/(*(ap+(N-1)*M+N-1));
 
            for(j=0;j<N;j++)
            {   
                *(ap+i*M+j)=(* (ap+i*M+j))+koeff*( *(ap+(N-1)*M+j));
                cout<< *(ap+i*M+j)<<endl;
            }
 
        }
            opred=( *(ap+(N-1)*M+(N-1))) * (pow(-1,p)) * opr(&ap[0], (N-1), M);
    }
}
    
    else 
    {
        opred = ( *(ap+(N-1)*M+(N-1)));
    }
    
    return opred;
}
0
IT_Exp
Эксперт
34794 / 4073 / 2104
Регистрация: 17.06.2006
Сообщений: 32,602
Блог
 07.01.2015, 00:21
Ответы с готовыми решениями:

Решение системы линейных уравнений
Помогите решить на Си

Решение системы линейных уравнений
Компилятор wxDev-c++. Ошибка:C:\Program Files\Dev-Cpp\çàäà÷è\ãëàâà1\Ïåðâàÿ çàäà÷à\Makefile.win Error 1 #include &lt;cstdlib&gt; ...

Решение системы линейных уравнений.
\begin{cases} &amp; \text{ } cos(ax+2)...... x&gt;2 \\ &amp; \text{ } tg |x-2a| .... .. x\leq 2 \end{cases} x а h по а = 0.2 ...

0
Надоела реклама? Зарегистрируйтесь и она исчезнет полностью.
BasicMan
Эксперт
29316 / 5623 / 2384
Регистрация: 17.02.2009
Сообщений: 30,364
Блог
 07.01.2015, 00:21
Помогаю со студенческими работами здесь

Решение системы линейных уравнений
Не могу сделать прогу на С++. Кто сделает - огромное спасибо. Решить систему линейных уравнений 3x-5y+2z=2, 3x-7y+2z=0, x+5y-4z=-2

Найти решение системы линейных уравнений
Даны действительные числа a1, b1, c1, a2, b2, c2. Выяснить, верно ли, что a1b2-a2b1 0.0001 и если верно, то найти решение системы...

Найти решение системы линейных уравнений
привет всем, пожалуста помогите составить программу или испроваить ошибки в моей условие задачи даны действительные числа a,b,c,d,e,f....

Численное решение системы линейных уравнений
Решить систему линейных уравнений: -сформировать матрицу А коэффициентов и матрицу В правых частей (матрицы должны быть размещены в...

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


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

Или воспользуйтесь поиском по форуму:
1
Ответ Создать тему
Новые блоги и статьи
Советы по крайней бережливости. Внимание, это ОЧЕНЬ длинный пост.
Programma_Boinc 28.12.2025
Советы по крайней бережливости. Внимание, это ОЧЕНЬ длинный пост. Налог на собак: https:/ / **********/ gallery/ V06K53e Финансовый отчет в Excel: https:/ / **********/ gallery/ bKBkQFf Пост отсюда. . .
Кто-нибудь знает, где можно бесплатно получить настольный компьютер или ноутбук? США.
Programma_Boinc 26.12.2025
Нашел на реддите интересную статью под названием Anyone know where to get a free Desktop or Laptop? Ниже её машинный перевод. После долгих разбирательств я наконец-то вернула себе. . .
Thinkpad X220 Tablet — это лучший бюджетный ноутбук для учёбы, точка.
Programma_Boinc 23.12.2025
Рецензия / Мнение/ Перевод Нашел на реддите интересную статью под названием The Thinkpad X220 Tablet is the best budget school laptop period . Ниже её машинный перевод. Thinkpad X220 Tablet —. . .
PhpStorm 2025.3: WSL Terminal всегда стартует в ~
and_y87 14.12.2025
PhpStorm 2025. 3: WSL Terminal всегда стартует в ~ (home), игнорируя директорию проекта Симптом: После обновления до PhpStorm 2025. 3 встроенный терминал WSL открывается в домашней директории. . .
Как объединить две одинаковые БД Access с разными данными
VikBal 11.12.2025
Помогите пожалуйста !! Как объединить 2 одинаковые БД Access с разными данными.
Новый ноутбук
volvo 07.12.2025
Всем привет. По скидке в "черную пятницу" взял себе новый ноутбук Lenovo ThinkBook 16 G7 на Амазоне: Ryzen 5 7533HS 64 Gb DDR5 1Tb NVMe 16" Full HD Display Win11 Pro
 Музыка, написанная Искусственным Интеллектом
volvo 04.12.2025
Всем привет. Некоторое время назад меня заинтересовало, что уже умеет ИИ в плане написания музыки для песен, и, собственно, исполнения этих самых песен. Стихов у нас много, уже вышли 4 книги, еще 3. . .
От async/await к виртуальным потокам в Python
IndentationError 23.11.2025
Армин Ронахер поставил под сомнение async/ await. Создатель Flask заявляет: цветные функции - провал, виртуальные потоки - решение. Не threading-динозавры, а новое поколение лёгких потоков. Откат?. . .
Поиск "дружественных имён" СОМ портов
Argus19 22.11.2025
Поиск "дружественных имён" СОМ портов На странице: https:/ / norseev. ru/ 2018/ 01/ 04/ comportlist_windows/ нашёл схожую тему. Там приведён код на С++, который показывает только имена СОМ портов, типа,. . .
Сколько Государство потратило денег на меня, обеспечивая инсулином.
Programma_Boinc 20.11.2025
Сколько Государство потратило денег на меня, обеспечивая инсулином. Вот решила сделать интересный приблизительный подсчет, сколько государство потратило на меня денег на покупку инсулинов. . . .
КиберФорум - форум программистов, компьютерный форум, программирование
Поддержать - Реклама - Условия использования - Обратная связь
Powered by vBulletin
Copyright ©2000 - 2026, CyberForum.ru