0 / 0 / 0
Регистрация: 01.05.2012
Сообщений: 21
1

Решение СЛАУ методом Гауса-Зейделя.

01.05.2012, 20:47. Показов 2179. Ответов 1
Метки нет (Все метки)

Студворк — интернет-сервис помощи студентам
Всем привет! Простите за оффтоп, но... Сначала хотелось бы сказать спасибо всем, кто здесь реально помогает, а не разводит полемику! И конечно же спасибо администраторам, форум реально очень сильно помог!)
А теперь ближе к телу. Необходимо решить СЛАУ методом Зейделя. По сути, СЛАУ имеет вид А*х = В. Где А - матрица коэффициентов при х размером n*n, В - строка свободных членов , длиною n, ну и х - строка неизвестных членов. Короче не буду расписывать все нюансы метода, суть в том что матрица динамическая, а значения матрицы А и строки В вводятся с клавиатуры.
Сам код я написал. Задача решается. Но! Проблема в том, что весь код написан в функции main. Мне же надо, чтобы значения вызывались из другой функции (где и находится сам метод) например double zeydel.
Кроме того, прежде чем преступить непосредственно к решению необходимо проверить выполняется ли условие преобладания диагонали в матрице А, т.е. элементы, находящиеся на главное диагонали должны быть больше, чем остальные элементы в строке. И это видимо тоже должно проверяться в отдельной функции, например doble proverka. Алгоритм проверки и как ее написать тоже не проблема.
Короче помогите пожалуйста разобраться с функцией return для данной задачи!

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
//Метод Гауса-Зейделя для решенcия СЛАУ
//------------------------------------------------------------------------------
#include <stdio.h>
#include <alloc.h>
#include <conio.h>
#include <math.h>
//------------------------------------------------------------------------------
 
int main ()
{
  int i,j;
  int n;
  double **A, *B, *U, *X, *Xp, *Xn, max_s, max_n, eps;
 
  printf("Enter the size of the matrix:\nn = ");
  scanf("%d",&n);
 
  printf("\nEnter the values of the marix A:\n");
 
  A=(double**)malloc(n*sizeof(double*));
 
  for(i=0;i<n;i++)
  {
        A[i]=(double*)malloc(n*sizeof(double));
 
        for(j=0;j<n;j++)
        {
                printf("A[%d][%d] = ", i,j);
                scanf("%lf", &(A[i][j]));
        }
  }
 
  printf("\nEnter the values of the marix B:\n");
 
  B=(double*)malloc(n*sizeof(double));
 
  for(i=0;i<n;i++)
  {
        printf("B[%d] = ", i);
        scanf("%lf", &(B[i]));
  }
 
  printf("\nEnter pogreshnost:\neps = ");
  scanf("%lf", &eps);
 
  clrscr();
 
  X=(double*)malloc(n*sizeof(double));
  Xn=(double*)malloc(n*sizeof(double));
  Xp=(double*)malloc(n*sizeof(double));
 
  for(i=0;i<n;i++)
  {
        X[i] = 0;
        Xn[i] = 0;
        Xp[i] = 0;
  }
 
  U=(double*)malloc(n*sizeof(double));
 
  for(i=0;i<n;i++)
  {
        for(j=0;j<n;j++)
        {
                if(i==j)
                {
                        U[i]=A[i][j];
                }
        }
  }
 
  for(i=0;i<n;i++)
  {
        for(j=0;j<n;j++)
        {
                if(i==j)
                {
                        A[i][j] = 0;
                }
                else
                {
                        A[i][j] = -(A[i][j])/U[i];
                }
        }
        B[i] = (B[i])/U[i];
  }
 
  free(U);
 
  do
  {
        for(i=0;i<n;i++)
        {
                Xn[i] = X[i];
                X[i] = 0;
        }
 
        for (i=0;i<n;i++)
        {
                for (j=0;j<n;j++)
                {
                        X[i] = A[i][j]*Xp[j] + X[i];
                }
                X[i] = X[i] + B[i];
                Xp[i] = X[i];
        }
 
        max_s = 0;
        max_n = 0;
 
        for(i=0;i<n;i++)
        {
                max_n = fabs(X[i] - Xn[i]);
                if(max_n > max_s)
                {
                        max_s = max_n;
                }
        }
 
        for(i=0;i<n;i++)
        {
        printf("%.3lf ", (X[i]));
        }
        printf("  %lf\n", max_s);
 
  }while(max_s > eps);
 
  free(X);
  free(A);
  free(B);
 
  getch();
  return 0;
}
//------------------------------------------------------------------------------
Добавлено через 7 минут
Как я понимаю должно получится что-то вроде:

C
1
2
3
4
5
6
7
8
9
10
11
12
13
14
double proverka (double **A, int n)
//проверка преобладания диагонали
if (...) return p = 1;
else reurn p = 0.
 
double zeydel (double **A, double *B, double *X, double eps, int n)
//решение СЛАУ
return X[i];
 
int main()
//введение с клавиатуры **A, *В, eps, n
//вывод результатов
1) СЛАУ не может быть решена
2) Решение
0
Programming
Эксперт
94731 / 64177 / 26122
Регистрация: 12.04.2006
Сообщений: 116,782
01.05.2012, 20:47
Ответы с готовыми решениями:

Решение СЛАУ методом Зейделя (Гаусса-Зейделя) с заданной матрицей
ребята, помогите написать на си, решение слау методом Зейделя(Гаусса-Зейделя) с заданной матрицей

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

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

Решение системы уравнений методом Зейделя
Привет. У меня есть код программы, которая решает систему. Можете посмотреть, или это метод...

1
0 / 0 / 0
Регистрация: 01.05.2012
Сообщений: 21
02.05.2012, 21:01  [ТС] 2
Неужели никто не может помочь!?
0
IT_Exp
Эксперт
87844 / 49110 / 22898
Регистрация: 17.06.2006
Сообщений: 92,604
02.05.2012, 21:01
Помогаю со студенческими работами здесь

Решение СЛАУ методом Гаусса
Задание 7. Алгоритм Гаусса. Реализовать алгоритм решения СЛАУ методом Гаусса. Ограничения:...

Решение СЛАУ методом LU-разложения
Написал лабу по по выч.мату. решил поделиться. Надеюсь кому то еще пригодиться. #include...

Решение СЛАУ методом Гаусса
помогите, пожалуйста! надо создать программу, которая считает систему линейных уравнений методом...

Решение СЛАУ методом Крамера
Нужно написать программу на СИ для решения СЛАУ с помощью метода Крамера. 2x1 - x2 + 3x3 + 2x4...


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

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

КиберФорум - форум программистов, компьютерный форум, программирование
Powered by vBulletin
Copyright ©2000 - 2023, CyberForum.ru