Минимизация методом парабол. (перевести из c++ в матлаб)

@Mixirist · Регистрация: 27.06.2014

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

Вот программа.

Кликните здесь для просмотра всего текста

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
#include<iostream.h>
#include<math.h>
#include<conio.h>
#include<stdio.h>
#include "stdafx.h"
double f(double x)
{double f;
f=pow(x, 3)-x+exp(-x);
return f;
}
void main()
{double a=0, b=1, c, ya, yb, yc, s, t, yt, x, y;
const double e=0.00001;
int N;
cout<<"Method parabol "<<endl<<endl;
do
{cout<<"Vv. c ";
cin>>c;
}while ((a>=c) || (b<=c) || (f(a)<=f(c)) || f(b)<=f(c));
ya=f(a);
yb=f(b);
yc=f(c);
N=3;//kol-vo experimentov
while (b-a>2*e)
  {
   s=c+0.5*( (b-c)*(b-c)*(ya-yc) - (c-a)*(c-a)*(yb-yc) )/( (b-c)*(ya-yc) + (c-a)*(yb-yc) );
   if (s==c) t=(a+c)/2;   //see here
     else t=s;
   yt=f(t); N++;
 
   if (t<c)
     {
     if (yt<yc) {b=c; yb=yc; c=t; yc=yt;}
       else if (yt>yc) {a=t; ya=yt;}
        else {a=t; ya=yt; b=c; yb=yc; c=(a+b)/2; yc=f(c); N++;}
      }
   else if (t>c)
      {
       if (yt<yc) {a=c; ya=yc; c=t; yc=yt;}
          else if(yt>yc) {b=t; yb=yt;}
         else {a=c; ya=yc; b=t; yb=yt; c=(a+b)/2; yc=f(c); N++; }
         }
   }
x=(a+b)/2; y=f(x);
 
printf("\nx = %.5f \n", x);
printf("y = %.5f \n \n", y);
cout<<"Kol-vo tochek: "<<N;
}

Буду признателен если все получится.

Эх, забыли про мою темку)))

Зосима · 30.06.2014, 11:59

держи

Matlab M
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
clear, clc
 
f = @(x) x.^3 - x + exp(-x);
 
a=0;
b=1;
ezplot(f,[a, b]), hold on
 
e=0.00001;
disp('Метод парабол')
c = (a+b)/2;
% c = input('Введите с: ');
% while (a>=c) || (b<=c) || (f(a)<=f(c)) || (f(b)<=f(c))
%     c = input('Введите с: ');
% end
 
ya=f(a);
yb=f(b);
yc=f(c);
N=3; % kol-vo experimentov
while b-a>2*e
    s=c+0.5*( (b-c)*(b-c)*(ya-yc) - (c-a)*(c-a)*(yb-yc) )/( (b-c)*(ya-yc) + (c-a)*(yb-yc) );
   if (s==c) 
       t=(a+c)/2;   % see here
   else
       t=s;
   end
   yt=f(t); 
   N = N+1;
  if (t<c)
    if (yt<yc) 
        b=c; yb=yc; c=t; yc=yt;
    elseif (yt>yc) 
        a=t; ya=yt;
    else
        a=t; ya=yt; b=c; yb=yc; c=(a+b)/2; yc=f(c); 
        N = N+1;
    end
      
  elseif (t>c)
      if (yt<yc) 
          a=c; ya=yc; c=t; yc=yt;
      elseif (yt>yc) 
          b=t; yb=yt;
      else
          a=c; ya=yc; b=t; yb=yt; c=(a+b)/2; yc=f(c);
              N=N+1; 
      end
  end
end
x=(a+b)/2; y=f(x);
 
fprintf('\nx = %.5f \n', x);
fprintf('y = %.5f \n \n', y);
fprintf('Kol-vo tochek: %i \n', N);
 
plot(x,y,'or'), grid on

@Mixirist · 02.07.2014, 08:33 **[ТС]**

Спасибо)))

Новые блоги и статьи Все статьи Все блоги /
Модульный подход на примере F# DevAlt 06.03.2026 В блоге дяди Боба наткнулся на такое определение: В этой книге («Подход, основанный на вариантах использования») Ивар утверждает, что архитектура программного обеспечения — это структуры,. . .	Управление камерой с помощью скрипта OrbitControls.js на Three.js: Вращение, зум и панорамирование 8Observer8 05.03.2026 Содержание блога Финальная демка в браузере работает на Desktop и мобильных браузерах. Итоговый код: orbit-controls-threejs-js. zip. Сканируйте QR-код на мобильном. Вращайте камеру одним пальцем,. . .	SDL3 для Web (WebAssembly): Синхронизация спрайтов SDL3 и тел Box2D 8Observer8 04.03.2026 Содержание блога Финальная демка в браузере. Итоговый код: finish-sync-physics-sprites-sdl3-c. zip На первой гифке отладочные линии отключены, а на второй включены:. . .	SDL3 для Web (WebAssembly): Идентификация объектов на Box2D v3 - использование userData и событий коллизий 8Observer8 02.03.2026 Содержание блога Финальная демка в браузере. Итоговый код: finish-collision-events-sdl3-c. zip Сканируйте QR-код на мобильном и вы увидите, что появится джойстик для управления главным героем. . . .
Реалии Hrethgir 01.03.2026 Нет, я не закончил до сих пор симулятор. Эта задача сложнее. Не получилось уйти в плавсостав, но оно и к лучшему, возможно. Точнее получалось - но сварщиком в палубную команду, а это значит, в моём. . .	Ритм жизни kumehtar 27.02.2026 Иногда приходится жить в ритме, где дел становится всё больше, а вовлечения в происходящее — всё меньше. Плотный график не даёт вниманию закрепиться ни на одном событии. Утро начинается с быстрых,. . .	SDL3 для Web (WebAssembly): Сборка библиотек: SDL3, Box2D, FreeType, SDL3_ttf, SDL3_mixer и SDL3_image из исходников с помощью CMake и Emscripten 8Observer8 27.02.2026 Недавно вышла версия 3. 4. 2 библиотеки SDL3. На странице официальной релиза доступны исходники, готовые DLL (для x86, x64, arm64), а также библиотеки для разработки под Android, MinGW и Visual Studio. . . .	SDL3 для Web (WebAssembly): Реализация движения на Box2D v3 - трение и коллизии с повёрнутыми стенами 8Observer8 20.02.2026 Содержание блога Box2D позволяет легко создать главного героя, который не проходит сквозь стены и перемещается с заданным трением о препятствия, которые можно располагать под углом, как верхнее. . .

@Mixirist 0 / 0 / 0 Регистрация: 27.06.2014 Сообщений: 7
	02.07.2014, 08:33 [ТС]
	Спасибо))) 0

Минимизация методом парабол. (перевести из c++ в матлаб)

Решение