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

LU - разложение (выдает NaN) - C++

Восстановить пароль Регистрация
 
ISunRiseI
0 / 0 / 0
Регистрация: 24.12.2012
Сообщений: 4
26.12.2012, 21:16     LU - разложение (выдает NaN) #1
Здравствуйте. Помогите исправить программу. Проблема заключается в том, что некоторые значение получаются NaN/
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
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
#include<iostream.h>
#include<conio.h>
#include<iomanip.h>
int main()
{
    int i,j,k,n;
    float a[10][10],b[10],x[10],l[10][10],z[10],u[10][10],sum;
    cout<<"Enter number of unknowns:";
    cin>>n;
    cout<<"Enter co-efficients "<<endl;
    for(i=0;i<n;i++)
    for(j=0;j<=n;j++)
    {
        if(j<n) cin>>a[i][j];
        else cin>>b[i];
    }
    for(i=0;i<n;i++)
    for(j=0;j<n;j++)
        if(i>j) u[i][j]=0;
        else if(i==j) l[i][j]=1;
        else l[i][j]=0;
    //--------------------------------------------------------------------
    //AX=B
    //Decomposition of A into L and U   
    for(i=0;i<n;i++)
    {
        for(j=0;j<n;j++)
        {
            sum=0;
            if(i<=j)
            {
                for(k=0;k<n;k++) 
                    if(k!=i)
                        sum=sum+l[i][k]*u[k][j];
                u[i][j]=a[i][j]-sum;
            }
            else 
            {
                for(k=0;k<n;k++)
                    if(k!=j)
                        sum=sum+l[i][k]*u[k][j];
                    l[i][j]=(a[i][j]-sum)/u[j][j];
            }
        }
    }
    cout.precision(3);
    cout<<"\n\n  L= ";
    for(i=0;i<n;i++)
    {
        for(j=0;j<n;j++)
             cout<<setw(6)<<l[i][j];
        cout<<endl;
        cout<<"     ";
    }
    cout<<"\n\n  U= ";
    for(i=0;i<n;i++)
    {
        for(j=0;j<n;j++)
             cout<<setw(6)<<u[i][j];
        cout<<endl;
        cout<<"     ";
    }
    //--------------------------------------------------------------------
    //Now LUX=B
    //i.e LZ=B
    //Obtaining Z by forward subtitiution 
    z[0]=b[0]/l[0][0];
    for(i=1;i<n;i++)
    {
        sum=0;
        for(j=0;j<i;j++)
        {
            sum=sum+z[j]*l[i][j];
            
        }
        z[i]=b[i]-sum;
    }
    
    cout<<"\n\n  Z= ";
    for(i=0;i<n;i++)
    {
        cout<<setw(6)<<z[i];
        cout<<endl;
        cout<<"     ";
    }
    //--------------------------------------------------------------------
    //Now UX=Z
    //Finding X by backward substitution
    x[n-1]=z[n-1]/u[n-1][n-1];
    for(i=n-2;i>=0;i--)
    {
        sum=0;
        for(j=n-1;j>i;j--)
            sum=sum+x[j]*u[i][j];
        x[i]=(z[i]-sum)/u[i][i];
 
    }
    cout.precision(6);
    cout<<"\nThe solution is :"<<endl;
    for(i=0;i<n;i++)
        cout<<"\n x"<<i+1<<" = "<<x[i];
    cout<<endl;
    return 0;
}
После регистрации реклама в сообщениях будет скрыта и будут доступны все возможности форума.
Kuzia domovenok
 Аватар для Kuzia domovenok
1882 / 1737 / 116
Регистрация: 25.03.2012
Сообщений: 5,907
Записей в блоге: 1
26.12.2012, 21:19     LU - разложение (выдает NaN) #2
Цитата Сообщение от ISunRiseI Посмотреть сообщение
l[i][j]=(a[i][j]-sum)/u[j][j];
строка 42 проверь, не равно ли какое-нибудь u[j][j] равно нулю возможно деление на ноль.
строка 95 тоже
ISunRiseI
0 / 0 / 0
Регистрация: 24.12.2012
Сообщений: 4
26.12.2012, 21:26  [ТС]     LU - разложение (выдает NaN) #3
Проверил, деления на ноль нету.
Kuzia domovenok
 Аватар для Kuzia domovenok
1882 / 1737 / 116
Регистрация: 25.03.2012
Сообщений: 5,907
Записей в блоге: 1
26.12.2012, 21:43     LU - разложение (выдает NaN) #4
покажи пример системы, который выдаёт ошибку.
ISunRiseI
0 / 0 / 0
Регистрация: 24.12.2012
Сообщений: 4
26.12.2012, 21:49  [ТС]     LU - разложение (выдает NaN) #5
Ну например вот эта http://www.cleverstudents.ru/theory/...ations/002.png
Kuzia domovenok
 Аватар для Kuzia domovenok
1882 / 1737 / 116
Регистрация: 25.03.2012
Сообщений: 5,907
Записей в блоге: 1
27.12.2012, 00:11     LU - разложение (выдает NaN) #6
Вот такое разложение работает(фрагмент)
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
   int i,j,k,n;
    float a[10][10], b[10],x[10],l[10][10],z[10],u[10][10],sum;
    //n=3;
    //float a[3][3]={{2,3,-1},{1,-2,1},{1,0,2}};
    cout<<"Enter number of unknowns:";
    cin>>n;
    cout<<"Enter co-efficients "<<endl;
    for(i=0;i<n;i++)
    for(j=0;j<=n;j++)
    {
        if(j<n) {
            cin>>a[i][j];
            l[i][j]=u[i][j]=0.0;
        }
        else cin>>b[i];
    }
    for(j=0;j<n;j++){
        u[0][j]=a[0][j];
        l[j][j]=1.0f;
    }
    for(j=1;j<n;j++)
        l[j][0]=a[j][0]/u[0][0];
    for(i=1;i<n;i++){
        for(j=i;j<n;j++){
            u[i][j]=a[i][j];
            for (k=0; k<i; k++)
                u[i][j]-=l[i][k]*u[k][j];
        }
        for(j=i+1; j<n;j++){
            l[j][i]=a[j][i];
            for (k=0; k<i; k++)
                l[j][i]-=l[j][k]*u[k][i];
            l[j][i]=l[j][i]/u[i][i];
        }
    }
    cout.precision(3);
    cout<<"\n\n  L= ";
    for(i=0;i<n;i++)
    {
        for(j=0;j<n;j++)
             cout<<'\t'<<l[i][j];
        cout<<endl;
        cout<<"     ";
    }
    cout<<"\n\n  U= ";
    for(i=0;i<n;i++)
    {
        for(j=0;j<n;j++)
             cout<<'\t'<<u[i][j];
        cout<<endl;
        cout<<"     ";
    }
Yandex
Объявления
27.12.2012, 00:11     LU - разложение (выдает NaN)
Ответ Создать тему
Опции темы

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