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

СЛАУ метод Зейделя - C++

Восстановить пароль Регистрация
 
Рейтинг: Рейтинг темы: голосов - 41, средняя оценка - 4.71
maloy-rom
0 / 0 / 0
Регистрация: 15.12.2009
Сообщений: 10
17.02.2010, 00:22     СЛАУ метод Зейделя #1
Решить СЛАР методом Зейделя вида АХ=В с точностью до 0,0001. Пример на рисунке ниже

что не могу не как реализовать нормально. Точнее недопонимаю самого метода

Язык С

Помогите буду очень благодарен
Миниатюры
СЛАУ метод Зейделя  
Similar
Эксперт
41792 / 34177 / 6122
Регистрация: 12.04.2006
Сообщений: 57,940
17.02.2010, 00:22     СЛАУ метод Зейделя
Посмотрите здесь:

Метод Зейделя C++
C++ Метод итераций и метод Зейделя
Решение СЛАУ методом Зейделя C++
C++ Метод Зейделя
C++ Метод Зейделя в С++
После регистрации реклама в сообщениях будет скрыта и будут доступны все возможности форума.
M128K145
Эксперт C++
 Аватар для M128K145
8272 / 3491 / 142
Регистрация: 03.07.2009
Сообщений: 10,707
17.02.2010, 00:32     СЛАУ метод Зейделя #2
Думаю, не смертельно, что на Java и сделано для трех уравнений
Java
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
public class Zeidel extends IOStream
{
    final private int n = 3;
    final private double[][] koef = new double[n][n];
    private double[] otvet = new double[n];
    private double eps = 0.0001; 
    void InputKoef()
    {
        PrintLine("Введите матрицу коэфициентов:");
        int i, j;
        for(i = 0; i < n; ++i)
            for(j = 0; j < n; ++j)
                koef[i][j] = ReadDouble();
        PrintLine("Введите ответы:");           
        for(i = 0; i < n; ++i)
            otvet[i] = ReadDouble();    
    }
    void Out()
    {
        int i, j;
        PrintLine("Система уравнений:");
        for(i = 0; i < n; ++i, PrintLine())
            for(j = 0; j < n; ++j)
                Print(koef[i][j] + "x" + (j + 1) + 
                      (j == n - 1 ? " = " + otvet[i] : " + "));
    }
    private double[] x = new double[n];
    private boolean Equal(double x0, double x1, double x2)
    {
        return
            (((koef[0][0] * x0 + koef[0][1] * x1 + koef[0][2] * x2) - otvet[0] <= eps)&&
            ((koef[1][0] * x0 + koef[1][1] * x1 + koef[1][2] * x2) - otvet[1] <= eps) &&
            ((koef[2][0] * x0 + koef[2][1] * x1 + koef[2][2] * x2) - otvet[2] <= eps));
    }
    
    double[] z = new double[n];
 
    private static double Round(double chislo, int sign, int fract)
    {
        return (double)(int)(chislo * Math.pow(10, sign + fract)) / Math.pow(10, fract);
    }
    void Calculate()
    {
        int i, j;
        for(i = 0; i < n; ++i)
        {
            otvet[i] /= koef[i][i];
            z[i] = koef[i][i];
            for(j = 0; j < n; ++j)
                koef[i][j] /= z[i];
        }
        x[1] = x[2] = 0;
        i ^= i;
        do
        {
            x[0] = otvet[0] - koef[0][1] * x[1] - koef[0][2] * x[2];
            x[1] = otvet[1] - koef[1][0] * x[0] - koef[1][2] * x[2];
            x[2] = otvet[2] - koef[2][0] * x[0] - koef[2][1] * x[1];
            PrintLine("\ni = " + ++i +"\tx1 = " + Round(x[0], 4, 5) + "\tx2 = " + Round(x[1], 4, 5) + "\tx3 = " + Round(x[2], 4, 5));
        }while(!Equal(x[0], x[1], x[2]));
        PrintLine("\nx1 = " + Round(x[0], 4, 5) + "\tx2 = " + Round(x[1], 4, 5) + "\tx3 = " + Round(x[2], 4, 5));
    }
    void Verification()
    {
        int i, j;
        PrintLine("\n\nСистема уравнений:");
        for(i = 0; i < n; ++i, PrintLine())
            for(j = 0; j < n; ++j)
                Print(koef[i][j]*z[i] + " * " + Round(x[j], 4, 5) + 
                      (j == n - 1 ? " = " + otvet[i]*z[i] : " + "));
    }
}
класс в котором вызывается
Java
1
2
3
4
5
6
7
8
9
10
11
12
public class Main extends IOStream
{
    public static void main(String[] args) 
    {
        jar = false;
        Zeidel z = new Zeidel();
        z.InputKoef();
        z.Out();
        z.Calculate();
        z.Verification();
    }
}
maloy-rom
0 / 0 / 0
Регистрация: 15.12.2009
Сообщений: 10
17.02.2010, 21:56  [ТС]     СЛАУ метод Зейделя #3
я с JAVA вообще не знаком

думаю что то типа такого но полностью реализовать не могу

C
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
for(i=0;i < N;i++)
{
x[i]=-b[i];
 
for(j=0;j < N;j++)
{
if(i!=j)
x[i]+=a[i][j]*x[j];
}
 
x[i]/=-a[i][i];
}
 
for(i=0;i < N;i++)
{
if(fabs(x[i]-xn[i]) > norma)
norma=fabs(x[i]-xn[i]);
xn[i]=x[i];
}
ДинисламХХХ
5 / 5 / 1
Регистрация: 29.11.2010
Сообщений: 99
01.12.2010, 19:04     СЛАУ метод Зейделя #4
maloy-rom, это же метод простых итераций??
Bekmurat
0 / 0 / 0
Регистрация: 14.02.2012
Сообщений: 7
10.03.2012, 10:32     СЛАУ метод Зейделя #5
Не нат Java а на С++

Добавлено через 1 минуту
А как будет полностю
Yandex
Объявления
10.03.2012, 10:32     СЛАУ метод Зейделя
Ответ Создать тему
Опции темы

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