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

системы нелинейныx уравнений - C++

Восстановить пароль Регистрация
 
anton-mix
0 / 0 / 0
Регистрация: 02.11.2011
Сообщений: 68
08.03.2012, 16:53     системы нелинейныx уравнений #1
задача такова
решить систему нелинейных уравнений методом Ньютона
F(x,y)=x^3-4*x+3-y^2
G(x,y)=x-(y-2)^2+1

Проблема состоит в том что не заполняется массив коэффициентов C[i]=W[i][j]*F[j]
где W[i][j]-обратная матрица матрицы Якоби
F[j]- численный вектор составленный из F(x,y) и G(x,y) в точке x0,y0

помогите пожалуйста с проблемой=)
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
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
#include<cmath>
#include<iostream>
using namespace std;
const int n=2;
//Считаем елементы матрицы w при нач приближении
double DF_DX(double x)
{
    return 3*x*x-4;
}
double DF_DY(double y)
{
    return -2*y;
}
double DG_DX(double x)
{
    return 1;
}
double DG_DY(double y)
{
    return -2*(y-2);
}
 
 
//Считаем детерминант
double determinant(double x, double y)
{
    return (DF_DX(x))*(DG_DY(y))-(DG_DX(x))*(DF_DY(y));
};
 
 
//Считаем элементы обратной матрицы
double W11(double x,double y)
{
    return DG_DY(y)/determinant(x,y);
}
double W12(double x,double y)
{
    return -DF_DY(y)/determinant(x,y);
}
double W21(double x,double y)
{
    return -DG_DX(x)/determinant(x,y);
}
double W22(double x,double y)
{
    return -DF_DX(x)/determinant(x,y);
}
 
double F1(double x, double y)
{
    return x*x*x-4*x+3-y*y;
}
double F2(double x, double y)
{
    return x-(y-2)*(y-2)+1;
}
 
double matrica(double x, double y)
{
    double W[n][n];
    double F[n];
    W[0][0]=W11(x,y);
    W[0][1]=W12(x,y);
    W[1][0]=W21(x,y);
    W[1][1]=W22(x,y);
    for(int i=0;i<2;i++)
    {
        for (int j=0;j<2;j++)
        {
            cout<<W[i][j]<<" \n";
        }
    }
    F[0]=F1(x,y);
    F[1]=F2(x,y);
 
    cout<<"znachenir function v x,y\n";
    for(int j=0;j<2;j++)
    {
        cout<<F[j]<<" \n";
    }
 
    double C[n];
    for (int i=0;i<2;i++)
    {
        
        for(int j=0;j<2;j++)
        {
            C[i]+=W[i][j]*F[j];
            cout<<C[i]<<" \n";
        }
    }
    /*
double x_new[n];
    for(int i=0;i<n;i++)
    {
        x_new[i]=x[i]-C[i];
    }
*/
    return 0;
}
 
int main()
{
    double x=2.9;
    double y=3.98;
    cout<<"koeff matrix W"<<endl;
    cout<<DF_DX(x)<<endl;
    cout<<DF_DY(y)<<endl;
    cout<<DG_DX(x)<<endl;
    cout<<DG_DY(y)<<endl;
    cout<<"determinant";
    cout<<determinant(x,y)<<endl;
    cout<<"koeff obr matrix W"<<endl;
    cout<<W11(x,y)<<endl;
    cout<<W12(x,y)<<endl;
    cout<<W21(x,y)<<endl;
    cout<<W22(x,y)<<endl;
    cout<<"koeff function matrix\n";
    matrica(x,y);
 
    return 0;
}
Similar
Эксперт
41792 / 34177 / 6122
Регистрация: 12.04.2006
Сообщений: 57,940
08.03.2012, 16:53     системы нелинейныx уравнений
Посмотрите здесь:

Коэффициенты системы линейных уравнений C++
C++ Решение системы линейных уравнений.
C++ Решение Системы уравнений
C++ Решение системы линейных уравнений
C++ Расчет системы уравнений
C++ Решение системы нелинейных уравнений
Найти решение системы уравнений C++
Решение системы алгебраических уравнений C++

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

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

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