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

С++ для начинающих

Войти
Регистрация
Восстановить пароль
 
Рейтинг: Рейтинг темы: голосов - 14, средняя оценка - 4.79
zidans
0 / 0 / 0
Регистрация: 30.10.2010
Сообщений: 45
#1

Матрицы: LU – разложение - C++

28.04.2012, 21:23. Просмотров 1844. Ответов 0
Метки нет (Все метки)

Сделана программа по заданию: Решить систему линейных алгебраических уравнений Ax=b методом LU – разложения (предварительно изменив порядок записи уравнений);

Помогите, пожалуйста вывести в ответе матрицы L и U.

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
#include "stdio.h"
#include "conio.h"
#include "iostream.h"
int const N = 4;
double a[N][N] = {   {25,     2,  -2,     1},
             {0,   17.3,   3,   0.9},
             {-4,  -0.9,   1,  18.7},
             {2,      5,  25,    -2}};
double b[N] = {236, -44.5, -22.1, -49};
double l[N][N], u[N][N];
double y[N], x[N];
int s, i, j, n = N;
double sum, tmp;
 
main()
{
  
  for (i = 0; i < N; i++)
  {
    for (j = 0; j < N; j++)
    {
      l[i][j] = 0;
      u[i][j] = 0;
    }
    y[i] = 0;
    x[i] = 0;
    u[i][i] = 1;
 }
  for (i = 0; i < n; i++)
    for (j = 0; j < n; j++)
    {
      sum = 0;
      if (i >= j)
      {
    for (s = 0; s <= j - 1; s++)
      sum += (l[i][s]) * (u[s][j]);
    l[i][j] = a[i][j] - sum;
      }
      else
      {
    for (s = 0; s <= i - 1; s++)
      sum += l[i][s] * u[s][j];
    u[i][j] = (a[i][j] - sum) / l[i][i];
      }
    }
  y[0] = b[0] / l[0][0];
  for (i = 1; i < n; i++)
  {
    sum = 0;
    for (j = 0; j < i; j++)
      sum += l[i][j] * y[j];
    y[i] = (b[i] - sum) / l[i][i];
 }
 
  x[n - 1] = y[n - 1];
  for (i = n - 2; i >= 0; i--)
  {
    sum = 0;
    for (j = n - 1; j >= i + 1; j--)
      sum += x[j] * u[i][j];
    x[i] = y[i] - sum;
  }
  for(i=0;i<n;i++)
    cout << x[i] << ' ' << endl;;
  while (!kbhit()) {;}
}
Similar
Эксперт
41792 / 34177 / 6122
Регистрация: 12.04.2006
Сообщений: 57,940
28.04.2012, 21:23     Матрицы: LU – разложение
Посмотрите здесь:

Разложение матрицы A на треугольные матрицы - C++
Создать разложения матрицы A размером N x N на треугольные матрицы, удовлетворяющие условию PA=LU, с выбором главного элемента (возможное...

LU разложение матрицы - C++
хотел написать программу для LU разложения матрицы. взял формулы с вики , не правильно находит элемент (3,3) в матрице U. (пробовал...

LU Разложение - C++
Есть код. #include &lt;iostream&gt; #include &lt;vector&gt; using namespace std; void LU(vector &lt;vector &lt;double&gt;&gt; A, vector &lt;vector...

LU-разложение - C++
Решить матрицу методом LU. матрицы задаются рандомно и выводятся на экран, если можно с коментами #include &quot;stdafx.h&quot; #include &quot;stdio.h&quot; ...

Разложение в ряд S(x) и Y(x) - C++
Всем привет! Необходимо вывести на экран таблицу значений функции Y(x) и ее разложения в ряд S(x) для x изменяющихся от xn до xk с...

Разложение в ряд - C++
#include &lt;stdio.h&gt; #include &lt;math.h&gt; int main () { double a,b,x,k; //k - to4noe printf(&quot;BBegute x\n&quot;); scanf(&quot;%lf&quot;,&amp;x); ...

разложение числа - C++
Как ,допустим, разложить число 1924 на 1 9 2 4. Даже идей нет

Разложение числа - C++
Всем привет! Есть некое число N и массив arr (k - размер массива). Нужно написать программу, которая выведет на экран все возможные...

Разложение числа - C++
Число можно разложить по 5 и по 3, то есть если это допустим 8 то выйдет 5 и 3, причем не должно быть остатка, допустим 22 можно разложить...

Разложение даты - C++
Как ввести дату в формате :ДД.ММ.ГГГГ Например:29.04.1995 И потом нужно обрабатывать день отдельно месяц отдельно и год отдельно . . . ...

разложение тангенса - C++
Бодрого утра....вот не могу справиться с ошибками...буду очень благодарна, кто поможет... #include &lt;iostream&gt; #include &lt;conio.h&gt; ...

Разложение в ряд - C++
Помогите пожалуйста Функция Разложение в ряд Область сходимости


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

Или воспользуйтесь поиском по форуму:
После регистрации реклама в сообщениях будет скрыта и будут доступны все возможности форума.
Ответ Создать тему
Опции темы

КиберФорум - форум программистов, компьютерный форум, программирование
Powered by vBulletin® Version 3.8.9
Copyright ©2000 - 2017, vBulletin Solutions, Inc.
Рейтинг@Mail.ru