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

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

Восстановить пароль Регистрация
 
Рейтинг: Рейтинг темы: голосов - 54, средняя оценка - 4.72
DVJ
1 / 1 / 0
Регистрация: 30.05.2011
Сообщений: 14
06.09.2011, 11:59     LU-разложение #1
Решить матрицу методом 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++
разложение тангенса C++
Разложение даты C++
Неполное LU разложение C++
LU разложение матрицы C++
После регистрации реклама в сообщениях будет скрыта и будут доступны все возможности форума.
sandye51
программист С++
 Аватар для sandye51
677 / 579 / 39
Регистрация: 19.12.2010
Сообщений: 2,016
06.09.2011, 12:09     LU-разложение #2
так решение написано, правда не проверял верное или нет
что еще надо?
DVJ
1 / 1 / 0
Регистрация: 30.05.2011
Сообщений: 14
06.09.2011, 13:22  [ТС]     LU-разложение #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. ну и по мелочам так
рассмешил

Yandex
Объявления
06.09.2011, 13:49     LU-разложение
Ответ Создать тему
Опции темы

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