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

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

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

LU-разложение - C++

06.09.2011, 11:59. Просмотров 7909. Ответов 3
Метки нет (Все метки)

Решить матрицу методом LU. матрицы задаются рандомно и выводятся на экран, если можно с коментами
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
#include "stdafx.h"
#include "stdio.h"
#include "conio.h"
#include "math.h"
#include "iostream"
//using namespace std;
 
int const N=5;
double a[N][N],l[N][N],u[N][N];
double b[N],y[N],x[N];
int s,i,j,n;
double sum,tmp;
 
int 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;
 }
 //==========================================================
 cin>>n;
 for(i=0;i<n;i++)
  for(j=0;j<n;j++)
   cin>>a[i][j]; 
 for(i=0;i<n;i++) 
   cin>>b[i];
//===========================================================
 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]<<' ';
 cout<<'\n';
 for(i=0;i<n;i++){
     for(j=0;j<n;j++)
         cout<<l[i][j]<<' ';
     cout<<endl;
     }
 for(i=0;i<n;i++){
     for(j=0;j<n;j++)
         cout<<u[i][j]<<' ';
     cout<<endl;
     }
 
     printf("%4lf ",x[i]);
 getch();
 return 0;
}
Similar
Эксперт
41792 / 34177 / 6122
Регистрация: 12.04.2006
Сообщений: 57,940
06.09.2011, 11:59
Здравствуйте! Я подобрал для вас темы с ответами на вопрос LU-разложение (C++):

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...

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

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

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

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

Надоела реклама? Зарегистрируйтесь и она исчезнет полностью.
sandye51
программист С++
682 / 584 / 39
Регистрация: 19.12.2010
Сообщений: 2,016
06.09.2011, 12:09 #2
так решение написано, правда не проверял верное или нет
что еще надо?
DVJ
1 / 1 / 0
Регистрация: 30.05.2011
Сообщений: 14
06.09.2011, 13:22  [ТС] #3
Дело в том что он не выводит матрицы, и матрицы задаются константой, делаю random вообще не чего не выходит

Добавлено через 8 минут
http://pcfu.ru/lu-razlozhenie/ тут вроде как все описано, но как то мутно, понять не чего не могу... Можно на другом языке переписать, знаю Pascal, c#, delphi, c++, asambler. ну и по мелочам так
sandye51
06.09.2011, 13:49     LU-разложение
  #4

Не по теме:

Цитата Сообщение от DVJ Посмотреть сообщение
знаю Pascal, c#, delphi, c++, asambler. ну и по мелочам так
рассмешил

MoreAnswers
Эксперт
37091 / 29110 / 5898
Регистрация: 17.06.2006
Сообщений: 43,301
06.09.2011, 13:49
Привет! Вот еще темы с ответами:

Разложение на слагаемые - C++
На входе у нас число (нат, пол) которое нужно разложить и ожидаймое количество слагаймых алгоритм решения таков..выделяем место для...

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

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

Разложение в ряд. - C++
Нужна помощь.Вычисление логарифма по основанию два вещественного числа...дан ряд n(x+1)=x-x^2/2+x^3/3-.....x^n/n - разложение натурального...


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

Или воспользуйтесь поиском по форуму:
Yandex
Объявления
06.09.2011, 13:49
Ответ Создать тему
Опции темы

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