Форум программистов, компьютерный форум CyberForum.ru

Решить уравнение методом Гаусса - C++

Восстановить пароль Регистрация
 
Salenser
 Аватар для Salenser
3 / 3 / 0
Регистрация: 09.12.2010
Сообщений: 21
25.03.2012, 21:57     Решить уравнение методом Гаусса #1
надо посчитать методом гаусса ур-е http://www.cyberforum.ru/cgi-bin/latex.cgi?A \cdot X = B

3-размерность

http://www.cyberforum.ru/cgi-bin/latex.cgi?<br />
A = \begin{pmatrix}<br />
2.74 & -1.18 & 3.17 \\<br />
1.12 & 0.83 & -2.16 \\ <br />
0.81 & 1.27 & 0.76<br />
\end{pmatrix}<br />
B = \begin{pmatrix}<br />
2.18 \\ -1.15 \\ 3.23<br />
\end{pmatrix}<br />
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
#include <iostream>
#include <fstream>
#include <stdlib.h>
#include <iostream>
#include <fstream>
#include <math.h>
#include <sstream>
#include <time.h>
#include <string.h>
#include <map>
using namespace std;
 
int main()
 
{ ifstream    in("input.txt"); //создание файла с которого будем считывать
  ofstream  out("output.txt") ; // создание файла с которого будем записывать
    int n ;
    int h;
 
   in>>n; размерность
 
    double **a=new double*[n];
for (int i=0; i<n;i++)
a[i]=new double[n];// инциализации матрицы А
    for( int i=0;i<=n-1 ;i++)
       for(int j=0;j<=n-1;j++)
       {
         in>>a[i][j];
    }
float b[n];
float x[n];
for (int i=0; i<=n-1; i++)
{
    in>>b[i];
 
 
}
for (int k=1;k<=n-1;k++) //прямой ход
    for (int i=k+1;i<=n;i++)
   {a[i][k]=a[i][k]/a[k][k];
 
        for (int j=k+1;j<=n;j++)
        a[i][j]=a[i][j]-a[k][j]*a[i][k];
 
    b[i]=b[i]-b[k]*a[i][k];
    }
     x[n]=b[n]/a[n][n];//обратный ход
  for (int i=n-1;i>=1;i--) {
    h=b[i];
    for (int j=i+1;j<=n;j++)
     h=h-a[i][j]*x[j];
      x[i]=h/a[i][i];
  }
    int s = 0;
    for (int i=0;i<=n;i++)
    for (int j=0;j<=n;j++)
    {cout<<a[i][j];
    cout<<"/t";
    s++;
    if ( s==n){cout<<"\n";
    s=0;}
    }
    for (int i=0;i<=n;i++)
    {
        cout<<x[i];
        cout<<"\n";
    }
    return 0;
}
Добавлено через 15 минут
когда я запускаю прогу у меня виндовс 7 сразу завершает программу
не пойму где ошибка , вроде проганял а без результатно.
После регистрации реклама в сообщениях будет скрыта и будут доступны все возможности форума.
Kuzia domovenok
 Аватар для Kuzia domovenok
1883 / 1738 / 116
Регистрация: 25.03.2012
Сообщений: 5,907
Записей в блоге: 1
25.03.2012, 22:47     Решить уравнение методом Гаусса #2
портится куча? Всё что выделено new[] должно быть освобождено delete[]
память под b[n] и x[n] тоже не выделена.
не понятно, почему в одних местах используется float, а в других double
ну и проверь в конце концов, успешно ли у тебя вся информация читается из файла и сохраняется в a и b?
Поставь точку останова и проверь, как прочиталось.
Операция "x[i]=h/a[i][i];" может создать проблему деления на ноль. Это надо учесть.
переменная h тоже не понятно почему целочисленная.
Dima249
10 / 10 / 0
Регистрация: 02.10.2011
Сообщений: 69
25.03.2012, 22:53     Решить уравнение методом Гаусса #3

Не по теме:

у меня где-то есть такая только криво работает щас найду и скину сюда



Добавлено через 1 минуту
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
#include <stdio.h>
#include <stdlib.h>
#include <iostream>
#include <time.h>
using namespace std;
int main(int argc, char *argv[])
{
    const int b=100;
    //Матрица системы
    int  A[b][b];
    double summa[b], A1[b][b], me, summa2[b];
    int l,o,t=0,n,m,i,j,k;
    cout<<"Vvtdite n: ";
    cin>>n;
    cout<<"Vvtdite m: ";
    cin>>m;
    for (i=0; i<n; i++)
    {
        for (j=0;j<m;j++) 
        {
            printf(" A[%i][%i]:",i+1,j+1);
            scanf("%i",&A[i][j]);
        }
    }
    for(i = 0; i < n; ++i) 
    {
        summa[i]=0;
        summa2[i]=0;
        for(j = 0; j < m; ++j) cout<<A[i][j]<<"\t"; 
        cout<<endl;
    }
 
    for(i = 0; i < n; ++i) 
    {     
        for(j = 0; j < m; ++j) 
            A1[i][j]=A[i][j];}
    //Число перестановок строк
    int rCnt=0;
    //Ведущий элемент
    double ldn;
    //Номер ведущего элемента
    int index;
    double temp;
 
    //Прямой ход
 
    for(j=0; j<m; j++)
    {
        //Выбор ведущего элемента столбца
        ldn=A1[j][j];
        index=j;
        for(i=j; i<n; i++)
            if(fabs(A1[i][j])>fabs(ldn))
            {
                ldn=A1[i][j];
                index=i;
            }
        //Если ведущий элемент - не на гл. диагонали
        //меняем строки местами
        if(index!=j)
        {
            ++rCnt;
            for(i=0; i<n; i++)
            {
                temp=A1[j][i];
                A1[j][i]=A1[index][i];
                A1[index][i]=temp;   
            }
        }
                
        for (i = 0; i<m; i++)
            A1[j][i] = A1[j][i]/ldn; //делим строку на ведущий элемент
 
        //Обнуление элементов ниже гл. диагонали
        for(i=j+1; i<n; ++i)
        {
            for(k=0; k<m; k++)
                A1[i][k]= fabs (A1[i][k]-A1[j][k]*A1[i][k]);
        }       
    }
 
    for(i = 0; i < n; i++)  
    {
        for(j = 0; j < m; j++) 
            cout<<A1[i][j]<<"\t"; 
        cout<<endl; 
    }
    for (i=0; i<n; i++)
    { 
        for (j=0; j<m; j++)
            summa2[i]=summa2[i]+A1[i][j];
        //cout<<summa[i]<<endl;
    }
 
    for (i=0; i<n; i++)
    {
        summa[i]=summa[i]/m;
        summa2[i]=summa2[i]/m;
        if (summa2[i]<summa[i]) 
        {
            t=t+1;
        }
    }
    system("PAUSE");       //для просмотра результата
    return 0;
}
тут ещё сумма я не помню зачем она) ну это задание из книжки лаба у меня такая была и тут только преобразование до треугольного вида
Kuzia domovenok
 Аватар для Kuzia domovenok
1883 / 1738 / 116
Регистрация: 25.03.2012
Сообщений: 5,907
Записей в блоге: 1
25.03.2012, 23:02     Решить уравнение методом Гаусса #4
Цитата Сообщение от Dima249 Посмотреть сообщение

Не по теме:

у меня где-то есть такая только криво работает щас найду и скину сюда



Добавлено через 1 минуту
тут ещё сумма я не помню зачем она) ну это задание из книжки лаба у меня такая была и тут только преобразование до треугольного вида
а при чём тут решение уравнений?
zewer
 Аватар для zewer
1019 / 710 / 71
Регистрация: 07.01.2011
Сообщений: 5,367
25.03.2012, 23:05     Решить уравнение методом Гаусса #5
а может запускай програму через Ctrl + F5 ???
Dima249
10 / 10 / 0
Регистрация: 02.10.2011
Сообщений: 69
26.03.2012, 07:18     Решить уравнение методом Гаусса #6
Цитата Сообщение от Kuzia domovenok Посмотреть сообщение
а при чём тут решение уравнений?
ну ему нужна как я понял программа которая будет решать уравнение системой гаусса в задании у него написано, но если просто составить матрицу треугольного вида, то моя программа подходит
-=ЮрА=-
Заблокирован
Автор FAQ
26.03.2012, 12:14     Решить уравнение методом Гаусса #7
Salenser, посмотри сюда
Методом Гаусса решить систему n линейных алгебраических уравнений
MoreAnswers
Эксперт
37091 / 29110 / 5898
Регистрация: 17.06.2006
Сообщений: 43,301
26.03.2012, 17:17     Решить уравнение методом Гаусса
Еще ссылки по теме:

C++ Дано уравнение ax2+bx+c=0. Решить уравнение, результат вывести на экран.
C++ Решить уравнение методом касательных с использованием функций
C++ Решение системы линейного уравнение методом Гаусса

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

Или воспользуйтесь поиском по форуму:
Salenser
 Аватар для Salenser
3 / 3 / 0
Регистрация: 09.12.2010
Сообщений: 21
26.03.2012, 17:17  [ТС]     Решить уравнение методом Гаусса #8
спасибо всем
Yandex
Объявления
26.03.2012, 17:17     Решить уравнение методом Гаусса
Ответ Создать тему
Опции темы

Текущее время: 02:24. Часовой пояс GMT +3.
КиберФорум - форум программистов, компьютерный форум, программирование
Powered by vBulletin® Version 3.8.9
Copyright ©2000 - 2016, vBulletin Solutions, Inc.
Рейтинг@Mail.ru