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

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

Войти
Регистрация
Восстановить пароль
 
anton-mix
0 / 0 / 0
Регистрация: 02.11.2011
Сообщений: 68
#1

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

08.03.2012, 16:53. Просмотров 297. Ответов 0
Метки нет (Все метки)

задача такова
решить систему нелинейных уравнений методом Ньютона
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++
a1x+b1y+c1=0 a2x+b2y+c2=0 Вывести решение системы (два числа – x и y) с шестью знаками после запятой. Если единственного решения не...

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

Решение Системы уравнений - C++
У меня задача - решать систему уравнений. матрица системы имеет порядок 3600. В будущем придется решать системы порядка 100 000. Мне дана...

Коэффициенты системы линейных уравнений - C++
Коэффициенты системы линейных уравнений заданы в виде прямоугольной матрицы. С помощью допустимых преобразований привести систему к...

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

Надоела реклама? Зарегистрируйтесь и она исчезнет полностью.
MoreAnswers
Эксперт
37091 / 29110 / 5898
Регистрация: 17.06.2006
Сообщений: 43,301
08.03.2012, 16:53
Привет! Вот еще темы с ответами:

Массивы. Системы линейных уравнений - C++
Нужна помощь с решение этой задачи с применением двумерного массива. Можно у кого что похожее есть\было? Не понимаю вообще как это...

Решение системы нелинейных уравнений - C++
найти решения уравнения графическим методом и методом поиска решений {cos(y+0,5)+x=0,8; { sinx-2y=1,6. P.S. только там скобка...

Решение системы линейных уравнений - C++
Помогите решить на Си

Решение системы линейных уравнений - C++
Не могу сделать прогу на С++. Кто сделает - огромное спасибо. Решить систему линейных уравнений 3x-5y+2z=2, 3x-7y+2z=0, x+5y-4z=-2


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

Или воспользуйтесь поиском по форуму:
Ответ Создать тему
Опции темы

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