Задача по методам оптимизации

@Starget-93 · Регистрация: 12.11.2012

Студворк — интернет-сервис помощи студентам

Всем доброго времени суток!
Необходимо запрограммировать метод тяжелого шарика.
Заранее спасибо!

@IrineK · 22.03.2013, 10:37

Оптимум какой функции искать?

@Starget-93 · 22.03.2013, 10:58 **[ТС]**

1 * x1^3 + 2 * x2^2 - 3 * x1 - 4 * x2

@IrineK · 22.03.2013, 13:23

Алгоритм и блок-схема на с.89-92

@Starget-93 · 22.03.2013, 13:41 **[ТС]**

Спасибо !

@IrineK · 23.03.2013, 04:38

Минимум - в точке (1; 1)

@Ded_Vasilij · 23.03.2013, 07:51

за книжечку спасибо

@kdk96 · 16.11.2015, 17:33

IrineK, Вы не могли бы скинуть текст программы, моя программа высчитывает приблизительно но почему то не так точно и всего за 4 итерации а не как у вас . Хотя писал код по методичке

Добавлено через 3 минуты
IrineK, вот мой код

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
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
153
154
155
156
157
158
159
160
161
162
163
#include "stdafx.h"
#include <cmath>
 
bool sravn(double x[], double y[], double eps, int n)
{
    int i;
    int k = 0;
    bool s = true;
    for (i = 0; i<n; i++)
        if (fabs(x[i] - y[i]) <= eps)
        {
            k++;
        }
    if (k == n) s = false;
    return s;
}
 
double pr(double (*f) (double*) , double eps, double *x, int i, int n)
{
    double f1 = 0, f2 = 0;
    double *x1 = new double[n];
    double *x2 = new double[n];
    for (int j = 0; j < n; j++)
    {
        x1[j] = x2[j] = x[j];
    }
    double dx = 0.1;
    x1[i] += dx;
    x2[i] -= dx;
    f1 = (f(x1) - f(x2)) / (2 * dx);
    dx = dx / 10;
    for (int j = 0; j < n; j++)
    {
        x1[j] = x2[j] = x[j];
    }
    x1[i] += dx;
    x2[i] -= dx;
    f2 = (f(x1) - f(x2)) / (2 * dx);
    while (fabs(f2 - f1) > eps)
    {
        f1 = f2;
        dx /= 10;
        for (int j = 0; j < n; j++)
        {
            x1[j] = x2[j] = x[j];
        }
        x1[i] += dx;
        x2[i] -= dx;
 
        f2 = (f(x1) - f(x2)) / (2 * dx);
        if (f2 == 0)
            break;
    }
    return f1;
}
 
//double *proizv(double *x, int n, double eps, double(*fu) (double*))
//{
//  double *f = new double[n];
//  for (int i = 0; i < n; i++)
//  {
//      f[i] = pr(fu, eps, x, i, n);
//  }
//  return f;
//}
 
double rosenbr(double *x)
{
    return (1 - x[0])*(1 - x[0]) + 100 * (x[1] - x[0] * x[0])*(x[1] - x[0] * x[0]);
}
 
double fo(double *x)
{
    return pow(x[0], 3) + 2 * pow(x[1], 2) - 3 * x[0] - 4 * x[1];
}
double fo1(double *x)
{
    return x[0] * x[0]-4;
}
double fo2(double *x)
{
    return 3 * x[0] * x[0] + x[0] * x[1] + 2 * x[1] * x[1] - x[0] - 4 * x[1];
}
double fo3(double *x)
{
    return -x[0] * x[0] - 5 * x[1] * x[1] - 3 * x[2] * x[2] + x[0] * x[1] - 2 * x[0] * x[2] + 2 * x[1] * x[2] + 11 * x[0] + 2 * x[1] + 18 * x[2] + 10;
}
 
double pr11(double x)
{
    return 3 * x*x - 3;
}
 
double pr22(double x)
{
    return 4 * x - 4;
}
int main()
{
    double(*f) (double*);
    f = &fo;
    double ff;
    double eps = 0.00000001;
    int n=2;
    int k=0, p=0, iter=0;
    double a=0.5, b=0.2;
    double *T = new double[n];
    double *x0 = new double[n];
    for (int i = 0; i < n; i++)
    {
        scanf("%lf", &x0[i]);
    }
    double *x1 = new double[n];
    for (int i = 0; i < n; i++)
    {
        scanf("%lf", &x1[i]);
    }
    double *x = new double[n];
    double f0, f1;
    f0 = f(x0);
    f1 = f(x1);
    //printf("%lf %lf", f0, f1);
    if (f1 > f0)
    {
        T = x0;
        x0 = x1;
        x1 = T;
    }
 
    do
    {
        if (k == 1)
        {
            x0 = x1;
            x1 = x;
        }
        k = 0;
        iter++;
        for (int i = 0; i < n; i++)
        {
            x[i] = x1[i] - a*pr(f, 0.000001, x1, i, n)+b*(x1[i]-x0[i]);
        }
        //x[0]=x1[0]-a*pr11(x1[0])+b*(x1[0] - x0[0]);
        //x[1] = x1[1] - a*pr11(x1[1]) + b*(x1[1] - x0[1]);
        ff = f(x);
        if (ff > f1 && p == 0 && b > 0.1)
            b -= 0.1;
        else
        {
            if (b <= 0.1) p = 1;
            if (ff > f1 && p == 1 && a >= eps) a /= 2;
            if (ff<f1) k = 1;
        }
    } while (sravn(x1,x,eps,n) && a>eps);
    printf("%lf %lf ", x[0], x[1]);
    printf("%lf   iter %d  a %lf  b %lf", f(x), iter, a ,b);
 
 
 
 
    system("pause");
    return 0;
}

Заранее спасибо!

Новые блоги и статьи Все статьи Все блоги /
Вывод данных через динамический список в справочнике Maks 01.04.2026 Реализация из решения ниже выполнена на примере нетипового справочника "Спецтехника" разработанного в конфигурации КА2. Задача: вывести данные из ТЧ нетипового документа. . .	Функция заполнения текстового поля в реквизите формы документа Maks 01.04.2026 Алгоритм из решения ниже реализован на нетиповом документе "ВыдачаОборудованияНаСпецтехнику" разработанного в конфигурации КА2, в дополнении к предыдущему решению. На форме документа создается. . .	К слову об оптимизации kumehtar 01.04.2026 Вспоминаю начало 2000-х, университет, когда я писал на Delphi. Тогда среди программистов на форумах активно обсуждали аккуратную работу с памятью: нужно было следить за переменными, вовремя. . .	Идея фильтра интернета (сервер = слой+фильтр). Hrethgir 31.03.2026 Суть идеи заключается в том, чтобы запустить свой сервер, о чём я если честно мечтал давно и давно приобрёл книгу как это сделать. Но не было причин его запускать. Очумелые учёные напечатали на. . .
Модель здравосоХранения 6. ESG-повестка и устойчивое развитие; углублённый анализ кадрового бренда anaschu 31.03.2026 В прикрепленном документе раздумья о том, как можно поменять модель в будущем	10 пpимет, которые всегда сбываются Maks 31.03.2026 1. Чтобы, наконец, пришла маршрутка, надо закурить. Если сигарета последняя, маршрутка придет еще до второй затяжки даже вопреки расписанию. 2. Нaдоели зима и снег? Не надо переезжать. Достаточно. . .	Перемещение выделенных строк ТЧ из одного документа в другой Maks 31.03.2026 Реализация из решения ниже выполнена на примере нетипового документа "ВыдачаОборудованияНаСпецтехнику" с единственной табличной частью "ОборудованиеИКомплектующие" разработанного в конфигурации КА2. . . .	Functional First Web Framework Suave DevAlt 30.03.2026 Sauve. IO Апнулись до NET10. Из зависимостей один пакет, работает одинаково хорошо как в режиме проекта так и в интерактивном режиме. из сложностей - чисто функциональный подход. Решил. . .

@Starget-93 0 / 0 / 0 Регистрация: 12.11.2012 Сообщений: 13

	Задача по методам оптимизации 22.03.2013, 09:04. Показов 2980. Ответов 7 Метки нет (Все метки) Всем доброго времени суток! Необходимо запрограммировать метод тяжелого шарика. Заранее спасибо! 0

@IrineK 2023 / 1641 / 425 Регистрация: 23.02.2011 Сообщений: 6,002 Записей в блоге: 25
	22.03.2013, 10:37
	Оптимум какой функции искать? 0

@Starget-93 0 / 0 / 0 Регистрация: 12.11.2012 Сообщений: 13
	22.03.2013, 10:58 [ТС]
	1 * x1^3 + 2 * x2^2 - 3 * x1 - 4 * x2 0

@Starget-93 0 / 0 / 0 Регистрация: 12.11.2012 Сообщений: 13
	22.03.2013, 13:41 [ТС]
	Спасибо ! 0

@IrineK 2023 / 1641 / 425 Регистрация: 23.02.2011 Сообщений: 6,002 Записей в блоге: 25
	23.03.2013, 04:38
	Минимум - в точке (1; 1) Миниатюры 0

@Ded_Vasilij 233 / 215 / 63 Регистрация: 01.09.2012 Сообщений: 2,103
	23.03.2013, 07:51
	за книжечку спасибо 0