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

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

Войти
Регистрация
Восстановить пароль
 
Nazz
WEB-developer
896 / 727 / 11
Регистрация: 12.03.2009
Сообщений: 2,802
Записей в блоге: 2
#1

Метод однократной интерполяции (ошибка в реализации) - C++

03.12.2010, 12:18. Просмотров 799. Ответов 1
Метки нет (Все метки)

имеется схема алгоритма, метода однократной интерполяции.
нужно реализовать ету схему используя с++.
вот пишу код:
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
#include <iostream.h>
#include <conio.h>
#include <string.h>
#include <math.h>
const e=0.0001;
const n=100;
 
void main(){
int a;
int b;
cout<<"Vvvedit znachennya zminnoi a: ";
cin>>a;
cout<<"\nVvedit znachnnya zminnoi b ";
cin>>b;
int k=0,c=0;
int xmin,x0,x1,x2,x3,x4,x;
float o,y0,y1,y2,y3,y4,S1,S2,f2,f1,y;
x1=a;
o=(a+b)/n;
 
//punct 1
p1:
x2=x1+o;
y1=4*pow(x1,2)+19*x1+14;
y2=4*pow(x2,2)+19*x2+14;
 
//perevirka ymovu
if (y1>=y2)
{
    goto p2;
}
 
 
//perevirka na povtorne vukonannya
if (k==0)
{
    if (y1<y2)
    {
        k=1;
        o=0-o;
        goto p1;
    }
}
if (k==1)
{
    if (y1<y2)
    {
        k=2;
        o=1/2*o;
        goto p1;
    }
}
if (k==2)
{
    goto p2;
}
 
 
//punct 2
p2:
o=2*o;
x3=x2+o;
y2=4*pow(x2,2)+19*x2+14;
y3=4*pow(x3,2)+19*x3+14;
if (y2>y3)
{
    goto p2;
}
if (y3>=y2)
{
    goto p3;
}
 
 
 
//punct 3
p3:
x4=x3-o;
y4=4*pow(x4,2)+19*x4+14;
x1=x1;
x2=x2;
x3=x3;
y1=y1;
y2=y2;
y3=y3;
 
 
//punct 5
if(y2>=y4)
{
    x1=x2;
    x2=x3;
    x3=x4;
}
 
 
//punct 6
if(y2<=y4)
{
    x1=x1;
    x2=x2;
    x3=x4;
}
 
 
//punct 7
S1=o*(y1-y3);
S2=2*(y1-2*y2+y3);
xmin=x2-S1/S2;
 
 
//punct 9
if (xmin<a)
{
    x=a;
}
 
if (xmin>b)
{
    x=b;
}
 
//punct 10
if (c==0)
{
    y=4*pow(x,2)+19*x+14;
}
if (c!=0)
    f1=y;
    y=4*pow(x,2)+19*x+14;
    f2=y;
    if ((f2-f1)<=e)
    {
        goto p4;
    }
}
c=c+1;
 
 
 
//punct 11
o=1/2*o;
x1=x;
goto p1;
 
 
 
 
p4:
cout<<"Kinets!";
 
}

получается зацыкливание после р2:
по логике я понимаю, что там и будет зацыкливание, но может я что-то не так реализовал(((...
подскажыте пожалуста?


текстовый алгоритм
делал не по схеме, а по текстовом описании, вот сейчас напишу его кратко по пунктах, так как и делал:
1.
х1=х0;
х2=х1+о;
f(x1);
f(x2);
естли f(x1)>f(x2) то до пункта 2
естли f(x2)>f(x1) то o=-o опять до пункта 1
естли опять f(x2)>f(x1) то o=1/2*o

2.
x3=x2+o;
естли f(x3)<f(x2) то o=2*o и опять етот же пункт, иначе до пункта 3.

3.
x4=x3-o;

4.
потом с 4-х точек выбираем наиболее отдалённую
естли f(x2)>=f(x4) то отбрасывается x1; точки переназываются
естли f(x4)>=f(x2) то отбрасывается x3; точки переназываются

5. после етого по етим точкам делается квадратическая интерполяция
s1=o*(f(x1)-f(x3));
s2=2*(f(x1)-2*f(x2)+f(x3));
xmin=x2-s1/s2;

6.
проверяем принадлежыт ли знайденое значение отрезку [a,b]

7. естли xmin<a => x=a;
естли xmin>b => x=b;

8. f(x)

9. ети вычисления продолжаются пока (f(x+1)-f(x))<=e;
естли условие не выполняется, то идём к первому пункту, и продожаем при чём x(1)=x и o=1/2*o.


воопще нужно минимизовать функцию методом однократной интерполяции
функция f=4x^2+19x+14
e=0.0001
промежуток: вводится пользователём(любой).
Миниатюры
Метод однократной интерполяции (ошибка в реализации)  
Similar
Эксперт
41792 / 34177 / 6122
Регистрация: 12.04.2006
Сообщений: 57,940
03.12.2010, 12:18     Метод однократной интерполяции (ошибка в реализации)
Посмотрите здесь:

метод линейной интерполяции - C++
Нужно решить уравнение sin(x)+log(x)=1 методом методом линейной интерполяции Помогите найти ошибку. Стока уже дней мучаюсь.. #...

Метод хорд - проверить корректность реализации кода - C++
Приветствую, форумчане. Хочу обратиться к вам за помощью, вот еле написала программу для метода хорд, методом шары если честно, так как...

Ошибка в реализации вектора - C++
здравствуйте, прошу помочь найти ошибку в простой реализации вектора. ошибка возникает когда вызывается pop_back(), однако совершенно не...

Ошибка в реализации связного списка - C++
Здравствуйте. Делаю заголовочный файл связного списка. В результате компиляции выдает ошибку: &quot;List is not a class template&quot;. Не понимаю в...

Ошибка при реализации стека - C++
Здравствуйте, помогите пожалуйста разобраться в ошибке.Пишу класс, в нём есть динамический символьный массив, вызываю функцию void без...

Ошибка при реализации калькулятора - C++
Собственно, пишу первую программу (банальный калькулятор), принцип его работы такой: Пользователь выбирает нужное ему действие с числами,...

Ошибка в реализации длинной арифметики - C++
Здравствуйте. Я скопировал с e-maxx'а и объединил всё в одну программу: #include &lt;stdio.h&gt; #include &lt;math.h&gt; #include &lt;cstdio&gt; ...

Ошибка в реализации задачи о ходе коня - C++
#include &quot;stdafx.h&quot; #include &lt;iostream&gt; using namespace std; int steps(int MAS, int N, int x, int y, int counter) { switch...

Ошибка в реализации хеш-функции SHA1 - C++
Здравствуйте, Решил написать простую реализацию, но результат вычислений оставляет желать лучшего Подскажите, пожалуйста, какой нюанс я...

Ошибка в реализации метода, проверяющего правильность данных. - C++
Вводиться ISBN номер. Метод реализован ckeck_ISBN(string isbn)! Нужно проверить правильность введенного номера по шаблону xx-xx-xx-c. xx...

Ошибка: функция уже имеет текст реализации - C++
Здравствуйте! Помогите, пожалуйста разобраться с ошибкой. При компиляции, компилятор выдаёт &quot;функция &quot;One::One(int,int)&quot; уже имеет...

Ошибка при использовании итератора в реализации класса - C++
На скрине всё понятно чо каво. Собственно, что я упустил? Что не дочитал? Проясните, пожалуйста. Из-за медленного интернета не могу...


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

Или воспользуйтесь поиском по форуму:
После регистрации реклама в сообщениях будет скрыта и будут доступны все возможности форума.
Nazz
WEB-developer
896 / 727 / 11
Регистрация: 12.03.2009
Сообщений: 2,802
Записей в блоге: 2
05.12.2010, 11:11  [ТС]     Метод однократной интерполяции (ошибка в реализации) #2
никто не подскажет как можна решыть ету проблему?
Ответ Создать тему
Опции темы

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