Метод простых итераций для решения СЛАУ (консоль -> форма)

@BAOlOEJlbLLlUK · Регистрация: 28.04.2013

Студворк — интернет-сервис помощи студентам

Доброго времени суток. Есть код реализации метода простых итераций для решения СЛАУ в консоли, нужно сделать так, чтобы матрица А заполнялась из одного Стринггрида, B из второго, n бралось из cspinedita, ну и eps из edita.Вот код:

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
#include <vcl.h>
#pragma hdrstop
#include <stdio.h>
#include <conio.h>
#include <math.h>
 
 
    const n = 5;
 
    float A[n][n] = {
                    { 79.2,   0,  35, 19.8,  24 },
                    { 39.6,  85,   0, 19.8,  25 },
                    { 19.8, -15,  45,    0,  10 },
                    { 49.5,  18,  20, 89.1,   0 },
                    {  9.9,  15,  20,-49.5,  95 },
                    };
 
    float  B[n]  = { 86, 55, 77, 5, -64};
 
    float  X[30][n] = {0};
    float  TEST[n];
    float  summ;
    float  eps;
    int    k = 0;
 
 
bool ExitFunc(float eps)
{
    if(k == 0) return true;
    for(int i=0; i<n; i++){
    if( fabs(X[k][i]-X[k-1][i]) > eps )
        return true;
    }return false;      
}
 
void main(void)
    {   
    printf("Enter eps:");
    scanf("%f",&eps);
    do 
    {
    for(int i=0; i<n; i++)        
        {
        summ = 0;
        for(int j=0; j<n; j++)
           if(i!=j)
              summ += A[i][j] * X[k][j];
        X[k+1][i] = (1/A[i][i]) * (B[i] - summ);            
        }
    k++;
    }while(ExitFunc(eps));
 
    //Проверка
    printf("\nResult:\n");
    for(int j=0; j<n; j++)
       printf("%4.8f   ",X[k][j]);
 
    printf("\nTest:\n");
    for(int i=0; i<n; i++)
        {
        for(int j=0; j<n; j++)
           TEST[i] +=(A[i][j] * X[k][j]);
        printf("%4.8f  ", TEST[i]);
        }
 
    printf("\nIterations:\n%i",k);
    getch();
    }

Помогите пожалуйста.Буду благодарен.

Добавлено через 8 минут
Попробовал реализовать, но выводит ошибки:

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
double X[30][n]={0};//тут: constant expression required
int    k = 0;
bool ExitFunc(float eps){// и тут: declaration syntax error
    if(k == 0){ return true;}
    for(int i=0; i<n; i++){
    if( fabs(X[k][i]-X[k-1][i]) > eps){
        return true;}
    }return false;
}
int i,j;
double** A;
double* B;
float  summ;
float  eps;
Edit1->Text="0,0001";
bool ExitFunc(float eps);
n=StrToInt(CSpinEdit4->Value);
B=new double[n];
//float  TEST[n];
for(int i=0;i<n;i++){
  A[i]=new double[n];
  }
for(int i=0; i<n ;i++){
    for (int j=0;j<n;j++){
        A[i][j]=StrToFloat(StringGrid8->Cells[i][j]);
        B[i]=StrToFloat(StringGrid9->Cells[0][j]);
        }
        }
 
eps=StrToFloat(Edit1->Text);
 
do
{
  for(int i=0; i<n; i++)
  {
    summ = 0;
    for(int j=0; j<n; j++)
    if(i!=j)
      summ += A[i][j] * X[k][j];
    X[k+1][i] = (1/A[i][i]) * (B[i] - summ);
    Memo4->Lines->Add(FloatToStr(k+1)+"-ÿ èòåðàöèÿ:");
    Memo4->Lines->Add(X[k+1][i]);
 
  }
  k++;
}while(ExitFunc(eps));
 
for(int j=0; j<n; j++){
       StringGrid10->Cells[0][j]=FloatToStr(X[k][j]);}

Добавлено через 28 минут
И все это происходит после нажатия кнопки

@BRcr · 19.04.2015, 21:42

Я не понял, вот этот последний код - он прописан в обработчике кнопки?

А в целом, чехарда какая-то. Размеры статического массива нельзя задавать не-константами, прототипы размещать после определения функций тоже ни к чему.

@BAOlOEJlbLLlUK · 20.04.2015, 00:28 **[ТС]**

Да, в кнопке, помогите, как правильно будет

@BRcr · 20.04.2015, 00:36

Для начала. Правильно будет - вынести определение функции ExitFunc за пределы обработчика клика кнопки. И задавать размеры массива X[][] константами. И объявление обозначенного массива тоже следует вынести либо в глобальную область видимости, либо в класс, где будет объявлен метод ExitFunc...

@BAOlOEJlbLLlUK · 20.04.2015, 12:35 **[ТС]**

но мне нужно значение n брать из edita

Добавлено через 37 секунд
т.е. из cspinedita

@BRcr · 20.04.2015, 13:04

Значит, статический массив отпадает.
Вот так можно создать и использовать массив произвольного размера:

C++
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
// объявляем, причем лучше в классе формы...
double * arr[ 30 ];
 
void __fastcall TForm1::btn_2Click( TObject * Sender )
{
    // создаем
    long N( sedit_1->Value );
 
    for ( int i( 0 ), i_limit( 30 ); i < i_limit; ++i )
    {
        arr[ i ] = new double[ N ];
    }
 
    // как-то используем массив
    // ...
 
    // и потом освобождаем его
    for ( int i( 0 ), i_limit( 30 ); i < i_limit; ++i )
    {
        delete[ ]arr[ i ];
    }
}

@BAOlOEJlbLLlUK · 20.04.2015, 13:40 **[ТС]**

после нажатия на кнопку выдает ошибку:
вот скрин

@BAOlOEJlbLLlUK · 20.04.2015, 14:35 **[ТС]**

вроде исправил написал

C++
1
2
3
4
A = new double * [n];
for(int i=0;i<n;i++){
 A[i]=new double[n];
  }

Добавлено через 5 минут
но потом после заполнения массивов и нажатия на кнопку новая ошибка:

@BRcr · 20.04.2015, 20:38

Хто такое массив "В"?
Лучше всего приводи полный код в его последнем виде. Иначе непонятно по одним только скринам, что не так.

Также изображения и прочие файлы следует загружать на форум. Расширенный режим редактора -> управление вложениями.

@_Dimon_ · 20.04.2015, 22:47

Сообщение от BAOlOEJlbLLlUK

но потом после заполнения массивов и нажатия на кнопку новая ошибка:

X[k+1][i]
что такое X и что такое k?

@BAOlOEJlbLLlUK · 21.04.2015, 23:51 **[ТС]**

B - это столбец свободных коэффициентов, X- это вектор решения, k- это итерация

Добавлено через 3 минуты
функция:

C++
1
2
3
4
5
6
7
8
9
10
int k=0;
double **X;
bool ExitFunc(float eps);
bool ExitFunc(float eps){
    if(k == 0){ return true;}
    for(int i=0; i<n; i++){
    if( fabs(X[k][i]-X[k-1][i]) > eps){
        return true;}
    }return false;
}

обработчик клика

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
int n;
n=StrToInt(CSpinEdit4->Value);
//double
//int    k = 0;
X = new double * [n];
for ( int i=0; i<n; ++i )
    {
        X[i] = new double[n];
    }
//int i,j;
double** A;
double* B;
float  summ;
float  eps;
//Edit1->Text="0,0001";
B=new double[n];
//float  TEST[n];
A = new double * [n];
for(int i=0;i<n;i++){
 A[i]=new double[n];
  }
for(int i=0; i<n ;i++){
    for (int j=0;j<n;j++){
        A[i][j]=StrToFloat(StringGrid8->Cells[i][j]);
        B[i]=StrToFloat(StringGrid9->Cells[0][j]);
        }
        }
do
{
  for(int i=0; i<n; i++)
  {
    summ = 0;
    for(int j=0; j<n; j++){
    if(i!=j){
      summ += A[i][j] * X[k][j];}}
    X[k+1][i] = (1/(A[i][i])*(B[i] - summ));
    Memo4->Lines->Add(FloatToStr(k+1)+"-я итерация:");
    Memo4->Lines->Add(X[k+1][i]);
 
  }
  k++;
}while(ExitFunc(eps));
 
for(int j=0; j<n; j++){
       StringGrid10->Cells[0][j]=FloatToStr(X[k][j]);}
for ( int i=0; i<n ; ++i )
    {
        delete[]X[i];
    } 
}

Новые блоги и статьи Все статьи Все блоги /
Уведомление о неверно выбранном значении справочника Maks 06.04.2026 Алгоритм из решения ниже реализован на примере нетипового документа "НарядПутевка", разработанного в конфигурации КА2. Задача: уведомлять пользователя, если в документе выбран неверный склад. . .	Установка Qt Creator для C и C++: ставим среду, CMake и MinGW без фреймворка Qt 8Observer8 05.04.2026 Среду разработки Qt Creator можно установить без фреймворка Qt. Есть отдельный репозиторий для этой среды: https:/ / github. com/ qt-creator/ qt-creator, где можно скачать установщик, на вкладке Releases:. . .	AkelPad-скрипты, структуры, и немного лирики.. testuser2 05.04.2026 Такая программа, как AkelPad существует уже давно, и также давно существуют скрипты под нее. Тем не менее, прога живет, периодически что-то не спеша дополняется, улучшается. Что меня в первую очередь. . .	Отображение реквизитов в документе по условию и контроль их заполнения Maks 04.04.2026 Алгоритм из решения ниже реализован на примере нетипового документа "ПланированиеСпецтехники", разработанного в конфигурации КА2. Данный документ берёт данные из другого нетипового документа. . .
Фото всей Земли с борта корабля Orion миссии Artemis II kumehtar 04.04.2026 Это первое подобное фото сделанное человеком за 50 лет. Снимок называют новым вариантом легендарной фотографии «The Blue Marble» 1972 года, сделанной с борта корабля «Аполлон-17». Новое фото. . .	Вывод диалогового окна перед закрытием, если документ не проведён Maks 04.04.2026 Алгоритм из решения ниже реализован на примере нетипового документа "СписаниеМатериалов", разработанного в конфигурации КА2. Задача: реализовать программный контроль на предмет проведения документа. . .	Программный контроль заполнения реквизитов табличной части документа Maks 02.04.2026 Алгоритм из решения ниже реализован на примере нетипового документа "СписаниеМатериалов", разработанного в конфигурации КА2. Задача: 1. Реализовать контроль заполнения реквизита. . .	wmic не является внутренней или внешней командой Maks 02.04.2026 Решение: DISM / Online / Add-Capability / CapabilityName:WMIC~~~~ Отсюда: https:/ / winitpro. ru/ index. php/ 2025/ 02/ 14/ komanda-wmic-ne-naydena/

@BRcr 4043 / 2333 / 292 Регистрация: 03.02.2011 Сообщений: 5,066 Записей в блоге: 10
	19.04.2015, 21:42
	Я не понял, вот этот последний код - он прописан в обработчике кнопки? А в целом, чехарда какая-то. Размеры статического массива нельзя задавать не-константами, прототипы размещать после определения функций тоже ни к чему. 0

@BAOlOEJlbLLlUK 0 / 0 / 0 Регистрация: 28.04.2013 Сообщений: 41
	20.04.2015, 00:28 [ТС]
	Да, в кнопке, помогите, как правильно будет 0

@BRcr 4043 / 2333 / 292 Регистрация: 03.02.2011 Сообщений: 5,066 Записей в блоге: 10
	20.04.2015, 00:36
	Для начала. Правильно будет - вынести определение функции ExitFunc за пределы обработчика клика кнопки. И задавать размеры массива X[][] константами. И объявление обозначенного массива тоже следует вынести либо в глобальную область видимости, либо в класс, где будет объявлен метод ExitFunc... 0

@BAOlOEJlbLLlUK 0 / 0 / 0 Регистрация: 28.04.2013 Сообщений: 41
	20.04.2015, 12:35 [ТС]
	но мне нужно значение n брать из edita Добавлено через 37 секунд т.е. из cspinedita 0

@BAOlOEJlbLLlUK 0 / 0 / 0 Регистрация: 28.04.2013 Сообщений: 41
	20.04.2015, 13:40 [ТС]
	после нажатия на кнопку выдает ошибку: вот скрин Миниатюры 0

@BRcr 4043 / 2333 / 292 Регистрация: 03.02.2011 Сообщений: 5,066 Записей в блоге: 10
	20.04.2015, 20:38
	Хто такое массив "В"? Лучше всего приводи полный код в его последнем виде. Иначе непонятно по одним только скринам, что не так. Также изображения и прочие файлы следует загружать на форум. Расширенный режим редактора -> управление вложениями. 0